/[MITgcm]/MITgcm/pkg/mom_common/mom_diagnostics_init.F
ViewVC logotype

Diff of /MITgcm/pkg/mom_common/mom_diagnostics_init.F

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

revision 1.1 by jmc, Sun Sep 4 19:22:48 2005 UTC revision 1.2 by jmc, Sun Sep 4 19:35:16 2005 UTC
# Line 47  C--    add diagnostics of simple derivat Line 47  C--    add diagnostics of simple derivat
47         diagCode  = 'SM      MR      '         diagCode  = 'SM      MR      '
48         diagUnits = 'm^2/s^2         '         diagUnits = 'm^2/s^2         '
49         CALL DIAGNOSTICS_ADD2LIST( diagNum,         CALL DIAGNOSTICS_ADD2LIST( diagNum,
50       I              diagName, diagCode, diagUnits, diagTitle, myThid )       I              diagName, diagCode, diagUnits, diagTitle, myThid )
51    
52  #ifdef ALLOW_MOM_VECINV  #ifdef ALLOW_MOM_VECINV
53         diagName  = 'momHDiv '         diagName  = 'momHDiv '
# Line 55  C--    add diagnostics of simple derivat Line 55  C--    add diagnostics of simple derivat
55         diagCode  = 'SM      MR      '         diagCode  = 'SM      MR      '
56         diagUnits = 's^-1            '         diagUnits = 's^-1            '
57         CALL DIAGNOSTICS_ADD2LIST( diagNum,         CALL DIAGNOSTICS_ADD2LIST( diagNum,
58       I              diagName, diagCode, diagUnits, diagTitle, myThid )       I              diagName, diagCode, diagUnits, diagTitle, myThid )
59    
60         diagName  = 'momVort3'         diagName  = 'momVort3'
61         diagTitle = '3rd component (vertical) of Vorticity'         diagTitle = '3rd component (vertical) of Vorticity'
62         diagCode  = 'SZ      MR      '         diagCode  = 'SZ      MR      '
63         diagUnits = 's^-1            '         diagUnits = 's^-1            '
64         CALL DIAGNOSTICS_ADD2LIST( diagNum,         CALL DIAGNOSTICS_ADD2LIST( diagNum,
65       I              diagName, diagCode, diagUnits, diagTitle, myThid )       I              diagName, diagCode, diagUnits, diagTitle, myThid )
66  #endif /* ALLOW_MOM_VECINV */  #endif /* ALLOW_MOM_VECINV */
67    
68  C--    add diagnostics of tendencies from each terms  C--    add diagnostics of tendencies from each terms
# Line 72  C--    add diagnostics of tendencies fro Line 72  C--    add diagnostics of tendencies fro
72         diagTitle = 'U momentum tendency from Dissipation'         diagTitle = 'U momentum tendency from Dissipation'
73         WRITE(diagCode,'(A,I3.3,A)') 'UU   ',diagNum+2,'MR      '         WRITE(diagCode,'(A,I3.3,A)') 'UU   ',diagNum+2,'MR      '
74         CALL DIAGNOSTICS_ADD2LIST( diagNum,         CALL DIAGNOSTICS_ADD2LIST( diagNum,
75       I              diagName, diagCode, diagUnits, diagTitle, myThid )       I              diagName, diagCode, diagUnits, diagTitle, myThid )
76         diagName  = 'Vm_Diss '         diagName  = 'Vm_Diss '
77         diagTitle = 'V momentum tendency from Dissipation'         diagTitle = 'V momentum tendency from Dissipation'
78         WRITE(diagCode,'(A,I3.3,A)') 'VV   ', diagNum, 'MR      '         WRITE(diagCode,'(A,I3.3,A)') 'VV   ', diagNum, 'MR      '
79         CALL DIAGNOSTICS_ADD2LIST( diagNum,         CALL DIAGNOSTICS_ADD2LIST( diagNum,
80       I              diagName, diagCode, diagUnits, diagTitle, myThid )       I              diagName, diagCode, diagUnits, diagTitle, myThid )
81    
82         diagName  = 'Um_Advec'         diagName  = 'Um_Advec'
83         diagTitle = 'U momentum tendency from Advection terms'         diagTitle = 'U momentum tendency from Advection terms'
84         WRITE(diagCode,'(A,I3.3,A)') 'UU   ',diagNum+2,'MR      '         WRITE(diagCode,'(A,I3.3,A)') 'UU   ',diagNum+2,'MR      '
85         CALL DIAGNOSTICS_ADD2LIST( diagNum,         CALL DIAGNOSTICS_ADD2LIST( diagNum,
86       I              diagName, diagCode, diagUnits, diagTitle, myThid )       I              diagName, diagCode, diagUnits, diagTitle, myThid )
87         diagName  = 'Vm_Advec'         diagName  = 'Vm_Advec'
88         diagTitle = 'V momentum tendency from Advection terms'         diagTitle = 'V momentum tendency from Advection terms'
89         WRITE(diagCode,'(A,I3.3,A)') 'VV   ', diagNum, 'MR      '         WRITE(diagCode,'(A,I3.3,A)') 'VV   ', diagNum, 'MR      '
90         CALL DIAGNOSTICS_ADD2LIST( diagNum,         CALL DIAGNOSTICS_ADD2LIST( diagNum,
91       I              diagName, diagCode, diagUnits, diagTitle, myThid )       I              diagName, diagCode, diagUnits, diagTitle, myThid )
92    
93         diagName  = 'Um_Cori '         diagName  = 'Um_Cori '
94         diagTitle = 'U momentum tendency from Coriolis term'         diagTitle = 'U momentum tendency from Coriolis term'
95         WRITE(diagCode,'(A,I3.3,A)') 'UU   ',diagNum+2,'MR      '         WRITE(diagCode,'(A,I3.3,A)') 'UU   ',diagNum+2,'MR      '
96         CALL DIAGNOSTICS_ADD2LIST( diagNum,         CALL DIAGNOSTICS_ADD2LIST( diagNum,
97       I              diagName, diagCode, diagUnits, diagTitle, myThid )       I              diagName, diagCode, diagUnits, diagTitle, myThid )
98         diagName  = 'Vm_Cori '         diagName  = 'Vm_Cori '
99         diagTitle = 'V momentum tendency from Coriolis term'         diagTitle = 'V momentum tendency from Coriolis term'
100         WRITE(diagCode,'(A,I3.3,A)') 'VV   ', diagNum, 'MR      '         WRITE(diagCode,'(A,I3.3,A)') 'VV   ', diagNum, 'MR      '
101         CALL DIAGNOSTICS_ADD2LIST( diagNum,         CALL DIAGNOSTICS_ADD2LIST( diagNum,
102       I              diagName, diagCode, diagUnits, diagTitle, myThid )       I              diagName, diagCode, diagUnits, diagTitle, myThid )
103    
104  #ifdef ALLOW_MOM_VECINV  #ifdef ALLOW_MOM_VECINV
105         diagName  = 'Um_AdvZ3'         diagName  = 'Um_AdvZ3'
106         diagTitle = 'U momentum tendency from Vorticity Advection'         diagTitle = 'U momentum tendency from Vorticity Advection'
107         WRITE(diagCode,'(A,I3.3,A)') 'UU   ',diagNum+2,'MR      '         WRITE(diagCode,'(A,I3.3,A)') 'UU   ',diagNum+2,'MR      '
108         CALL DIAGNOSTICS_ADD2LIST( diagNum,         CALL DIAGNOSTICS_ADD2LIST( diagNum,
109       I              diagName, diagCode, diagUnits, diagTitle, myThid )       I              diagName, diagCode, diagUnits, diagTitle, myThid )
110         diagName  = 'Vm_AdvZ3'         diagName  = 'Vm_AdvZ3'
111         diagTitle = 'V momentum tendency from Vorticity Advection'         diagTitle = 'V momentum tendency from Vorticity Advection'
112         WRITE(diagCode,'(A,I3.3,A)') 'VV   ', diagNum, 'MR      '         WRITE(diagCode,'(A,I3.3,A)') 'VV   ', diagNum, 'MR      '
113         CALL DIAGNOSTICS_ADD2LIST( diagNum,         CALL DIAGNOSTICS_ADD2LIST( diagNum,
114       I              diagName, diagCode, diagUnits, diagTitle, myThid )       I              diagName, diagCode, diagUnits, diagTitle, myThid )
115    
116         diagName  = 'Um_AdvRe'         diagName  = 'Um_AdvRe'
117         diagTitle = 'U momentum tendency from vertical Advection'         diagTitle = 'U momentum tendency from vertical Advection'
118       &           //' (Explicit part)'       &           //' (Explicit part)'
119         WRITE(diagCode,'(A,I3.3,A)') 'UU   ',diagNum+2,'MR      '         WRITE(diagCode,'(A,I3.3,A)') 'UU   ',diagNum+2,'MR      '
120         CALL DIAGNOSTICS_ADD2LIST( diagNum,         CALL DIAGNOSTICS_ADD2LIST( diagNum,
121       I              diagName, diagCode, diagUnits, diagTitle, myThid )       I              diagName, diagCode, diagUnits, diagTitle, myThid )
122         diagName  = 'Vm_AdvRe'         diagName  = 'Vm_AdvRe'
123         diagTitle = 'V momentum tendency from vertical Advection'         diagTitle = 'V momentum tendency from vertical Advection'
124       &           //' (Explicit part)'       &           //' (Explicit part)'
125         WRITE(diagCode,'(A,I3.3,A)') 'VV   ', diagNum, 'MR      '         WRITE(diagCode,'(A,I3.3,A)') 'VV   ', diagNum, 'MR      '
126         CALL DIAGNOSTICS_ADD2LIST( diagNum,         CALL DIAGNOSTICS_ADD2LIST( diagNum,
127       I              diagName, diagCode, diagUnits, diagTitle, myThid )       I              diagName, diagCode, diagUnits, diagTitle, myThid )
128  #endif /* ALLOW_MOM_VECINV */  #endif /* ALLOW_MOM_VECINV */
129    
130  C--    add diagnostics of advective & viscous flux :  C--    add diagnostics of advective & viscous flux :
# Line 140  C-     Advective flux of Zonal momentum: Line 140  C-     Advective flux of Zonal momentum:
140         diagTitle = 'Zonal      Advective Flux of U momentum'         diagTitle = 'Zonal      Advective Flux of U momentum'
141         WRITE(diagCode,'(A,I3.3,A)') 'UM   ',diagNum+2,'MR      '         WRITE(diagCode,'(A,I3.3,A)') 'UM   ',diagNum+2,'MR      '
142         CALL DIAGNOSTICS_ADD2LIST( diagNum,         CALL DIAGNOSTICS_ADD2LIST( diagNum,
143       I              diagName, diagCode, diagUnits, diagTitle, myThid )       I              diagName, diagCode, diagUnits, diagTitle, myThid )
144         diagName  = 'ADVy_Um '         diagName  = 'ADVy_Um '
145         diagTitle = 'Meridional Advective Flux of U momentum'         diagTitle = 'Meridional Advective Flux of U momentum'
146         WRITE(diagCode,'(A,I3.3,A)') 'VZ   ', diagNum, 'MR      '         WRITE(diagCode,'(A,I3.3,A)') 'VZ   ', diagNum, 'MR      '
147         CALL DIAGNOSTICS_ADD2LIST( diagNum,         CALL DIAGNOSTICS_ADD2LIST( diagNum,
148       I              diagName, diagCode, diagUnits, diagTitle, myThid )       I              diagName, diagCode, diagUnits, diagTitle, myThid )
149         diagName  = 'ADVrE_Um'         diagName  = 'ADVrE_Um'
150         diagTitle = 'Vertical   Advective Flux of U momentum'         diagTitle = 'Vertical   Advective Flux of U momentum'
151       &           //' (Explicit part)'       &           //' (Explicit part)'
152         diagCode  = 'WU      LR      '         diagCode  = 'WU      LR      '
153         CALL DIAGNOSTICS_ADD2LIST( diagNum,         CALL DIAGNOSTICS_ADD2LIST( diagNum,
154       I              diagName, diagCode, diagUnits, diagTitle, myThid )       I              diagName, diagCode, diagUnits, diagTitle, myThid )
155  #endif /* ALLOW_MOM_FLUXFORM */  #endif /* ALLOW_MOM_FLUXFORM */
156  c      diagName  = 'ADVrI_Um'  c      diagName  = 'ADVrI_Um'
157  c      diagTitle = 'Vertical   Advective Flux of U momentum'  c      diagTitle = 'Vertical   Advective Flux of U momentum'
158  c    &           //' (Implicit part)'  c    &           //' (Implicit part)'
159  c      diagCode  = 'WU      LR      '  c      diagCode  = 'WU      LR      '
160  c      CALL DIAGNOSTICS_ADD2LIST( diagNum,  c      CALL DIAGNOSTICS_ADD2LIST( diagNum,
161  c    I              diagName, diagCode, diagUnits, diagTitle, myThid )  c    I              diagName, diagCode, diagUnits, diagTitle, myThid )
162    
163  C-     Advective flux of Meridional momentum:  C-     Advective flux of Meridional momentum:
164  #ifdef ALLOW_MOM_FLUXFORM  #ifdef ALLOW_MOM_FLUXFORM
# Line 166  C-     Advective flux of Meridional mome Line 166  C-     Advective flux of Meridional mome
166         diagTitle = 'Zonal      Advective Flux of V momentum'         diagTitle = 'Zonal      Advective Flux of V momentum'
167         WRITE(diagCode,'(A,I3.3,A)') 'UZ   ',diagNum+2,'MR      '         WRITE(diagCode,'(A,I3.3,A)') 'UZ   ',diagNum+2,'MR      '
168         CALL DIAGNOSTICS_ADD2LIST( diagNum,         CALL DIAGNOSTICS_ADD2LIST( diagNum,
169       I              diagName, diagCode, diagUnits, diagTitle, myThid )       I              diagName, diagCode, diagUnits, diagTitle, myThid )
170         diagName  = 'ADVy_Vm '         diagName  = 'ADVy_Vm '
171         diagTitle = 'Meridional Advective Flux of V momentum'         diagTitle = 'Meridional Advective Flux of V momentum'
172         WRITE(diagCode,'(A,I3.3,A)') 'VM   ', diagNum, 'MR      '         WRITE(diagCode,'(A,I3.3,A)') 'VM   ', diagNum, 'MR      '
173         CALL DIAGNOSTICS_ADD2LIST( diagNum,         CALL DIAGNOSTICS_ADD2LIST( diagNum,
174       I              diagName, diagCode, diagUnits, diagTitle, myThid )       I              diagName, diagCode, diagUnits, diagTitle, myThid )
175         diagName  = 'ADVrE_Vm'         diagName  = 'ADVrE_Vm'
176         diagTitle = 'Vertical   Advective Flux of V momentum'         diagTitle = 'Vertical   Advective Flux of V momentum'
177       &           //' (Explicit part)'       &           //' (Explicit part)'
178         diagCode  = 'WV      LR      '         diagCode  = 'WV      LR      '
179         CALL DIAGNOSTICS_ADD2LIST( diagNum,         CALL DIAGNOSTICS_ADD2LIST( diagNum,
180       I              diagName, diagCode, diagUnits, diagTitle, myThid )       I              diagName, diagCode, diagUnits, diagTitle, myThid )
181  #endif /* ALLOW_MOM_FLUXFORM */  #endif /* ALLOW_MOM_FLUXFORM */
182  c      diagName  = 'ADVrI_Vm'  c      diagName  = 'ADVrI_Vm'
183  c      diagTitle = 'Vertical   Advective Flux of V momentum'  c      diagTitle = 'Vertical   Advective Flux of V momentum'
184  c    &           //' (Implicit part)'  c    &           //' (Implicit part)'
185  c      diagCode  = 'WV      LR      '  c      diagCode  = 'WV      LR      '
186  c      CALL DIAGNOSTICS_ADD2LIST( diagNum,  c      CALL DIAGNOSTICS_ADD2LIST( diagNum,
187  c    I              diagName, diagCode, diagUnits, diagTitle, myThid )  c    I              diagName, diagCode, diagUnits, diagTitle, myThid )
188    
189  c-     Viscous flux of Zonal momentum :  c-     Viscous flux of Zonal momentum :
190  #ifdef ALLOW_MOM_FLUXFORM  #ifdef ALLOW_MOM_FLUXFORM
# Line 192  c-     Viscous flux of Zonal momentum : Line 192  c-     Viscous flux of Zonal momentum :
192         diagTitle = 'Zonal      Viscous Flux of U momentum'         diagTitle = 'Zonal      Viscous Flux of U momentum'
193         WRITE(diagCode,'(A,I3.3,A)') 'UM   ',diagNum+2,'MR      '         WRITE(diagCode,'(A,I3.3,A)') 'UM   ',diagNum+2,'MR      '
194         CALL DIAGNOSTICS_ADD2LIST( diagNum,         CALL DIAGNOSTICS_ADD2LIST( diagNum,
195       I              diagName, diagCode, diagUnits, diagTitle, myThid )       I              diagName, diagCode, diagUnits, diagTitle, myThid )
196         diagName  = 'VISCy_Um'         diagName  = 'VISCy_Um'
197         diagTitle = 'Meridional Viscous Flux of U momentum'         diagTitle = 'Meridional Viscous Flux of U momentum'
198         WRITE(diagCode,'(A,I3.3,A)') 'VZ   ', diagNum, 'MR      '         WRITE(diagCode,'(A,I3.3,A)') 'VZ   ', diagNum, 'MR      '
199         CALL DIAGNOSTICS_ADD2LIST( diagNum,         CALL DIAGNOSTICS_ADD2LIST( diagNum,
200       I              diagName, diagCode, diagUnits, diagTitle, myThid )       I              diagName, diagCode, diagUnits, diagTitle, myThid )
201         diagName  = 'VISrE_Um'         diagName  = 'VISrE_Um'
202         diagTitle = 'Vertical   Viscous Flux of U momentum'         diagTitle = 'Vertical   Viscous Flux of U momentum'
203       &           //' (Explicit part)'       &           //' (Explicit part)'
204         diagCode  = 'WU      LR      '         diagCode  = 'WU      LR      '
205         CALL DIAGNOSTICS_ADD2LIST( diagNum,         CALL DIAGNOSTICS_ADD2LIST( diagNum,
206       I              diagName, diagCode, diagUnits, diagTitle, myThid )       I              diagName, diagCode, diagUnits, diagTitle, myThid )
207  #endif /* ALLOW_MOM_FLUXFORM */  #endif /* ALLOW_MOM_FLUXFORM */
208         diagName  = 'VISrI_Um'         diagName  = 'VISrI_Um'
209         diagTitle = 'Vertical   Viscous Flux of U momentum'         diagTitle = 'Vertical   Viscous Flux of U momentum'
210       &           //' (Implicit part)'       &           //' (Implicit part)'
211         diagCode  = 'WU      LR      '         diagCode  = 'WU      LR      '
212         CALL DIAGNOSTICS_ADD2LIST( diagNum,         CALL DIAGNOSTICS_ADD2LIST( diagNum,
213       I              diagName, diagCode, diagUnits, diagTitle, myThid )       I              diagName, diagCode, diagUnits, diagTitle, myThid )
214    
215  c-     Viscous flux of Meridional momentum :  c-     Viscous flux of Meridional momentum :
216  #ifdef ALLOW_MOM_FLUXFORM  #ifdef ALLOW_MOM_FLUXFORM
# Line 218  c-     Viscous flux of Meridional moment Line 218  c-     Viscous flux of Meridional moment
218         diagTitle = 'Zonal      Viscous Flux of V momentum'         diagTitle = 'Zonal      Viscous Flux of V momentum'
219         WRITE(diagCode,'(A,I3.3,A)') 'UZ   ',diagNum+2,'MR      '         WRITE(diagCode,'(A,I3.3,A)') 'UZ   ',diagNum+2,'MR      '
220         CALL DIAGNOSTICS_ADD2LIST( diagNum,         CALL DIAGNOSTICS_ADD2LIST( diagNum,
221       I              diagName, diagCode, diagUnits, diagTitle, myThid )       I              diagName, diagCode, diagUnits, diagTitle, myThid )
222         diagName  = 'VISCy_Vm'         diagName  = 'VISCy_Vm'
223         diagTitle = 'Meridional Viscous Flux of V momentum'         diagTitle = 'Meridional Viscous Flux of V momentum'
224         WRITE(diagCode,'(A,I3.3,A)') 'VM   ', diagNum, 'MR      '         WRITE(diagCode,'(A,I3.3,A)') 'VM   ', diagNum, 'MR      '
225         CALL DIAGNOSTICS_ADD2LIST( diagNum,         CALL DIAGNOSTICS_ADD2LIST( diagNum,
226       I              diagName, diagCode, diagUnits, diagTitle, myThid )       I              diagName, diagCode, diagUnits, diagTitle, myThid )
227         diagName  = 'VISrE_Vm'         diagName  = 'VISrE_Vm'
228         diagTitle = 'Vertical   Viscous Flux of V momentum'         diagTitle = 'Vertical   Viscous Flux of V momentum'
229       &           //' (Explicit part)'       &           //' (Explicit part)'
230         diagCode  = 'WV      LR      '         diagCode  = 'WV      LR      '
231         CALL DIAGNOSTICS_ADD2LIST( diagNum,         CALL DIAGNOSTICS_ADD2LIST( diagNum,
232       I              diagName, diagCode, diagUnits, diagTitle, myThid )       I              diagName, diagCode, diagUnits, diagTitle, myThid )
233  #endif /* ALLOW_MOM_FLUXFORM */  #endif /* ALLOW_MOM_FLUXFORM */
234         diagName  = 'VISrI_Vm'         diagName  = 'VISrI_Vm'
235         diagTitle = 'Vertical   Viscous Flux of V momentum'         diagTitle = 'Vertical   Viscous Flux of V momentum'
236       &           //' (Implicit part)'       &           //' (Implicit part)'
237         diagCode  = 'WV      LR      '         diagCode  = 'WV      LR      '
238         CALL DIAGNOSTICS_ADD2LIST( diagNum,         CALL DIAGNOSTICS_ADD2LIST( diagNum,
239       I              diagName, diagCode, diagUnits, diagTitle, myThid )       I              diagName, diagCode, diagUnits, diagTitle, myThid )
240    
241  C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|  C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|
242    
Legend:
Removed from v.1.1  
changed lines
  Added in v.1.2
  ViewVC Help
Powered by ViewVC 1.1.22