Changeset 32289cd for init_energy.f
- Timestamp:
- 11/19/09 11:29:41 (14 years ago)
- Branches:
- master
- Children:
- 38d77eb
- Parents:
- 6650a56
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
init_energy.f
r6650a56 r32289cd 5 5 ! 6 6 ! Copyright 2003-2005 Frank Eisenmenger, U.H.E. Hansmann, 7 ! Shura Hayryan, Chin-Ku 7 ! Shura Hayryan, Chin-Ku 8 8 ! Copyright 2007 Frank Eisenmenger, U.H.E. Hansmann, 9 9 ! Jan H. Meinke, Sandipan Mohanty … … 17 17 ! PURPOSE: initialize energy parameters 18 18 ! 0 => ECEPP2 or ECEPP3 depending on the value of sh2 19 ! 1 => FLEX 19 ! 1 => FLEX 20 20 ! 2 => Lund force field 21 21 ! 3 => ECEPP with Abagyan corrections … … 29 29 include 'INCL.H' 30 30 31 integer ll, iendst, its 32 31 33 character libdir*(*),tesfil*80 32 34 33 35 if (ientyp.eq.1) then 34 36 flex = .true. 35 else 37 else 36 38 flex = .false. 37 39 end if 38 40 39 41 lunlib=10 40 42 ll=iendst(libdir) … … 56 58 endif 57 59 58 if (ientyp .eq. 2) then 60 if (ientyp .eq. 2) then 59 61 reslib=libdir(1:ll)//'lib.lun' 60 62 endif … … 78 80 79 81 else 80 82 81 83 if (itysol.ne.0) then 82 84 write(*,'(a)') ' init_energy> undefined solvent type !' … … 110 112 111 113 include 'INCL.H' 114 115 double precision hbc, hba, ri, ai, aei, aic, a, aj, c, rij 116 117 integer i, j, iac, ido, jac, jdo 112 118 113 119 dimension hbc(mxhbdo,mxhbac),hba(mxhbdo,mxhbac) … … 308 314 309 315 include 'INCL.H' 316 317 integer i 310 318 311 319 ! Atom types ------------------------------------------------------------ … … 517 525 & 'ile ','leu ','lys ','lys+','met ','phe ','cpro','pro ','cpru', 518 526 & 'prou','pron','pro+','ser ','thr ','trp ','tyr ','val ' / 519 527 520 528 data (onltcd(i),i=1,nrsty)/ ! One-letter codes for amino acid types 521 529 & 'A', 'R', 'R', 'N', 'D', 'D', 'C', 'C', 'Q', … … 527 535 ! coefficients for their solvation free energy (kcal/molxA**2) 528 536 529 ! Method: 530 531 ! itysol=1 : OONS --> T.Ooi, et al, 537 ! Method: 538 539 ! itysol=1 : OONS --> T.Ooi, et al, 532 540 ! Proc. Natl. Acad. Sci. USA 8 (1987) 3086-3090. 533 ! itysol=2 : JRF --> J.Vila, et al, 541 ! itysol=2 : JRF --> J.Vila, et al, 534 542 ! PROTEINS: Struct Funct Genet 10(1991) 199-218. 535 ! itysol=3 : WE92 --> L.Wesson, D.Eisenberg, 543 ! itysol=3 : WE92 --> L.Wesson, D.Eisenberg, 536 544 ! Protein Science 1 (1992) 227-235. 537 545 ! itysol=4 : SCH1 --> D.Eisenberg, et al, … … 539 547 ! itysol=5 : SCH2 --> A.H.Juffer, et al, 540 548 ! Proteine Science 4 (1995) 2499-2509. 541 ! itysol=6 : SCH3 --> L.Wesson, D.Eisenberg, 549 ! itysol=6 : SCH3 --> L.Wesson, D.Eisenberg, 542 550 ! Protein Science 1 (1992) 227-235. 543 551 ! itysol=7 : SCH4 --> C.A. Schiffer, et al, 544 552 ! Mol. Simul. 10(1993) 121-149. 545 ! itysol=8 : EM86 --> D.Eisenberg, A.D. Mclachlan, 553 ! itysol=8 : EM86 --> D.Eisenberg, A.D. Mclachlan, 546 554 ! Nature 319 (1986) 199-203. 547 555 ! itysol=9 : BM --> B. Freyberg, et al,
Note:
See TracChangeset
for help on using the changeset viewer.