integer ndih_constr,idih_constr(maxdih_constr)
      integer ndih_nconstr,idih_nconstr(maxdih_constr)
      double precision phi0(maxdih_constr),drange(maxdih_constr),
     & ftors(maxdih_constr)
      common /torcnstr/ phi0,drange,ftors,ndih_constr,idih_constr,
     &                  ndih_nconstr,idih_nconstr