--- MITgcm/model/src/write_grid.F	2012/06/01 03:19:53	1.35
+++ MITgcm/model/src/write_grid.F	2013/07/26 14:22:22	1.38
@@ -1,9 +1,9 @@
-C $Header: /home/ubuntu/mnt/e9_copy/MITgcm/model/src/write_grid.F,v 1.35 2012/06/01 03:19:53 jmc Exp $
+C $Header: /home/ubuntu/mnt/e9_copy/MITgcm/model/src/write_grid.F,v 1.38 2013/07/26 14:22:22 jmc Exp $
 C $Name:  $
 
-#include "AD_CONFIG.h"
 #include "PACKAGES_CONFIG.h"
 #include "CPP_OPTIONS.h"
+#include "AD_CONFIG.h"
 
 C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|
 CBOP
@@ -93,6 +93,8 @@
         IF ( usingCurvilinearGrid .OR. rotateGrid ) THEN
           CALL WRITE_FLD_XY_RS( 'AngleCS',' ',angleCosC,-1,myThid)
           CALL WRITE_FLD_XY_RS( 'AngleSN',' ',angleSinC,-1,myThid)
+          CALL WRITE_FLD_XY_RS( 'U2zonDir',' ',u2zonDir,-1,myThid)
+          CALL WRITE_FLD_XY_RS( 'V2zonDir',' ',v2zonDir,-1,myThid)
         ENDIF
 C     Write 3D geometry arrays
         CALL WRITE_FLD_XY_RS( 'Depth',' ',tmpfld,-1,myThid)
@@ -121,7 +123,7 @@
 C     Write vertical grid arrays
         CALL WRITE_GLVEC_RS( 'RC', ' ', rC,  Nr, -1, myThid )
         CALL WRITE_GLVEC_RS( 'RF', ' ', rF, 1+Nr,-1, myThid )
-        CALL WRITE_GLVEC_RS( 'DRC',' ', drC, Nr, -1, myThid )
+        CALL WRITE_GLVEC_RS( 'DRC',' ', drC,1+Nr,-1, myThid )
         CALL WRITE_GLVEC_RS( 'DRF',' ', drF, Nr, -1, myThid )
         IF ( selectSigmaCoord.NE.0 ) THEN
          CALL WRITE_GLVEC_RS( 'AHybSigF',' ',aHybSigmF,1+Nr,-1,myThid )