--- MITgcm/pkg/diagnostics/diagnostics_main_init.F 2005/02/10 15:01:34 1.9 +++ MITgcm/pkg/diagnostics/diagnostics_main_init.F 2005/08/16 21:29:38 1.20 @@ -1,4 +1,4 @@ -C $Header: /home/ubuntu/mnt/e9_copy/MITgcm/pkg/diagnostics/diagnostics_main_init.F,v 1.9 2005/02/10 15:01:34 jmc Exp $ +C $Header: /home/ubuntu/mnt/e9_copy/MITgcm/pkg/diagnostics/diagnostics_main_init.F,v 1.20 2005/08/16 21:29:38 molod Exp $ C $Name: $ #include "DIAG_OPTIONS.h" @@ -44,6 +44,8 @@ C & R & levels = Nr \\ C & L & levels = MAX(Nr,NrPhys) \\ C & M & levels = MAX(Nr,NrPhys) - 1 \\ +C & G & levels = Ground_level Number \\ +C & I & levels = sea-Ice_level Number \\ C \end{tabular} C \end{center} @@ -51,7 +53,6 @@ IMPLICIT NONE #include "SIZE.h" #include "EEPARAMS.h" -c #include "EESUPPORT.h" #include "PARAMS.h" C !INPUT PARAMETERS: @@ -84,6 +85,11 @@ IF ( fluidIsAir ) THEN tUnit4c= 'K ' sUnit4c= 'g/kg' +#ifdef ALLOW_FIZHI + IF(useFIZHI)THEN + sUnit4c= 'kg/kg' + ENDIF +#endif /* ALLOW_FIZHI */ ELSE tUnit4c= 'degC' sUnit4c= 'psu ' @@ -120,8 +126,36 @@ CALL DIAGNOSTICS_ADD2LIST( diagNum, I diagName, diagCode, diagUnits, diagTitle, myThid ) + diagName = 'SST ' + diagTitle = 'Sea Surface Temperature (degC,K)' + diagUnits = DIAGS_MK_UNITS( tUnit4c, myThid ) + diagCode = 'SM M1 ' + CALL DIAGNOSTICS_ADD2LIST( diagNum, + I diagName, diagCode, diagUnits, diagTitle, myThid ) + diagName = 'SALT ' - diagTitle = 'Salt or Water Vapor mixing ratio (g/kg)' + diagTitle = 'Salt or Water Vapor mixing ratio' + diagUnits = DIAGS_MK_UNITS( sUnit4c, myThid ) + diagCode = 'SM MR ' + CALL DIAGNOSTICS_ADD2LIST( diagNum, + I diagName, diagCode, diagUnits, diagTitle, myThid ) + + diagName = 'RELHUM ' + diagTitle = 'Relative Humidity ' + diagUnits = 'percent ' + diagCode = 'SM MR ' + CALL DIAGNOSTICS_ADD2LIST( diagNum, + I diagName, diagCode, diagUnits, diagTitle, myThid ) + + diagName = 'SSS ' + diagTitle = 'Sea Surface Salinity ' + diagUnits = DIAGS_MK_UNITS( sUnit4c, myThid ) + diagCode = 'SM M1 ' + CALL DIAGNOSTICS_ADD2LIST( diagNum, + I diagName, diagCode, diagUnits, diagTitle, myThid ) + + diagName = 'SALTanom' + diagTitle = 'Salt anomaly (=SALT-35)' diagUnits = DIAGS_MK_UNITS( sUnit4c, myThid ) diagCode = 'SM MR ' CALL DIAGNOSTICS_ADD2LIST( diagNum, @@ -141,6 +175,20 @@ CALL DIAGNOSTICS_ADD2LIST( diagNum, I diagName, diagCode, diagUnits, diagTitle, myThid ) + diagName = 'UVEL_k2 ' + diagTitle = 'Zonal Component of Velocity at level 2 (m/s)' + diagUnits = 'm/s ' + WRITE(diagCode,'(A,I3.3,A)') 'UU ',diagNum+2,'M1 ' + CALL DIAGNOSTICS_ADD2LIST( diagNum, + I diagName, diagCode, diagUnits, diagTitle, myThid ) + + diagName = 'VVEL_k2 ' + diagTitle = 'Meridional Component of Velocity at level 2 (m/s)' + diagUnits = 'm/s ' + WRITE(diagCode,'(A,I3.3,A)') 'VV ', diagNum ,'M1 ' + CALL DIAGNOSTICS_ADD2LIST( diagNum, + I diagName, diagCode, diagUnits, diagTitle, myThid ) + diagName = 'WVEL ' diagTitle = 'Vertical Component of Velocity (r_units/s)' diagUnits = DIAGS_MK_UNITS( rUnit2c//'/s', myThid ) @@ -162,6 +210,13 @@ CALL DIAGNOSTICS_ADD2LIST( diagNum, I diagName, diagCode, diagUnits, diagTitle, myThid ) + diagName = 'SALTSQan' + diagTitle = 'Square of Salt anomaly (=(SALT-35)^2 (g^2/kg^2)' + diagUnits = DIAGS_MK_UNITS( '('//sUnit4c//')^2', myThid ) + diagCode = 'SM MR ' + CALL DIAGNOSTICS_ADD2LIST( diagNum, + I diagName, diagCode, diagUnits, diagTitle, myThid ) + diagName = 'UVELSQ ' diagTitle = 'Square of Zonal Comp of Velocity (m^2/s^2)' diagUnits = 'm^2/s^2 ' @@ -228,7 +283,7 @@ diagName = 'WVELMASS' diagTitle = 'Vertical Mass-Weighted Comp of Velocity (m/s)' diagUnits = DIAGS_MK_UNITS( rUnit2c//'/s', myThid ) - diagCode = 'WM MR ' + diagCode = 'WM LR ' CALL DIAGNOSTICS_ADD2LIST( diagNum, I diagName, diagCode, diagUnits, diagTitle, myThid ) @@ -249,28 +304,28 @@ diagName = 'WTHMASS ' diagTitle = 'Vertical Mass-Weight Transp of Pot Temp (K.m/s)' diagUnits = DIAGS_MK_UNITS(tUnit4c//'.'//rUnit2c//'/s', myThid ) - diagCode = 'WM MR ' + diagCode = 'WM LR ' CALL DIAGNOSTICS_ADD2LIST( diagNum, I diagName, diagCode, diagUnits, diagTitle, myThid ) diagName = 'USLTMASS' - diagTitle = 'Zonal Mass-Weight Transp of Salt (g/kg.m/s)' + diagTitle = 'Zonal Mass-Weight Transp of Salt ' diagUnits = DIAGS_MK_UNITS(sUnit4c//'.m/s', myThid ) WRITE(diagCode,'(A,I3.3,A)') 'UU ',diagNum+2,'MR ' CALL DIAGNOSTICS_ADD2LIST( diagNum, I diagName, diagCode, diagUnits, diagTitle, myThid ) diagName = 'VSLTMASS' - diagTitle = 'Meridional Mass-Weight Transp of Salt (g/kg.m/s)' + diagTitle = 'Meridional Mass-Weight Transp of Salt' diagUnits = DIAGS_MK_UNITS(sUnit4c//'.m/s', myThid ) WRITE(diagCode,'(A,I3.3,A)') 'VV ', diagNum ,'MR ' CALL DIAGNOSTICS_ADD2LIST( diagNum, I diagName, diagCode, diagUnits, diagTitle, myThid ) diagName = 'WSLTMASS' - diagTitle = 'Vertical Mass-Weight Transp of Salt (g/kg.m/s)' + diagTitle = 'Vertical Mass-Weight Transp of Salt ' diagUnits = DIAGS_MK_UNITS(sUnit4c//'.'//rUnit2c//'/s', myThid ) - diagCode = 'WM MR ' + diagCode = 'WM LR ' CALL DIAGNOSTICS_ADD2LIST( diagNum, I diagName, diagCode, diagUnits, diagTitle, myThid ) @@ -291,28 +346,28 @@ diagName = 'WVELTH ' diagTitle = 'Vertical Transp of Pot Temp (K.m/s)' diagUnits = DIAGS_MK_UNITS(tUnit4c//'.'//rUnit2c//'/s', myThid ) - diagCode = 'WM MR ' + diagCode = 'WM LR ' CALL DIAGNOSTICS_ADD2LIST( diagNum, I diagName, diagCode, diagUnits, diagTitle, myThid ) diagName = 'UVELSLT ' - diagTitle = 'Zonal Transp of Salt (g/kg.m/s)' + diagTitle = 'Zonal Transp of Salt ' diagUnits = DIAGS_MK_UNITS( sUnit4c//'.m/s', myThid ) WRITE(diagCode,'(A,I3.3,A)') 'UU ',diagNum+2,'MR ' CALL DIAGNOSTICS_ADD2LIST( diagNum, I diagName, diagCode, diagUnits, diagTitle, myThid ) diagName = 'VVELSLT ' - diagTitle = 'Meridional Transp of Salt (g/kg.m/s)' + diagTitle = 'Meridional Transp of Salt ' diagUnits = DIAGS_MK_UNITS( sUnit4c//'.m/s', myThid ) WRITE(diagCode,'(A,I3.3,A)') 'VV ', diagNum ,'MR ' CALL DIAGNOSTICS_ADD2LIST( diagNum, I diagName, diagCode, diagUnits, diagTitle, myThid ) diagName = 'WVELSLT ' - diagTitle = 'Vertical Transp of Salt (g/kg.m/s)' + diagTitle = 'Vertical Transp of Salt ' diagUnits = DIAGS_MK_UNITS(sUnit4c//'.'//rUnit2c//'/s', myThid ) - diagCode = 'WM MR ' + diagCode = 'WM LR ' CALL DIAGNOSTICS_ADD2LIST( diagNum, I diagName, diagCode, diagUnits, diagTitle, myThid ) @@ -325,6 +380,34 @@ CALL DIAGNOSTICS_ADD2LIST( diagNum, I diagName, diagCode, diagUnits, diagTitle, myThid ) + diagName = 'RHOANOSQ' + diagTitle = 'Square of Density Anomaly (=(Rho-rhoConst)^2)' + diagUnits = 'kg^2/m^6 ' + diagCode = 'SM MR ' + CALL DIAGNOSTICS_ADD2LIST( diagNum, + I diagName, diagCode, diagUnits, diagTitle, myThid ) + + diagName = 'URHOMASS' + diagTitle = 'Zonal Transport of Density' + diagUnits = 'kg/m^2/s ' + WRITE(diagCode,'(A,I3.3,A)') 'UU ',diagNum+2,'MR ' + CALL DIAGNOSTICS_ADD2LIST( diagNum, + I diagName, diagCode, diagUnits, diagTitle, myThid ) + + diagName = 'VRHOMASS' + diagTitle = 'Meridional Transport of Density' + diagUnits = 'kg/m^2/s ' + WRITE(diagCode,'(A,I3.3,A)') 'VV ',diagNum,'MR ' + CALL DIAGNOSTICS_ADD2LIST( diagNum, + I diagName, diagCode, diagUnits, diagTitle, myThid ) + + diagName = 'WRHOMASS' + diagTitle = 'Vertical Transport of Potential Density' + diagUnits = 'kg/m^2/s ' + diagCode = 'WM LR ' + CALL DIAGNOSTICS_ADD2LIST( diagNum, + I diagName, diagCode, diagUnits, diagTitle, myThid ) + diagName = 'PHIHYD ' diagTitle = 'Hydrostatic (ocean) pressure / (atmos) geo-Potential' diagUnits = 'm^2/s^2 ' @@ -332,6 +415,13 @@ CALL DIAGNOSTICS_ADD2LIST( diagNum, I diagName, diagCode, diagUnits, diagTitle, myThid ) + diagName = 'PHIHYDSQ' + diagTitle = 'Square of Hyd. (ocean) press / (atmos) geoPotential' + diagUnits = 'm^4/s^4 ' + diagCode = 'SM MR ' + CALL DIAGNOSTICS_ADD2LIST( diagNum, + I diagName, diagCode, diagUnits, diagTitle, myThid ) + diagName = 'PHIBOT ' diagTitle = 'ocean bottom pressure / top. atmos geo-Potential' diagUnits = 'm^2/s^2 ' @@ -419,7 +509,30 @@ CALL DIAGNOSTICS_ADD2LIST( diagNum, I diagName, diagCode, diagUnits, diagTitle, myThid ) +#ifdef SHORTWAVE_HEATING + diagName = 'SWFLUX ' + diagTitle = 'net upward SW radiation, >0 increases theta ' + diagUnits = 'W/m^2 ' + diagCode = 'SM U1 ' + CALL DIAGNOSTICS_ADD2LIST( diagNum, + I diagName, diagCode, diagUnits, diagTitle, myThid ) +#endif + C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----| + diagName = 'PRESSURE' + diagTitle = 'Atmospheric Pressure (Pa) ' + diagUnits = DIAGS_MK_UNITS( rUnit2c, myThid ) + diagCode = 'SM MR ' + CALL DIAGNOSTICS_ADD2LIST( diagNum, + I diagName, diagCode, diagUnits, diagTitle, myThid ) + + diagName = 'RSURF ' + diagTitle = 'Surface (Pressure, Height) (Pa,m)' + diagUnits = DIAGS_MK_UNITS( rUnit2c, myThid ) + diagCode = 'SM M1 ' + CALL DIAGNOSTICS_ADD2LIST( diagNum, + I diagName, diagCode, diagUnits, diagTitle, myThid ) + RETURN END