-C Calculate patrial derivative for theta angle
-#ifdef NEWCORR
- geel_loc_ij=a22*gmuij1(1)
- & +a23*gmuij1(2)
- & +a32*gmuij1(3)
- & +a33*gmuij1(4)
-c write(iout,*) "derivative over thatai"
-c write(iout,*) a22*gmuij1(1), a23*gmuij1(2) ,a32*gmuij1(3),
-c & a33*gmuij1(4)
- gloc(nphi+i,icg)=gloc(nphi+i,icg)+
- & geel_loc_ij*wel_loc
-c write(iout,*) "derivative over thatai-1"
-c write(iout,*) a22*gmuij2(1), a23*gmuij2(2) ,a32*gmuij2(3),
-c & a33*gmuij2(4)
- geel_loc_ij=a22*gmuij2(1)+a23*gmuij2(2)+a32*gmuij2(3)
- & +a33*gmuij2(4)
- gloc(nphi+i-1,icg)=gloc(nphi+i-1,icg)+
- & geel_loc_ij*wel_loc
- geel_loc_ji=a22*gmuji1(1)+a23*gmuji1(2)+a32*gmuji1(3)
- & +a33*gmuji1(4)
-c write(iout,*) "derivative over thataj"
-c write(iout,*) a22*gmuji1(1), a23*gmuji1(2) ,a32*gmuji1(3),
-c & a33*gmuji1(4)