- Timestamp:
- 11/19/09 11:29:48 (14 years ago)
- Branches:
- master
- Children:
- 7137e5d
- Parents:
- 32289cd
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
canon.f
r32289cd r38d77eb 69 69 eol = energy() 70 70 e_min = eol 71 write ( *,'(a,e12.5,/)') 'energy of start configuration:',eol71 write (logString, '(a,e12.5,/)') 'energy of start configuration:',eol 72 72 73 73 ! Write start configuration in pdb-format into file … … 79 79 call metropolis(eol,acz,can_weight) 80 80 end do 81 write (*,*) 'Energy after equilibration:',eol81 write (logString, *) 'Energy after equilibration:',eol 82 82 83 83 !======================Simulation in canonical ensemble … … 86 86 call metropolis(eol,acz,can_weight) 87 87 if (eol.lt.e_min) then 88 write ( *,*) "New minimum energy:", eol, "t=", nsw89 write ( *,*) eyel,eyhb,eyvr,eysl88 write (logString, *) "New minimum energy:", eol, "t=", nsw 89 write (logString, *) eyel,eyhb,eyvr,eysl 90 90 e_min = eol 91 91 call outpdb(0, "minen.pdb") … … 116 116 117 117 acz = acz/dble(nsw*nvr) 118 write (*,*) 'acceptance rate:',acz119 write (*,*)118 write (logString, *) 'acceptance rate:',acz 119 write (logString, *) 120 120 ! ------------ Output Dihedreals of final configuration 121 write (*,*) 'last energy',eol121 write (logString, *) 'last energy',eol 122 122 call outvar(0,'lastconf.var') 123 123 ! Output final conformation as pdb-file
Note:
See TracChangeset
for help on using the changeset viewer.