--- MITgcm/model/src/ini_dynvars.F	2001/07/13 14:26:57	1.3
+++ MITgcm/model/src/ini_dynvars.F	2003/03/07 23:10:21	1.5.4.2
@@ -1,19 +1,25 @@
-C $Header: /home/ubuntu/mnt/e9_copy/MITgcm/model/src/ini_dynvars.F,v 1.3 2001/07/13 14:26:57 heimbach Exp $
+C $Header: /home/ubuntu/mnt/e9_copy/MITgcm/model/src/ini_dynvars.F,v 1.5.4.2 2003/03/07 23:10:21 heimbach Exp $
 C $Name:  $
 
 #include "CPP_OPTIONS.h"
 
+CBOP
+C     !ROUTINE: INI_DYNVARS
+C     !INTERFACE:
       SUBROUTINE INI_DYNVARS( myThid )
-C     /==========================================================\
-C     | SUBROUTINE INI_DYNVARS                                   |
-C     | o Initialise to zero all DYNVARS.h arrays                | 
-C     \==========================================================/
-      IMPLICIT NONE
-
-C     Sets all the State variables to zero.
-C     Sets all the Gs (arrays used for multi-level time-stepping)
-C     to zero.
+C     !DESCRIPTION: \bv
+C     *==========================================================*
+C     | SUBROUTINE INI_DYNVARS                                    
+C     | o Initialise to zero all DYNVARS.h arrays                  
+C     *==========================================================*
+C     | Sets all the State variables to zero.
+C     | Sets all the Gs (arrays used for multi-level time-stepping)
+C     | to zero.
+C     *==========================================================*
+C     \ev
 
+C     !USES:
+      IMPLICIT NONE
 C     === Global variables ===
 #include "SIZE.h"
 #include "EEPARAMS.h"
@@ -23,16 +29,18 @@
 #include "GW.h"
 #endif
 
+C     !INPUT/OUTPUT PARAMETERS:
 C     == Routine arguments ==
 C     myThid -  Number of this instance of INI_UVEL 
       INTEGER myThid
-CEndOfInterface
 
+C     !LOCAL VARIABLES:
 C     == Local variables ==
 C     bi,bj  - Loop counters
 C     I,J,K
       INTEGER bi, bj
       INTEGER  I,  J, K
+CEOP
 
 C--   Over all tiles
       DO bj = myByLo(myThid), myByHi(myThid)
@@ -72,6 +80,7 @@
            gWnm1(I,J,K,bi,bj)=0. _d 0
 #endif
 
+           totPhiHyd(I,J,K,bi,bj)=0. _d 0
           ENDDO
          ENDDO
         ENDDO
@@ -80,7 +89,11 @@
         DO J=1-Oly,sNy+Oly
          DO I=1-Olx,sNx+Olx
           etaN(I,J,bi,bj)=0. _d 0
+          etaH(I,J,bi,bj)=0. _d 0
+#ifdef INCLUDE_CD_CODE
           etaNm1(I,J,bi,bj)=0. _d 0
+#endif
+          phiHydLow(I,J,bi,bj)=0. _d 0
          ENDDO
         ENDDO