/[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.3 by heimbach, Fri Jul 13 14:26:57 2001 UTC revision 1.12 by edhill, Thu Oct 30 12:00:41 2003 UTC
# Line 1  Line 1 
1  C $Header$  C $Header$
2  C $Name$  C $Name$
3    
4    #include "PACKAGES_CONFIG.h"
5  #include "CPP_OPTIONS.h"  #include "CPP_OPTIONS.h"
6    
7    CBOP
8    C     !ROUTINE: INI_DYNVARS
9    C     !INTERFACE:
10        SUBROUTINE INI_DYNVARS( myThid )        SUBROUTINE INI_DYNVARS( myThid )
11  C     /==========================================================\  C     !DESCRIPTION: \bv
12  C     | SUBROUTINE INI_DYNVARS                                   |  C     *==========================================================*
13  C     | o Initialise to zero all DYNVARS.h arrays                |  C     | SUBROUTINE INI_DYNVARS                                    
14  C     \==========================================================/  C     | o Initialise to zero all DYNVARS.h arrays                  
15        IMPLICIT NONE  C     *==========================================================*
16    C     | Sets all the State variables to zero.
17  C     Sets all the State variables to zero.  C     | Sets all the Gs (arrays used for multi-level time-stepping)
18  C     Sets all the Gs (arrays used for multi-level time-stepping)  C     | to zero.
19  C     to zero.  C     *==========================================================*
20    C     \ev
21    
22    C     !USES:
23          IMPLICIT NONE
24  C     === Global variables ===  C     === Global variables ===
25  #include "SIZE.h"  #include "SIZE.h"
26  #include "EEPARAMS.h"  #include "EEPARAMS.h"
# Line 22  C     === Global variables === Line 29  C     === Global variables ===
29  #ifdef ALLOW_NONHYDROSTATIC  #ifdef ALLOW_NONHYDROSTATIC
30  #include "GW.h"  #include "GW.h"
31  #endif  #endif
32    #ifdef ALLOW_CD_CODE
33    #include "CD_CODE_VARS.h"
34    #endif
35    
36    C     !INPUT/OUTPUT PARAMETERS:
37  C     == Routine arguments ==  C     == Routine arguments ==
38  C     myThid -  Number of this instance of INI_UVEL  C     myThid -  Number of this instance of INI_UVEL
39        INTEGER myThid        INTEGER myThid
 CEndOfInterface  
40    
41    C     !LOCAL VARIABLES:
42  C     == Local variables ==  C     == Local variables ==
43  C     bi,bj  - Loop counters  C     bi,bj  - Loop counters
44  C     I,J,K  C     I,J,K
45        INTEGER bi, bj        INTEGER bi, bj
46        INTEGER  I,  J, K        INTEGER  I,  J, K
47    CEOP
48    
49  C--   Over all tiles  C--   Over all tiles
50        DO bj = myByLo(myThid), myByHi(myThid)        DO bj = myByLo(myThid), myByHi(myThid)
# Line 58  C-      3D arrays Line 70  C-      3D arrays
70             gTnm1(I,J,K,bi,bj)=0. _d 0             gTnm1(I,J,K,bi,bj)=0. _d 0
71             gSnm1(I,J,K,bi,bj)=0. _d 0             gSnm1(I,J,K,bi,bj)=0. _d 0
72    
73  #ifdef INCLUDE_CD_CODE  #ifdef ALLOW_CD_CODE
            guCD(I,J,K,bi,bj)=0. _d 0  
            gvCD(I,J,K,bi,bj)=0. _d 0  
74             uNM1(I,J,K,bi,bj)=0. _d 0             uNM1(I,J,K,bi,bj)=0. _d 0
75             vNM1(I,J,K,bi,bj)=0. _d 0             vNM1(I,J,K,bi,bj)=0. _d 0
76             uVeld(I,J,K,bi,bj)=0. _d 0             uVeld(I,J,K,bi,bj)=0. _d 0
# Line 68  C-      3D arrays Line 78  C-      3D arrays
78  #endif  #endif
79    
80  #ifdef ALLOW_NONHYDROSTATIC  #ifdef ALLOW_NONHYDROSTATIC
81    ceh3 needs an IF ( useNONHYDROSTATIC ) THEN
82             gW(I,J,K,bi,bj)=0. _d 0             gW(I,J,K,bi,bj)=0. _d 0
83             gWnm1(I,J,K,bi,bj)=0. _d 0             gWnm1(I,J,K,bi,bj)=0. _d 0
84  #endif  #endif
85    
86               totPhiHyd(I,J,K,bi,bj)=0. _d 0
87            ENDDO            ENDDO
88           ENDDO           ENDDO
89          ENDDO          ENDDO
# Line 80  C-      2D arrays Line 92  C-      2D arrays
92          DO J=1-Oly,sNy+Oly          DO J=1-Oly,sNy+Oly
93           DO I=1-Olx,sNx+Olx           DO I=1-Olx,sNx+Olx
94            etaN(I,J,bi,bj)=0. _d 0            etaN(I,J,bi,bj)=0. _d 0
95              etaH(I,J,bi,bj)=0. _d 0
96    #ifdef ALLOW_CD_CODE
97            etaNm1(I,J,bi,bj)=0. _d 0            etaNm1(I,J,bi,bj)=0. _d 0
98    #endif
99              phiHydLow(I,J,bi,bj)=0. _d 0
100           ENDDO           ENDDO
101          ENDDO          ENDDO
102    
Legend:
Removed from v.1.3  
changed lines
  Added in v.1.12
  ViewVC Help
Powered by ViewVC 1.1.22