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

Diff of /MITgcm/pkg/gmredi/gmredi_diags.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 37  C     == Local variables == Line 37  C     == Local variables ==
37        CHARACTER*(MAX_LEN_MBUF) suff        CHARACTER*(MAX_LEN_MBUF) suff
38    
39    
40  #ifdef ALLOW_TIMEAVE  #if ( defined (GM_NON_UNITY_DIAGONAL) || defined (GM_EXTRA_DIAGONAL) )
41    
42          IF ( DIFFERENT_MULTIPLE(diagFreq,myTime,myTime-deltaTClock)
43         &   ) THEN
44            WRITE(suff,'(I10.10)') myIter
45    #ifdef GM_NON_UNITY_DIAGONAL
46            CALL WRITE_FLD_XYZ_RL( 'GM_Kux.',suff,Kux,myIter,myThid)
47            CALL WRITE_FLD_XYZ_RL( 'GM_Kvy.',suff,Kvy,myIter,myThid)
48    #endif
49    #ifdef GM_EXTRA_DIAGONAL
50           IF (GM_ExtraDiag) THEN
51            CALL WRITE_FLD_XYZ_RL( 'GM_Kuz.',suff,Kuz,myIter,myThid)
52            CALL WRITE_FLD_XYZ_RL( 'GM_Kvz.',suff,Kvz,myIter,myThid)
53           ENDIF
54    #endif
55          ENDIF
56    
57    #endif /* GM_NON_UNITY_DIAGONAL || GM_EXTRA_DIAGONAL */
58    
59    
60    #ifdef ALLOW_TIMEAVE
61    
62  C     Dump files and restart average computation if needed  C     Dump files and restart average computation if needed
63        IF ( DIFFERENT_MULTIPLE(taveFreq,myTime,myTime-deltaTClock)        IF ( DIFFERENT_MULTIPLE(taveFreq,myTime,myTime-deltaTClock)
64       &     ) THEN       &     ) THEN
# Line 54  C      Normalize by integrated time Line 73  C      Normalize by integrated time
73           CALL TIMEAVE_NORMALIZ(GM_Kwz_T,GM_timeave, Nr,           CALL TIMEAVE_NORMALIZ(GM_Kwz_T,GM_timeave, Nr,
74       &                 bi,bj,myThid)       &                 bi,bj,myThid)
75  #ifdef GM_VISBECK_VARIABLE_K  #ifdef GM_VISBECK_VARIABLE_K
76           CALL TIMEAVE_NORMALIZ(Visbeck_K_T,GM_timeave, 1,          IF (GM_Visbeck_alpha.NE.0.)
77         &   CALL TIMEAVE_NORMALIZ(Visbeck_K_T,GM_timeave, 1,
78         &                 bi,bj,myThid)
79    #endif
80    #ifdef GM_BOLUS_ADVEC
81            IF (GM_AdvForm) THEN
82             CALL TIMEAVE_NORMALIZ(GM_PsiXtave,GM_timeave, Nr,
83       &                 bi,bj,myThid)       &                 bi,bj,myThid)
84             CALL TIMEAVE_NORMALIZ(GM_PsiYtave,GM_timeave, Nr,
85         &                 bi,bj,myThid)
86            ENDIF
87  #endif  #endif
88          ENDDO          ENDDO
89         ENDDO         ENDDO
# Line 66  C      Normalize by integrated time Line 94  C      Normalize by integrated time
94         CALL WRITE_FLD_XYZ_RL('GM_Kwy-T.',suff,GM_Kwy_T,myIter,myThid)         CALL WRITE_FLD_XYZ_RL('GM_Kwy-T.',suff,GM_Kwy_T,myIter,myThid)
95         CALL WRITE_FLD_XYZ_RL('GM_Kwz-T.',suff,GM_Kwz_T,myIter,myThid)         CALL WRITE_FLD_XYZ_RL('GM_Kwz-T.',suff,GM_Kwz_T,myIter,myThid)
96  #ifdef GM_VISBECK_VARIABLE_K  #ifdef GM_VISBECK_VARIABLE_K
97         CALL WRITE_FLD_XY_RL('Visbeck_K-T.',suff,Visbeck_K_T,         IF (GM_Visbeck_alpha.NE.0.)
98         &  CALL WRITE_FLD_XY_RL('Visbeck_K-T.',suff,Visbeck_K_T,
99       &                        myIter,myThid)       &                        myIter,myThid)
100  #endif  #endif
101    #ifdef GM_BOLUS_ADVEC
102           IF (GM_AdvForm) THEN
103            CALL WRITE_FLD_XYZ_RL('GM_PsiXtave.',suff,GM_PsiXtave,
104         &                        myIter,myThid)
105            CALL WRITE_FLD_XYZ_RL('GM_PsiYtave.',suff,GM_PsiYtave,
106         &                        myIter,myThid)
107           ENDIF
108    #endif
109                    
110  C      Reset averages to zero  C      Reset averages to zero
111         DO bj = myByLo(myThid), myByHi(myThid)         DO bj = myByLo(myThid), myByHi(myThid)
# Line 76  C      Reset averages to zero Line 113  C      Reset averages to zero
113           CALL TIMEAVE_RESET(GM_Kwx_T,Nr,bi,bj,myThid)           CALL TIMEAVE_RESET(GM_Kwx_T,Nr,bi,bj,myThid)
114           CALL TIMEAVE_RESET(GM_Kwy_T,Nr,bi,bj,myThid)           CALL TIMEAVE_RESET(GM_Kwy_T,Nr,bi,bj,myThid)
115           CALL TIMEAVE_RESET(GM_Kwz_T,Nr,bi,bj,myThid)           CALL TIMEAVE_RESET(GM_Kwz_T,Nr,bi,bj,myThid)
116    #ifdef GM_VISBECK_VARIABLE_K
117             IF (GM_Visbeck_alpha.NE.0.)
118         &    CALL TIMEAVE_RESET(Visbeck_K_T,1,bi,bj,myThid)
119    #endif
120    #ifdef GM_BOLUS_ADVEC
121             IF (GM_AdvForm) THEN
122              CALL TIMEAVE_RESET(GM_PsiXtave,Nr,bi,bj,myThid)
123              CALL TIMEAVE_RESET(GM_PsiYtave,Nr,bi,bj,myThid)
124             ENDIF
125    #endif
126           DO k=1,Nr           DO k=1,Nr
127            GM_TimeAve(k,bi,bj)=0.            GM_TimeAve(k,bi,bj)=0.
128           ENDDO           ENDDO
 #ifdef GM_VISBECK_VARIABLE_K  
          CALL TIMEAVE_RESET(Visbeck_K_T,1,bi,bj,myThid)  
 #endif  
129          ENDDO          ENDDO
130         ENDDO         ENDDO
131                    
Legend:
Removed from v.1.5  
changed lines
  Added in v.1.6
  ViewVC Help
Powered by ViewVC 1.1.22