- Timestamp:
- 11/19/09 11:29:16 (14 years ago)
- Branches:
- master
- Children:
- 2d5d396
- Parents:
- ffd2d46
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
canon.f
rffd2d46 rcb47b9c 40 40 integer nmes 41 41 ! temp: Temperature of simulation 42 double precision temp 42 double precision temp, e_min 43 43 ! 44 44 ! common/bet/beta … … 63 63 64 64 eol = energy() 65 e_min = eol 65 66 write (*,'(a,e12.5,/)') 'energy of start configuration:',eol 66 67 … … 79 80 do nsw=0,nswp 80 81 call metropolis(eol,acz,can_weight) 82 if (eol.lt.e_min) then 83 write (*,*) "New minimum energy:", eol, "t=", nsw 84 write (*,*) eyel,eyhb,eyvr,eysl 85 e_min = eol 86 call outpdb(0, "minen.pdb") 87 endif 81 88 ! 82 89 if(mod(nsw,nmes).eq.0) then … … 96 103 end do 97 104 ! Write down information on actual conformation 98 write(13,'(i5,2f12.3,5i7)') nsw, eol, rgy, 99 & nhel,mhel,nbet,mbet,mhb 105 write(13,'(i7,2f15.3,5i7,4f15.3)') nsw, eol, rgy, 106 & nhel,mhel,nbet,mbet,mhb, 107 & eyel,eyhb,eyvr,eysl 100 108 end if 101 109 !
Note:
See TracChangeset
for help on using the changeset viewer.