Ignore:
Timestamp:
2022-01-17T19:18:33+01:00 (2 years ago)
Author:
agnes.ducharne
Message:

Changes to make the IMPOSE_SOILT mode functional. This mode is not anymore dependent on IMPOSE_VEG, so we can impose soil
properties even when we read a vegetation map. A new output variable is added to export "ksref" before being vertically modified.

With IMPOSE_SOILT, we can either impose a certain USDA texture (via SOIL_FRACTIONS) and this propagates to all the soil parameters, but we can as well impose specific parameters in isolation (with a default txture as Loam). This committ was tested with various configs, and it also works if we use a restart file, and we can even change the soil texture map in this case.

Location:
branches/ORCHIDEE_2_2/ORCHIDEE/src_xml
Files:
2 edited

Legend:

Unmodified
Added
Removed

  • branches/ORCHIDEE_2_2/ORCHIDEE/src_xml/field_def_orchidee.xml

    r7265 r7432  
    349349    <field id="fldcapacity" name="fldcapacity" long_name="Field capacity (in % of soil volume)" grid_ref="grid_nslm" unit="%" operation="once" > mcfc*100 </field> 
    350350    <field id="ksat" name="ksat" long_name="Saturated Hydraulic conductivity (at each hydrologic node)" unit="10^(-6) m/s" grid_ref="grid_nslm" operation="once" freq_offset="0ts"> this*1000/86400 </field> <!-- ym check this value. jg why ?  --> 
     351    <field id="ksref" name="ksref" long_name="Reference saturated hydraulic conductivity (before imposing a vertical profile)" unit="mm/d" operation="once" freq_offset="0ts" />  
    351352    <field id="alphavg" name="alphavg" long_name="Van Genuchten alpha parameter (at each hydrologic node)" unit="m(-1)" grid_ref="grid_nslm" operation="once" freq_offset="0ts"/>  
    352353    <field id="nvg" name="nvg" long_name="Van Genuchten n parameter (at each hydrologic node)" unit="-" grid_ref="grid_nslm" operation="once" freq_offset="0ts"/> 

  • branches/ORCHIDEE_2_2/ORCHIDEE/src_xml/file_def_orchidee.xml

    r7265 r7432  
    5656    <field field_ref="mrsow" level="1"/> 
    5757    <field field_ref="ksat" grid_ref="grid_nslm_out" level="1"/> 
     58    <field field_ref="ksref" level="1"/> 
    5859    <field field_ref="precip_rain" level="0" name="rain" unit="mm/d" > @precip_rain*86400  </field> 
    5960    <field field_ref="precip_snow" level="0" name="snowf" unit="mm/d" > @precip_snow*86400  </field> 
Note: See TracChangeset for help on using the changeset viewer.