Changeset bd2278d for EXAMPLES/parallel_tempering_s.f
- Timestamp:
- 09/05/08 11:49:42 (16 years ago)
- Branches:
- master
- Children:
- fafe4d6
- Parents:
- 2ebb8b6
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
EXAMPLES/parallel_tempering_s.f
r2ebb8b6 rbd2278d 20 20 integer switch 21 21 22 c=================================================== Energy setup22 ! =================================================== Energy setup 23 23 24 cDirectory for SMMP libraries25 cChange the following directory path to where you want to put SMMP26 clibraries of residues.24 ! Directory for SMMP libraries 25 ! Change the following directory path to where you want to put SMMP 26 ! libraries of residues. 27 27 libdir='../SMMP/' 28 28 29 cThe switch in the following line is now not used.29 ! The switch in the following line is now not used. 30 30 flex=.false. ! .true. for Flex / .false. for ECEPP 31 31 32 cChoose energy type with the following switch instead ...32 ! Choose energy type with the following switch instead ... 33 33 ientyp = 0 34 c0 => ECEPP2 or ECEPP3 depending on the value of sh235 c1 => FLEX36 c2 => Lund force field37 c3 => ECEPP with Abagyan corrections38 c 34 ! 0 => ECEPP2 or ECEPP3 depending on the value of sh2 35 ! 1 => FLEX 36 ! 2 => Lund force field 37 ! 3 => ECEPP with Abagyan corrections 38 ! 39 39 40 40 sh2=.false. ! .true. for ECEPP/2; .false. for ECEPP3 … … 48 48 call init_energy(libdir) 49 49 50 c================================================= Structure setup50 ! ================================================= Structure setup 51 51 52 52 grpn = 'nh2' ! N-terminal group … … 59 59 ntlml = 0 60 60 write (*,*) 'Solvent: ', itysol 61 cInitialize random number generator.61 ! Initialize random number generator. 62 62 call sgrnd(31433) 63 63 … … 69 69 70 70 call init_molecule(iabin,grpn,grpc,seqfile,varfile) 71 cDecide if and when to use BGS, and initialize Lund data structures71 ! Decide if and when to use BGS, and initialize Lund data structures 72 72 bgsprob=0.75 ! Prob for BGS, given that it is possible 73 cupchswitch= 0 => No BGS 1 => BGS with probability bgsprob74 c2 => temperature dependent choice73 ! upchswitch= 0 => No BGS 1 => BGS with probability bgsprob 74 ! 2 => temperature dependent choice 75 75 upchswitch=1 76 76 rndord=.true. … … 80 80 81 81 82 c======================================== Add your task down here82 ! ======================================== Add your task down here 83 83 num_rep = 5 84 84 nequi = 100 … … 87 87 newsta = .true. 88 88 switch = 1 89 cparallel tempering on a single CPU89 ! parallel tempering on a single CPU 90 90 eol = energy() 91 91 write (*,*) "Energy before randomization:", eol … … 94 94 write (*,*) "Final energy:", eol 95 95 96 c======================================== End of main96 ! ======================================== End of main 97 97 end
Note:
See TracChangeset
for help on using the changeset viewer.