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,wcatprot,&
- wcatcat,wscbase,wpepbase,wscpho,wpeppho,wdihc,wcatnucl,wcat_tran,wcat_ang
+ wcatcat,wscbase,wpepbase,wscpho,wpeppho,wdihc,wcatnucl,wcat_tran,wcat_ang,&
+ wlip_prot,wmartini
#ifdef CLUSTER
real(kind=8) :: scalscp
#endif
real(kind=8),dimension(:,:,:),allocatable :: alphasurcat,&
alphisocat,wqdipcat,dtailcat,wstatecat
real(kind=8),dimension(:,:,:,:),allocatable :: dheadcat
+
+
+ real(kind=8),dimension(:,:),allocatable :: alphapolmart, alphapolmart2, &
+ epsheadmart,sig0headmart,sigiso1mart,sigiso2mart,rborn1mart,rborn2mart,&
+ sigmap1mart,sigmap2mart,chis1mart,chis2mart,wquadmart,chipp1mart,chipp2mart,&
+ epsintabmart,debaykapmart,chi1mart,chi2mart,sigmamart, nstatemart, epsmart,&
+ aa_aq_mart,bb_aq_mart
+
+ real(kind=8),dimension(:,:,:),allocatable :: alphasurmart,&
+ alphisomart,wqdipmart,dtailmart,wstatemart
+ real(kind=8),dimension(:,:,:,:),allocatable :: dheadmart
+
+
integer,dimension(60000) :: contlistscpi_f,contlistscpj_f
! integer :: ifirstrun,ilist_scp_first
! real(kind=8),dimension(:,:),allocatable :: alphapol,epshead,&
newcontlistcatscangfk,&
newcontlistcatscangti,newcontlistcatscangtj,&
newcontlistcatscangtk,newcontlistcatscangtl,&
- newcontlistcatcatnormi,newcontlistcatcatnormj
+ newcontlistcatcatnormi,newcontlistcatcatnormj,&
+ newcontlistmartsci,newcontlistmartscj,&
+ newcontlistmartpi,newcontlistmartpj
+
g_listcatscang_start,g_listcatscang_end,&
g_listcatscangf_start,g_listcatscangf_end,&
g_listcatscangt_start,g_listcatscangt_end,&
- g_listcatcatnorm_start,g_listcatcatnorm_end,g_ilist_catcatnorm
+ g_listcatcatnorm_start,g_listcatcatnorm_end,g_ilist_catcatnorm,&
+ g_ilist_martp,g_ilist_martsc,&
+ g_listmartsc_start,g_listmartsc_end,&
+ g_listmartp_start,g_listmartp_end
+
! MARTINI FORCE FIELD