call int_bounds(ndih_constr,idihconstr_start,idihconstr_end)
endif
if (ntheta_constr.eq.0) then
- idihconstr_start=1
- idihconstr_end=0
+ ithetaconstr_start=1
+ ithetaconstr_end=0
else
call int_bounds
& (ntheta_constr,ithetaconstr_start,ithetaconstr_end)
endif
+C print *,ntheta_constr,ithetaconstr_start,ithetaconstr_end
c nsumgrad=(nres-nnt)*(nres-nnt+1)/2
c nlen=nres-nnt+1
nsumgrad=(nres-nnt)*(nres-nnt+1)/2