Changeset 10824

Timestamp:

2019-04-02T15:32:49+02:00 (6 years ago)

Author:

clem

Message:

output cfl diagnostics (courant numbers) in netcdf => add cfl_cu, cfl_cv, cfl_cw

Location:

NEMO/trunk

Files:

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

NEMO/trunk/cfgs/SHARED/field_def_nemo-oce.xml

@@ -108 +108 @@
         <field id="tpt_dep"      long_name="T-point depth"                  standard_name="depth_below_geoid"   unit="m"   grid_ref="grid_T_3D" />
         <field id="e3tdef"       long_name="T-cell thickness deformation"                                       unit="%"   grid_ref="grid_T_3D" />
+
+        <!-- variables available with ln_diacfl=.true. -->
+        <field id="cfl_cu"       long_name="u-courant number"   unit="#" />
+        <field id="cfl_cv"       long_name="v-courant number"   unit="#" />
+        <field id="cfl_cw"       long_name="w-courant number"   unit="#" />
+
       </field_group> <!-- grid_T -->
 

NEMO/trunk/src/OCE/DIA/diacfl.F90

@@ -17 +17 @@
    USE lbclnk          ! ocean lateral boundary condition (or mpp link)
    USE in_out_manager  ! I/O manager
+   USE iom             ! 
    USE timing          ! Performance output
 
@@ -78 +79 @@
       END DO
       !
+      ! write outputs
+      IF( iom_use('cfl_cu') )   CALL iom_put( 'cfl_cu', MAXVAL( zCu_cfl, dim=3 ) )
+      IF( iom_use('cfl_cv') )   CALL iom_put( 'cfl_cv', MAXVAL( zCv_cfl, dim=3 ) )
+      IF( iom_use('cfl_cw') )   CALL iom_put( 'cfl_cw', MAXVAL( zCw_cfl, dim=3 ) )
+
       !                    ! calculate maximum values and locations
       IF( lk_mpp ) THEN