Changeset bd2278d for enyflx.f


Ignore:
Timestamp:
09/05/08 11:49:42 (16 years ago)
Author:
baerbaer <baerbaer@…>
Branches:
master
Children:
fafe4d6
Parents:
2ebb8b6
Message:

Reformatting comments and continuation marks.

Fortran 90 and higher use ! to mark comments no matter where they are in the
code. The only valid continuation marker is &.
I also added the SMMP.kdevelop.filelist to the repository to make it easier
to use kdevelop.

git-svn-id: svn+ssh://svn.berlios.de/svnroot/repos/smmp/trunk@12 26dc1dd8-5c4e-0410-9ffe-d298b4865968

File:
1 edited

Legend:

Unmodified
Added
Removed
  • enyflx.f

    r2ebb8b6 rbd2278d  
    1 c **************************************************************
    2 c
    3 c This file contains the subroutines: enyflx
    4 c
    5 c Copyright 2003       Frank Eisenmenger, U.H.E. Hansmann,
    6 c                      Shura Hayryan, Chin-Ku
    7 c Copyright 2007       Frank Eisenmenger, U.H.E. Hansmann,
    8 c                      Jan H. Meinke, Sandipan Mohanty
    9 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! **************************************************************
    1111
    1212
    1313      real*8 function enyflx(nml)
    1414
    15 c .......................................................................
    16 c
    17 c  PURPOSE: Calculate internal energy of molecule 'nml' with FLEX dataset
    18 c
    19 c  CALLS: none
    20 c
    21 c .......................................................................
     15! .......................................................................
     16!
     17!  PURPOSE: Calculate internal energy of molecule 'nml' with FLEX dataset
     18!
     19!  CALLS: none
     20!
     21! .......................................................................
    2222
    2323      include 'INCL.H'
     
    2626      if (ntlvr.eq.0) then
    2727        write (*,'(a,i4)')
    28      #           ' enyflx> No variables defined in molecule #',nml
     28     &           ' enyflx> No variables defined in molecule #',nml
    2929        return
    3030      endif
     
    9292                rij=sqrt(rij2)
    9393                if(epsd) then
    94 c --------------------------------- distance dependent dielectric constant
     94! --------------------------------- distance dependent dielectric constant
    9595                sr=slp_f*rij
    9696                ep=plt-(sr*sr+2.0*sr+2.0)*(plt-1.0)*exp(-sr)/2.0
     
    119119
    120120                  cth=(xij*px+yij*py+zij*pz)/(rij*
    121      #                 sqrt(px*px+py*py+pz*pz))
     121     &                 sqrt(px*px+py*py+pz*pz))
    122122
    123123                  if (cth.gt.0.0) then
    124124                    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 )
    127127                  else                             ! No Hydrogen Bond
    128128                    eyvw=eyvw + evw
     
    153153              rij=sqrt(rij2)
    154154              if(epsd) then
    155 c --------------------------------- distance dependent dielectric constant
     155! --------------------------------- distance dependent dielectric constant
    156156              sr=slp_f*rij
    157157              ep=plt-(sr*sr+2.0*sr+2.0)*(plt-1.)*exp(-sr)/2.0
     
    180180
    181181                cth=(xij*px+yij*py+zij*pz)/(rij*
    182      #               sqrt(px*px+py*py+pz*pz))
     182     &               sqrt(px*px+py*py+pz*pz))
    183183
    184184                if (cth.gt.0.0) then
    185185                  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 )
    188188                else                             ! No Hydrogen Bond
    189189                  eyvw=eyvw + evw
Note: See TracChangeset for help on using the changeset viewer.