double precision ww,ww0,ww_low,ww_up,ww_orig,ww_oorig,x_orig,
     &  epp_low,epp_up,rpp_low,rpp_up,elpp6_low,elpp6_up,elpp3_low,
     &  elpp3_up,b_low,b_up,bnew1_low,bnew1_up,bnew2_low,bnew2_up,
     &  ccnew_low,ccnew_up,ddnew_low,ddnew_up,e0new_low,e0new_up,
     &  eenew_low,eenew_up,epscp_low,epscp_up,rscp_low,rscp_up,
     &  bnew1tor_low,bnew1tor_up,bnew2tor_low,bnew2tor_up,
     &  ccnewtor_low,ccnewtor_up,ddnewtor_low,ddnewtor_up,e0newtor_low,
     &  e0newtor_up,eenewtor_low,eenewtor_up,e0,wello_PMF,wturn_PMF,
     &  wello_PMF_low,wello_PMF_up,wturn_PMF_low,wturn_PMF_up,
     &  e0_low,e0_up,
     &  x_up,x_low,xm,epss_low,epss_up,epsp_low,epsp_up,escal,
     &  escalprim,escalbis,weitor,weitor_low,weitor_up,sigma_low,
     &  sigma_up,v1bend_low,v1bend_up,wpmf
      integer imask,mask_elec,mask_fourier,mask_bnew1,mask_bnew2,
     &  mask_bnewtor1,mask_bnewtor2,
     &  mask_ccnew,mask_ddnew,mask_e0new,mask_eenew,mod_fourier,
     &  mod_fouriertor,mask_fouriertor,
     &  mask_bnew1tor,mask_bnew2tor,mask_ccnewtor,mask_ddnewtor,
     &  mask_e0newtor,mask_eenewtor,
     &  mask_scp,iw,nsingle_sc,npair_sc,ityp_ssc,ityp_psc,mask_tor,
     &  mask_ang,mask_sigma,mask_e0,mask_ello_PMF,mask_turn_PMF,
     &  n_weight_PMF
      logical mod_other_params,mod_elec,mod_scp,mod_side,mod_sccor,
     &  mod_tor,mod_side_other,mod_ang,torsion_pmf,turn_pmf,eello_pmf
      common /chujec/ ww(max_ene),ww0(max_ene),ww_low(max_ene),
     &  ww_up(max_ene),ww_orig(max_ene),ww_oorig(max_ene),wpmf,
     &  weitor(0:ntyp,-ntyp:ntyp,-ntyp:ntyp,2),
     &  weitor_low(0:ntyp,-ntyp:ntyp,-ntyp:ntyp,2),
     &  weitor_up(0:ntyp,-ntyp:ntyp,-ntyp:ntyp,2),x_orig(max_paropt),
     &  epp_low(2,2),epp_up(2,2),rpp_low(2,2),rpp_up(2,2),
     &  elpp6_low(2,2),elpp6_up(2,2),elpp3_low(2,2),elpp3_up(2,2),
     &  b_low(13,0:ntyp),b_up(13,0:ntyp),
     &  bnew1_low(3,2,0:ntyp),bnew1_up(3,2,0:ntyp),
     &  bnew2_low(3,2,0:ntyp),bnew2_up(3,2,0:ntyp),
     &  ccnew_low(3,2,0:ntyp),ccnew_up(3,2,0:ntyp),
     &  ddnew_low(3,2,0:ntyp),ddnew_up(3,2,0:ntyp),
     &  e0new_low(3,0:ntyp),e0new_up(3,0:ntyp),
     &  eenew_low(2,2,2,0:ntyp),eenew_up(2,2,2,0:ntyp),
     &  bnew1tor_low(3,2,0:ntyp),bnew1tor_up(3,2,0:ntyp),
     &  bnew2tor_low(3,2,0:ntyp),bnew2tor_up(3,2,0:ntyp),
     &  ccnewtor_low(3,2,0:ntyp),ccnewtor_up(3,2,0:ntyp),
     &  ddnewtor_low(3,2,0:ntyp),ddnewtor_up(3,2,0:ntyp),
     &  e0newtor_low(3,0:ntyp),e0newtor_up(3,0:ntyp),
     &  eenewtor_low(2,2,2,0:ntyp),eenewtor_up(2,2,2,0:ntyp),
     &  e0(0:ntyp,-ntyp:ntyp),wello_PMF,wturn_PMF,
     &  e0_low(0:ntyp,-ntyp:ntyp),e0_up(0:ntyp,-ntyp:ntyp),
     &  wello_PMF_low,wello_PMF_up,wturn_PMF_low,wturn_PMF_up,
     &  v1bend_low(maxang_kcc,0:ntyp),v1bend_up(maxang_kcc,0:ntyp),
     &  x_up(max_paropt),x_low(max_paropt),
     &  epscp_low(0:ntyp,2),epscp_up(0:ntyp,2),rscp_low(0:ntyp,2),
     &  rscp_up(0:ntyp,2),epss_low(ntyp),epss_up(ntyp),epsp_low(nntyp),
     &  epsp_up(nntyp),sigma_low(5,ntyp),sigma_up(5,ntyp),
     &  xm(max_paropt),
     &  escal(max_ene,maxT,maxbatch,maxprot),
     &  escalprim(max_ene,maxT,maxbatch,maxprot),
     &  escalbis(max_ene,maxT,maxbatch,maxprot),
     &  imask(max_ene),nsingle_sc,npair_sc,ityp_ssc(ntyp),
     &  ityp_psc(2,nntyp),mask_elec(2,2,4),
     &  mask_tor(0:ntyp,-ntyp:ntyp,-ntyp:ntyp,2),
     &  mask_ang(0:ntyp),
     &  mask_fourier(13,0:ntyp),mask_bnew1(3,2,0:ntyp),
     &  mask_bnew2(3,2,0:ntyp),mask_bnewtor1(3,2,0:ntyp),
     &  mask_bnewtor2(3,2,0:ntyp),
     &  mask_ccnew(3,2,0:ntyp),mask_ddnew(3,2,0:ntyp),
     &  mask_e0new(3,0:ntyp),
     &  mask_eenew(2,2,2,0:ntyp),mask_sigma(5,ntyp),
     &  mask_fouriertor(13,0:ntyp),mask_bnew1tor(3,2,0:ntyp),
     &  mask_bnew2tor(3,2,0:ntyp),
     &  mask_ccnewtor(3,2,0:ntyp),mask_ddnewtor(3,2,0:ntyp),
     &  mask_e0newtor(3,0:ntyp),
     &  mask_eenewtor(2,2,2,0:ntyp),
     &  mask_scp(0:ntyp,2,2),mod_other_params,mod_fourier(0:ntyp),
     &  mask_e0(0:ntyp,-ntyp:ntyp),mask_ello_PMF,mask_turn_PMF,
     &  mod_fouriertor(0:ntyp),
     &  mod_elec,mod_scp,mod_side,mod_side_other,mod_sccor,mod_tor,
     &  mod_ang,iw(max_ene),torsion_pmf,turn_pmf,eello_pmf,n_weight_PMF