Changeset 5702 for branches/UKMO

Timestamp:

2015-08-18T18:21:29+02:00 (9 years ago)

Author:

davestorkey

Message:

Add in new diagnostics.

Location:

branches/UKMO/smagorinsky_diagnostics/NEMOGCM

Files:

• CONFIG/SHARED/field_def.xml (modified) (2 diffs)
• NEMO/OPA_SRC/LDF/ldfdyn_smag.F90 (modified) (2 diffs)
• NEMO/OPA_SRC/LDF/ldftra_smag.F90 (modified) (1 diff)

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

@@ -375 +375 @@
          <field id="uoces"        long_name="ocean transport along i-axis times salinity (CRS)"                                                  unit="1e-3*m/s"   grid_ref="grid_U_3D" />
 
+         <!-- variables available with key_dynldf_smag -->
+         <field id="ahmu_lap_smag"      long_name="ocean laplacian momentum diffusivity at U-points"   unit="m4/s"   grid_ref="grid_U_3D" />
+         <field id="ahmu_bilap_smag"    long_name="ocean bilaplacian momentum diffusivity at U-points"   unit="m4/s"   grid_ref="grid_U_3D" />
+
+         <!-- variables available with key_traldf_smag -->
+         <field id="ahtu_lap_smag"      long_name="ocean laplacian tracer diffusivity at U-points"   unit="m4/s"   grid_ref="grid_U_3D" />
+
          <!-- variables available with MLE -->
          <field id="psiu_mle"     long_name="MLE streamfunction along i-axis"   unit="m3/s"   grid_ref="grid_U_3D" />
@@ -409 +416 @@
          <field id="vocet"        long_name="ocean transport along j-axis times temperature (CRS)"                                               unit="degC*m/s"   grid_ref="grid_V_3D" />
          <field id="voces"        long_name="ocean transport along j-axis times salinity (CRS)"                                                  unit="1e-3*m/s"   grid_ref="grid_V_3D" />
+
+         <!-- variables available with key_dynldf_smag -->
+         <field id="ahmv_lap_smag"      long_name="ocean laplacian momentum diffusivity at V-points"   unit="m4/s"   grid_ref="grid_V_3D" />
+         <field id="ahmv_bilap_smag"    long_name="ocean bilaplacian momentum diffusivity at V-points"   unit="m4/s"   grid_ref="grid_V_3D" />
+
+         <!-- variables available with key_traldf_smag -->
+         <field id="ahtv_lap_smag"      long_name="ocean laplacian tracer diffusivity at V-points"   unit="m4/s"   grid_ref="grid_V_3D" />
 
          <!-- variables available with MLE -->

branches/UKMO/smagorinsky_diagnostics/NEMOGCM/NEMO/OPA_SRC/LDF/ldfdyn_smag.F90

@@ -181 +181 @@
       CALL lbc_lnk( ahm1, 'T', 1. )   ! Lateral boundary conditions on ( ahtt )
       CALL lbc_lnk( ahm2, 'F', 1. )   ! Lateral boundary conditions on ( ahtt )
+      CALL iom_put( 'ahmu_lap_smag', ahm1(:,:,:) ) 
+      CALL iom_put( 'ahmv_lap_smag', ahm2(:,:,:) )
 
       ENDIF    ! ln_dynldf_lap
@@ -253 +255 @@
       CALL lbc_lnk( ahm3, 'U', 1. )   ! Lateral boundary conditions
       CALL lbc_lnk( ahm4, 'V', 1. )
+      CALL iom_put( 'ahmu_bilap_smag', ahm3(:,:,:) ) 
+      CALL iom_put( 'ahmv_bilap_smag', ahm4(:,:,:) )
+
    ENDIF
 

branches/UKMO/smagorinsky_diagnostics/NEMOGCM/NEMO/OPA_SRC/LDF/ldftra_smag.F90

@@ -174 +174 @@
         CALL lbc_lnk( ahtu, 'U', 1. )   ! Lateral boundary conditions
         CALL lbc_lnk( ahtv, 'V', 1. )
+        CALL iom_put( 'ahtu_lap_smag', aht1(:,:,:) ) 
+        CALL iom_put( 'ahtv_lap_smag', aht2(:,:,:) )
 
 IF(  kt == nit000 ) THEN