/[MITgcm]/MITgcm/model/src/ini_dynvars.F
ViewVC logotype

Diff of /MITgcm/model/src/ini_dynvars.F

Parent Directory Parent Directory | Revision Log Revision Log | View Revision Graph Revision Graph | View Patch Patch

revision 1.12 by edhill, Thu Oct 30 12:00:41 2003 UTC revision 1.20 by jmc, Sun Aug 24 21:32:17 2008 UTC
# Line 1  Line 1 
1  C $Header$  C $Header$
2  C $Name$  C $Name$
3    
 #include "PACKAGES_CONFIG.h"  
4  #include "CPP_OPTIONS.h"  #include "CPP_OPTIONS.h"
5    
6  CBOP  CBOP
# Line 10  C     !INTERFACE: Line 9  C     !INTERFACE:
9        SUBROUTINE INI_DYNVARS( myThid )        SUBROUTINE INI_DYNVARS( myThid )
10  C     !DESCRIPTION: \bv  C     !DESCRIPTION: \bv
11  C     *==========================================================*  C     *==========================================================*
12  C     | SUBROUTINE INI_DYNVARS                                      C     | SUBROUTINE INI_DYNVARS
13  C     | o Initialise to zero all DYNVARS.h arrays                    C     | o Initialise to zero all DYNVARS.h arrays
14  C     *==========================================================*  C     *==========================================================*
15  C     | Sets all the State variables to zero.  C     | Sets all the State variables to zero.
16  C     | Sets all the Gs (arrays used for multi-level time-stepping)  C     | Sets all the Gs (arrays used for multi-level time-stepping)
# Line 26  C     === Global variables === Line 25  C     === Global variables ===
25  #include "EEPARAMS.h"  #include "EEPARAMS.h"
26  #include "PARAMS.h"  #include "PARAMS.h"
27  #include "DYNVARS.h"  #include "DYNVARS.h"
 #ifdef ALLOW_NONHYDROSTATIC  
 #include "GW.h"  
 #endif  
 #ifdef ALLOW_CD_CODE  
 #include "CD_CODE_VARS.h"  
 #endif  
28    
29  C     !INPUT/OUTPUT PARAMETERS:  C     !INPUT/OUTPUT PARAMETERS:
30  C     == Routine arguments ==  C     == Routine arguments ==
31  C     myThid -  Number of this instance of INI_UVEL  C     myThid  :: my Thread Id. number
32        INTEGER myThid        INTEGER myThid
33    
34  C     !LOCAL VARIABLES:  C     !LOCAL VARIABLES:
35  C     == Local variables ==  C     == Local variables ==
36  C     bi,bj  - Loop counters  C     bi,bj   :: tile indices
37  C     I,J,K  C     I,J,K   :: loop index
38        INTEGER bi, bj        INTEGER bi, bj
39        INTEGER  I,  J, K        INTEGER I, J, K
40  CEOP  CEOP
41    
42  C--   Over all tiles  C--   Over all tiles
# Line 65  C-      3D arrays Line 58  C-      3D arrays
58             gV(I,J,K,bi,bj)=0. _d 0             gV(I,J,K,bi,bj)=0. _d 0
59             gT(I,J,K,bi,bj)=0. _d 0             gT(I,J,K,bi,bj)=0. _d 0
60             gS(I,J,K,bi,bj)=0. _d 0             gS(I,J,K,bi,bj)=0. _d 0
61             gUnm1(I,J,K,bi,bj)=0. _d 0  #ifdef ALLOW_ADAMSBASHFORTH_3
62             gVnm1(I,J,K,bi,bj)=0. _d 0             guNm(I,J,K,bi,bj,1)=0. _d 0
63             gTnm1(I,J,K,bi,bj)=0. _d 0             gvNm(I,J,K,bi,bj,1)=0. _d 0
64             gSnm1(I,J,K,bi,bj)=0. _d 0             gtNm(I,J,K,bi,bj,1)=0. _d 0
65               gsNm(I,J,K,bi,bj,1)=0. _d 0
66  #ifdef ALLOW_CD_CODE             guNm(I,J,K,bi,bj,2)=0. _d 0
67             uNM1(I,J,K,bi,bj)=0. _d 0             gvNm(I,J,K,bi,bj,2)=0. _d 0
68             vNM1(I,J,K,bi,bj)=0. _d 0             gtNm(I,J,K,bi,bj,2)=0. _d 0
69             uVeld(I,J,K,bi,bj)=0. _d 0             gsNm(I,J,K,bi,bj,2)=0. _d 0
70             vVeld(I,J,K,bi,bj)=0. _d 0  #else
71  #endif             guNm1(I,J,K,bi,bj) =0. _d 0
72               gvNm1(I,J,K,bi,bj) =0. _d 0
73  #ifdef ALLOW_NONHYDROSTATIC             gtNm1(I,J,K,bi,bj) =0. _d 0
74  ceh3 needs an IF ( useNONHYDROSTATIC ) THEN             gsNm1(I,J,K,bi,bj) =0. _d 0
            gW(I,J,K,bi,bj)=0. _d 0  
            gWnm1(I,J,K,bi,bj)=0. _d 0  
75  #endif  #endif
   
76             totPhiHyd(I,J,K,bi,bj)=0. _d 0             totPhiHyd(I,J,K,bi,bj)=0. _d 0
77               IVDConvCount(i,j,k,bi,bj) = 0. _d 0
78    #ifdef ALLOW_ADDFLUID
79               addMass(I,J,K,bi,bj) = 0. _d 0
80    #endif
81            ENDDO            ENDDO
82           ENDDO           ENDDO
83          ENDDO          ENDDO
# Line 93  C-      2D arrays Line 87  C-      2D arrays
87           DO I=1-Olx,sNx+Olx           DO I=1-Olx,sNx+Olx
88            etaN(I,J,bi,bj)=0. _d 0            etaN(I,J,bi,bj)=0. _d 0
89            etaH(I,J,bi,bj)=0. _d 0            etaH(I,J,bi,bj)=0. _d 0
 #ifdef ALLOW_CD_CODE  
           etaNm1(I,J,bi,bj)=0. _d 0  
 #endif  
90            phiHydLow(I,J,bi,bj)=0. _d 0            phiHydLow(I,J,bi,bj)=0. _d 0
91              hMixLayer(I,J,bi,bj)=0. _d 0
92           ENDDO           ENDDO
93          ENDDO          ENDDO
94    
Legend:
Removed from v.1.12  
changed lines
  Added in v.1.20
  ViewVC Help
Powered by ViewVC 1.1.22