double precision Kh(MaxQ,MaxR,MaxT_h,max_parm),
& q0(MaxQ,MaxR,MaxT_h,max_parm),delta,deltrms,deltrgy,fimin,
& f(maxR,maxT_h,max_parm),beta_h(MaxT_h,max_parm)
double precision Kh(MaxQ,MaxR,MaxT_h,max_parm),
& q0(MaxQ,MaxR,MaxT_h,max_parm),delta,deltrms,deltrgy,fimin,
& f(maxR,maxT_h,max_parm),beta_h(MaxT_h,max_parm)