if (ndih_constr.gt.0) then
allocate(idih_constr(ndih_constr),idih_nconstr(ndih_constr)) !(maxdih_constr)
allocate(phi0(ndih_constr),drange(ndih_constr)) !(maxdih_constr)
if (ndih_constr.gt.0) then
allocate(idih_constr(ndih_constr),idih_nconstr(ndih_constr)) !(maxdih_constr)
allocate(phi0(ndih_constr),drange(ndih_constr)) !(maxdih_constr)