Changeset 13114 for branches

Timestamp:

2020-06-16T17:48:15+02:00 (4 years ago)

Author:

mattmartin

Message:

Include diagnostic for SSH tendencies and their absolute values.

Location:

branches/UKMO/dev_r5518_GO6_package_FOAMv14/NEMOGCM

Files:

• CONFIG/SHARED/field_def.xml (modified) (1 diff)
• NEMO/OPA_SRC/ASM/bias.F90 (modified) (1 diff)
• NEMO/OPA_SRC/DIA/diawri.F90 (modified) (1 diff)

branches/UKMO/dev_r5518_GO6_package_FOAMv14/NEMOGCM/CONFIG/SHARED/field_def.xml

@@ -65 +65 @@
          <field id="ssh2"         long_name="square of sea surface height"   standard_name="square_of_sea_surface_height_above_geoid"   unit="m2" > ssh * ssh </field >
          <field id="sshmax"       long_name="max of sea surface height"   field_ref="ssh"   operation="maximum"                                  />
-
+         <field id="sshdt"        long_name="sea surface height tendency"    standard_name="tendency_of_sea_surface_height_above_sea_level" unit="m/s" />          
+         <field id="sshdtabs"     long_name="absolute sea surface height tendency"    standard_name="absolute_tendency_of_sea_surface_height_above_sea_level" unit="m/s" />          
+         
          <field id="mldkz5"       long_name="Turbocline depth (Kz = 5e-4)"                       standard_name="ocean_mixed_layer_thickness_defined_by_vertical_tracer_diffusivity"                unit="m"          />
          <field id="mldr10_1"     long_name="Mixed Layer Depth (dsigma = 0.01 wrt 10m)"          standard_name="ocean_mixed_layer_thickness_defined_by_sigma_theta"                                unit="m"          />

branches/UKMO/dev_r5518_GO6_package_FOAMv14/NEMOGCM/NEMO/OPA_SRC/ASM/bias.F90

@@ -283 +283 @@
          WRITE(numout,*) '     bias T and S apply             ln_bias_ts_app = ',ln_bias_ts_app
          WRITE(numout,*) '     bias pressure correctn apply   ln_bias_pc_app = ',ln_bias_pc_app
-         WRITE(numout,*) '     bias pressure correctn apply   ln_bias_pc_app = ',ln_bias_pc_app
          WRITE(numout,*) '     lat ramp for bias correction   nn_lat_ramp    = ',nn_lat_ramp
          WRITE(numout,*) '     time step for writing bias fld nn_bias_itwrt  = ',nn_bias_itwrt

branches/UKMO/dev_r5518_GO6_package_FOAMv14/NEMOGCM/NEMO/OPA_SRC/DIA/diawri.F90

@@ -161 +161 @@
 
 
-      CALL iom_put( "ssh" , sshn )                 ! sea surface height
-      
+      CALL iom_put( "ssh" , sshn )                    ! sea surface height
+      CALL iom_put( "sshdt" , ( sshn - sshb ) / rdt ) ! sea surface height tendency (per second)
+      CALL iom_put( "sshdtabs" , ABS ( ( sshn - sshb ) / rdt ) ) ! absolute value of sea surface height tendency (per second)
+            
       CALL iom_put( "toce", tsn(:,:,:,jp_tem) )    ! 3D temperature
       CALL theta2t ! in-situ temperature conversion