-c write (iout,*) "nphi+i",nphi+i," gloc",gloc(nphi+i,icg)
-c write (iout,*) "i",i," gradc",(gradc(j,i,icg),j=1,3)
-c write (iout,*) "i",i," gradx",(gradx(j,i,icg),j=1,3)
+ write (iout,*) "nphi+i",nphi+i," gloc",gloc(nphi+i,icg)
+ write (iout,*) "i",i," gradc",(gradc(j,i,icg),j=1,3)
+ write (iout,*) "i",i," gradx",(gradx(j,i,icg),j=1,3)