/[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.50 by jmc, Sat Jul 19 16:46:25 2014 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    #ifdef NONLIN_FRSURF
563          diagName  = 'PHIHYDcR'
564          diagTitle = DIAGS_MK_TITLE( 'Hydrostatic '
565         I                       //pTitle//' Anomaly @ const r', myThid )
566          diagUnits = 'm^2/s^2         '
567          diagCode  = 'SMR     MR      '
568          CALL DIAGNOSTICS_ADDTOLIST( diagNum,
569         I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
570    #endif
571    
572  #ifdef ALLOW_NONHYDROSTATIC  #ifdef ALLOW_NONHYDROSTATIC
573        diagName  = 'PHI_NH  '        diagName  = 'PHI_NH  '
574        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 759  c     diagTitle = 'Free-Surface r-Positi
759    
760        diagName  = 'TOTSTEND'        diagName  = 'TOTSTEND'
761        diagTitle = DIAGS_MK_TITLE('Tendency of '//sTitle, myThid )        diagTitle = DIAGS_MK_TITLE('Tendency of '//sTitle, myThid )
762        diagUnits = DIAGS_MK_UNITS( sUnit4c//'/day', myThid )        diagUnits = DIAGS_MK_UNITS( sUnit5c//'/day', myThid )
763          diagCode  = 'SMR     MR      '
764          CALL DIAGNOSTICS_ADDTOLIST( diagNum,
765         I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
766    
767          diagName  = 'MoistCor'
768          diagTitle = 'Heating correction due to moist thermodynamics'
769          diagUnits = 'W/m^2           '
770          diagCode  = 'SM      MR      '
771          CALL DIAGNOSTICS_ADDTOLIST( diagNum,
772         I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
773    
774    #ifdef ALLOW_GENERIC_ADVDIFF
775          diagName  = 'gT_Forc '
776          diagTitle = 'Potential Temp. forcing tendency'
777          diagUnits = DIAGS_MK_UNITS( tUnit4c//'/s', myThid )
778          diagCode  = 'SMR     MR      '
779          CALL DIAGNOSTICS_ADDTOLIST( diagNum,
780         I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
781    
782          diagName  = 'gS_Forc '
783          diagTitle = DIAGS_MK_TITLE(
784         &            sTitle//'forcing tendency', myThid )
785          diagUnits = DIAGS_MK_UNITS( sUnit5c//'/s', myThid )
786          diagCode  = 'SMR     MR      '
787          CALL DIAGNOSTICS_ADDTOLIST( diagNum,
788         I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
789    
790          diagName  = 'AB_gT   '
791          diagTitle = 'Potential Temp. tendency from Adams-Bashforth'
792          diagUnits = DIAGS_MK_UNITS( tUnit4c//'/s', myThid )
793          diagCode  = 'SMR     MR      '
794          CALL DIAGNOSTICS_ADDTOLIST( diagNum,
795         I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
796    
797          diagName  = 'AB_gS   '
798          diagTitle = DIAGS_MK_TITLE(
799         &            sTitle//'tendency from Adams-Bashforth', myThid )
800          diagUnits = DIAGS_MK_UNITS( sUnit5c//'/s', myThid )
801        diagCode  = 'SMR     MR      '        diagCode  = 'SMR     MR      '
802        CALL DIAGNOSTICS_ADDTOLIST( diagNum,        CALL DIAGNOSTICS_ADDTOLIST( diagNum,
803       I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )       I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
804    
805          diagName  = 'gTinAB  '
806          diagTitle = 'Potential Temp. tendency going in Adams-Bashforth'
807          diagUnits = DIAGS_MK_UNITS( tUnit4c//'/s', myThid )
808          diagCode  = 'SMR     MR      '
809          CALL DIAGNOSTICS_ADDTOLIST( diagNum,
810         I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
811    
812          diagName  = 'gSinAB  '
813          diagTitle = DIAGS_MK_TITLE(
814         &            sTitle//'tendency going in Adams-Bashforth', myThid )
815          diagUnits = DIAGS_MK_UNITS( sUnit5c//'/s', myThid )
816          diagCode  = 'SMR     MR      '
817          CALL DIAGNOSTICS_ADDTOLIST( diagNum,
818         I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
819    #endif /* ALLOW_GENERIC_ADVDIFF */
820    
821          diagName  = 'AB_gU   '
822          diagTitle = 'U momentum tendency from Adams-Bashforth'
823          diagUnits = 'm/s^2           '
824          diagCode  = 'UUR     MR      '
825          diagMate  = diagNum + 2
826          CALL DIAGNOSTICS_ADDTOLIST( diagNum,
827         I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
828    
829          diagName  = 'AB_gV   '
830          diagTitle = 'V momentum tendency from Adams-Bashforth'
831          diagUnits = 'm/s^2           '
832          diagCode  = 'VVR     MR      '
833          diagMate  = diagNum
834          CALL DIAGNOSTICS_ADDTOLIST( diagNum,
835         I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
836    
837    #ifdef ALLOW_NONHYDROSTATIC
838          diagName  = 'AB_gW   '
839          diagTitle = 'W momentum tendency from Adams-Bashforth'
840          diagUnits = DIAGS_MK_UNITS( rUnit2c//'/s^2', myThid )
841          diagCode  = 'WM      LR      '
842          CALL DIAGNOSTICS_ADDTOLIST( diagNum,
843         I          diagName, diagCode, diagUnits, diagTitle, 0, myThid )
844    #endif /* ALLOW_NONHYDROSTATIC */
845    
846    #ifdef ALLOW_EDDYPSI
847          diagName  = 'TAUXEDDY'
848          diagTitle = 'Zonal Eddy Stress'
849          diagUnits = 'N/m**2          '
850          diagCode  = 'UU P    MR      '
851          diagMate  = diagNum
852          CALL DIAGNOSTICS_ADDTOLIST( diagNum,
853         I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
854    
855          diagName  = 'TAUYEDDY'
856          diagTitle = 'Meridional Eddy Stress'
857          diagUnits = 'N/m**2          '
858          diagCode  = 'VV P    MR      '
859          diagMate  = diagNum
860          CALL DIAGNOSTICS_ADDTOLIST( diagNum,
861         I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
862    
863    # ifdef ALLOW_GMREDI
864          diagName  = 'UVELMEAN'
865          diagTitle = 'Zonal Component of Mean Velocity (m/s)'
866          diagUnits = 'm/s             '
867          diagCode  = 'UUR     MR      '
868          diagMate  = diagNum + 2
869          CALL DIAGNOSTICS_ADDTOLIST( diagNum,
870         I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
871    
872          diagName  = 'VVELMEAN'
873          diagTitle = 'Meridional Component of Mean Velocity (m/s)'
874          diagUnits = 'm/s             '
875          diagCode  = 'VVR     MR      '
876          diagMate  = diagNum
877          CALL DIAGNOSTICS_ADDTOLIST( diagNum,
878         I   diagName, diagCode, diagUnits, diagTitle, diagMate, myThid )
879    # endif /* ALLOW_GMREDI */
880    #endif /* ALLOW_EDDYPSI */
881    
882        RETURN        RETURN
883        END        END
Legend:
Removed from v.1.38  
changed lines
  Added in v.1.50
  ViewVC Help
Powered by ViewVC 1.1.22