- Timestamp:
- 09/05/08 11:49:42 (16 years ago)
- Branches:
- master
- Children:
- fafe4d6
- Parents:
- 2ebb8b6
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
main.f
r2ebb8b6 rbd2278d 1 c**************************************************************2 cThis file contains the: main (SINGLE PROCESSOR JOBS ONLY,3 CFOR PARALLEL JOBS USE pmain)4 c 5 cCopyright 2003-2005 Frank Eisenmenger, U.H.E. Hansmann,6 cShura Hayryan, Chin-Ku7 cCopyright 2007 Frank Eisenmenger, U.H.E. Hansmann,8 cJan H. Meinke, Sandipan Mohanty9 c 10 CCALLS: init_energy,init_molecule11 CCALLS TASK SUBROUTINE: anneal,canon,elp,minim,mulcan_par,12 cmulcan_sim,partem_s, or regul13 CCAN ALSO CALL MEASUREMENT ROUTINES: cnteny,contacts,helix,hbond,14 Coutpdb,outvar,rgyr,15 Crmsinit and rsmdfun,zimmer16 c$Id: main.f 334 2007-08-07 09:23:59Z meinke $17 c**************************************************************1 ! ************************************************************** 2 ! This file contains the: main (SINGLE PROCESSOR JOBS ONLY, 3 ! FOR PARALLEL JOBS USE pmain) 4 ! 5 ! Copyright 2003-2005 Frank Eisenmenger, U.H.E. Hansmann, 6 ! Shura Hayryan, Chin-Ku 7 ! Copyright 2007 Frank Eisenmenger, U.H.E. Hansmann, 8 ! Jan H. Meinke, Sandipan Mohanty 9 ! 10 ! CALLS: init_energy,init_molecule 11 ! CALLS TASK SUBROUTINE: anneal,canon,elp,minim,mulcan_par, 12 ! mulcan_sim,partem_s, or regul 13 ! CAN ALSO CALL MEASUREMENT ROUTINES: cnteny,contacts,helix,hbond, 14 ! outpdb,outvar,rgyr, 15 ! rmsinit and rsmdfun,zimmer 16 ! $Id: main.f 334 2007-08-07 09:23:59Z meinke $ 17 ! ************************************************************** 18 18 19 19 program main … … 26 26 logical lrand,bgsposs 27 27 28 c=================================================== Energy setup28 ! =================================================== Energy setup 29 29 30 cDirectory for SMMP libraries31 cChange the following directory path to where you want to put SMMP32 clibraries of residues.30 ! Directory for SMMP libraries 31 ! Change the following directory path to where you want to put SMMP 32 ! libraries of residues. 33 33 libdir='./SMMP/' 34 34 35 cThe switch in the following line is now not used.35 ! The switch in the following line is now not used. 36 36 flex=.false. ! .true. for Flex / .false. for ECEPP 37 37 38 cChoose energy type with the following switch instead ...38 ! Choose energy type with the following switch instead ... 39 39 ientyp = 0 40 c0 => ECEPP2 or ECEPP3 depending on the value of sh241 c1 => FLEX42 c2 => Lund force field43 c3 => ECEPP with Abagyan corrections44 c 40 ! 0 => ECEPP2 or ECEPP3 depending on the value of sh2 41 ! 1 => FLEX 42 ! 2 => Lund force field 43 ! 3 => ECEPP with Abagyan corrections 44 ! 45 45 46 46 sh2=.false. ! .true. for ECEPP/2; .false. for ECEPP3 … … 54 54 call init_energy(libdir) 55 55 56 c================================================= Structure setup56 ! ================================================= Structure setup 57 57 58 58 grpn = 'nh2' ! N-terminal group … … 67 67 ntlml = 0 68 68 write (*,*) 'Solvent: ', itysol 69 cInitialize random number generator.69 ! Initialize random number generator. 70 70 call sgrnd(31433) 71 71 … … 78 78 call init_molecule(iabin,grpn,grpc,seqfile,varfile) 79 79 80 cDecide if and when to use BGS, and initialize Lund data structures80 ! Decide if and when to use BGS, and initialize Lund data structures 81 81 bgsprob=0.75 ! Prob for BGS, given that it is possible 82 cupchswitch= 0 => No BGS 1 => BGS with probability bgsprob83 c2 => temperature dependent choice82 ! upchswitch= 0 => No BGS 1 => BGS with probability bgsprob 83 ! 2 => temperature dependent choice 84 84 upchswitch=1 85 85 rndord=.true. … … 89 89 90 90 91 c======================================== Add your task down here91 ! ======================================== Add your task down here 92 92 93 93 imin = 1 ! Quasi-Newton … … 95 95 eps = 1.0d-7 ! requested precision 96 96 call minim(imin, maxit, eps) 97 cTo do a canonical Monte Carlo simulation uncomment the lines below97 ! To do a canonical Monte Carlo simulation uncomment the lines below 98 98 ! nequi = 100 99 99 ! nsweep = 50000 … … 101 101 ! temp = 300.0 102 102 ! lrand = .true. 103 cCanonical Monte Carlo103 ! Canonical Monte Carlo 104 104 ! call canon(nequi, nsweep, nmes, temp, lrand) 105 105 106 cFor simulated annealing uncomment the lines below106 ! For simulated annealing uncomment the lines below 107 107 ! tmin = 200.0 108 108 ! tmax = 500.0 109 109 ! call anneal(nequi, nsweep, nmes, tmax, tmin, lrand); 110 c======================================== End of main110 ! ======================================== End of main 111 111 end
Note:
See TracChangeset
for help on using the changeset viewer.