/[MITgcm]/MITgcm/pkg/diagnostics/diagnostics_main_init.F
ViewVC logotype

Diff of /MITgcm/pkg/diagnostics/diagnostics_main_init.F

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

revision 1.38 by jmc, Tue Jun 14 00:18:37 2011 UTC revision 1.52 by jmc, Sat May 28 23:24:47 2016 UTC
# Line 76  C     sTitle     :: "salt" title Line 76  C     sTitle     :: "salt" title
76        CHARACTER*16   diagUnits        CHARACTER*16   diagUnits
77        CHARACTER*(80) diagTitle        CHARACTER*(80) diagTitle
78        CHARACTER*2    rUnit2c        CHARACTER*2    rUnit2c
79        CHARACTER*4    tUnit4c, sUnit4c        CHARACTER*4    tUnit4c
80          CHARACTER*5    sUnit5c
81        CHARACTER*(10) rTitle, eTitle, fTitle        CHARACTER*(10) rTitle, eTitle, fTitle
82        CHARACTER*(20) pTitle, sTitle        CHARACTER*(20) pTitle, sTitle
83    
# Line 102  C--------------------------------------- Line 103  C---------------------------------------
103        ENDIF        ENDIF
104        IF ( fluidIsAir ) THEN        IF ( fluidIsAir ) THEN
105          tUnit4c= 'K   '          tUnit4c= 'K   '
106          sUnit4c= 'g/kg'          sUnit5c= 'kg/kg'
107          sTitle = ' Specific Humidity  '          sTitle = ' Specific Humidity  '
108  #ifdef ALLOW_FIZHI          IF (useAIM) sUnit5c= 'g/kg '
109          IF (useFIZHI) sUnit4c= 'kg/kg'        ELSEIF ( eosType.EQ.'TEOS10' ) THEN
110  #endif /* ALLOW_FIZHI */          tUnit4c= 'degC'
111            sUnit5c= 'g/kg '
112    c       tTitle = 'Conservative Temp.   '
113            sTitle = ' Absolute Salinity  '
114        ELSE        ELSE
115          tUnit4c= 'degC'          tUnit4c= 'degC'
116          sUnit4c= 'psu '          sUnit5c= 'psu  '
117    c       tTitle = 'Potential Temperature'
118          sTitle = ' Salinity           '          sTitle = ' Salinity           '
119        ENDIF        ENDIF
120  C-    free-surface (eTitle) and fixed-boundary (fTitle) position:  C-    free-surface (eTitle) and fixed-boundary (fTitle) position:
# Line 127  C-    free-surface (eTitle) and fixed-bo Line 132  C-    free-surface (eTitle) and fixed-bo
132  C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|  C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|
133  C-    state variables of the main code (and related quadratic var):  C-    state variables of the main code (and related quadratic var):
134    
   
135        diagName  = 'ETAN    '        diagName  = 'ETAN    '
136        diagTitle = DIAGS_MK_TITLE( eTitle//rTitle//' Anomaly', myThid )        diagTitle = DIAGS_MK_TITLE( eTitle//rTitle//' Anomaly', myThid )
137  c     IF ( fluidIsWater .AND. usingZCoords )  c     IF ( fluidIsWater .AND. usingZCoords )
# Line 169  c    I          diagName, diagCode, diag Line 173  c    I          diagName, diagCode, diag
173    
174        diagName  = 'SALT    '        diagName  = 'SALT    '
175        diagTitle = DIAGS_MK_TITLE( sTitle,  myThid )        diagTitle = DIAGS_MK_TITLE( sTitle,  myThid )
176        diagUnits = DIAGS_MK_UNITS( sUnit4c, myThid )        diagUnits = DIAGS_MK_UNITS( sUnit5c, myThid )
177        diagCode  = 'SMR     MR      '        diagCode  = 'SMR     MR      '
178        CALL DIAGNOSTICS_ADDTOLIST( diagNum,        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
179       I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )       I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
# Line 183  c    I          diagName, diagCode, diag Line 187  c    I          diagName, diagCode, diag
187    
188  c     diagName  = 'SSS     '  c     diagName  = 'SSS     '
189  c     diagTitle = 'Sea Surface Salinity '  c     diagTitle = 'Sea Surface Salinity '
190  c     diagUnits = DIAGS_MK_UNITS( sUnit4c, myThid )  c     diagUnits = DIAGS_MK_UNITS( sUnit5c, myThid )
191  c     diagCode  = 'SM      M1      '  c     diagCode  = 'SM      M1      '
192  c     CALL DIAGNOSTICS_ADDTOLIST( diagNum,  c     CALL DIAGNOSTICS_ADDTOLIST( diagNum,
193  c    I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )  c    I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
# Line 191  c    I          diagName, diagCode, diag Line 195  c    I          diagName, diagCode, diag
195        IF ( fluidIsWater ) THEN        IF ( fluidIsWater ) THEN
196        diagName  = 'SALTanom'        diagName  = 'SALTanom'
197        diagTitle = 'Salt anomaly (=SALT-35; g/kg)'        diagTitle = 'Salt anomaly (=SALT-35; g/kg)'
198        diagUnits = DIAGS_MK_UNITS( sUnit4c, myThid )        diagUnits = DIAGS_MK_UNITS( sUnit5c, myThid )
199        diagCode  = 'SMR     MR      '        diagCode  = 'SMR     MR      '
200        CALL DIAGNOSTICS_ADDTOLIST( diagNum,        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
201       I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )       I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
# Line 229  c    I          diagName, diagCode, diag Line 233  c    I          diagName, diagCode, diag
233    
234        diagName  = 'SALTSQ  '        diagName  = 'SALTSQ  '
235        diagTitle = DIAGS_MK_TITLE( 'Square of '//sTitle, myThid )        diagTitle = DIAGS_MK_TITLE( 'Square of '//sTitle, myThid )
236        diagUnits = DIAGS_MK_UNITS( '('//sUnit4c//')^2', myThid )        diagUnits = DIAGS_MK_UNITS( '('//sUnit5c//')^2', myThid )
237        diagCode  = 'SMRP    MR      '        diagCode  = 'SMRP    MR      '
238        CALL DIAGNOSTICS_ADDTOLIST( diagNum,        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
239       I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )       I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
# Line 237  c    I          diagName, diagCode, diag Line 241  c    I          diagName, diagCode, diag
241        IF ( fluidIsWater ) THEN        IF ( fluidIsWater ) THEN
242        diagName  = 'SALTSQan'        diagName  = 'SALTSQan'
243        diagTitle = 'Square of Salt anomaly (=(SALT-35)^2 (g^2/kg^2)'        diagTitle = 'Square of Salt anomaly (=(SALT-35)^2 (g^2/kg^2)'
244        diagUnits = DIAGS_MK_UNITS( '('//sUnit4c//')^2', myThid )        diagUnits = DIAGS_MK_UNITS( '('//sUnit5c//')^2', myThid )
245        diagCode  = 'SMRP    MR      '        diagCode  = 'SMRP    MR      '
246        CALL DIAGNOSTICS_ADDTOLIST( diagNum,        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
247       I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )       I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
# Line 343  C-    use 'UVELMASS' as mate. Line 347  C-    use 'UVELMASS' as mate.
347        CALL DIAGNOSTICS_ADDTOLIST( diagNum,        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
348       I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )       I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
349    
350  c     diagName  = 'PsiVEL  '        diagName  = 'PsiVEL  '
351  c     diagTitle = 'Horizontal Velocity Stream-Function'        diagTitle = 'Horizontal Velocity Stream-Function'
352  c     diagUnits = 'm^2/s           '        diagUnits = DIAGS_MK_UNITS( rUnit2c//'.m^2/s', myThid )
353  c     diagUnits = DIAGS_MK_UNITS( rUnit2c//'.m^2/s', myThid )        diagCode  = 'SZ  P   MR      '
 c     diagCode  = 'SZ  P   MR      '  
354  C-    use 'PhiVEL' as mate.  C-    use 'PhiVEL' as mate.
355  c     diagMate  = diagNum        diagMate  = diagNum
356  c     CALL DIAGNOSTICS_ADDTOLIST( diagNum,        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
357  c    I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )       I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
358    
359        diagName  = 'UTHMASS '        diagName  = 'UTHMASS '
360        diagTitle = 'Zonal Mass-Weight Transp of Pot Temp'        diagTitle = 'Zonal Mass-Weight Transp of Pot Temp'
# Line 379  c    I   diagName, diagCode, diagUnits, Line 382  c    I   diagName, diagCode, diagUnits,
382        diagName  = 'USLTMASS'        diagName  = 'USLTMASS'
383        diagTitle = DIAGS_MK_TITLE( 'Zonal Mass-Weight Transp of '        diagTitle = DIAGS_MK_TITLE( 'Zonal Mass-Weight Transp of '
384       I                           //sTitle, myThid )       I                           //sTitle, myThid )
385        diagUnits = DIAGS_MK_UNITS(sUnit4c//'.m/s', myThid )        diagUnits = DIAGS_MK_UNITS(sUnit5c//'.m/s', myThid )
386        diagCode  = 'UUr     MR      '        diagCode  = 'UUr     MR      '
387        diagMate  = diagNum + 2        diagMate  = diagNum + 2
388        CALL DIAGNOSTICS_ADDTOLIST( diagNum,        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
# Line 388  c    I   diagName, diagCode, diagUnits, Line 391  c    I   diagName, diagCode, diagUnits,
391        diagName  = 'VSLTMASS'        diagName  = 'VSLTMASS'
392        diagTitle = DIAGS_MK_TITLE( 'Meridional Mass-Weight Transp of '        diagTitle = DIAGS_MK_TITLE( 'Meridional Mass-Weight Transp of '
393       I                           //sTitle, myThid )       I                           //sTitle, myThid )
394        diagUnits = DIAGS_MK_UNITS(sUnit4c//'.m/s', myThid )        diagUnits = DIAGS_MK_UNITS(sUnit5c//'.m/s', myThid )
395        diagCode  = 'VVr     MR      '        diagCode  = 'VVr     MR      '
396        diagMate  = diagNum        diagMate  = diagNum
397        CALL DIAGNOSTICS_ADDTOLIST( diagNum,        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
# Line 397  c    I   diagName, diagCode, diagUnits, Line 400  c    I   diagName, diagCode, diagUnits,
400        diagName  = 'WSLTMASS'        diagName  = 'WSLTMASS'
401        diagTitle = DIAGS_MK_TITLE( 'Vertical Mass-Weight Transp of '        diagTitle = DIAGS_MK_TITLE( 'Vertical Mass-Weight Transp of '
402       I                           //sTitle, myThid )       I                           //sTitle, myThid )
403        diagUnits = DIAGS_MK_UNITS(sUnit4c//'.'//rUnit2c//'/s', myThid )        diagUnits = DIAGS_MK_UNITS(sUnit5c//'.'//rUnit2c//'/s', myThid )
404        diagCode  = 'WM      LR      '        diagCode  = 'WM      LR      '
405        CALL DIAGNOSTICS_ADDTOLIST( diagNum,        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
406       I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )       I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
# Line 428  c    I   diagName, diagCode, diagUnits, Line 431  c    I   diagName, diagCode, diagUnits,
431        diagName  = 'UVELSLT '        diagName  = 'UVELSLT '
432        diagTitle = DIAGS_MK_TITLE( 'Zonal Transport of '        diagTitle = DIAGS_MK_TITLE( 'Zonal Transport of '
433       I                          //sTitle, myThid )       I                          //sTitle, myThid )
434        diagUnits = DIAGS_MK_UNITS( sUnit4c//'.m/s', myThid )        diagUnits = DIAGS_MK_UNITS( sUnit5c//'.m/s', myThid )
435        diagCode  = 'UUR     MR      '        diagCode  = 'UUR     MR      '
436        diagMate  = diagNum + 2        diagMate  = diagNum + 2
437        CALL DIAGNOSTICS_ADDTOLIST( diagNum,        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
# Line 437  c    I   diagName, diagCode, diagUnits, Line 440  c    I   diagName, diagCode, diagUnits,
440        diagName  = 'VVELSLT '        diagName  = 'VVELSLT '
441        diagTitle = DIAGS_MK_TITLE( 'Meridional Transport of '        diagTitle = DIAGS_MK_TITLE( 'Meridional Transport of '
442       I                          //sTitle, myThid )       I                          //sTitle, myThid )
443        diagUnits = DIAGS_MK_UNITS( sUnit4c//'.m/s', myThid )        diagUnits = DIAGS_MK_UNITS( sUnit5c//'.m/s', myThid )
444        diagCode  = 'VVR     MR      '        diagCode  = 'VVR     MR      '
445        diagMate  = diagNum        diagMate  = diagNum
446        CALL DIAGNOSTICS_ADDTOLIST( diagNum,        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
# Line 446  c    I   diagName, diagCode, diagUnits, Line 449  c    I   diagName, diagCode, diagUnits,
449        diagName  = 'WVELSLT '        diagName  = 'WVELSLT '
450        diagTitle = DIAGS_MK_TITLE( 'Vertical Transport of '        diagTitle = DIAGS_MK_TITLE( 'Vertical Transport of '
451       I                          //sTitle, myThid )       I                          //sTitle, myThid )
452        diagUnits = DIAGS_MK_UNITS(sUnit4c//'.'//rUnit2c//'/s', myThid )        diagUnits = DIAGS_MK_UNITS(sUnit5c//'.'//rUnit2c//'/s', myThid )
453        diagCode  = 'WM      LR      '        diagCode  = 'WM      LR      '
454        CALL DIAGNOSTICS_ADDTOLIST( diagNum,        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
455       I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )       I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
# Line 461  c    I   diagName, diagCode, diagUnits, Line 464  c    I   diagName, diagCode, diagUnits,
464       I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )       I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
465    
466        diagName  = 'VVELPHI '        diagName  = 'VVELPHI '
467        diagTitle = DIAGS_MK_TITLE( 'Meridional Mass-Weight Transp of '        diagTitle = DIAGS_MK_TITLE( 'Merid. Mass-Weight Transp of '
468       I                 //pTitle//' Anomaly', myThid )       I                 //pTitle//' Anomaly', myThid )
469        diagUnits = 'm^3/s^3         '        diagUnits = 'm^3/s^3         '
470        diagCode  = 'VVr     MR      '        diagCode  = 'VVr     MR      '
# Line 480  C---+----1----+----2----+----3----+----4 Line 483  C---+----1----+----2----+----3----+----4
483    
484        diagName  = 'RHOANOSQ'        diagName  = 'RHOANOSQ'
485        diagTitle = 'Square of Density Anomaly (=(Rho-rhoConst)^2)'        diagTitle = 'Square of Density Anomaly (=(Rho-rhoConst)^2)'
486        diagUnits = 'kg^2/m^6          '        diagUnits = 'kg^2/m^6        '
487        diagCode  = 'SMRP    MR      '        diagCode  = 'SMRP    MR      '
488        CALL DIAGNOSTICS_ADDTOLIST( diagNum,        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
489       I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )       I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
# Line 502  C---+----1----+----2----+----3----+----4 Line 505  C---+----1----+----2----+----3----+----4
505       I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )       I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
506    
507        diagName  = 'WRHOMASS'        diagName  = 'WRHOMASS'
508        diagTitle = 'Vertical Transport of Potential Density'        diagTitle = 'Vertical Transport of Density'
509          diagUnits = 'kg/m^2/s        '
510          diagCode  = 'WM      LR      '
511          CALL DIAGNOSTICS_ADDTOLIST( diagNum,
512         I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
513    
514          diagName  = 'WdRHO_P '
515          diagTitle = 'Vertical velocity times delta^k(Rho)_at-const-P'
516          diagUnits = 'kg/m^2/s        '
517          diagCode  = 'WM      LR      '
518          CALL DIAGNOSTICS_ADDTOLIST( diagNum,
519         I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
520    
521          diagName  = 'WdRHOdP '
522          diagTitle = 'Vertical velocity times delta^k(Rho)_at-const-T,S'
523        diagUnits = 'kg/m^2/s        '        diagUnits = 'kg/m^2/s        '
524        diagCode  = 'WM      LR      '        diagCode  = 'WM      LR      '
525        CALL DIAGNOSTICS_ADDTOLIST( diagNum,        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
# Line 542  c     diagTitle = 'Square of ocean botto Line 559  c     diagTitle = 'Square of ocean botto
559        CALL DIAGNOSTICS_ADDTOLIST( diagNum,        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
560       I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )       I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
561    
562          diagName  = 'PHI_SURF'
563          diagTitle = DIAGS_MK_TITLE('Surface Dynamical '//pTitle,myThid)
564          diagUnits = 'm^2/s^2         '
565          diagCode  = 'SM      M1      '
566          CALL DIAGNOSTICS_ADDTOLIST( diagNum,
567         I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
568    
569    #ifdef NONLIN_FRSURF
570          diagName  = 'PHIHYDcR'
571          diagTitle = DIAGS_MK_TITLE( 'Hydrostatic '
572         I                       //pTitle//' Anomaly @ const r', myThid )
573          diagUnits = 'm^2/s^2         '
574          diagCode  = 'SMR     MR      '
575          CALL DIAGNOSTICS_ADDTOLIST( diagNum,
576         I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
577    #endif
578    
579  #ifdef ALLOW_NONHYDROSTATIC  #ifdef ALLOW_NONHYDROSTATIC
580        diagName  = 'PHI_NH  '        diagName  = 'PHI_NH  '
581        diagTitle = DIAGS_MK_TITLE( 'Non-Hydrostatic '//pTitle, myThid )        diagTitle = DIAGS_MK_TITLE( 'Non-Hydrostatic '//pTitle, myThid )
# Line 732  c     diagTitle = 'Free-Surface r-Positi Line 766  c     diagTitle = 'Free-Surface r-Positi
766    
767        diagName  = 'TOTSTEND'        diagName  = 'TOTSTEND'
768        diagTitle = DIAGS_MK_TITLE('Tendency of '//sTitle, myThid )        diagTitle = DIAGS_MK_TITLE('Tendency of '//sTitle, myThid )
769        diagUnits = DIAGS_MK_UNITS( sUnit4c//'/day', myThid )        diagUnits = DIAGS_MK_UNITS( sUnit5c//'/day', myThid )
770          diagCode  = 'SMR     MR      '
771          CALL DIAGNOSTICS_ADDTOLIST( diagNum,
772         I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
773    
774          diagName  = 'MoistCor'
775          diagTitle = 'Heating correction due to moist thermodynamics'
776          diagUnits = 'W/m^2           '
777          diagCode  = 'SM      MR      '
778          CALL DIAGNOSTICS_ADDTOLIST( diagNum,
779         I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
780    
781    #ifdef ALLOW_FRICTION_HEATING
782          diagName  = 'HeatDiss'
783          diagTitle = 'Heating from frictional dissipation'
784          diagUnits = 'W/m^2           '
785          diagCode  = 'SM      MR      '
786          CALL DIAGNOSTICS_ADDTOLIST( diagNum,
787         I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
788    #endif /* ALLOW_FRICTION_HEATING */
789    
790    #ifdef ALLOW_GENERIC_ADVDIFF
791          diagName  = 'gT_Forc '
792          diagTitle = 'Potential Temp. forcing tendency'
793          diagUnits = DIAGS_MK_UNITS( tUnit4c//'/s', myThid )
794          diagCode  = 'SMR     MR      '
795          CALL DIAGNOSTICS_ADDTOLIST( diagNum,
796         I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
797    
798          diagName  = 'gS_Forc '
799          diagTitle = DIAGS_MK_TITLE(
800         &            sTitle//'forcing tendency', myThid )
801          diagUnits = DIAGS_MK_UNITS( sUnit5c//'/s', myThid )
802          diagCode  = 'SMR     MR      '
803          CALL DIAGNOSTICS_ADDTOLIST( diagNum,
804         I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
805    
806          diagName  = 'AB_gT   '
807          diagTitle = 'Potential Temp. tendency from Adams-Bashforth'
808          diagUnits = DIAGS_MK_UNITS( tUnit4c//'/s', myThid )
809        diagCode  = 'SMR     MR      '        diagCode  = 'SMR     MR      '
810        CALL DIAGNOSTICS_ADDTOLIST( diagNum,        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
811       I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )       I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
812    
813          diagName  = 'AB_gS   '
814          diagTitle = DIAGS_MK_TITLE(
815         &            sTitle//'tendency from Adams-Bashforth', myThid )
816          diagUnits = DIAGS_MK_UNITS( sUnit5c//'/s', myThid )
817          diagCode  = 'SMR     MR      '
818          CALL DIAGNOSTICS_ADDTOLIST( diagNum,
819         I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
820    
821          diagName  = 'gTinAB  '
822          diagTitle = 'Potential Temp. tendency going in Adams-Bashforth'
823          diagUnits = DIAGS_MK_UNITS( tUnit4c//'/s', myThid )
824          diagCode  = 'SMR     MR      '
825          CALL DIAGNOSTICS_ADDTOLIST( diagNum,
826         I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
827    
828          diagName  = 'gSinAB  '
829          diagTitle = DIAGS_MK_TITLE(
830         &            sTitle//'tendency going in Adams-Bashforth', myThid )
831          diagUnits = DIAGS_MK_UNITS( sUnit5c//'/s', myThid )
832          diagCode  = 'SMR     MR      '
833          CALL DIAGNOSTICS_ADDTOLIST( diagNum,
834         I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
835    #endif /* ALLOW_GENERIC_ADVDIFF */
836    
837          diagName  = 'AB_gU   '
838          diagTitle = 'U momentum tendency from Adams-Bashforth'
839          diagUnits = 'm/s^2           '
840          diagCode  = 'UUR     MR      '
841          diagMate  = diagNum + 2
842          CALL DIAGNOSTICS_ADDTOLIST( diagNum,
843         I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
844    
845          diagName  = 'AB_gV   '
846          diagTitle = 'V momentum tendency from Adams-Bashforth'
847          diagUnits = 'm/s^2           '
848          diagCode  = 'VVR     MR      '
849          diagMate  = diagNum
850          CALL DIAGNOSTICS_ADDTOLIST( diagNum,
851         I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
852    
853    #ifdef ALLOW_NONHYDROSTATIC
854          diagName  = 'AB_gW   '
855          diagTitle = 'W momentum tendency from Adams-Bashforth'
856          diagUnits = DIAGS_MK_UNITS( rUnit2c//'/s^2', myThid )
857          diagCode  = 'WM      LR      '
858          CALL DIAGNOSTICS_ADDTOLIST( diagNum,
859         I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
860    #endif /* ALLOW_NONHYDROSTATIC */
861    
862    #ifdef ALLOW_EDDYPSI
863          diagName  = 'TAUXEDDY'
864          diagTitle = 'Zonal Eddy Stress'
865          diagUnits = 'N/m^2           '
866          diagCode  = 'UU      LR      '
867          diagMate  = diagNum + 2
868          CALL DIAGNOSTICS_ADDTOLIST( diagNum,
869         I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
870    
871          diagName  = 'TAUYEDDY'
872          diagTitle = 'Meridional Eddy Stress'
873          diagUnits = 'N/m^2           '
874          diagCode  = 'VV      LR      '
875          diagMate  = diagNum
876          CALL DIAGNOSTICS_ADDTOLIST( diagNum,
877         I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
878    
879    # ifdef ALLOW_GMREDI
880          diagName  = 'U_EulerM'
881          diagTitle = 'Zonal Eulerian-Mean Velocity (m/s)'
882          diagUnits = 'm/s             '
883          diagCode  = 'UUR     MR      '
884          diagMate  = diagNum + 2
885          CALL DIAGNOSTICS_ADDTOLIST( diagNum,
886         I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
887    
888          diagName  = 'V_EulerM'
889          diagTitle = 'Meridional Eulerian-Mean Velocity (m/s)'
890          diagUnits = 'm/s             '
891          diagCode  = 'VVR     MR      '
892          diagMate  = diagNum
893          CALL DIAGNOSTICS_ADDTOLIST( diagNum,
894         I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
895    # endif /* ALLOW_GMREDI */
896    #endif /* ALLOW_EDDYPSI */
897    
898        RETURN        RETURN
899        END        END
Legend:
Removed from v.1.38  
changed lines
  Added in v.1.52
  ViewVC Help
Powered by ViewVC 1.1.22