- Timestamp:
- 11/19/09 11:29:48 (14 years ago)
- Branches:
- master
- Children:
- 7137e5d
- Parents:
- 32289cd
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
anneal.f
r32289cd r38d77eb 84 84 85 85 eol=energy() 86 write (*,'(a,e12.5,/)') 'energy of start configuration: ',eol 86 write (logString, '(a,e12.5,/)') 87 & 'energy of start configuration: ',eol 87 88 88 89 ! Write start configuration in pdb-format into file … … 94 95 call metropolis(eol,acz,can_weight) 95 96 end do 96 write(*,*) 'Energy after equilibration:',eol 97 write (logString, *) 98 & 'Energy after equilibration:', eol 97 99 98 100 !======================Simulation by simulated annealing … … 131 133 132 134 acz = acz/dble(nsw*nvr) 133 write (*,*) 'acceptance rate:',acz134 write (*,*)135 write (logString, *) 'acceptance rate:',acz 136 write (logString, *) 135 137 ! ------------ Output Dihedreals of final configuration 136 write (*,*) 'last energy',eol138 write (logString, *) 'last energy',eol 137 139 call outvar(0,' ') 138 140 ! Output final conformation as pdb-file 139 141 call outpdb(0,"final.pdb") 140 write (*,*)142 write (logString, *) 141 143 142 144 ! ------------ Output Dihedreals of conformation with lowest energy 143 write (*,*) 'lowest energy ever found:',nemin,ymin145 write (logString, *) 'lowest energy ever found:',nemin,ymin 144 146 close(14) 145 147 ! =====================
Note:
See TracChangeset
for help on using the changeset viewer.