|
Last change
on this file since bd2278d was bd2278d, checked in by baerbaer <baerbaer@…>, 16 years ago |
|
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
|
-
Property mode
set to
100644
|
|
File size:
878 bytes
|
| Line | |
|---|
| 1 | ! This file contains enyreg
|
|---|
| 2 | !
|
|---|
| 3 | ! Copyright 2007 Frank Eisenmenger, U.H.E. Hansmann,
|
|---|
| 4 | ! Jan H. Meinke, Sandipan Mohanty
|
|---|
| 5 | ! *******************************
|
|---|
| 6 | real*8 function enyreg(nml)
|
|---|
| 7 |
|
|---|
| 8 | ! ----------------------------------------------------
|
|---|
| 9 | !
|
|---|
| 10 | ! PURPOSE: sum( ( R_i - R^ref_j )**2 )
|
|---|
| 11 | !
|
|---|
| 12 | ! with: R_i - atom position i in SMMP structure
|
|---|
| 13 | ! R^ref_j - corresponding atom j in PDB str.
|
|---|
| 14 | !
|
|---|
| 15 | ! CALLS: none
|
|---|
| 16 | !
|
|---|
| 17 | ! ----------------------------------------------------
|
|---|
| 18 |
|
|---|
| 19 | include 'INCL.H'
|
|---|
| 20 | include 'INCP.H'
|
|---|
| 21 |
|
|---|
| 22 |
|
|---|
| 23 | eny = 0.d0
|
|---|
| 24 |
|
|---|
| 25 | do i=iatrs1(irsml1(nml)),iatrs2(irsml2(nml))
|
|---|
| 26 |
|
|---|
| 27 | j=ixatp(i)
|
|---|
| 28 | if (j.gt.0) then ! corresp. atom in ref. structure
|
|---|
| 29 |
|
|---|
| 30 | eny = eny + (xat(i)-xatp(j))**2+(yat(i)-yatp(j))**2+
|
|---|
| 31 | & (zat(i)-zatp(j))**2
|
|---|
| 32 | endif
|
|---|
| 33 | enddo ! atoms
|
|---|
| 34 |
|
|---|
| 35 | enyreg = eny
|
|---|
| 36 |
|
|---|
| 37 | return
|
|---|
| 38 | end
|
|---|
| 39 |
|
|---|
Note:
See
TracBrowser
for help on using the repository browser.
