Changeset bd2278d for dihedr.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
  • dihedr.f

    r2ebb8b6 rbd2278d  
    1 c **************************************************************
    2 c
    3 c This file contains the subroutines: dihedr,valang
    4 c
    5 c Copyright 2003-2005  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: dihedr,valang
     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! **************************************************************
    1111
    1212      real*8 function dihedr(i1,i2,i3,i4)
    1313
    14 c .............................................
    15 c  PURPOSE: return dihedral angle (i1,i2,i3,i4)
    16 c           [in rad.]
    17 c
    18 c  INPUT:   i1,i2,i3,i4 - indices of four atoms
    19 c
    20 c  CALLS:   none
    21 c .............................................
     14! .............................................
     15!  PURPOSE: return dihedral angle (i1,i2,i3,i4)
     16!           [in rad.]
     17!
     18!  INPUT:   i1,i2,i3,i4 - indices of four atoms
     19!
     20!  CALLS:   none
     21! .............................................
    2222
    2323      include 'INCL.H'
     
    4949        dihedr=acos(a)
    5050        if (ux1*(uy2*z2-uz2*y2)+uy1*(uz2*x2-ux2*z2)+
    51      #      uz1*(ux2*y2-uy2*x2).lt.zero) dihedr =-dihedr
     51     &      uz1*(ux2*y2-uy2*x2).lt.zero) dihedr =-dihedr
    5252        return
    5353      else
    5454        write (*,'(a,4i5)')' dihedr> Error in coordinates of atoms #: '
    55      #                     ,i1,i2,i3,i4
     55     &                     ,i1,i2,i3,i4
    5656
    5757        write (*,*) 'stored coordinates are xvals :',
    58      #       xat(i1),xat(i2),xat(i3),xat(i4)
     58     &       xat(i1),xat(i2),xat(i3),xat(i4)
    5959        write (*,*) 'yvals:', yat(i1),yat(i2),yat(i3),yat(i4)
    6060        write (*,*) 'zvals:', zat(i1),zat(i2),zat(i3),zat(i4)
     
    6464
    6565      end
    66 c ************************************
     66! ************************************
    6767      real*8 function valang(i1,i2,i3)
    6868
    69 c .........................................
    70 c  PURPOSE: return valence angle (i1,i2,i3)
    71 c           [in rad.] with 'i2' as vertex
    72 c
    73 c  INPUT:   i1,i2,i3 - indices of 3 atoms
    74 c
    75 c  CALLS:   none
    76 c .............................................
     69! .........................................
     70!  PURPOSE: return valence angle (i1,i2,i3)
     71!           [in rad.] with 'i2' as vertex
     72!
     73!  INPUT:   i1,i2,i3 - indices of 3 atoms
     74!
     75!  CALLS:   none
     76! .............................................
    7777
    7878      include 'INCL.H'
     
    101101      else
    102102        write (*,'(a,3i5)')' valang> Error in coordinates of atoms #: '
    103      #                     ,i1,i2,i3
     103     &                     ,i1,i2,i3
    104104        write (*,*) 'stored coordinates are xvals :',
    105      #       xat(i1),xat(i2),xat(i3)
     105     &       xat(i1),xat(i2),xat(i3)
    106106        write (*,*) 'yvals:', yat(i1),yat(i2),yat(i3)
    107107        write (*,*) 'zvals:', zat(i1),zat(i2),zat(i3)
Note: See TracChangeset for help on using the changeset viewer.