C & rad2deg*theta_constr0(i), rad2deg*theta_drange(i),
C & rad2deg*difi,0.25d0*for_thet_constr(i)*difi**4,
C & gloc(itheta+nphi-2,icg)
C & rad2deg*theta_constr0(i), rad2deg*theta_drange(i),
C & rad2deg*difi,0.25d0*for_thet_constr(i)*difi**4,
C & gloc(itheta+nphi-2,icg)