integer ndih_constr,idih_constr(maxdih_constr),ntheta_constr,
     & itheta_constr(maxdih_constr)
      integer ndih_nconstr,idih_nconstr(maxdih_constr)
      integer idihconstr_start,idihconstr_end,ithetaconstr_start,
     & ithetaconstr_end
      double precision phi0(maxdih_constr),drange(maxdih_constr),
     & ftors(maxdih_constr),theta_constr0(maxdih_constr),
     & theta_drange(maxdih_constr),for_thet_constr(maxdih_constr)
      common /torcnstr/ phi0,drange,ftors,ndih_constr,idih_constr,
     &  ndih_nconstr,idih_nconstr,idihconstr_start,idihconstr_end