newest dihed constr added

[unres4.git] / source / unres / data / energy_data.f90

diff --git a/source/unres/data/energy_data.f90 b/source/unres/data/energy_data.f90
index baf7f55..a26a0aa 100644 (file)
--- a/source/unres/data/energy_data.f90
+++ b/source/unres/data/energy_data.f90
@@ -3,7 +3,7 @@
       use names
 !-----------------------------------------------------------------------------
 ! Max. number of energy intervals
-      integer,parameter :: max_ene=21 !10
+      integer,parameter :: max_ene=49 !10
 !-----------------------------------------------------------------------------
 ! Maximum number of terms in SC bond-stretching potential
       integer,parameter :: maxbondterm=3
@@ -66,7 +66,8 @@
        wscloc,wcorr,wcorr4,wcorr5,wcorr6,wsccor,wel_loc,wturn3,wturn4,&
        wturn6,wvdwpp,wliptran,wshield,lipscale,wtube, &
        wbond_nucl,wang_nucl,wcorr_nucl,wcorr3_nucl,welpp,wtor_nucl,&
-       wtor_d_nucl,welsb,wsbloc,wvdwsb,welpsb,wvdwpp_nucl,wvdwpsb
+       wtor_d_nucl,welsb,wsbloc,wvdwsb,welpsb,wvdwpp_nucl,wvdwpsb,wcatprot,&
+       wcatcat,wscbase,wpepbase,wscpho,wpeppho
 #ifdef CLUSTER
       real(kind=8) :: scalscp
 #endif
@@ -273,9 +274,9 @@
       integer,dimension(:),allocatable :: idih_constr,idih_nconstr,itheta_constr !(maxdih_constr)
       integer :: idihconstr_start,idihconstr_end, &
        ithetaconstr_start,ithetaconstr_end
-      real(kind=8) :: ftors
+!      real(kind=8) :: ftors
       real(kind=8),dimension(:),allocatable :: drange,theta_constr0,theta_drange !(maxdih_constr)
-      real(kind=8),dimension(:),allocatable :: phi0 !(maxdih_constr)
+      real(kind=8),dimension(:),allocatable :: phi0,ftors !(maxdih_constr)
       real(kind=8),dimension(:),allocatable :: for_thet_constr !(maxdih_constr)
 
 !-----------------------------------------------------------------------------
@@ -336,8 +337,36 @@
 ! AFM
        real(kind=8) :: distafminit,forceAFMconst,velAFMconst
       integer :: afmend,afmbeg
+      real(kind=8),dimension(:,:), allocatable :: catprm
 
+         real(kind=8),dimension(:,:), allocatable ::  eps_scbase, &
+        sigma_scbase,                         &
+        sigmap1_scbase,sigmap2_scbase,        &
+        dhead_scbasei, dhead_scbasej, epshead_scbase,&
+        sig0head_scbase,  rborn_scbasei,rborn_scbasej,alphapol_scbase,epsintab_scbase,&
+        aa_scbase,bb_scbase
+         real(kind=8),dimension(:,:,:), allocatable :: alphasur_scbase, &
+        wdipdip_scbase,wqdip_scbase,chi_scbase,chipp_scbase,chis_scbase
 
+        real(kind=8),dimension(:), allocatable ::  eps_pepbase, &
+        sigma_pepbase,                         &
+        sigmap1_pepbase,sigmap2_pepbase,&
+        aa_pepbase,bb_pepbase
 
+         real(kind=8),dimension(:,:), allocatable :: alphasur_pepbase, &
+        wdipdip_pepbase,chi_pepbase,chipp_pepbase,chis_pepbase
 
+        real(kind=8),dimension(:), allocatable ::  eps_scpho, &
+        sigma_scpho,                         &
+        sigmap1_scpho,sigmap2_scpho,&
+        aa_scpho,bb_scpho,wqq_scpho,epsintab_scpho,alphapol_scpho,&
+        rborn_scphoi,rborn_scphoj,dhead_scphoi,alphi_scpho
+
+         real(kind=8),dimension(:,:), allocatable :: alphasur_scpho, &
+        chi_scpho,chipp_scpho,chis_scpho,              &
+        wqdip_scpho
+         real(kind=8) ,dimension(4) :: alphasur_peppho
+         real(kind=8) ,dimension(2) :: wqdip_peppho
+         real(kind=8) :: eps_peppho,sigma_peppho,sigmap1_peppho,sigmap2_peppho, &
+         aa_peppho,bb_peppho
       end module energy_data