double precision PMFtab(0:maxHdim,maxT_h,maxN),delta_q
      integer tmin(0:maxT_h,maxN),tmax(maxT_h,maxN)
      common /PMF/ PMFtab,delta_q,tmin,tmax