Changeset 3003 for LMDZ5/trunk/DefLists

Timestamp:

Oct 3, 2017, 4:45:26 PM (7 years ago)

Author:

Laurent Fairhead

Message:

Modifications to the code and xml files to output Ap and B, the coefficients
of the hybrid coordinates as requested by the CMIP6 DataRequest?
LF (with guidance from A. Caubel and S. Senesi)

Location:

LMDZ5/trunk/DefLists

Files:

• context_lmdz.xml (modified) (2 diffs)
• field_def_lmdz.xml (modified) (1 diff)
• file_def_histmth_lmdz.xml (modified) (1 diff)

LMDZ5/trunk/DefLists/context_lmdz.xml

@@ -48 +48 @@
     <axis id="Bhyb" standard_name="Bhyb comp of Hyb Cord" unit="">
     </axis>
+    <axis id="Ahyb_inter" standard_name="A comp of Hyb Cord at interface" unit="Pa">
+    </axis> 
+    <axis id="Bhyb_inter" standard_name="B comp of Hyb Cord at interface" unit="">
+    </axis>
     <axis id="Alt" standard_name="Height approx for scale heigh of 8km at levels" unit="km">
     </axis>
     <axis id="plev" standard_name="model_level_number" unit="Pa">
+    </axis>
+    <axis id="klev"  prec="8" long_name = "hybrid sigma pressure coordinate"
+        standard_name ="atmosphere_hybrid_sigma_pressure_coordinate" unit="1">
+    <!-- <variable name="formula" type="string"> p = ap + b*ps </variable> -->
+    <!-- <variable name="formula_terms" type="string"> ap: ap b: b ps: ps </variable> -->
+    </axis>
+
+    <axis id="bnds" standard_name="bounds" unit="1" >
     </axis>
 
@@ -84 +96 @@
     <grid id="grid_scalar">
     </grid>
+
+   <grid id="klev_bnds"> <axis axis_ref="klev" /> <axis axis_ref="bnds" /> </grid>
 
      <grid id="grid_glo">

LMDZ5/trunk/DefLists/field_def_lmdz.xml

@@ -7 +7 @@
 <!--    </field_group>  -->
     
+    <field_group id="coord_hyb">
+        <field id="Ahyb"      axis_ref="klev" />
+        <field id="Ahyb_inter" grid_ref="klev_bnds" />
+        <field id="Bhyb"       axis_ref="klev"  />
+        <field id="Bhyb_inter"  grid_ref="klev_bnds" />
+    </field_group>
+
     <field_group id="fields_2D" grid_ref="grid_glo">
         <field id="phis"    long_name="Surface geop.height"     unit="m2/s2" />

LMDZ5/trunk/DefLists/file_def_histmth_lmdz.xml

@@ -9 +9 @@
 <!--              <field field_ref="Alt" level="1" />  -->
 <!--            </field_group>  -->
+
+            <field_group operation="once">
+              <field field_ref="Ahyb" level="1" />
+              <field field_ref="Ahyb_inter" level="1" />
+              <field field_ref="Bhyb" level="1" />
+              <field field_ref="Bhyb_inter" level="1" />
+            </field_group>
 
             <field_group operation="instant">