- Timestamp:
- 11/19/09 11:29:48 (14 years ago)
- Branches:
- master
- Children:
- 7137e5d
- Parents:
- 32289cd
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
redvar.f
r32289cd r38d77eb 101 101 ntlvr=ivrml1(ntlml)+nvrml(ntlml)-1 102 102 if (ntlvr.eq.0) then 103 write (*,*) 'redvar> No variables defined in molecule(s)'103 write(logString,*)'redvar> No variables defined in molecule(s)' 104 104 return 105 105 endif … … 108 108 io=iopfil(lunvar,varfil,'old','formatted') 109 109 if (io.eq.0) then 110 write ( *,'(a,/,a,i3,2a)')110 write (logString, '(a,/,a,i3,2a)') 111 111 & ' redvar> ERROR opening file to set variables:', 112 112 & ' LUN=',lunvar,' FILE=',varfil(1:iendst(varfil)) … … 136 136 vr=blat(iatvr(i)) 137 137 else 138 write( *,*) 'redvar> unknown variable type: ',it,' !'138 write(logString,*) 'redvar> unknown variable type: ',it,' !' 139 139 stop 140 140 endif … … 159 159 if (iredin(lincmd,nml).le.0.or. 160 160 & nml.le.0.or.nml.gt.ntlml) then 161 write (*,*) 'redvar> ','Incorrect molecule number >', 161 write (logString, *) 162 & 'redvar> ','Incorrect molecule number >', 162 163 & lincmd(1:l),'< Must be in range [1,', 163 164 & ntlml,'] !' … … 310 311 311 312 if (inum.le.0.or.inum.gt.ntlml) then 312 write ( *,*) ' # 1: ',inum313 write (logString, *) ' # 1: ',inum 313 314 goto 104 314 315 endif … … 324 325 325 326 if (k.lt.nfi.or.k.gt.irsml2(j)) then 326 write ( *,*) ' # 2: ',inum327 write (logString, *) ' # 2: ',inum 327 328 goto 104 328 329 endif … … 343 344 if (l.lt.nfi.or. 344 345 & l.gt.nfi+nvrrs(k)-1) then 345 write ( *,*) ' # 3: ',inum346 write (logString, *) ' # 3: ',inum 346 347 goto 104 347 348 endif … … 427 428 428 429 if (ide.eq.1.and.ihy.eq.ide) then 429 write ( *,*) ' # 4: ',ide430 write (logString, *) ' # 4: ',ide 430 431 goto 104 431 432 endif … … 441 442 if (iredin(linh,ibz).le.0.or.ibz.le.0) 442 443 & then 443 write ( *,*) ' # 5 '444 write (logString, *) ' # 5 ' 444 445 goto 104 445 446 endif … … 456 457 if (iredin(linh,iez).le.0.or.iez.le.0.or. 457 458 & iez.lt.ibz) then 458 write ( *,*) ' # 6 '459 write (logString, *) ' # 6 ' 459 460 goto 104 460 461 endif … … 524 525 endif 525 526 enddo 526 if (.not.did) write ( *,'(3a)')527 if (.not.did) write (logString, '(3a)') 527 528 & ' redvar> No variables affected by command >', 528 529 & lincmd(1:ice),'<' … … 531 532 ll1=ibegst(linfld(nfld)) 532 533 ll2=iendst(linfld(nfld)) 533 write ( *,*) 'll1,ll2, linfld(nfld): ',ll1,ll2,534 write (logString, *) 'll1,ll2, linfld(nfld): ',ll1,ll2, 534 535 & '>',linfld(nfld)(ll1:ll2),'<' 535 536 … … 559 560 enddo 560 561 561 write (*,'(3a,/,1x,5(a,2x),a)') ' redvar> ',nmml(i)(1:ie), 562 write (logString, '(3a,/,1x,5(a,2x),a)') ' redvar> ', 563 & nmml(i)(1:ie), 562 564 & ' with global parameters:', 563 565 & (strg(k)(ibegst(strg(k)):),k=1,6) … … 579 581 did=.true. 580 582 if (fxvr(iv)) then 581 write (*,'(3a,i4,1x,4a,f10.3,a)') ' redvar> ',583 write(logString,'(3a,i4,1x,4a,f10.3,a)') ' redvar> ', 582 584 & nmml(i)(1:ie),': residue ',j-jb,seq(j), 583 585 & ': ',nmvr(iv),' set ',vlvrx(iv),' Fixed' 584 586 else 585 write ( *,'(3a,i4,1x,4a,f10.3)') ' redvar> ',587 write (logString, '(3a,i4,1x,4a,f10.3)') ' redvar> ', 586 588 & nmml(i)(1:ie),': residue ',j-jb,seq(j), 587 589 & ': ',nmvr(iv),' set ',vlvrx(iv) … … 597 599 enddo 598 600 if (did) then 599 if (in.gt.0) write ( *,'(3a,i5,a)')601 if (in.gt.0) write (logString, '(3a,i5,a)') 600 602 & ' redvar> Molecule ',nmml(i)(1:ie),': ',in, 601 603 & ' variable(s) remain unchanged' 602 604 call setvar(iml,vlvrx) 603 605 else 604 write ( *,'(3a)') ' redvar> Molecule ',606 write (logString, '(3a)') ' redvar> Molecule ', 605 607 & nmml(i)(1:ie),': No internal variables changed' 606 608 endif … … 610 612 return 611 613 ! ____________________________________________________________ Errors 612 100 write ( *,'(3a)') ' redvar> Cannot interpret command >',614 100 write (logString, '(3a)') ' redvar> Cannot interpret command >', 613 615 & lincmd(1:ice),'<' 614 616 close(lunvar) 615 617 stop 616 101 write (*,'(a,i5,a)') ' redvar> Command number ',ncmd,' reached' 618 101 write (logString, '(a,i5,a)') ' redvar> Command number ',ncmd, 619 & ' reached' 617 620 close(lunvar) 618 621 stop 619 102 write ( *,'(3a)') ' redvar> Cannot read value from >',622 102 write (logString, '(3a)') ' redvar> Cannot read value from >', 620 623 & lincmd(1:ice),'<' 621 624 close(lunvar) 622 625 stop 623 103 write ( *,'(a,i3,3a)') ' redvar> Cannot read >',maxide,626 103 write (logString, '(a,i3,3a)') ' redvar> Cannot read >',maxide, 624 627 & ' identifiers from >',linfld(ii)(1:ife),'<' 625 628 close(lunvar) 626 629 stop 627 104 write ( *,'(5a)') ' redvar> Error in identifier >',630 104 write (logString, '(5a)') ' redvar> Error in identifier >', 628 631 & linide(1:ide),'< of command >',lincmd(1:ice),'<' 629 632 close(lunvar) 630 633 stop 631 105 write (*,'(a,/,a,/,2a,/)') ' redvar> line with global paramters:', 634 105 write (logString, '(a,/,a,/,2a,/)') 635 & ' redvar> line with global paramters:', 632 636 & line(1:ile),' must contain 6 floating', 633 637 & ' point numbers separated by commas !' … … 635 639 stop 636 640 637 106 write (*,'(a,/,a,/,2a,/)') ' redvar> line with global paramters:', 641 106 write (logString, '(a,/,a,/,2a,/)') 642 & ' redvar> line with global paramters:', 638 643 & line(1:ile),' angles must be inside ', 639 644 &'ranges [-180,180], [-90,90], and [-180,180] Deg., respectively !'
Note:
See TracChangeset
for help on using the changeset viewer.