- Timestamp:
- 09/05/08 11:49:42 (16 years ago)
- Branches:
- master
- Children:
- fafe4d6
- Parents:
- 2ebb8b6
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
enyflx.f
r2ebb8b6 rbd2278d 1 c**************************************************************2 c 3 cThis file contains the subroutines: enyflx4 c 5 cCopyright 2003 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 c**************************************************************1 ! ************************************************************** 2 ! 3 ! This file contains the subroutines: enyflx 4 ! 5 ! Copyright 2003 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 ! ************************************************************** 11 11 12 12 13 13 real*8 function enyflx(nml) 14 14 15 c.......................................................................16 c 17 cPURPOSE: Calculate internal energy of molecule 'nml' with FLEX dataset18 c 19 cCALLS: none20 c 21 c.......................................................................15 ! ....................................................................... 16 ! 17 ! PURPOSE: Calculate internal energy of molecule 'nml' with FLEX dataset 18 ! 19 ! CALLS: none 20 ! 21 ! ....................................................................... 22 22 23 23 include 'INCL.H' … … 26 26 if (ntlvr.eq.0) then 27 27 write (*,'(a,i4)') 28 #' enyflx> No variables defined in molecule #',nml28 & ' enyflx> No variables defined in molecule #',nml 29 29 return 30 30 endif … … 92 92 rij=sqrt(rij2) 93 93 if(epsd) then 94 c--------------------------------- distance dependent dielectric constant94 ! --------------------------------- distance dependent dielectric constant 95 95 sr=slp_f*rij 96 96 ep=plt-(sr*sr+2.0*sr+2.0)*(plt-1.0)*exp(-sr)/2.0 … … 119 119 120 120 cth=(xij*px+yij*py+zij*pz)/(rij* 121 #sqrt(px*px+py*py+pz*pz))121 & sqrt(px*px+py*py+pz*pz)) 122 122 123 123 if (cth.gt.0.0) then 124 124 eyhb=eyhb+ evw + cth*( 125 #(ahb(ity,jty)-aij(ity,jty))/rij12-126 #(chb(ity,jty)-cij(ity,jty))/rij6 )125 & (ahb(ity,jty)-aij(ity,jty))/rij12- 126 & (chb(ity,jty)-cij(ity,jty))/rij6 ) 127 127 else ! No Hydrogen Bond 128 128 eyvw=eyvw + evw … … 153 153 rij=sqrt(rij2) 154 154 if(epsd) then 155 c--------------------------------- distance dependent dielectric constant155 ! --------------------------------- distance dependent dielectric constant 156 156 sr=slp_f*rij 157 157 ep=plt-(sr*sr+2.0*sr+2.0)*(plt-1.)*exp(-sr)/2.0 … … 180 180 181 181 cth=(xij*px+yij*py+zij*pz)/(rij* 182 #sqrt(px*px+py*py+pz*pz))182 & sqrt(px*px+py*py+pz*pz)) 183 183 184 184 if (cth.gt.0.0) then 185 185 eyhb=eyhb+ evw + cth*( 186 #(ahb(ity,jty)-a14(ity,jty))/rij12-187 #(chb(ity,jty)-cij(ity,jty))/rij6 )186 & (ahb(ity,jty)-a14(ity,jty))/rij12- 187 & (chb(ity,jty)-cij(ity,jty))/rij6 ) 188 188 else ! No Hydrogen Bond 189 189 eyvw=eyvw + evw
Note:
See TracChangeset
for help on using the changeset viewer.