- Timestamp:
- 11/19/09 11:29:41 (14 years ago)
- Branches:
- master
- Children:
- 38d77eb
- Parents:
- 6650a56
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
rmsdfun.f
r6650a56 r32289cd 5 5 ! 6 6 ! Copyright 2003-2005 Frank Eisenmenger, U.H.E. Hansmann, 7 ! Shura Hayryan, Chin-Ku 7 ! Shura Hayryan, Chin-Ku 8 8 ! Copyright 2007 Frank Eisenmenger, U.H.E. Hansmann, 9 9 ! Jan H. Meinke, Sandipan Mohanty … … 26 26 ! 27 27 ! Input 28 double precision xrf, yrf, zrf, rm, av1, av2, rssd 29 30 integer nml, ir1, ir2, ixat, isl 31 28 32 dimension ixat(mxat),xrf(mxatp),yrf(mxatp),zrf(mxatp) 29 33 ! Local … … 57 61 ! 58 62 ! NB uncomment last lines in 'fitmol' to return coordinates 59 ! in 'x2' after fitting the ref. str. onto SMMP structure 63 ! in 'x2' after fitting the ref. str. onto SMMP structure 60 64 ! ---------------------------------------------------------------- 61 65 62 66 include 'INCL.H' 63 67 include 'INCP.H' 68 69 double precision x1, x2, xrf, yrf, zrf, rm, av1, av2, rmsd 70 71 integer nml, nr, na, n, im, ir, ia, ir1, ir2, ix, ixat, isl 72 64 73 65 74 !-------------------------------------------------------- input … … 147 156 ! ....................................................... 148 157 !f2py intent(out) rmsd 149 158 150 159 include 'INCL.H' 151 160 ! implicit real*8 (a-h,o-z) 152 161 ! implicit integer*4 (i-n) 153 162 154 163 ! ------------------------------------------- input/output 164 double precision dn, a1, a2, x1, x2, q, dm, dp, dxm, dym, dzm, dxp 165 double precision dyp, dzp, e, v, em, rmsd, rm 166 167 integer n, i, j, ndim4, im 168 155 169 dimension x1(3,mxat),x2(3,mxat) 156 170 ! -------------------------------------------------- local … … 270 284 !f2py intent(out) d 271 285 !f2py intent(out) v 286 integer nmax 287 272 288 parameter (NMAX=500) 273 289 274 290 integer n,nrot,i,ip,iq,j 275 291 … … 398 414 ! 399 415 !------------------------------------------------------------------------------ 400 ! Reads in pdb-file 'string' into INCP.H and initalizes 416 ! Reads in pdb-file 'string' into INCP.H and initalizes 401 417 ! the files that 'rmdsopt' needs to calculate the rmsd 402 418 ! of a configuration with the pdb-configuration … … 409 425 include 'INCP.H' 410 426 427 integer i, nml, ier 428 411 429 character string*(*) 412 430 413 431 if(string.eq.'smmp') then 414 432 ! … … 428 446 ! 429 447 call pdbread(string,ier) 430 if(ier.ne.0) stop 448 if(ier.ne.0) stop 431 449 call atixpdb(nml) 432 450 !
Note:
See TracChangeset
for help on using the changeset viewer.