X-Git-Url: http://mmka.chem.univ.gda.pl/gitweb/?a=blobdiff_plain;f=source%2Funres%2Fsrc_MD-M%2FCOMMON.MD;h=131ea7cd8461fc69af96364a81765c134ebc1ab4;hb=f4a7e366d8b9513779319977c2954e72d85293af;hp=7b437fe44a33e51ca8169fef6b23085389b5142a;hpb=a45ce81aae2c857c30b821e39e45a7a08356de4d;p=unres.git

diff --git a/source/unres/src_MD-M/COMMON.MD b/source/unres/src_MD-M/COMMON.MD
index 7b437fe..131ea7c 100644
--- a/source/unres/src_MD-M/COMMON.MD
+++ b/source/unres/src_MD-M/COMMON.MD
@@ -13,6 +13,22 @@
      & d_a(3,0:MAXRES2),d_a_work(6*MAXRES),kinetic_force(MAXRES6),
      & Gsqrp(MAXRES2,MAXRES2),Gsqrm(MAXRES2,MAXRES2),
      & vtot(MAXRES2),Gvec(maxres2,maxres2),Geigen(maxres2)
+
+       real*8 odl(max_template,maxdim),sigma_odl(max_template,maxdim),
+     &    dih(max_template,maxres),sigma_dih(max_template,maxres),
+     &    sigma_odlir(max_template,maxdim)
+c
+c    Specification of new variables used in  subroutine e_modeller
+c    modified by FP (Nov.,2014)
+       real*8 xxtpl(max_template,maxres),yytpl(max_template,maxres),
+     &        zztpl(max_template,maxres),thetatpl(max_template,maxres),
+     &        sigma_theta(max_template,maxres),
+     &        sigma_d(max_template,maxres)
+c
+
+       integer ires_homo(maxdim),
+     & jres_homo(maxdim),idomain(max_template,maxres)
+
        double precision v_ini,d_time,d_time0,t_bath,tau_bath,
      & EK,potE,potEcomp(0:n_ene+4),totE,totT,amax,kinetic_T,dvmax,damax,
      & edriftmax,
@@ -28,21 +44,33 @@
       integer n_timestep,ntwx,ntwe,lang,count_reset_moment,
      & count_reset_vel,reset_fricmat,nfrag,npair,nfrag_back,
      & ifrag_back(3,maxfrag_back,maxprocs/20),ntime_split,ntime_split0,
-     & maxtime_split
+     & maxtime_split,lim_odl,lim_dih,link_start_homo,link_end_homo,
+     & idihconstr_start_homo,idihconstr_end_homo
       logical large,print_compon,tbf,rest,reset_moment,reset_vel,
-     & surfarea,rattle,usampl,mdpdb,RESPA
+     & surfarea,rattle,usampl,mdpdb,RESPA,preminim,
+     & l_homo(max_template,maxdim)
       integer igmult_start,igmult_end,my_ng_count,ng_start,ng_counts,
      & nginv_start,nginv_counts,myginv_ng_count
       common /back_constr/ uconst_back,utheta,ugamma,uscdiff,
      & dutheta,dugamma,duscdiff,duscdiffx,
      & wfrag_back,nfrag_back,ifrag_back
+       common /homrestr/ odl,dih,sigma_dih,sigma_odl,
+     & lim_odl,lim_dih,ires_homo,jres_homo,link_start_homo,
+     & link_end_homo,idihconstr_start_homo,idihconstr_end_homo,
+     & idomain,l_homo
+c
+c    FP (30/10/2014,04/03/2015)
+c      
+       common /homrestr_double/
+     & xxtpl,yytpl,zztpl,thetatpl,sigma_theta,sigma_d,sigma_odlir
+c
       common /qmeas/ qfrag,qpair,qinfrag,qinpair,wfrag,wpair,eq_time,
      & Ucdfrag,Ucdpair,dUdconst,dUdxconst,dqwol,dxqwol,Uconst,
      & iset,mset,nset,usampl,ifrag,ipair,npair,nfrag
       common /mdpar/ v_ini,d_time,d_time0,scal_fric,
      & t_bath,tau_bath,dvmax,damax,n_timestep,mdpdb,
      & ntime_split,ntime_split0,maxtime_split,
-     & ntwx,ntwe,large,print_compon,tbf,rest
+     & ntwx,ntwe,large,print_compon,tbf,rest,preminim
       common /MDcalc/ totT,totE,potE,potEcomp,EK,amax,edriftmax,
      & kinetic_T
       common /lagrange/ d_t,d_t_old,d_t_new,d_t_work,
@@ -54,7 +82,7 @@
      & rattle,RESPA
       common /inertia/ IP,ISC,MP,MSC
       double precision scal_fric,rwat,etawat,gamp,
-     & gamsc(ntyp),stdfp,stdfsc(ntyp),stdforcp(MAXRES),
+     & gamsc(ntyp1),stdfp,stdfsc(ntyp),stdforcp(MAXRES),
      & stdforcsc(MAXRES),pstok,restok(ntyp+1),cPoise,Rb
       common /langevin/ pstok,restok,gamp,gamsc,
      & stdfp,stdfsc,stdforcp,stdforcsc,rwat,etawat,cPoise,Rb,surfarea,