/[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.2 by adcroft, Tue May 29 14:01:37 2001 UTC revision 1.5.4.2 by heimbach, Fri Mar 7 23:10:21 2003 UTC
# Line 3  C $Name$ Line 3  C $Name$
3    
4  #include "CPP_OPTIONS.h"  #include "CPP_OPTIONS.h"
5    
6    CBOP
7    C     !ROUTINE: INI_DYNVARS
8    C     !INTERFACE:
9        SUBROUTINE INI_DYNVARS( myThid )        SUBROUTINE INI_DYNVARS( myThid )
10  C     /==========================================================\  C     !DESCRIPTION: \bv
11  C     | SUBROUTINE INI_DYNVARS                                   |  C     *==========================================================*
12  C     | o Initialise to zero all DYNVARS.h arrays                |  C     | SUBROUTINE INI_DYNVARS                                    
13  C     \==========================================================/  C     | o Initialise to zero all DYNVARS.h arrays                  
14        IMPLICIT NONE  C     *==========================================================*
15    C     | Sets all the State variables to zero.
16  C     Sets all the State variables to zero.  C     | Sets all the Gs (arrays used for multi-level time-stepping)
17  C     Sets all the Gs (arrays used for multi-level time-stepping)  C     | to zero.
18  C     to zero.  C     *==========================================================*
19    C     \ev
20    
21    C     !USES:
22          IMPLICIT NONE
23  C     === Global variables ===  C     === Global variables ===
24  #include "SIZE.h"  #include "SIZE.h"
25  #include "EEPARAMS.h"  #include "EEPARAMS.h"
# Line 23  C     === Global variables === Line 29  C     === Global variables ===
29  #include "GW.h"  #include "GW.h"
30  #endif  #endif
31    
32    C     !INPUT/OUTPUT PARAMETERS:
33  C     == Routine arguments ==  C     == Routine arguments ==
34  C     myThid -  Number of this instance of INI_UVEL  C     myThid -  Number of this instance of INI_UVEL
35        INTEGER myThid        INTEGER myThid
 CEndOfInterface  
36    
37    C     !LOCAL VARIABLES:
38  C     == Local variables ==  C     == Local variables ==
39  C     bi,bj  - Loop counters  C     bi,bj  - Loop counters
40  C     I,J,K  C     I,J,K
41        INTEGER bi, bj        INTEGER bi, bj
42        INTEGER  I,  J, K        INTEGER  I,  J, K
43    CEOP
44    
45  C--   Over all tiles  C--   Over all tiles
46        DO bj = myByLo(myThid), myByHi(myThid)        DO bj = myByLo(myThid), myByHi(myThid)
# Line 43  C-      3D arrays Line 51  C-      3D arrays
51           DO J=1-Oly,sNy+Oly           DO J=1-Oly,sNy+Oly
52            DO I=1-Olx,sNx+Olx            DO I=1-Olx,sNx+Olx
53    
54             uVel(I,J,K,bi,bj)=0.             uVel(I,J,K,bi,bj)=0. _d 0
55             vVel(I,J,K,bi,bj)=0.             vVel(I,J,K,bi,bj)=0. _d 0
56             wVel(I,J,K,bi,bj)=0.             wVel(I,J,K,bi,bj)=0. _d 0
57             theta(I,J,K,bi,bj)=0.             theta(I,J,K,bi,bj)=0. _d 0
58             salt(I,J,K,bi,bj)=0.             salt(I,J,K,bi,bj)=0. _d 0
59    
60             gU(I,J,K,bi,bj)=0.             gU(I,J,K,bi,bj)=0. _d 0
61             gV(I,J,K,bi,bj)=0.             gV(I,J,K,bi,bj)=0. _d 0
62             gT(I,J,K,bi,bj)=0.             gT(I,J,K,bi,bj)=0. _d 0
63             gS(I,J,K,bi,bj)=0.             gS(I,J,K,bi,bj)=0. _d 0
64             gUnm1(I,J,K,bi,bj)=0.             gUnm1(I,J,K,bi,bj)=0. _d 0
65             gVnm1(I,J,K,bi,bj)=0.             gVnm1(I,J,K,bi,bj)=0. _d 0
66             gTnm1(I,J,K,bi,bj)=0.             gTnm1(I,J,K,bi,bj)=0. _d 0
67             gSnm1(I,J,K,bi,bj)=0.             gSnm1(I,J,K,bi,bj)=0. _d 0
68    
69  #ifdef INCLUDE_CD_CODE  #ifdef INCLUDE_CD_CODE
70             guCD(I,J,K,bi,bj)=0.             guCD(I,J,K,bi,bj)=0. _d 0
71             gvCD(I,J,K,bi,bj)=0.             gvCD(I,J,K,bi,bj)=0. _d 0
72             uNM1(I,J,K,bi,bj)=0.             uNM1(I,J,K,bi,bj)=0. _d 0
73             vNM1(I,J,K,bi,bj)=0.             vNM1(I,J,K,bi,bj)=0. _d 0
74             uVeld(I,J,K,bi,bj)=0.             uVeld(I,J,K,bi,bj)=0. _d 0
75             vVeld(I,J,K,bi,bj)=0.             vVeld(I,J,K,bi,bj)=0. _d 0
76  #endif  #endif
77    
78  #ifdef ALLOW_NONHYDROSTATIC  #ifdef ALLOW_NONHYDROSTATIC
79             gW(I,J,K,bi,bj)=0.             gW(I,J,K,bi,bj)=0. _d 0
80             gWnm1(I,J,K,bi,bj)=0.             gWnm1(I,J,K,bi,bj)=0. _d 0
81  #endif  #endif
82    
83               totPhiHyd(I,J,K,bi,bj)=0. _d 0
84            ENDDO            ENDDO
85           ENDDO           ENDDO
86          ENDDO          ENDDO
# Line 79  C-      3D arrays Line 88  C-      3D arrays
88  C-      2D arrays  C-      2D arrays
89          DO J=1-Oly,sNy+Oly          DO J=1-Oly,sNy+Oly
90           DO I=1-Olx,sNx+Olx           DO I=1-Olx,sNx+Olx
91            etaN(I,J,bi,bj)=0.            etaN(I,J,bi,bj)=0. _d 0
92            etaNm1(I,J,bi,bj)=0.            etaH(I,J,bi,bj)=0. _d 0
93    #ifdef INCLUDE_CD_CODE
94              etaNm1(I,J,bi,bj)=0. _d 0
95    #endif
96              phiHydLow(I,J,bi,bj)=0. _d 0
97           ENDDO           ENDDO
98          ENDDO          ENDDO
99    
Legend:
Removed from v.1.2  
changed lines
  Added in v.1.5.4.2
  ViewVC Help
Powered by ViewVC 1.1.22