source/unres/data/names.F90

   1       module names
   2 !-----------------------------------------------------------------------------
   3 !-----------------------------------------------------------------------------
   4 ! Number of AA types (at present only natural AA's will be handled
   5       integer,parameter :: ntyp=27,ntyp1=ntyp+1,ntyp11=2*ntyp1+1
   6       integer,dimension(5) :: ntyp_molec=(/27,5,0,0,5/),ntyp1_molec=(/28,6,0,0,6/)
   7       integer,parameter ::maxmolec=5
   8
   9 !-----------------------------------------------------------------------------
  10 ! common.names
  11 !      common /names/
  12 !el      character(len=3),dimension(:),allocatable :: restyp    !(-ntyp1:ntyp1)
  13 !el      character(len=1),dimension(:),allocatable :: onelet    !(-ntyp1:ntyp1)
  14 !-----------------------------------------------------------------------------
  15 !-----------------------------------------------------------------------------
  16 !      block data nazwy
  17 !el      allocate(restyp(-ntyp1:ntyp1)) !(-ntyp1:ntyp1)
  18         character(len=3),dimension(-ntyp1:ntyp1,maxmolec) :: restyp = reshape(&
  19         (/'DD ','DDX','DDY','DDZ','DAU','DAI','DDB','DSM','DPR','DLY', &
  20           'DAR','DHI','DAS',&
  21        'DGL','DSG','DGN','DSN','DTH',&
  22        'DYY','DAL','DTY','DTR','DVA','DLE','DIL','DPN','MED','DCY','ZER',&
  23        'CYS','MET','PHE','ILE','LEU','VAL','TRP','TYR','ALA','GLY','THR',&
  24        'SER','GLN','ASN','GLU','ASP','HIS','ARG','LYS','PRO','SME','DBZ',&
  25        'AIB','ABU','PTR','TPO','SEP','D  ',&
  26        '   ','   ','   ','   ','   ','   ','   ','   ','   ','   ', &
  27        '   ','   ','   ',&
  28        '   ','   ','   ','   ','   ','   ','   ','   ','   ','   ',&
  29        '   ','   ','   ','   ','   ','   ',&
  30        'A  ','G  ','C  ','T  ','U  ','X  ','   ','   ','   ',&
  31        '   ','   ','   ','   ','   ','   ','   ','   ','   ',&
  32        '   ','   ','   ','   ','   ','   ','   ','   ','   ','   ', &
  33        '   ','   ','   ',&
  34        '   ','   ','   ','   ','   ','   ','   ','   ','   ','   ',&
  35        '   ','   ','   ','   ','   ','   ',&
  36        '   ','   ','   ','   ','   ','   ','   ','   ','   ',&
  37        '   ','   ','   ','   ','   ','   ','   ','   ','   ',&
  38        '   ','   ','   ','   ','   ','   ','   ','   ','   ','   ', &
  39        '   ','   ','   ',&
  40        '   ','   ','   ','   ','   ','   ','   ','   ','   ','   ',&
  41        '   ','   ','   ','   ','   ','   ',&
  42        '   ','   ','   ','   ','   ','   ','   ','   ','   ',&
  43        '   ','   ','   ','   ','   ','   ','   ','   ','   ',&
  44        '   ','   ','   ','   ','   ','   ','   ','   ','   ','   ', &
  45        '   ','   ','   ',&
  46        '   ','   ','   ','   ','   ','   ','   ','   ','   ','   ',&
  47        '   ','   ','   ','   ','   ','   ',&
  48        'NA+','MG2','K+ ','CA2','CL-','   ','   ','   ',&
  49        '   ','   ','   ','   ','   ','   ','   ','   ','   ','   ',&
  50        '   ','   ','   ','   ','   ','   ','   ','   ','   ','   ',&
  51        '   ','   ','   ','   ','   ','   ','   ','   ','   ','   ',&
  52        '   ','   ','   ','   ','   ','   ','   ','   ','   ','   ',&
  53        '   ','   ','   ','   ','   ','   ','   ','   ','   ','   ',&
  54        '   ','   ','   ','   ','   ','   ','   ','   ','   ','   '&
  55        /),(/ntyp11,maxmolec/))
  56 !el      allocate(onelet(-ntyp1:ntyp1))         !(-ntyp1:ntyp1)
  57         character(len=1),dimension(-ntyp1:ntyp1) :: onelet = reshape(&
  58         (/'z','z','z','z','z','z','z','z',&
  59         'p','k','r','h','d','e','n','q','s',&
  60         't','g','a','y','w','v','l','i','f','m','c','x',&
  61         'C','M','F','I','L','V','W','Y','A','G','T',&
  62         'S','Q','N','E','D','H','R','K','P',&
  63         'z','z','z','z','z','z','z','X'/),(/ntyp11/))
  64 !        character(len=1),dimension(ntyp1_nucl) :: restyp_nucl = &
  65 !        (/'A','G','C','T','U','X'/)
  66
  67 !-----------------------------------------------------------------------------
  68 !-----------------------------------------------------------------------------
  69 ! Number of energy components
  70       integer,parameter :: n_ene=49
  71       integer :: n_ene2=2*n_ene
  72 !-----------------------------------------------------------------------------
  73 ! common.names
  74 !#ifndef WHAM_RUN
  75 !      common /namterm/
  76 !      character(len=10),dimension(n_ene) :: ename = &
  77 !      (/"EVDW SC-SC","EVDW2 SC-p","EES p-p   ","ECORR4    ","ECORR5    ",&
  78 !        "ECORR6    ","EELLO     ","ETURN3    ","ETURN4    ","ETURN6    ",&
  79 !        "EBE bend  ","ESC SCloc ","ETORS     ","ETORSD    ","EHPB      ","EVDWPP    ",&
  80 !        "ESTR      ","EVDW2_14  ","UCONST    ","          ","ESCCOR    "/)
  81 !      character(len=10),dimension(n_ene) :: wname = &
  82 !      (/"WSC       ","WSCP      ","WELEC     ","WCORR     ","WCORR5    ","WCORR6    ","WEL_LOC   ",&
  83 !        "WTURN3    ","WTURN4    ","WTURN6    ","WANG      ","WSCLOC    ","WTOR      ","WTORD     ",&
  84 !        "WSTRAIN   ","WVDWPP    ","WBOND     ","SCAL14    ","          ","          ","WSCCOR    "/)
  85 !      integer :: nprint_ene = 20
  86 !      integer,dimension(n_ene) :: print_order = &
  87 !         (/1,2,3,11,12,13,14,4,5,6,7,8,9,10,19,18,15,17,16,21,0/)
  88 !#else
  89       character(len=10),dimension(n_ene) :: ename = &
  90       (/"EVDW SC-SC","EVDW2 SC-p","EES p-p   ","ECORR4    ","ECORR5    ",&
  91         "ECORR6    ","EELLO     ","ETURN3    ","ETURN4    ","ETURN6    ",&
  92         "EBE bend  ","ESC SCloc ","ETORS     ","ETORSD    ","EHPB      ","EVDWPP    ",&
  93         "EVDW2_14  ","ESTR      ","ESCCOR    ","EDIHC     ","EVDW_T    ",&
  94         "ELT       ","          ","          ","ETUBE     ",&
  95         "EVDWPP    ","EESPP     ","EVDWPSB   ","EESPSB    ","EVDWSB    ",&
  96         "EESSB     ","ESTR      ","EBE       ","ESBLOC    ","ETORS     ",&
  97         "ETORSD    ","ECORR     ","ECORR3    ","NULL      ","NULL      ",&
  98         "ECATPROT  ","ECATCAT   ","NULL      ","NULL      ","NULL      ",&
  99         "ESCBASE   ","EPEPBASE  ","ESCPHO    ","EPEPPHO   "/)
 100
 101       character(len=10),dimension(n_ene) :: wname = &
 102       (/"WSC       ","WSCP      ","WELEC     ","WCORR     ","WCORR5    ","WCORR6    ","WEL_LOC   ",&
 103         "WTURN3    ","WTURN4    ","WTURN6    ","WANG      ","WSCLOC    ","WTOR      ","WTORD     ",&
 104         "WHPB      ","WVDWPP    ","WSCP14    ","WBOND     ","WSCCOR    ","WDIHC     ","WSC       ",&
 105         "WLT       ","          ","          ","WTUBE     " ,&
 106         "WVDWPP    ","WELPP     ","WVDWPSB   ","WELPSB    ","WVDWSB    ",&
 107         "WELSB     ","WBOND     ","WANG      ","WSBLOC    ","WTOR      ",&
 108         "WTORD     ","WCORR     ","WCORR3    ","WNULL     ","WNULL     ",&
 109         "WCATPROT  ","WCATCAT   ","WNULL     ","WNULL     ","WNULL     ",&
 110         "WSCBASE   ","WPEPBASE  ","WSCPHO    ","WPEPPHO   "/)
 111
 112       integer :: nprint_ene = 21
 113       integer,dimension(n_ene) :: print_order = &
 114          (/1,2,3,18,11,12,13,14,4,5,6,7,8,9,10,19,16,15,17,20,21,22,23,24,25,&
 115          26,27,28,29,30,31,32,33,34,35,36,37,38,39,40,41,42,43,44,45,46,47,&
 116          48,49/)
 117
 118       character(len=1), dimension(2) :: sugartyp = (/'D',' '/)
 119 !#endif
 120 !-----------------------------------------------------------------------------
 121 !-----------------------------------------------------------------------------
 122       end module names