& f(maxR,maxT_h,max_parm),beta_h(MaxT_h,max_parm)
integer nR(maxT_h,max_parm),snk(MaxR,MaxT_h,max_parm,MaxSlice),
& nT_h(max_parm),maxit,totraj(maxR,max_parm),nRR(maxT_h,max_parm)
logical replica(max_parm),umbrella(max_parm),read_iset(max_parm)
& f(maxR,maxT_h,max_parm),beta_h(MaxT_h,max_parm)
integer nR(maxT_h,max_parm),snk(MaxR,MaxT_h,max_parm,MaxSlice),
& nT_h(max_parm),maxit,totraj(maxR,max_parm),nRR(maxT_h,max_parm)
logical replica(max_parm),umbrella(max_parm),read_iset(max_parm)