source/unres/data/MD_data.F90

   1       module MD_data
   2 !-----------------------------------------------------------------------------
   3 #ifndef LANG0
   4 ! commom.langevin
   5 !      common /langforc/
   6       real(kind=8),dimension(:,:),allocatable :: friction,stochforc !(3,0:MAXRES2)
   7       real(kind=8),dimension(:),allocatable :: fric_work,stoch_work,&
   8        fricgam !(MAXRES6)
   9       real(kind=8),dimension(:,:),allocatable :: fricmat,fricvec,&
  10        pfric_mat,vfric_mat,afric_mat,prand_mat,vrand_mat1,&
  11        vrand_mat2 !(MAXRES2,MAXRES2)
  12       real(kind=8),dimension(:,:,:),allocatable :: pfric0_mat,&
  13        afric0_mat,vfric0_mat,prand0_mat,vrand0_mat1,vrand0_mat2 !(MAXRES2,MAXRES2,0:maxflag_stoch)
  14       logical,dimension(:),allocatable :: flag_stoch !(0:maxflag_stoch)
  15 !      common /langmat/
  16       real(kind=8),dimension(:,:),allocatable :: mt1,mt2,mt3 !(maxres2,maxres2)
  17 !-----------------------------------------------------------------------------
  18 #else
  19 ! commom.langevin.lang0
  20 !      common /langforc/
  21       real(kind=8),dimension(:,:),allocatable :: friction,stochforc !(3,0:MAXRES2)
  22       real(kind=8),dimension(:,:),allocatable :: fricmat,fricvec !(MAXRES2,MAXRES2)
  23       real(kind=8),dimension(:),allocatable :: fric_work,stoch_work,&
  24        fricgam !(MAXRES6)
  25       logical,dimension(:),allocatable :: flag_stoch !(0:maxflag_stoch)
  26       real(kind=8) :: vrand_mat1,vrand_mat2,prand_mat,vfric_mat,&
  27        afric_mat,pfric_mat,pfric0_mat,afric0_mat,vfric0_mat,prand0_mat,&
  28        vrand0_mat1,vrand0_mat2
  29 !      common /langmat/
  30       integer :: mt1,mt2,mt3
  31 #endif
  32 !-----------------------------------------------------------------------------
  33 ! commom.hairpin in CSA module
  34 !-----------------------------------------------------------------------------
  35 ! common.mce in MCM_MD module
  36 !-----------------------------------------------------------------------------
  37 ! common.MD
  38 !      common /mdgrad/ in module.energy
  39 !      common /back_constr/ in module.energy
  40 !      common /qmeas/ in module.energy
  41 !      common /mdpar/
  42       real(kind=8) :: v_ini,d_time,d_time0,scal_fric,&
  43        t_bath,tau_bath,dvmax,damax
  44       integer :: n_timestep,ntime_split,ntime_split0,maxtime_split,&
  45        ntwx,ntwe
  46       logical :: mdpdb,large,print_compon,tbf,rest
  47 !      common /MDcalc/
  48       real(kind=8) :: totT,totE,potE,EK,amax,edriftmax,kinetic_T
  49       real(kind=8),dimension(:),allocatable :: potEcomp !(0:n_ene+4)
  50 !      common /lagrange/
  51       real(kind=8),dimension(:,:),allocatable :: d_t,d_a,d_t_old !(3,0:MAXRES2)
  52       real(kind=8),dimension(:),allocatable :: d_a_work !(6*MAXRES)
  53       real(kind=8),dimension(:,:),allocatable :: Gmat,Ginv,A,&
  54        Gsqrp,Gsqrm,Gvec !(maxres2,maxres2)
  55       real(kind=8),dimension(:),allocatable :: Geigen !(maxres2)
  56       real(kind=8),dimension(:),allocatable ::vtot !(maxres2)
  57       logical :: reset_moment,reset_vel,rattle,RESPA
  58       integer :: dimen,dimen1,dimen3
  59       integer :: lang,count_reset_moment,count_reset_vel
  60 !      common /inertia/
  61       real(kind=8) :: IP,mp
  62       real(kind=8),dimension(:),allocatable :: ISC,msc !(ntyp+1)
  63 !      common /langevin/
  64       real(kind=8) :: rwat,etawat,stdfp,pstok,gamp!,Rb
  65       real(kind=8) :: cPoise=2.9361d0, Rb=0.001986d0
  66       real(kind=8),dimension(:),allocatable :: gamsc !(ntyp1)
  67       real(kind=8),dimension(:),allocatable :: stdfsc !(ntyp)
  68
  69       real(kind=8),dimension(:),allocatable :: restok !(ntyp+1)
  70       logical :: surfarea
  71       integer :: reset_fricmat
  72 !      common /mdpmpi/
  73       integer :: igmult_start,igmult_end,my_ng_count,myginv_ng_count
  74       integer,dimension(:),allocatable :: ng_start,ng_counts,&
  75        nginv_counts !(0:MaxProcs-1)
  76       integer,dimension(:),allocatable :: nginv_start !(0:MaxProcs)
  77 !-----------------------------------------------------------------------------
  78 ! common.muca
  79 !      common /double_muca/
  80       real(kind=8) :: elow,ehigh,factor,hbin,factor_min
  81       real(kind=8),dimension(:),allocatable :: emuca,nemuca,&
  82        nemuca2,hist !(4*maxres)
  83 !      common /integer_muca/
  84       integer :: nmuca,imtime,muca_smooth
  85 !      common /mucarem/
  86       real(kind=8),dimension(:),allocatable :: elowi,ehighi !(maxprocs)
  87 !-----------------------------------------------------------------------------
  88 ! Maximum number of timesteps for which stochastic MD matrices can be stored
  89       integer,parameter :: maxflag_stoch=0
  90 !-----------------------------------------------------------------------------
  91 !      common /przechowalnia/ subroutines: setup_MD_matrices
  92       real(kind=8),dimension(:,:),allocatable :: Gcopy !(maxres2,maxres2), maxres2=2*maxres
  93 !-----------------------------------------------------------------------------
  94 !      common /przechowalnia/ subroutines: setup_fricmat,setup_MD_matrices
  95       real(kind=8),dimension(:),allocatable :: Ghalf
  96 !-----------------------------------------------------------------------------
  97 !      COMMON /BANII/ D
  98       real(kind=8),DIMENSION(:),allocatable :: D_ban !(MAXRES6) maxres6=6*maxres
  99 !-----------------------------------------------------------------------------
 100       end module MD_data