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
& ' ivec_start',ivec_start,' ivec_end',ivec_end,
& ' iset_start',iset_start,' iset_end',iset_end,
& ' idihconstr_start',idihconstr_start,' idihconstr_end',
- & idihconstr_end
+ & idihconstr_end,
& ' ithetaconstr_start',ithetaconstr_start,' ithetaconstr_end',
& ithetaconstr_end