/[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.19 by dimitri, Wed Nov 28 00:18:17 2007 UTC revision 1.27 by jmc, Tue Jan 20 20:46:55 2015 UTC
# Line 9  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 28  C     === Global variables === Line 28  C     === Global variables ===
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 44  C--   Over all tiles Line 44  C--   Over all tiles
44         DO bi = myBxLo(myThid), myBxHi(myThid)         DO bi = myBxLo(myThid), myBxHi(myThid)
45    
46  C-      3D arrays  C-      3D arrays
47          DO K=1,Nr          DO k=1,Nr
48           DO J=1-Oly,sNy+Oly           DO j=1-OLy,sNy+OLy
49            DO I=1-Olx,sNx+Olx            DO i=1-OLx,sNx+OLx
50    
51             uVel(I,J,K,bi,bj)=0. _d 0             uVel (i,j,k,bi,bj) = 0. _d 0
52             vVel(I,J,K,bi,bj)=0. _d 0             vVel (i,j,k,bi,bj) = 0. _d 0
53             wVel(I,J,K,bi,bj)=0. _d 0             wVel (i,j,k,bi,bj) = 0. _d 0
54             theta(I,J,K,bi,bj)=0. _d 0             theta(i,j,k,bi,bj) = 0. _d 0
55             salt(I,J,K,bi,bj)=0. _d 0             salt (i,j,k,bi,bj) = 0. _d 0
56    
57             gU(I,J,K,bi,bj)=0. _d 0             gU(i,j,k,bi,bj) = 0. _d 0
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  #ifdef USE_OLD_EXTERNAL_FORCING
60             gS(I,J,K,bi,bj)=0. _d 0             gT(i,j,k,bi,bj) = 0. _d 0
61               gS(i,j,k,bi,bj) = 0. _d 0
62    #endif
63  #ifdef ALLOW_ADAMSBASHFORTH_3  #ifdef ALLOW_ADAMSBASHFORTH_3
64             guNm(I,J,K,bi,bj,1)=0. _d 0             guNm(i,j,k,bi,bj,1) = 0. _d 0
65             gvNm(I,J,K,bi,bj,1)=0. _d 0             gvNm(i,j,k,bi,bj,1) = 0. _d 0
66             gtNm(I,J,K,bi,bj,1)=0. _d 0             gtNm(i,j,k,bi,bj,1) = 0. _d 0
67             gsNm(I,J,K,bi,bj,1)=0. _d 0             gsNm(i,j,k,bi,bj,1) = 0. _d 0
68             guNm(I,J,K,bi,bj,2)=0. _d 0             guNm(i,j,k,bi,bj,2) = 0. _d 0
69             gvNm(I,J,K,bi,bj,2)=0. _d 0             gvNm(i,j,k,bi,bj,2) = 0. _d 0
70             gtNm(I,J,K,bi,bj,2)=0. _d 0             gtNm(i,j,k,bi,bj,2) = 0. _d 0
71             gsNm(I,J,K,bi,bj,2)=0. _d 0             gsNm(i,j,k,bi,bj,2) = 0. _d 0
72  #else  #else
73             guNm1(I,J,K,bi,bj) =0. _d 0             guNm1(i,j,k,bi,bj) = 0. _d 0
74             gvNm1(I,J,K,bi,bj) =0. _d 0             gvNm1(i,j,k,bi,bj) = 0. _d 0
75             gtNm1(I,J,K,bi,bj) =0. _d 0             gtNm1(i,j,k,bi,bj) = 0. _d 0
76             gsNm1(I,J,K,bi,bj) =0. _d 0             gsNm1(i,j,k,bi,bj) = 0. _d 0
77  #endif  #endif
78               totPhiHyd (i,j,k,bi,bj) = 0. _d 0
79             totPhiHyd(I,J,K,bi,bj)=0. _d 0             rhoInSitu (i,j,k,bi,bj) = 0. _d 0
80             IVDConvCount(i,j,k,bi,bj) = 0. _d 0             IVDConvCount(i,j,k,bi,bj) = 0. _d 0
81            ENDDO            ENDDO
82           ENDDO           ENDDO
83          ENDDO          ENDDO
84    
85  C-      2D arrays  C-      2D arrays
86          DO J=1-Oly,sNy+Oly          DO J=1-OLy,sNy+OLy
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
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            hMixLayer(i,j,bi,bj) = 0. _d 0
92           ENDDO           ENDDO
93          ENDDO          ENDDO
94    
Legend:
Removed from v.1.19  
changed lines
  Added in v.1.27
  ViewVC Help
Powered by ViewVC 1.1.22