double precision zz,y,yc,diff,h,w,thet_point,pmf_thet,wthet_PMF integer npoint,np,npel,nthet_point common /deriv/ zz(8,maxpoint),y(4,maxpoint), & yc(4,maxpoint),diff(4,maxpoint),h(maxvar,maxvar), & w(4,maxpoint),npoint,np(-2:2,-2:2),npel(-2:2,-2:2) common /PMFtheta/ thet_point(180,0:ntyp),pmf_thet(180,0:ntyp), & wthet_PMF,nthet_point(0:ntyp)