Opened 9 months ago

Closed 8 months ago

#928 closed defect (fixed)

Interpolating of soils_param.nc file usda/zobler and soilcolor.nc

Reported by: jgipsl Owned by: somebody
Priority: major Milestone: ORCHIDEE 4.2
Component: Anthropogenic processes Version: trunc
Keywords: Cc:

Description (last modified by jgipsl)

Interpolating soils_param_usda.nc file with XIOS was not validated before.

Done in [8026], changes were done to allow reading soiltext variable from one file and soilcolor from another file. This option was already working with standard interpolation method. Among these changes, there are also a missing initialization of sandfraction and siltfraction when interpolating with XIOS, in slowproc.

Still to be done:

  • To read soils_param_usda.nc with XIOS, a manual change in context_input_orchidee.xml is needed. This is because the soils_param.nc and soils_param_usda.nc, the latitudes are not ordered in the same way and the coordinate variables for x and y are not in the file. A permanent modification with adapted input files so that XIOS can read directly the grid from the file has to be done. Currently, following line needs to be changed:
    <domain id="dom_file_soiltext" type="rectilinear">
      <!-- for soils_param_usda.nc <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>
    
  • Reading of soilcolor variable (from soils_param.nc) needs to be validated when read by XIOS. This was not done when validating the interpolations in ticket #812 because in version 2, this variable was not read (alb_bg.nc was read instead).

Attachments (1)

script_add_coordinate_variables.sh (1.3 KB) - added by jgipsl 8 months ago.

Download all attachments as: .zip

Change History (8)

comment:1 Changed 9 months ago by jgipsl

  • Description modified (diff)

comment:2 Changed 9 months ago by xnwang

Last edited 9 months ago by xnwang (previous) (diff)

Changed 8 months ago by jgipsl

comment:3 Changed 8 months ago by jgipsl

The coordinate variables x(x) and y(y) have been added in soils_param.nc and soils_param_usdatop.nc. nco was used with following script (done at irene-skl): script_add_coordinate_variables.sh

The new files are suffixed _coord. They are stored in the shared repository:
${R_IN}/SRF/SOIL/soils_param_zobler_coord.nc
${R_IN}/SRF/SOIL/soils_param_usdatop_coord.nc

Using these new files, the bounds for longitudes and latitudes set in context_input_orchidee.xml is not longer taken into account. This line, will therefor be kept in the context_input_orchidee.xml file to ensure that using the old file, no changes have been done. Following line will be kept:

<generate_rectilinear_domain bounds_lon_start="-180" bounds_lon_end="180" bounds_lat_start="90" bounds_lat_end="-90"/>

But remember, using the old file soils_param_usdatop.nc, the current context_input_orchidee.xml is not correct for the above line. Therefore, using the new files should be done.

Note that the new files were also changed to be in netCDF-4 format before they were put to the shared repository (ncks -k ...).

Last edited 8 months ago by jgipsl (previous) (diff)

comment:4 Changed 8 months ago by jgipsl

[8061] : Now use the new files as default in offline experiments

rev 6520 in IPSLCM7 : use the new files in the coupled IPSLCM7 configuration

comment:5 Changed 8 months ago by jgipsl

The issue in this ticket is now solved for the default soils_param files : soils_param_zobler.nc and soils_param_usdatop.nc. Maybe the other soils_param files needs to be adapted as well. Following files in IGCM/SRF/SOIL must be considered:

SRF/SOIL/soils_param_fao2.nc
SRF/SOIL/soils_param_faotop.nc
SRF/SOIL/soils_param_fao_nc.nc
SRF/SOIL/soils_param_SoilGrids_5km.nc
SRF/SOIL/soils_param_SoilGrids_halfdeg.nc
SRF/SOIL/soils_param_sp_mip.nc
SRF/SOIL/soils_param_usda_with_clay_oxisol.nc
Last edited 8 months ago by jgipsl (previous) (diff)

comment:6 Changed 8 months ago by jgipsl

The files above have been done as well, using the same script. For the files which were in netCDF-4 format, the script didn't work. I transformed them into format classic (using ncks -3 ) before. They are transformed back in to netCDF-4 in the end (ncks -4). All files are added to to same folder in the sharded acount IGCM/SRF/SOIL, and the README is adapted.

comment:7 Changed 8 months ago by jgipsl

  • Milestone set to ORCHIDEE 4.2
  • Resolution set to fixed
  • Status changed from new to closed
Note: See TracTickets for help on using tickets.