Changeset 13828 for NEMO/branches/UKMO

Timestamp:

2020-11-19T11:19:28+01:00 (3 years ago)

Author:

davestorkey

Message:

UKMO/NEMO_4.0.1_density_diag : code changes.

Location:

NEMO/branches/UKMO/NEMO_4.0.1_density_diag

Files:

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

NEMO/branches/UKMO/NEMO_4.0.1_density_diag/cfgs/SHARED/field_def_nemo-oce.xml

@@ -70 +70 @@
         <field id="bn2"          long_name="squared Brunt-Vaisala frequency"                                           unit="s-1"    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="rho_insitu"   long_name="in situ density"              standard_name="sea_water_density"       unit="kg/m3"  grid_ref="grid_T_3D" />
+        <field id="rho_insitu_e3t"  long_name="in situ density (thickness weighted)"                      unit="kg"     grid_ref="grid_T_3D" > rho_insitu * e3t </field >
 
         <!-- Energy - horizontal divergence -->

NEMO/branches/UKMO/NEMO_4.0.1_density_diag/src/OCE/DIA/diawri.F90

@@ -168 +168 @@
          CALL iom_put( "sbs", z2d )                ! bottom salinity
       ENDIF
+
+      IF ( iom_use( "rho_insitu" ) ) CALL iom_put( "rho_insitu", rhd(:,:,:) * rau0 + rau0 )  ! 3D in situ density
 
       IF ( iom_use("taubot") ) THEN                ! bottom stress