Changeset 3809

Timestamp:

2013-02-13T19:49:05+01:00 (11 years ago)

Author:

cetlod

Message:

2013/dev_r3411_CNRS4_IOCRS : add new diagnostics

Location:

branches/2013/dev_r3411_CNRS4_IOCRS/NEMOGCM

Files:

• CONFIG/ORCA2_LIM_CRS/EXP00/iodef.xml (modified) (8 diffs)
• CONFIG/ORCA2_LIM_CRS/MY_SRC/diawri.F90 (added)
• CONFIG/ORCA2_LIM_CRS/MY_SRC/oce.F90 (added)
• NEMO/OPA_SRC/CRS/crs_dom.F90 (modified) (2 diffs)
• NEMO/OPA_SRC/CRS/crsdiawri.F90 (modified) (1 diff)

branches/2013/dev_r3411_CNRS4_IOCRS/NEMOGCM/CONFIG/ORCA2_LIM_CRS/EXP00/iodef.xml

@@ -18 +18 @@
       
       <group id="grid_T" axis_ref="none" grid_ref="grid_T">
-   <field id="toce"         description="temperature"                               unit="degC" axis_ref="deptht"   />
+        <field id="toce"         description="temperature"                               unit="degC" axis_ref="deptht"   />
          <field id="soce"         description="salinity"                                  unit="psu"  axis_ref="deptht"   />
    <field id="sst"          description="sea surface temperature"                   unit="degC"                     />
@@ -31 +31 @@
    <field id="mldr10_1"     description="Mixed Layer Depth 0.01 ref.10m"            unit="m"                        />
          <field id="rhop"         description="potential density (sigma0)"                unit="kg/m3" axis_ref="deptht"  />
+   <field id="eken"         description="kinetic energy"                            unit="m2/s2" axis_ref="deptht"   />
    <!-- next variables available with key_diahth -->
    <field id="mlddzt"       description="Thermocline Depth (max dT/dz)"             unit="m"                        />
@@ -294 +295 @@
    <field id="sss_crs"    description="sea surface salinity"    unit="psu"                    />
    <field id="hdiv_crs"   description="horizontal divergence"   unit="s-1" axis_ref="deptht"  />
-
+        <field id="eken_crs"   description="kinetic energy"          unit="m2/s2" axis_ref="deptht" />
       </group>
 
@@ -341 +342 @@
       
       <group id="1d" output_freq="86400"  output_level="10" enabled=".FALSE.">                      <!-- 1d files -->
+
+          <file id="1d_grid_T" name="auto" description="ocean T grid coarsened variables" >
+             <field ref="toce"    name="votemper" />
+             <field ref="soce"    name="vosaline" />
+             <field ref="sst"     name="sosstsst" />
+             <field ref="sss"     name="sosaline" />
+             <field ref="ssh"     name="sossheig" />
+             <field ref="eken"         name="energkin"  />
+          </file>
+
           <file id="1d_gcrs_T" name="auto" description="ocean T grid coarsened variables" >
              <field ref="toce_crs"    name="votemper" />
@@ -348 +359 @@
              <field ref="ssh_crs"     name="sossheig" />
              <field ref="hdiv_crs"    name="vohdiver" />
+             <field ref="eken_crs"    name="energkin" />
           </file>
+      
 <!-- 
 ............................................................................................................
@@ -363 +376 @@
          <field ref="toce" name="votemper"  />
          <field ref="uoce" name="vozocrtx"  />     <!-- include a U-grid variable in the list => switch to T-grid -->
+              <field ref="eken"         name="energkin"  />
+
        </group>
      </file>
@@ -398 +413 @@
       <group id="3d" output_freq="259200" output_level="10" enabled=".FALSE.">                      <!-- 3d files -->
       </group>
-      
+
       <group id="5d" output_freq="432000" output_level="10" enabled=".TRUE.">                      <!-- 5d files -->
-
-   <file id="5d_grid_T" name="auto" description="ocean T grid variables" >
-     <field ref="toce"         name="votemper"  />
+        <file id="5d_grid_T" name="auto" description="ocean T grid coarsened variables" >
+     <field ref="eken"         name="energkin"  />
+          <field ref="toce"         name="votemper"  />
      <field ref="soce"         name="vosaline"  />
      <field ref="sst"          name="sosstsst"  />
@@ -426 +441 @@
 
         <file id="5d_gcrs_T" name="auto" description="ocean T grid coarsened variables" >
+           <field ref="eken_crs"    name="energkin" />
            <field ref="toce_crs"    name="votemper" />
            <field ref="soce_crs"    name="vosaline" />

branches/2013/dev_r3411_CNRS4_IOCRS/NEMOGCM/NEMO/OPA_SRC/CRS/crs_dom.F90

@@ -107 +107 @@
       ! Physical and dynamical ocean fields for output or passing to TOP, time-mean fields
       REAL(wp), DIMENSION(:,:,:,:), ALLOCATABLE      :: tsn_crs
-      REAL(wp), DIMENSION(:,:,:)  , ALLOCATABLE      :: un_crs, vn_crs, wn_crs
+      REAL(wp), DIMENSION(:,:,:)  , ALLOCATABLE      :: un_crs, vn_crs, wn_crs, rke_crs
       REAL(wp), DIMENSION(:,:,:)  , ALLOCATABLE      :: hdivn_crs    
       REAL(wp), DIMENSION(:,:)    , ALLOCATABLE      :: sshn_crs    
@@ -191 +191 @@
 
       ALLOCATE( un_crs(jpi_crs,jpj_crs,jpk) , vn_crs(jpi_crs,jpj_crs,jpk) , &
-         &      wn_crs(jpi_crs,jpj_crs,jpk) , hdivn_crs(jpi_crs,jpj_crs,jpk), STAT=ierr(11))
+         &      wn_crs(jpi_crs,jpj_crs,jpk) , hdivn_crs(jpi_crs,jpj_crs,jpk),&
+         &      rke_crs(jpi_crs,jpj_crs,jpk),                                STAT=ierr(11))
 
       ALLOCATE( sshn_crs(jpi_crs,jpj_crs),  emp_crs(jpi_crs,jpj_crs)    , &

branches/2013/dev_r3411_CNRS4_IOCRS/NEMOGCM/NEMO/OPA_SRC/CRS/crsdiawri.F90

@@ -189 +189 @@
          &         p_fse3=zfse3v, p_pfield=zs, p_cfield3d=zs_crs )
       CALL crs_iom_put( "voces_crs" , pv_r3d=zs_crs )   ! vS
+      
+      !  Kinetic energy
+      CALL crsfun( p_e1e2t=e1e2t, cd_type='T', cd_op='VOL', p_cmask=tmask_crs, p_ptmask=tmask, &
+         &         p_pfield3d_1=zfse3t, p_pfield3d_2=rke, p_cfield3d=rke_crs )
+      rke_crs(:,:,:) = rke_crs(:,:,:) * tmask_crs(:,:,:)
+      CALL crs_iom_put( "eken_crs", pv_r3d=rke_crs )
+   
       !