- Timestamp:
- 09/05/08 11:49:42 (16 years ago)
- Branches:
- master
- Children:
- fafe4d6
- Parents:
- 2ebb8b6
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
rgyr.f
r2ebb8b6 rbd2278d 1 c**************************************************************2 c 3 cThis file contains the subroutines: rgyr4 c 5 cCopyright 2003-2005 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: rgyr 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 ! ************************************************************** 11 11 12 12 13 13 subroutine rgyr(nml, rgy, ee) 14 14 15 CCALCULATES THE RADIUS-OF-GYRATION AND THE END-TO-END DISTANCE16 CFOR A GIVEN PROTEIN CONFORMATION17 CIf nml == 0, calculate the radius of gyration for all molecules18 C 19 Crgy = radius-of-gyration20 Cee = end-to-end distance21 C 22 CREQUIREMENTS: c_alfa has to be called BEFORE call of this subroutine23 C 24 CCALLS: NONE25 C 15 ! CALCULATES THE RADIUS-OF-GYRATION AND THE END-TO-END DISTANCE 16 ! FOR A GIVEN PROTEIN CONFORMATION 17 ! If nml == 0, calculate the radius of gyration for all molecules 18 ! 19 ! rgy = radius-of-gyration 20 ! ee = end-to-end distance 21 ! 22 ! REQUIREMENTS: c_alfa has to be called BEFORE call of this subroutine 23 ! 24 ! CALLS: NONE 25 ! 26 26 include 'INCL.H' 27 cf2py intent(in) nml28 cf2py intent(out) rgy29 cf2py intent(out) ee27 !f2py intent(in) nml 28 !f2py intent(out) rgy 29 !f2py intent(out) ee 30 30 integer typ 31 31 if (nml.eq.0) then … … 42 42 if (nat.le.0) then 43 43 write (*,'(a,i4)') 44 #' rgyr> No atoms found for molecule #',nml44 & ' rgyr> No atoms found for molecule #',nml 45 45 return 46 46 endif … … 128 128 129 129 ee = sqrt((xat(i2)-xat(i1))**2+(yat(i2)-yat(i1))**2 130 #+(zat(i2)-zat(i1))**2)130 & +(zat(i2)-zat(i1))**2) 131 131 132 132 return
Note:
See TracChangeset
for help on using the changeset viewer.