source: branches/ORCHIDEE_2_2/ORCHIDEE/src_xml/context_input_orchidee.xml @ 8579

Last change on this file since 8579 was 8579, checked in by antoine.bierjon, 3 weeks ago

Adapted read and interpolation of irrigation input files with XIOS_interpolation.
Added a default value for AEI_SW interpolation (in case of missing value, with and without XIOS), which is chosen to be the average of the GMIA5 dataset from FAO.
Related to ticket #857

  • Property svn:keywords set to Date Revision HeadURL
File size: 25.3 KB
Line 
1<!-- =========================================================================================================== -->
2<!-- ORCHIDEE context for input files                                                                            -->
3<!-- context_input_orchidee.xml : Configuration file for ORCHIDEE for reading of input files using XIOS          -->
4<!-- This file must be included from iodef.xml                                                                   -->
5<!-- All specifications specific for the reading and interpolating of input files are regrouped in this file.    -->
6<!-- Find below sections for file_definition, field_definition, domain_definition, grid_definition and           -->
7<!-- variable_definition.                                                                                        -->   
8<!-- Recent change:  added block to read params_sp_mip.nc (Salma Tafasca and Agnes Ducharne, November 2020)      -->
9<!--                                                                                                             --> 
10<!-- $HeadURL$          -->
11<!-- $Date$                                                      -->
12<!-- $Revision$                                                                                           -->
13<!-- =========================================================================================================== -->
14<context id="orchidee">
15
16
17  <!-- ========================================================================================================= -->
18  <!-- Definition of input files                                                                                 -->
19  <!-- The input files are activated or desactivated from the source code depending on simulation set up.        -->
20  <!-- The variables are here on the input source grid.                                                          -->
21  <!-- Different treatment are done accordig to the file.                                                        -->
22  <!-- Values higher than 1e10 are for some files considered as missing values and the field is then set to 0.   -->
23  <!-- A normalization are sometimes done with the fraction of the grid cell.                                    -->
24  <!-- Sometimes, if only missing values are found for a grid cell, then a default value is set.                 -->
25  <!-- ========================================================================================================= -->
26 
27  <file_definition>
28
29    <file id="soils_param_file" mode="read" output_freq="1y" enabled="true">
30      <field id="soilcolor" name="soilcolor" domain_ref="dom_file_soil" operation="instant" freq_offset="1ts" enabled="false"/>
31      <field id="soiltext" name="soiltext" domain_ref="dom_file_soil" operation="instant" freq_offset="1ts"/>
32    </file>
33
34    <!-- 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
35                  This file is read in the initialization phase only if there is no restart file.
36         For soilbulk and soilph, set values lower than min_sechiba as masked values because missing values are zero in the file.
37        Renormalization will be done later using the interpolated mask.  -->
38    <file id="soilbulk_file" mode="read" output_freq="1y" enabled="true">
39      <field id="soilbulk_mask_read" name="soilbulk" domain_ref="dom_file_soilbulk" operation="instant" freq_offset="1ts" > this &gt; $min_sechiba ? 1 : 0 </field>
40      <field id="soilbulk_read" name="soilbulk" domain_ref="dom_file_soilbulk" operation="instant" freq_offset="1ts" > this &gt; $min_sechiba ? this : 0 </field>
41    </file>
42
43    <!-- Default file name: params_sp_mip.nc for both all soil hydraulic parameters. -->
44    <!--file id="soilparam_file" mode="read" output_freq="1y" enabled="true"-->
45      <!--field id="soilks" name="soilks" domain_ref="dom_file_soil" operation="instant" freq_offset="1ts" > this &gt; $min_sechiba ? 1e10 : 0 </field-->
46      <!--field id="soilnvan" name="soilnvan" domain_ref="dom_file_soil" operation="instant" freq_offset="1ts" > this &gt; $min_sechiba ? 1e10 : 0 </field-->
47      <!--field id="soilavan" name="soilavan" domain_ref="dom_file_soil" operation="instant" freq_offset="1ts" > this &gt; $min_sechiba ? 1e10 : 0 </field-->
48      <!--field id="soilmcr" name="soilmcr" domain_ref="dom_file_soil" operation="instant" freq_offset="1ts" > this &gt; $min_sechiba ? 1e10 : 0 </field-->
49      <!--field id="soilmcs" name="soilmcs" domain_ref="dom_file_soil" operation="instant" freq_offset="1ts" > this &gt; $min_sechiba ? 1e10 : 0 </field-->
50      <!--field id="soilmcfc" name="soilmcfc" domain_ref="dom_file_soil" operation="instant" freq_offset="1ts" > this &gt; $min_sechiba ? 1e10 : 0 </field-->
51      <!--field id="soilmcw" name="soilmcw" domain_ref="dom_file_soil" operation="instant" freq_offset="1ts" > this &gt; $min_sechiba ? 1e10 : 0 </field-->
52    <!--/file-->
53   
54    <file  id="PFTmap_file"  mode="read" output_freq="1y">
55      <field id="maxvegetfrac_read"  name="maxvegetfrac" domain_ref="dom_file_pft" axis_ref="nvm" operation="instant"  freq_offset="1ts"> this &lt; 1e10 ? this : 0 </field>
56      <field id="maxvegetfrac_mask"  name="maxvegetfrac" domain_ref="dom_file_pft" axis_ref="nvm" operation="instant"  freq_offset="1ts"> this &lt; 1e10 ? 1 : 0 </field>
57    </file>
58   
59    <file id="topography_file"  mode="read" output_freq="1y">
60      <!--      <field id="pente"  domain_ref="dom_file_topo"  operation="instant"  freq_offset="1ts" > (this &gt; 1e10) ? 1 : this </field>     ==> In reference, missing value are use for computing, this probably wrong -->
61      <field id="pente_read"       name="pente" domain_ref="dom_file_topo"  operation="instant"  freq_offset="1ts" > ((this/$slope_noreinf) &lt; 1)? this/$slope_noreinf : 1 </field>
62      <field id="mask_pente_read"  name="pente" domain_ref="dom_file_topo"  operation="instant"  freq_offset="1ts" > (this &lt; $min_sechiba)? 0 : 1                         </field>  <!-- create mask -->
63    </file>
64   
65    <file id="albedo_file" mode="read" output_freq="1y">
66      <field id="mask_alb_read"    name="mask"       domain_ref="dom_file_albedo"  operation="instant"  freq_offset="1ts" />
67      <field id="bg_alb_vis_read"  name="bg_alb_vis" domain_ref="dom_file_albedo"  operation="instant"  freq_offset="1ts" >  mask_alb_read &gt; 0.99 ? this : $missing_value  </field>
68      <field id="bg_alb_nir_read"  name="bg_alb_nir" domain_ref="dom_file_albedo"  operation="instant"  freq_offset="1ts" >  mask_alb_read &gt; 0.99 ? this : $missing_value  </field>
69    </file>
70   
71    <file id="lai_file" mode="read" output_freq="1y" time_counter_name="toto" cyclic="true">
72      <field id="lai_read"       name="LAI" grid_ref="lai" operation="instant"  freq_offset="1ts" > this &gt; 20 ? 0 : this </field>
73      <field id="mask_lai_read"  name="LAI" grid_ref="lai" operation="instant"  freq_offset="1ts" > this &gt; 20 ? 0 : 1    </field>  <!-- create mask -->
74    </file>
75
76    <file id="fertilizer_file" mode="read" output_freq="1y" cyclic="true">
77      <field id="N_qt_OTHER_year_read"  name="N_qt_OTHER_year" domain_ref="dom_file_fertilizer" operation="instant"  freq_offset="1ts"/>
78      <field id="N_qt_WRICE_year_read"  name="N_qt_WRICE_year" domain_ref="dom_file_fertilizer" operation="instant"  freq_offset="1ts"/>
79    </file>
80   
81    <file id="bbg_clim_file" mode="read" output_freq="1y" cyclic="true">
82      <field id="flx_co2_bbg_year_read"  name="flx_co2_bbg_year" domain_ref="dom_file_bbg_clim" operation="instant"  freq_offset="1ts"/>
83    </file>
84   
85    <file id="reftemp_file" mode="read" output_freq="1y" cyclic="true">
86      <field id="reftemp_read"  name="temperature" domain_ref="dom_file_reftemp" operation="instant"  freq_offset="1ts"/>
87    </file>
88   
89    <file id="woodharvest_file" mode="read" output_freq="1y" cyclic="true">
90      <field id="woodharvest_read"  name="woodharvest" domain_ref="dom_file_woodharvest" operation="instant"  freq_offset="1ts"/>
91    </file>
92   
93    <file id="irrigmap_file" mode="read" output_freq="1y" cyclic="true">
94      <field id="irrigmap_read"  name="irrig" domain_ref="dom_file_irrigmap" operation="instant"  freq_offset="1ts" > (this &gt; 101) ? $missing_value : this </field>
95    </file>
96   
97    <file id="aei_sw_file" mode="read" output_freq="1y" cyclic="true">
98      <field id="aei_sw_read"  name="aeisw_pct" domain_ref="dom_file_aei_sw" operation="instant"  freq_offset="1ts" > (this &gt; 101) ? $missing_value : this </field>
99    </file>
100   
101       
102    <file id="nudge_moistc" name="nudge_moistc" type="one_file" mode="read" output_freq="1d">
103      <field  id="moistc_in" name="moistc" operation="instant" grid_ref="grid_file_moistc" freq_offset="1ts" > (this &gt; 1e10) ? 0 : this </field>
104      <field  id="mask_moistc_in" name="moistc" operation="instant" grid_ref="grid_file_moistc" freq_offset="1ts" > (this &gt; 1e10) ? 0 : 1 </field>
105    </file>
106   
107    <file id="nudge_snow" name="nudge_snow" type="one_file" mode="read" output_freq="1d">
108      <field  id="snowdz_in" name="snowdz" operation="instant" grid_ref="grid_file_snow" freq_offset="1ts" > (this &gt; 1e10) ? 0 : this </field>
109      <field  id="snowrho_in" name="snowrho" operation="instant" grid_ref="grid_file_snow" freq_offset="1ts" > (this &gt; 1e10) ? 0 : this </field>
110      <field  id="snowtemp_in" name="snowtemp" operation="instant" grid_ref="grid_file_snow" freq_offset="1ts" > (this &gt; 1e10) ? 0 : this </field>
111      <field  id="mask_snow_in" name="snowdz" operation="instant" grid_ref="grid_file_snow" freq_offset="1ts" > (this &gt; 1e10) ? 0 : 1 </field>
112    </file>
113   
114  </file_definition>
115
116
117  <!-- ========================================================================================================= -->
118  <!-- Definition of all field coming from files but now interpolated to the model grid.                         -->
119  <!-- ========================================================================================================= -->
120
121  <field_definition>
122
123    <field_group enabled="false" read_access="true">
124
125      <field_group id="soil_color" enabled="false">
126        <field id="soilcolor0" field_ref="soilcolor" domain_ref="dom_from_soil" > soilcolor == 0 </field>
127        <field id="soilcolor1" field_ref="soilcolor" domain_ref="dom_from_soil" > soilcolor == 1 </field>
128        <field id="soilcolor2" field_ref="soilcolor" domain_ref="dom_from_soil" > soilcolor == 2 </field>
129        <field id="soilcolor3" field_ref="soilcolor" domain_ref="dom_from_soil" > soilcolor == 3 </field>
130        <field id="soilcolor4" field_ref="soilcolor" domain_ref="dom_from_soil" > soilcolor == 4 </field>
131        <field id="soilcolor5" field_ref="soilcolor" domain_ref="dom_from_soil" > soilcolor == 5 </field>
132        <field id="soilcolor6" field_ref="soilcolor" domain_ref="dom_from_soil" > soilcolor == 6 </field>
133        <field id="soilcolor7" field_ref="soilcolor" domain_ref="dom_from_soil" > soilcolor == 7 </field>
134        <field id="soilcolor8" field_ref="soilcolor" domain_ref="dom_from_soil" > soilcolor == 8 </field>
135        <field id="soilcolor9" field_ref="soilcolor" domain_ref="dom_from_soil" > soilcolor == 9 </field>
136      </field_group>
137
138      <field_group id="soil_text" field_ref="soiltext" domain_ref="dom_from_soil" enabled="true" read_access="true">
139        <field id="soiltext1"  > soiltext == 1 </field>
140        <field id="soiltext2"  > soiltext == 2 </field>
141        <field id="soiltext3"  > soiltext == 3 </field>
142        <field id="soiltext4"  > soiltext == 4 </field>
143        <field id="soiltext5"  > soiltext == 5 </field>
144        <field id="soiltext6"  > soiltext == 6 </field>
145        <field id="soiltext7"  > soiltext == 7 </field>
146        <field id="soiltext8"  > soiltext == 8 </field>
147        <field id="soiltext9"  > soiltext == 9 </field>
148        <field id="soiltext10" > soiltext == 10 </field>
149        <field id="soiltext11" > soiltext == 11 </field>
150        <field id="soiltext12" > soiltext == 12 </field>
151        <field id="soiltext13" > soiltext == 13 </field>
152      </field_group>
153
154      <!-- If soilbulk_mask .gt. min_sechiba, then renormalize by dividing with soilbulk_mask else set bulk_default.
155      soilbulk_mask interpolated to the model grid corresponds to the fraction of the grid cell covered by non-masked
156      values in the input file. Where soilbulk_mask=0, no data where found on the source grid. -->
157      <field id="soilbulk_tmp"  field_ref="soilbulk_read" domain_ref="dom_from_soilbulk" />
158      <field id="soilbulk_mask" field_ref="soilbulk_mask_read" domain_ref="dom_from_soilbulk" read_access="true" enabled="true" />
159      <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>
160     
161      <field id="frac_veget"      field_ref="maxvegetfrac_read" domain_ref="dom_from_pft" read_access="true" enabled="true" />
162      <field id="frac_veget_frac" field_ref="maxvegetfrac_mask" domain_ref="dom_from_pft" read_access="true" enabled="true" />
163
164      <field id="slope_tmp"           field_ref="pente_read"       domain_ref="dom_from_topo_read"                                  />
165      <field id="frac_slope_interp"   field_ref="mask_pente_read"  domain_ref="dom_from_topo_read" read_access="true" enabled="true"/>
166      <field id="reinf_slope_interp"                               domain_ref="dom_from_topo" read_access="true" enabled="true"> (frac_slope_interp &gt; $min_sechiba ) ? 1-slope_tmp/(frac_slope_interp+1e-100) : $slope_default </field>
167
168      <field id="bg_alb_vis_tmp"   field_ref="bg_alb_vis_read"  domain_ref="dom_from_albedo"  />
169      <field id="bg_alb_vis_interp" field_ref="bg_alb_vis_tmp"  read_access="true" enabled="true" > this == $missing_value ? $albbg_vis_default : this </field>
170
171      <field id="bg_alb_nir_tmp"    field_ref="bg_alb_nir_read"  domain_ref="dom_from_albedo"  />
172      <field id="bg_alb_nir_interp" field_ref="bg_alb_nir_tmp"  read_access="true"  enabled="true" >  this == $missing_value ? $albbg_nir_default : this </field>
173     
174      <field id="lai_tmp"          field_ref="lai_read"      grid_ref="grid_from_lai" />
175      <field id="frac_lai_interp"  field_ref="mask_lai_read" grid_ref="grid_from_lai" read_access="true" enabled="true" />
176      <field id="lai_interp"       field_ref="lai_tmp"                           read_access="true" enabled="true"  >  lai_tmp/(frac_lai_interp+1e-100) </field>  <!-- renormalize -->
177     
178      <field id="N_qt_WRICE_year_interp"  field_ref="N_qt_WRICE_year_read"  domain_ref="dom_from_fertilizer" read_access="true" enabled="true" />
179      <field id="N_qt_OTHER_year_interp"  field_ref="N_qt_OTHER_year_read"  domain_ref="dom_from_fertilizer" read_access="true" enabled="true" />
180     
181      <field id="flx_co2_bbg_year_interp"  field_ref="flx_co2_bbg_year_read"  domain_ref="dom_from_bbg_clim" read_access="true" enabled="true" />
182
183      <field id="reftemp_interp"  field_ref="reftemp_read"  domain_ref="dom_from_reftemp" read_access="true" enabled="true" />
184
185      <field id="woodharvest_interp"  field_ref="woodharvest_read"  domain_ref="dom_from_woodharvest" read_access="true" enabled="true" />
186
187      <field id="irrigmap_tmp"   field_ref="irrigmap_read"  domain_ref="dom_from_irrigmap"  />
188      <field id="irrigmap_interp"  field_ref="irrigmap_tmp" read_access="true" enabled="true" > this == $missing_value ? 0 : this </field>
189
190      <field id="aei_sw_tmp"   field_ref="aei_sw_read"  domain_ref="dom_from_aei_sw"  />
191      <field id="aei_sw_interp"  field_ref="aei_sw_tmp" read_access="true" enabled="true" > this == $missing_value ? $aei_sw_default : this </field>
192     
193    </field_group>
194
195    <field_group id="orchidee_input_mc" grid_ref="grid_nslm_nstm_read">
196      <!-- Interpolate the mask. This gives now the fraction of the source area underlaying the model grid cell -->
197      <field id="mask_moistc_interp" field_ref="mask_moistc_in" read_access="true" />
198      <!-- Interpolate the variable. Note that where the source file had missing values, the variables was set to 0. Therefor normalization is needed.     -->
199      <!-- Normalize the interpolated values with the mask which corresponds to the fraction of underlaying valid source grid cells.                       -->
200      <!-- The default value 0.2 is used if the mask is 0. This is the case if no information at all was found on the source grid for the model grid cell. -->
201      <field id="moistc_tmp" field_ref="moistc_in"  />
202      <field id="moistc_interp" field_ref="moistc_tmp" read_access="true" > (mask_moistc_interp &lt; 1e-4) ? 0.2 : this/(mask_moistc_interp+1e-100) </field>
203    </field_group>
204
205    <field_group id="orchidee_input_snow" grid_ref="grid_nsnow2" >
206      <field id="mask_snow_interp" field_ref="mask_snow_in" read_access="true" />
207      <field id="snowdz_tmp" field_ref="snowdz_in" />
208      <!-- Set snowdz=0 where no information found on source grid, for all other points normlize with the fracion of underlaying source grid cells      -->
209      <field id="snowdz_interp" field_ref="snowdz_tmp" read_access="true" > (mask_snow_interp &lt; 1e-4) ? 0 : this/(mask_snow_interp+1e-100) </field>
210      <field id="snowrho_tmp" field_ref="snowrho_in" />
211      <!-- Set snowrho=50 where no information found on source grid, for all other points normlize with the fracion of underlaying source grid cells    -->
212      <field id="snowrho_interp" field_ref="snowrho_tmp" read_access="true" > (mask_snow_interp &lt; 1e-4) ? 50 : this/(mask_snow_interp+1e-100) </field>
213      <field id="snowtemp_tmp" field_ref="snowtemp_in" />
214      <!-- Set snowtemp=273 where no information found on source grid, for all other points normlize with the fracion of underlaying source grid cells  -->
215      <field id="snowtemp_interp" field_ref="snowtemp_tmp" read_access="true" > (mask_snow_interp &lt; 1e-4) ? 273 : this/(mask_snow_interp+1e-100) </field>
216    </field_group>
217 
218  </field_definition>
219
220
221
222  <!-- =========================================================================================================== -->
223  <!-- Definition of a horizontal domain for each input file                                                       -->
224  <!-- 2 different domains are specified for each input file X: dom_file_X and dom_from_X                          -->
225  <!--                                                                                                             -->
226  <!-- dom_file_X is the horizontal original grid in the input file. In this domain, specifications to help XIOS   -->
227  <!-- read the grid can be added such as generate_rectilinear_domain                                              -->
228  <!--                                                                                                             -->
229  <!-- dom_from_X describes how to make the interpolation from the file to the model grid. This domain contains    -->
230  <!-- the interpolation methode specific to the file.                                                             -->
231  <!--                                                                                                             -->
232  <!-- If the variable contains more dimensions than the horizontal domain and 1 axis, 2 grids must be specified:  -->
233  <!-- One for the orignial grid in the input file and one for the interpolation to model grid. These two grids    -->
234  <!-- uses 2 domains specified as above.                                                                          -->
235  <!-- We here suppose that the interpolation is only done on the horizontal domain                                -->
236  <!-- =========================================================================================================== -->
237
238  <domain_definition>
239
240    <domain id="dom_file_moistc" type="rectilinear">
241      <generate_rectilinear_domain />
242    </domain>
243   
244    <domain id="dom_file_snow" type="rectilinear">
245      <generate_rectilinear_domain />
246    </domain>
247   
248    <domain id="dom_from_soil" domain_ref="domain_landpoints">
249      <interpolate_domain order="1" renormalize="false"/>
250    </domain>
251
252    <domain id="dom_from_soilbulk" domain_ref="domain_landpoints">
253    <interpolate_domain order="1" renormalize="false" />
254    </domain>
255
256    <domain id="dom_from_pft" domain_ref="domain_landpoints">
257      <interpolate_domain order="1" renormalize="false" />
258    </domain>
259   
260    <domain id="dom_from_topo" domain_ref="domain_landpoints">
261    </domain>
262   
263    <domain id="dom_from_topo_read" domain_ref="domain_landpoints">
264      <interpolate_domain order="1" renormalize="false" />
265    </domain>
266
267    <domain id="dom_from_albedo" domain_ref="domain_landpoints">
268      <interpolate_domain order="1" renormalize="true" detect_missing_value="true"/>
269    </domain>
270   
271    <domain id="dom_from_lai" domain_ref="domain_landpoints">
272      <interpolate_domain order="1" renormalize="false" />
273    </domain>
274   
275    <domain id="dom_from_bbg_clim" domain_ref="domain_landpoints">
276      <interpolate_domain order="1" renormalize="false" />
277    </domain>
278   
279    <domain id="dom_from_fertilizer" domain_ref="domain_landpoints">
280      <interpolate_domain order="1" renormalize="false" />
281    </domain>
282   
283    <domain id="dom_from_reftemp" domain_ref="domain_landpoints">
284      <interpolate_domain order="1" renormalize="false" />
285    </domain>
286   
287    <domain id="dom_from_woodharvest" domain_ref="domain_landpoints">
288      <interpolate_domain order="1" renormalize="false" />
289    </domain>
290   
291    <domain id="dom_from_irrigmap" domain_ref="domain_landpoints">
292      <interpolate_domain order="1" renormalize="true" detect_missing_value="true" />
293    </domain>
294   
295    <domain id="dom_from_aei_sw" domain_ref="domain_landpoints">
296      <interpolate_domain order="1" renormalize="true" detect_missing_value="true" />
297    </domain>
298   
299    <domain id="dom_file_soil" type="rectilinear">
300      <!-- generate_rectilinear_domain bounds_lon_start="-180" bounds_lon_end="180" bounds_lat_start="-90" bounds_lat_end="90"/-->
301      <generate_rectilinear_domain bounds_lon_start="-180" bounds_lon_end="180" bounds_lat_start="90" bounds_lat_end="-90"/>
302    </domain>
303
304    <domain id="dom_file_soilbulk" type="rectilinear">
305    <generate_rectilinear_domain bounds_lon_start="-180" bounds_lon_end="180" bounds_lat_start="90" bounds_lat_end="-90"/>/>
306    </domain>
307   
308    <domain id="dom_file_pft" type="rectilinear">
309      <generate_rectilinear_domain/>
310    </domain>
311   
312    <domain id="dom_file_topo" type="rectilinear">
313      <generate_rectilinear_domain/>
314    </domain>
315   
316    <domain id="dom_file_albedo" type="rectilinear">
317      <generate_rectilinear_domain />
318    </domain>
319   
320    <domain id="dom_file_lai" type="rectilinear">
321      <generate_rectilinear_domain/>
322    </domain>
323   
324    <domain id="dom_file_bbg_clim" type="rectilinear">
325      <generate_rectilinear_domain/>
326    </domain>
327   
328    <domain id="dom_file_fertilizer" type="rectilinear">
329      <generate_rectilinear_domain/>
330    </domain>
331
332    <domain id="dom_file_reftemp" type="rectilinear">
333      <generate_rectilinear_domain bounds_lon_start="-180" bounds_lon_end="180" bounds_lat_start="90" bounds_lat_end="-90"/>
334    </domain>
335
336    <domain id="dom_file_woodharvest" type="rectilinear">
337      <generate_rectilinear_domain/>
338    </domain>
339
340    <domain id="dom_file_irrigmap" type="rectilinear">
341      <generate_rectilinear_domain/>
342    </domain>
343
344    <domain id="dom_file_aei_sw" type="rectilinear">
345      <generate_rectilinear_domain/>
346    </domain>
347   
348  </domain_definition>
349
350
351  <!-- =========================================================================================================== -->
352  <!-- Definition of input grids and interpolated grids                                                            -->
353  <!-- =========================================================================================================== -->
354
355  <grid_definition>
356
357    <grid id="grid_file_moistc">
358      <domain domain_ref="dom_file_moistc"/>
359      <axis axis_ref="nslm" />
360      <axis axis_ref="nstm" />
361    </grid> 
362
363    <grid id="grid_file_snow">
364      <domain domain_ref="dom_file_snow"/>
365      <axis axis_ref="nsnow" />
366    </grid> 
367   
368    <grid id="grid_from_lai">
369      <domain domain_ref="dom_from_lai" />
370      <axis axis_ref="nvm" />
371      <axis axis_ref="lai_time" />
372    </grid>
373   
374    <grid id="lai">
375      <domain domain_ref="dom_file_lai" />
376      <axis axis_ref="nvm" />
377      <axis axis_ref="lai_time" />
378    </grid>
379   
380    <grid id="grid_nsnow2">
381      <domain domain_ref="domain_landpoints" >
382        <interpolate_domain order="1"/>
383      </domain>
384      <axis axis_ref="nsnow" />
385    </grid>
386
387  </grid_definition>
388 
389 
390  <!-- =========================================================================================================== -->
391  <!-- Variable denfinition needed for interpolation of input files                                                -->
392  <!-- The values for below variables are set from the fortran code.                                               -->
393  <!-- =========================================================================================================== -->
394
395  <variable_definition>
396   
397    <variable id="slope_default" type="double">xxx</variable>
398    <variable id="slope_noreinf" type="double">xxx</variable>
399    <variable id="min_sechiba" type="double">xxx</variable>
400    <variable id="albbg_vis_default" type="double">xxx</variable>
401    <variable id="albbg_nir_default" type="double">xxx</variable>
402    <variable id="bulk_default" type="double">xxx</variable>
403    <variable id="aei_sw_default" type="double">xxx</variable>
404   
405  </variable_definition>
406 
407
408</context>
Note: See TracBrowser for help on using the repository browser.