Changeset 15699

Timestamp:

2022-02-08T11:37:37+01:00 (2 years ago)

Author:

davestorkey

Message:

Fix bug in potential density diagnostic (duplicate iom_put call)
and also add thickness-weighted version.

Location:

NEMO/branches/UKMO/NEMO_4.0.4_GO8_package

Files:

• cfgs/SHARED/field_def_nemo-oce.xml (modified) (1 diff)
• src/OCE/DIA/diaprod.F90 (modified) (1 diff)

NEMO/branches/UKMO/NEMO_4.0.4_GO8_package/cfgs/SHARED/field_def_nemo-oce.xml

@@ -157 +157 @@
         <field id="beta"         long_name="haline contraction"                                                        unit="1e3"    grid_ref="grid_T_3D" />
         <field id="rhop"         long_name="potential density (sigma0)"        standard_name="sea_water_sigma_theta"   unit="kg/m3"  grid_ref="grid_T_3D" />
+        <field id="rhop_e3t"     long_name="potential density (sigma0) thickness weighted"        standard_name="sea_water_sigma_theta"   unit="kg/m3"  grid_ref="grid_T_3D"> rhop * e3t </field>
         <field id="N_2d"         long_name="Depth-mean N"                                                              unit="m/s" />
         <field id="modslp"       long_name="sqrt( slpi^2 + slpj^2 )"                                                   unit="1"      grid_ref="grid_T_3D" />

NEMO/branches/UKMO/NEMO_4.0.4_GO8_package/src/OCE/DIA/diaprod.F90

@@ -65 +65 @@
       !
 
-      IF( iom_use("urhop") .OR. iom_use("vrhop") .OR. iom_use("wrhop") &
-#if ! defined key_diaar5
-     &  .OR. iom_use("rhop") &
-#endif
-     & ) THEN
+      IF( iom_use("urhop") .OR. iom_use("vrhop") .OR. iom_use("wrhop") ) THEN
          CALL eos( tsn, z3d, zrhop )                 ! now in situ and potential density
          zrhop(:,:,:) = zrhop(:,:,:)-1000.e0         ! reference potential density to 1000 to avoid precision issues in rhop2 calculation
          zrhop(:,:,jpk) = 0._wp
-#if ! defined key_diaar5
-         CALL iom_put( 'rhop', zrhop )
-#else
-         ! If key_diaar5 set then there is already an iom_put call to output rhop.
-         ! Really should be a standard diagnostics option?
-#endif
       ENDIF