--- MITgcm/pkg/diagnostics/diagnostics_main_init.F	2005/02/10 15:01:34	1.9
+++ MITgcm/pkg/diagnostics/diagnostics_main_init.F	2005/02/15 20:57:38	1.10
@@ -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.10 2005/02/15 20:57:38 jmc Exp $
 C $Name:  $
 
 #include "DIAG_OPTIONS.h"
@@ -228,7 +228,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,7 +249,7 @@
       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 )
 
@@ -270,7 +270,7 @@
       diagName  = 'WSLTMASS'
       diagTitle = 'Vertical Mass-Weight Transp of Salt (g/kg.m/s)'
       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,7 +291,7 @@
       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 )
 
@@ -312,7 +312,7 @@
       diagName  = 'WVELSLT '
       diagTitle = 'Vertical Transp of Salt (g/kg.m/s)'
       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 +325,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         '