/[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.10 by jmc, Tue Apr 22 22:03:17 2008 UTC revision 1.11 by jmc, Mon Jan 11 19:49:28 2010 UTC
# Line 7  CBOP Line 7  CBOP
7  C     !ROUTINE: MOM_DIAGNOSTICS_INIT  C     !ROUTINE: MOM_DIAGNOSTICS_INIT
8  C     !INTERFACE:  C     !INTERFACE:
9        SUBROUTINE MOM_DIAGNOSTICS_INIT( myThid )        SUBROUTINE MOM_DIAGNOSTICS_INIT( myThid )
10    
11  C     !DESCRIPTION:  C     !DESCRIPTION:
12  C     Routine to initialize Momentum diagnostics  C     Routine to initialize Momentum diagnostics
13  C      (from mom_common, mom_fluxform & mom_vecinv packages)  C      (from mom_common, mom_fluxform & mom_vecinv packages)
# Line 20  C     === Global variables === Line 21  C     === Global variables ===
21    
22  C     !INPUT/OUTPUT PARAMETERS:  C     !INPUT/OUTPUT PARAMETERS:
23  C     === Routine arguments ===  C     === Routine arguments ===
24  C     myThid -  Number of this instance of GAD_INIT  C     myThid :: my Thread Id. number
25        INTEGER myThid        INTEGER myThid
26  CEOP  CEOP
27    
# Line 28  CEOP Line 29  CEOP
29  #ifdef ALLOW_MOM_COMMON  #ifdef ALLOW_MOM_COMMON
30  C     !LOCAL VARIABLES:  C     !LOCAL VARIABLES:
31  C     === Local variables ===  C     === Local variables ===
32  C     msgBuf      - Informational/error meesage buffer  C     msgBuf      :: Informational/error message buffer
33  c     CHARACTER*(MAX_LEN_MBUF) msgBuf  c     CHARACTER*(MAX_LEN_MBUF) msgBuf
34    
35        INTEGER       diagNum        INTEGER       diagNum
# Line 229  C--   add diagnostics for viscosity coef Line 230  C--   add diagnostics for viscosity coef
230  C--   add diagnostics of simple derivative quantities:  C--   add diagnostics of simple derivative quantities:
231        diagName  = 'momKE   '        diagName  = 'momKE   '
232        diagTitle = 'Kinetic Energy (in momentum Eq.)'        diagTitle = 'Kinetic Energy (in momentum Eq.)'
233        diagCode  = 'SM      MR      '        diagCode  = 'SMR     MR      '
234        diagUnits = 'm^2/s^2         '        diagUnits = 'm^2/s^2         '
235        CALL DIAGNOSTICS_ADDTOLIST( diagNum,        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
236       I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )       I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
# Line 237  C--   add diagnostics of simple derivati Line 238  C--   add diagnostics of simple derivati
238  #ifdef ALLOW_MOM_VECINV  #ifdef ALLOW_MOM_VECINV
239        diagName  = 'momHDiv '        diagName  = 'momHDiv '
240        diagTitle = 'Horizontal Divergence (in momentum Eq.)'        diagTitle = 'Horizontal Divergence (in momentum Eq.)'
241        diagCode  = 'SM      MR      '        diagCode  = 'SMR     MR      '
242        diagUnits = 's^-1            '        diagUnits = 's^-1            '
243        CALL DIAGNOSTICS_ADDTOLIST( diagNum,        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
244       I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )       I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
245    
246        diagName  = 'momVort3'        diagName  = 'momVort3'
247        diagTitle = '3rd component (vertical) of Vorticity'        diagTitle = '3rd component (vertical) of Vorticity'
248        diagCode  = 'SZ      MR      '        diagCode  = 'SZR     MR      '
249        diagUnits = 's^-1            '        diagUnits = 's^-1            '
250        CALL DIAGNOSTICS_ADDTOLIST( diagNum,        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
251       I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )       I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
252    
253        diagName  = 'Strain  '        diagName  = 'Strain  '
254        diagTitle = 'Horizontal Strain of Horizontal Velocities'        diagTitle = 'Horizontal Strain of Horizontal Velocities'
255        diagCode  = 'SZ      MR      '        diagCode  = 'SZR     MR      '
256        diagUnits = 's^-1            '        diagUnits = 's^-1            '
257        CALL DIAGNOSTICS_ADDTOLIST( diagNum,        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
258       I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )       I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
259    
260        diagName  = 'Tension '        diagName  = 'Tension '
261        diagTitle = 'Horizontal Tension of Horizontal Velocities'        diagTitle = 'Horizontal Tension of Horizontal Velocities'
262        diagCode  = 'SM      MR      '        diagCode  = 'SMR     MR      '
263        diagUnits = 's^-1            '        diagUnits = 's^-1            '
264        CALL DIAGNOSTICS_ADDTOLIST( diagNum,        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
265       I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )       I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
# Line 269  C--   add diagnostics of tendencies from Line 270  C--   add diagnostics of tendencies from
270    
271        diagName  = 'UBotDrag'        diagName  = 'UBotDrag'
272        diagTitle = 'U momentum tendency from Bottom Drag'        diagTitle = 'U momentum tendency from Bottom Drag'
273        diagCode  = 'UU      MR      '        diagCode  = 'UUR     MR      '
274        diagMate  = diagNum + 2        diagMate  = diagNum + 2
275        CALL DIAGNOSTICS_ADDTOLIST( diagNum,        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
276       I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )       I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
277        diagName  = 'VBotDrag'        diagName  = 'VBotDrag'
278        diagTitle = 'V momentum tendency from Bottom Drag'        diagTitle = 'V momentum tendency from Bottom Drag'
279        diagCode  = 'VV      MR      '        diagCode  = 'VVR     MR      '
280        diagMate  = diagNum        diagMate  = diagNum
281        CALL DIAGNOSTICS_ADDTOLIST( diagNum,        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
282       I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )       I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
283        diagName  = 'USidDrag'        diagName  = 'USidDrag'
284        diagTitle = 'U momentum tendency from Side Drag'        diagTitle = 'U momentum tendency from Side Drag'
285        diagCode  = 'UU      MR      '        diagCode  = 'UUR     MR      '
286        diagMate  = diagNum + 2        diagMate  = diagNum + 2
287        CALL DIAGNOSTICS_ADDTOLIST( diagNum,        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
288       I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )       I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
289        diagName  = 'VSidDrag'        diagName  = 'VSidDrag'
290        diagTitle = 'V momentum tendency from Side Drag'        diagTitle = 'V momentum tendency from Side Drag'
291        diagCode  = 'VV      MR      '        diagCode  = 'VVR     MR      '
292        diagMate  = diagNum        diagMate  = diagNum
293        CALL DIAGNOSTICS_ADDTOLIST( diagNum,        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
294       I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )       I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
295    
296        diagName  = 'Um_Diss '        diagName  = 'Um_Diss '
297        diagTitle = 'U momentum tendency from Dissipation'        diagTitle = 'U momentum tendency from Dissipation'
298        diagCode  = 'UU      MR      '        diagCode  = 'UUR     MR      '
299        diagMate  = diagNum + 2        diagMate  = diagNum + 2
300        CALL DIAGNOSTICS_ADDTOLIST( diagNum,        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
301       I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )       I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
302        diagName  = 'Vm_Diss '        diagName  = 'Vm_Diss '
303        diagTitle = 'V momentum tendency from Dissipation'        diagTitle = 'V momentum tendency from Dissipation'
304        diagCode  = 'VV      MR      '        diagCode  = 'VVR     MR      '
305        diagMate  = diagNum        diagMate  = diagNum
306        CALL DIAGNOSTICS_ADDTOLIST( diagNum,        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
307       I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )       I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
308    
309        diagName  = 'Um_Advec'        diagName  = 'Um_Advec'
310        diagTitle = 'U momentum tendency from Advection terms'        diagTitle = 'U momentum tendency from Advection terms'
311        diagCode  = 'UU      MR      '        diagCode  = 'UUR     MR      '
312        diagMate  = diagNum + 2        diagMate  = diagNum + 2
313        CALL DIAGNOSTICS_ADDTOLIST( diagNum,        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
314       I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )       I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
315        diagName  = 'Vm_Advec'        diagName  = 'Vm_Advec'
316        diagTitle = 'V momentum tendency from Advection terms'        diagTitle = 'V momentum tendency from Advection terms'
317        diagCode  = 'VV      MR      '        diagCode  = 'VVR     MR      '
318        diagMate  = diagNum        diagMate  = diagNum
319        CALL DIAGNOSTICS_ADDTOLIST( diagNum,        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
320       I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )       I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
321    
322        diagName  = 'Um_Cori '        diagName  = 'Um_Cori '
323        diagTitle = 'U momentum tendency from Coriolis term'        diagTitle = 'U momentum tendency from Coriolis term'
324        diagCode  = 'UU      MR      '        diagCode  = 'UUR     MR      '
325        diagMate  = diagNum + 2        diagMate  = diagNum + 2
326        CALL DIAGNOSTICS_ADDTOLIST( diagNum,        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
327       I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )       I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
328        diagName  = 'Vm_Cori '        diagName  = 'Vm_Cori '
329        diagTitle = 'V momentum tendency from Coriolis term'        diagTitle = 'V momentum tendency from Coriolis term'
330        diagCode  = 'VV      MR      '        diagCode  = 'VVR     MR      '
331        diagMate  = diagNum        diagMate  = diagNum
332        CALL DIAGNOSTICS_ADDTOLIST( diagNum,        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
333       I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )       I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
334    
335        diagName  = 'Um_Ext  '        diagName  = 'Um_Ext  '
336        diagTitle = 'U momentum tendency from external forcing'        diagTitle = 'U momentum tendency from external forcing'
337        diagCode  = 'UU      MR      '        diagCode  = 'UUR     MR      '
338        diagMate  = diagNum + 2        diagMate  = diagNum + 2
339        CALL DIAGNOSTICS_ADDTOLIST( diagNum,        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
340       I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )       I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
341        diagName  = 'Vm_Ext  '        diagName  = 'Vm_Ext  '
342        diagTitle = 'V momentum tendency from external forcing'        diagTitle = 'V momentum tendency from external forcing'
343        diagCode  = 'VV      MR      '        diagCode  = 'VVR     MR      '
344        diagMate  = diagNum        diagMate  = diagNum
345        CALL DIAGNOSTICS_ADDTOLIST( diagNum,        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
346       I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )       I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
# Line 347  C--   add diagnostics of tendencies from Line 348  C--   add diagnostics of tendencies from
348  #ifdef ALLOW_MOM_VECINV  #ifdef ALLOW_MOM_VECINV
349        diagName  = 'Um_AdvZ3'        diagName  = 'Um_AdvZ3'
350        diagTitle = 'U momentum tendency from Vorticity Advection'        diagTitle = 'U momentum tendency from Vorticity Advection'
351        diagCode  = 'UU      MR      '        diagCode  = 'UUR     MR      '
352        diagMate  = diagNum + 2        diagMate  = diagNum + 2
353        CALL DIAGNOSTICS_ADDTOLIST( diagNum,        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
354       I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )       I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
355        diagName  = 'Vm_AdvZ3'        diagName  = 'Vm_AdvZ3'
356        diagTitle = 'V momentum tendency from Vorticity Advection'        diagTitle = 'V momentum tendency from Vorticity Advection'
357        diagCode  = 'VV      MR      '        diagCode  = 'VVR     MR      '
358        diagMate  = diagNum        diagMate  = diagNum
359        CALL DIAGNOSTICS_ADDTOLIST( diagNum,        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
360       I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )       I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
# Line 361  C--   add diagnostics of tendencies from Line 362  C--   add diagnostics of tendencies from
362        diagName  = 'Um_AdvRe'        diagName  = 'Um_AdvRe'
363        diagTitle = 'U momentum tendency from vertical Advection'        diagTitle = 'U momentum tendency from vertical Advection'
364       &          //' (Explicit part)'       &          //' (Explicit part)'
365        diagCode  = 'UU      MR      '        diagCode  = 'UUR     MR      '
366        diagMate  = diagNum + 2        diagMate  = diagNum + 2
367        CALL DIAGNOSTICS_ADDTOLIST( diagNum,        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
368       I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )       I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
369        diagName  = 'Vm_AdvRe'        diagName  = 'Vm_AdvRe'
370        diagTitle = 'V momentum tendency from vertical Advection'        diagTitle = 'V momentum tendency from vertical Advection'
371       &          //' (Explicit part)'       &          //' (Explicit part)'
372        diagCode  = 'VV      MR      '        diagCode  = 'VVR     MR      '
373        diagMate  = diagNum        diagMate  = diagNum
374        CALL DIAGNOSTICS_ADDTOLIST( diagNum,        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
375       I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )       I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
# Line 384  C-    vertical momentum tendencies Line 385  C-    vertical momentum tendencies
385        ENDIF        ENDIF
386        diagName  = 'Wm_Diss '        diagName  = 'Wm_Diss '
387        diagTitle = 'W momentum tendency from Dissipation'        diagTitle = 'W momentum tendency from Dissipation'
388        diagCode  = 'WM      LR      '        diagCode  = 'WMr     LR      '
389        CALL DIAGNOSTICS_ADDTOLIST( diagNum,        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
390       I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )       I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
391    
392        diagName  = 'Wm_Advec'        diagName  = 'Wm_Advec'
393        diagTitle = 'W momentum tendency from Advection terms'        diagTitle = 'W momentum tendency from Advection terms'
394        diagCode  = 'WM      LR      '        diagCode  = 'WMr     LR      '
395        CALL DIAGNOSTICS_ADDTOLIST( diagNum,        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
396       I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )       I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
397    
398        diagName  = 'WSidDrag'        diagName  = 'WSidDrag'
399        diagTitle = 'Vertical momentum tendency from Side Drag'        diagTitle = 'Vertical momentum tendency from Side Drag'
400        diagCode  = 'WM      LR      '        diagCode  = 'WMr     LR      '
401        CALL DIAGNOSTICS_ADDTOLIST( diagNum,        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
402       I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )       I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
403  #endif /* ALLOW_NONHYDROSTATIC */  #endif /* ALLOW_NONHYDROSTATIC */
Legend:
Removed from v.1.10  
changed lines
  Added in v.1.11
  ViewVC Help
Powered by ViewVC 1.1.22