/[MITgcm]/MITgcm/pkg/gmredi/gmredi_init.F
ViewVC logotype

Diff of /MITgcm/pkg/gmredi/gmredi_init.F

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

revision 1.5 by jmc, Tue Mar 6 18:01:07 2001 UTC revision 1.6 by jmc, Sun Dec 16 18:54:49 2001 UTC
# Line 31  C     msgBuf      - Informational/error Line 31  C     msgBuf      - Informational/error
31  C     iUnit       - Work variable for IO unit number  C     iUnit       - Work variable for IO unit number
32        CHARACTER*(MAX_LEN_MBUF) msgBuf        CHARACTER*(MAX_LEN_MBUF) msgBuf
33        INTEGER iUnit        INTEGER iUnit
34        INTEGER k,bi,bj        INTEGER i,j,k,bi,bj
35    
 #ifdef ALLOW_TIMEAVE  
 C     Initialize averages to zero  
36        DO bj = myByLo(myThid), myByHi(myThid)        DO bj = myByLo(myThid), myByHi(myThid)
37         DO bi = myBxLo(myThid), myBxHi(myThid)         DO bi = myBxLo(myThid), myBxHi(myThid)
38    
39    C     Initialize arrays in common blocks :
40            DO k=1,Nr
41             DO j=1-Oly,sNy+OLy
42              DO i=1-Olx,sNx+Olx
43               Kwx(i,j,k,bi,bj) = 0.
44               Kwy(i,j,k,bi,bj) = 0.
45               Kwz(i,j,k,bi,bj) = 0.
46    #ifdef GM_EXTRA_DIAGONAL
47               Kuz(i,j,k,bi,bj) = 0.
48               Kvz(i,j,k,bi,bj) = 0.
49    #endif
50    #ifdef GM_NON_UNITY_DIAGONAL
51               Kux(i,j,k,bi,bj) = 0.
52               Kvy(i,j,k,bi,bj) = 0.
53    #endif
54    #ifdef GM_BOLUS_ADVEC
55               GM_PsiX(i,j,k,bi,bj) = 0.
56               GM_PsiY(i,j,k,bi,bj) = 0.
57    #endif
58    #ifdef GM_VISBECK_VARIABLE_K
59              VisbeckK(i,j,bi,bj) = 0.
60    #endif
61              ENDDO
62             ENDDO
63            ENDDO
64    
65    #ifdef ALLOW_TIMEAVE
66    C     Initialize averages to zero
67          CALL TIMEAVE_RESET(GM_Kwx_T,Nr, bi,bj,myThid)          CALL TIMEAVE_RESET(GM_Kwx_T,Nr, bi,bj,myThid)
68          CALL TIMEAVE_RESET(GM_Kwy_T,Nr, bi,bj,myThid)          CALL TIMEAVE_RESET(GM_Kwy_T,Nr, bi,bj,myThid)
69          CALL TIMEAVE_RESET(GM_Kwz_T,Nr, bi,bj,myThid)          CALL TIMEAVE_RESET(GM_Kwz_T,Nr, bi,bj,myThid)
# Line 46  C     Initialize averages to zero Line 73  C     Initialize averages to zero
73  #ifdef GM_VISBECK_VARIABLE_K  #ifdef GM_VISBECK_VARIABLE_K
74          CALL TIMEAVE_RESET(Visbeck_K_T, 1, bi,bj,myThid)          CALL TIMEAVE_RESET(Visbeck_K_T, 1, bi,bj,myThid)
75  #endif  #endif
76         ENDDO  #ifdef GM_BOLUS_ADVEC
77        ENDDO          CALL TIMEAVE_RESET(GM_PsiXtave,Nr, bi,bj,myThid)
78            CALL TIMEAVE_RESET(GM_PsiYtave,Nr, bi,bj,myThid)
79    #endif
80  #endif /* ALLOW_TIMEAVE */  #endif /* ALLOW_TIMEAVE */
81    
82    C- end bi,bj loops
83           ENDDO
84          ENDDO
85  #endif /* ALLOW_GMREDI */  #endif /* ALLOW_GMREDI */
86    
87        return        return
Legend:
Removed from v.1.5  
changed lines
  Added in v.1.6
  ViewVC Help
Powered by ViewVC 1.1.22