diff --git a/src/gray.f b/src/gray.f
index 7f1ee7e..47dd34f 100644
--- a/src/gray.f
+++ b/src/gray.f
@@ -910,7 +910,11 @@ c     versus psi, rhop, rhot
       end if
 
       if (iequil.eq.1)   call surf_anal
-
+c
+c     print psi,X,Y on the R,z grid
+c
+      call print_2d
+c
       if (iequil.ne.2.or.ipass.lt.0) then
 c     set simple limiter as two cylindrical walls at rwallm and r00
         nlim=5
@@ -1642,6 +1646,34 @@ c
       end
 c
 c
+c
+      subroutine print_2d
+      implicit real*8 (a-h,o-z)
+      parameter(nnw=501,nnh=501)
+      dimension rv(nnw),zv(nnh),psin(nnw,nnh)
+      dimension btotal(nnw,nnh)
+c
+      common/cpsin/rv,zv,psin
+      common/psival/psinv
+      common/eqnn/nr,nz,npp,nintp
+      common/parpl/brr,bphi,bzz,ajphi
+      common/xgxg/xg
+      common/ygyg/yg
+c
+      do j=1,nr
+        rjcm=rv(j)*1.d2
+        do k=1,nz
+          zkcm=zv(k)*1.d2
+          call plas_deriv(rjcm,0.d0,zkcm)
+          write(91,111) rv(j),zv(k),psinv,brr,bphi,bzz,xg,yg
+        enddo
+        write(91,*) ' '
+      enddo
+      return
+ 111  format(12(1x,e12.5))
+      end
+c
+c
 c
       subroutine print_prof
       implicit real*8 (a-h,o-z)