Changeset 8418 for branches/ORCHIDEE_2_2/ORCHIDEE/src_xml/context_input_orchidee.xml

Timestamp:

2024-02-12T19:49:09+01:00 (5 months ago)

Author:

bertrand.guenet

Message:

The Moyano function describing the soil moisture effect on OM decomposition is added

File:

• branches/ORCHIDEE_2_2/ORCHIDEE/src_xml/context_input_orchidee.xml (modified) (5 diffs)

branches/ORCHIDEE_2_2/ORCHIDEE/src_xml/context_input_orchidee.xml

@@ -30 +30 @@
       <field id="soilcolor" name="soilcolor" domain_ref="dom_file_soil" operation="instant" freq_offset="1ts" enabled="false"/>
       <field id="soiltext" name="soiltext" domain_ref="dom_file_soil" operation="instant" freq_offset="1ts"/>
+    </file>
+
+    <!-- Default file name: soil_bulk_and_ph.nc for both soilbulk and soilph. The variables can be read from different files, specified in run.def
+                  This file is read in the initialization phase only if there is no restart file.
+         For soilbulk and soilph, set values lower than min_sechiba as masked values because missing values are zero in the file.
+        Renormalization will be done later using the interpolated mask.  -->
+    <file id="soilbulk_file" mode="read" output_freq="1y" enabled="true">
+      <field id="soilbulk_mask_read" name="soilbulk" domain_ref="dom_file_soilbulk" operation="instant" freq_offset="1ts" > this &gt; $min_sechiba ? 1 : 0 </field>
+      <field id="soilbulk_read" name="soilbulk" domain_ref="dom_file_soilbulk" operation="instant" freq_offset="1ts" > this &gt; $min_sechiba ? this : 0 </field>
     </file>
 
@@ -134 +143 @@
         <field id="soiltext13" > soiltext == 13 </field>
       </field_group>
+
+      <!-- If soilbulk_mask .gt. min_sechiba, then renormalize by dividing with soilbulk_mask else set bulk_default.
+      soilbulk_mask interpolated to the model grid corresponds to the fraction of the grid cell covered by non-masked
+      values in the input file. Where soilbulk_mask=0, no data where found on the source grid. -->
+      <field id="soilbulk_tmp"  field_ref="soilbulk_read" domain_ref="dom_from_soilbulk" />
+      <field id="soilbulk_mask" field_ref="soilbulk_mask_read" domain_ref="dom_from_soilbulk" read_access="true" enabled="true" />
+      <field id="soilbulk"      field_ref="soilbulk_tmp"  domain_ref="dom_from_soilbulk" read_access="true" enabled="true"  > soilbulk_mask &gt; $min_sechiba ? soilbulk_tmp/(soilbulk_mask+1e-100) : $bulk_default </field>
       
       <field id="frac_veget"      field_ref="maxvegetfrac_read" domain_ref="dom_from_pft" read_access="true" enabled="true" />
@@ -220 +236 @@
     </domain>
 
+    <domain id="dom_from_soilbulk" domain_ref="domain_landpoints">
+    <interpolate_domain order="1" renormalize="false" />
+    </domain>
+
     <domain id="dom_from_pft" domain_ref="domain_landpoints">
       <interpolate_domain order="1" renormalize="false" />
@@ -258 +278 @@
       <!-- generate_rectilinear_domain bounds_lon_start="-180" bounds_lon_end="180" bounds_lat_start="-90" bounds_lat_end="90"/-->
       <generate_rectilinear_domain bounds_lon_start="-180" bounds_lon_end="180" bounds_lat_start="90" bounds_lat_end="-90"/>
+    </domain>
+
+    <domain id="dom_file_soilbulk" type="rectilinear">
+    <generate_rectilinear_domain bounds_lon_start="-180" bounds_lon_end="180" bounds_lat_start="90" bounds_lat_end="-90"/>/>
     </domain>
     
@@ -346 +370 @@
     <variable id="albbg_vis_default" type="double">xxx</variable>
     <variable id="albbg_nir_default" type="double">xxx</variable>
+    <variable id="bulk_default" type="double">xxx</variable>
     
   </variable_definition>