source/cluster/wham/src/COMMON.DFA

   1 C =======
   2 C COMMON.DFA
   3 C =======
   4 C 2010/12/20 By Juyong Lee
   5 C
   6 c parameter
   7 C [ 8 * ( Nres - 8 ) ] distance restraints
   8 C [ 2 * ( Nres - 8 ) ] angle restraints
   9 C [ Nres ]             neighbor restraints
  10 C Total : ~ 11 * Nres restraints
  11 C
  12 C
  13       INTEGER IDFAMAX,IDFAMX2,IDFACMD,IDMAXMIN, MAXN
  14       PARAMETER(IDFAMAX=4000,IDFAMX2=1000,IDFACMD=500,IDMAXMIN=500)
  15       PARAMETER(MAXN=4)
  16       real*8 wwdist,wwangle,wwnei
  17       parameter(wwdist=1.0d0,wwangle=1.0d0,wwnei=1.0d0)
  18
  19 C IDFAMAX  - maximum number of DFA restraint including distance, angle and
  20 C            number of neighbors ( Max of assign statement )
  21 C IDFAMX2  - maximum number of atoms which are targets of restraints
  22 C IDFACMD  - maximum number of 'DFA' command call
  23 C IDMAXMIN - Maximum number of minima of dist, angle and neighbor info. from fragments
  24 C MAXN     - Maximum Number of shell, currently 4
  25 C MAXRES   - Maximum number of CAs
  26
  27 CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCc
  28 C INTEGER
  29 C DFANUM  - Number of ALL DFA restrants
  30 c IDFA[DIS, PHI, THE, NEI] - NUMBER of restraints
  31 c IDISNUM - number of minima for a distance restraint
  32 c IPHINUM - number of minima for a phi angle restraint
  33 c ITHENUM - number of minima for a theta angle restraint
  34 c INEINUM - number of minima for a number of neighbors restraint
  35
  36 c IDISLIS - atom number of two atoms for distance restraint
  37 c IPHILIS - atom numbers of four atoms for angle restraint
  38 c ITHELIS - atom numbers of four atoms for angle restraint
  39 c INEILIS - atom number of center of neighbor calculation
  40 c JNEILIS - atom number of target of neighboring calculation
  41 c JNEINUM - number of target atoms of neighboring term
  42 C KSHELL  - SHELL number
  43
  44 C ishiftca - index shift for CA atoms in UNRES (1 if the 1st aa != GLY)
  45 C ilastca  - index of the last CA atom in UNRES (nres-1 if last aa != GLY)
  46
  47 C     old only for CHARMM
  48 C STOAGDF - Store assign information ( How many assign within one command )
  49 C NMAP    - mapping between dfanum and ndis, nphi, nthe, nnei
  50
  51       INTEGER IDFADIS,IDFAPHI,IDFATHE,IDFANEI,
  52      &        IDISLIS,IPHILIS,ITHELIS,INEILIS,
  53      &        IDISNUM,IPHINUM,ITHENUM,INEINUM,
  54      &        FNEI,DFACMD, DFANUM,
  55      &        NCA,ICAIDX,
  56      &        STOAGDF, NMAP, IDFACAT, KDISNUM, KSHELL
  57      &        ishiftca,ilastca
  58       COMMON /IDFA/ DFACMD, DFANUM,
  59      &              IDFADIS, IDFAPHI, IDFANEI, IDFATHE,
  60      &              IDISNUM(IDFAMAX), IPHINUM(IDFAMAX),
  61      &              ITHENUM(IDFAMAX), INEINUM(IDFAMAX),
  62      &              FNEI(IDFAMAX,IDMAXMIN), IDISLIS(2,IDFAMAX),
  63      &              IPHILIS(5,IDFAMAX), ITHELIS(5,IDFAMAX),
  64      &              INEILIS(IDFAMAX),
  65      &               KSHELL(IDFAMAX),
  66      &              IDFACAT(IDFACMD),
  67      &              KDISNUM(IDFAMAX),
  68      &              NCA, ICAIDX(MAXRES)
  69       COMMON /IDFA2/ ishiftca,ilastca
  70
  71 CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC
  72 C
  73 C REAL VARIABLES
  74 C
  75 c SCC[DIST, PHI, THE] - weight of each calculations
  76 c FDIST  - distance minima
  77 C FPHI   - phi minima
  78 c FTHE   - theta minima
  79 C DFAEXP  : calculate expential function in advance
  80 C
  81       REAL*8 SCCDIST, SCCPHI, SCCTHE, SCCNEI, FDIST, FPHI1, FPHI2,
  82      &       FTHE1, FTHE2,
  83      &       DIS_INC, PHI_INC, THE_INC, NEI_INC, BETA_INC,
  84      &       WSHET, EDFABET,
  85      &       CK, SCK, S1, S2
  86 c    &       ,DFAEXP
  87
  88       COMMON /RDFA/ SCCDIST(IDFAMAX,IDMAXMIN),FDIST(IDFAMAX,IDMAXMIN),
  89      &             SCCPHI(IDFAMAX,IDMAXMIN), SCCTHE(IDFAMAX,IDMAXMIN),
  90      &             SCCNEI(IDFAMAX,IDMAXMIN),
  91      &             FPHI1(IDFAMAX,IDMAXMIN), FPHI2(IDFAMAX,IDMAXMIN),
  92      &             FTHE1(IDFAMAX,IDMAXMIN), FTHE2(IDFAMAX,IDMAXMIN),
  93      &             DIS_INC, PHI_INC, THE_INC, NEI_INC, BETA_INC,
  94      &             WSHET(MAXRES,MAXRES), EDFABET,
  95      &             CK(4),SCK(4),S1(4),S2(4)
  96 c    &             ,DFAEXP(15001),
  97
  98       DATA CK/1.0D0,1.58740105197D0,2.08008382305D0,2.51984209979D0/
  99       DATA SCK/1.0D0,1.25992104989D0,1.44224957031D0,1.58740105197D0/
 100       DATA S1/3.75D0,5.75D0,7.75D0,9.75D0/
 101       DATA S2/4.25D0,6.25D0,8.25D0,10.25D0/