+ if (ntheta_constr.eq.0) then
+ 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
& ' ivec_start',ivec_start,' ivec_end',ivec_end,
& ' iset_start',iset_start,' iset_end',iset_end,
& ' idihconstr_start',idihconstr_start,' idihconstr_end',
& ' ivec_start',ivec_start,' ivec_end',ivec_end,
& ' iset_start',iset_start,' iset_end',iset_end,
& ' idihconstr_start',idihconstr_start,' idihconstr_end',
write (*,*) 'Processor:',fg_rank,myrank,' igrad_start',
& igrad_start,' igrad_end',igrad_end,' ngrad_start',ngrad_start,
& ' ngrad_end',ngrad_end
write (*,*) 'Processor:',fg_rank,myrank,' igrad_start',
& igrad_start,' igrad_end',igrad_end,' ngrad_start',ngrad_start,
& ' ngrad_end',ngrad_end