double precision efree(maxT,maxdimnat+2,maxprot),
     & efree0(maxT,maxprot),
     & elowest(maxT,maxdimnat+2,maxprot),
     & elowest_ent(maxT,maxdimnat+2,maxprot),
     & emean_ftot(maxT,maxprot),
     & ebis_ftot(maxT,maxprot),
     & eave_pftot(max_ene,maxT,maxprot),
     & eave_nat_pftot(max_ene,maxnatlike,maxT,maxprot),
     & eave_pfprim(max_ene,maxT,maxprot),
     & eave_pfprimtot(max_ene,maxT,maxprot),
     & eave_pfbistot(max_ene,maxT,maxprot),
     & esquare_ftot(maxT,maxprot),
     & emix_pftot(max_ene,maxT,maxprot),
     & emix_pfprimtot(max_ene,maxT,maxprot),
     & emix_pfbistot(max_ene,maxT,maxprot),
     & emixsq_pftot(max_ene,maxT,maxprot), 
     & eneps(3,nntyp,maxstr_proc,maxprot),
     & eneps_temp(3,nntyp),
     & etor_temp(0:2*maxterm,0:3,-ntyp+1:ntyp-1,-ntyp+1:ntyp-1,2),
     & ebend_temp_kcc(maxang_kcc,0:ntyp),
     & etor_temp_kcc(maxval_kcc,maxval_kcc,0:2*maxterm,
     & -ntyp+1:ntyp-1,-ntyp+1:ntyp-1),
     & eneps_mix_ftot(nntyp,maxT,maxprot),
     & eneps_mix_fbistot(nntyp,maxT,maxprot),
     & eneps_mixsq_ftot(nntyp,maxT,maxprot),
     & enepsave_ftot(nntyp,maxT,maxprot),
     & enepsave_nat_ftot(nntyp,maxT,maxnatlike,maxprot),
     & enepslowest_f(nntyp,maxT,maxprot),
     & entor(maxtor_var,maxstr_proc,maxprot),
     & entor_mix_fprimtot(maxtor_var,maxT,maxprot),
     & entor_mix_ftot(maxtor_var,maxT,maxprot),
     & entor_mix_fbistot(maxtor_var,maxT,maxprot),
     & entor_mixsq_ftot(maxtor_var,maxT,maxprot),
     & entorave_ftot(maxtor_var,maxT,maxprot),
     & entorave_nat_ftot(maxtor_var,maxT,maxnatlike,maxprot),
     & entorave_fprimtot(maxtor_var,maxT,maxprot),
     & entorave_fbistot(maxtor_var,maxT,maxprot),
     & entorlowest_f(maxtor_var,maxT,maxprot),
     & nuepsave_f(nntyp,maxdimnat,maxT,maxnatlike,maxprot),
     & nutorave_f(maxtor_var,maxdimnat,maxT,maxnatlike,maxprot),
     & rder(max_paropt,maxdimnat),
     & reps(nntyp,maxdimnat),rtor(maxtor_var,maxdimnat),
     & nuave(maxdimnat,maxT,maxnatlike,maxprot),
     & nu_pf(max_ene,maxdimnat,maxT,maxnatlike,maxprot),
     & zvtot(maxT,maxprot),
     & heatbase(maxprot),
     & zvdertot(max_paropt,maxT,maxprot),
     & zvepstot(nntyp,maxT,maxprot),
     & zvtortot(maxtor_var,maxT,maxprot),
     & sumlik(maxT,maxprot),
     & sumlikder(max_paropt,maxT,maxprot),
     & sumlikeps(nntyp,maxT,maxprot),
     & sumliktor(maxtor_var,maxT,maxprot),
     & sumqder(max_paropt,maxT,maxprot),
     & sumqeps(nntyp,maxT,maxprot),
     & sumqtor(maxtor_var,maxT,maxprot)
      integer ind_lowest(maxT,maxnatlike+2,maxprot)
      common /weightder/ eneps,eneps_temp,etor_temp,ebend_temp_kcc,
     & etor_temp_kcc,enepsave_ftot,eneps_mix_ftot,eneps_mix_fbistot,
     & eneps_mixsq_ftot,entor,
     & entorave_ftot,entorave_fprimtot,entorave_fbistot,
     & entor_mix_ftot,entor_mix_fprimtot,entor_mix_fbistot,
     & entor_mixsq_ftot,
     & efree,efree0,elowest,elowest_ent,
     & emean_ftot,ebis_ftot,
     & eave_pftot,eave_pfprimtot,
     & eave_pfbistot,emix_pftot,
     & emix_pfprimtot,emix_pfbistot,
     & esquare_ftot,emixsq_pftot,nuave,
     & nu_pf,nuepsave_f,
     & nutorave_f,
     & rder,reps,zvtot,heatbase,zvdertot,zvepstot,
     & rtor,zvtortot,
     & sumqder,sumqeps,sumqtor,sumlik,sumlikder,sumlikeps,sumliktor,
     & enepslowest_f,
     & entorlowest_f,ind_lowest