Branch ORCHIDEE-DYNAMICO

This branch is created in order to describe the work done to couple ORCHIDEE to the new atmospheric dynamical core developed by LMD and leaded by Thomas Dubos called DYNAMICO ​http://forge.ipsl.jussieu.fr/dynamico.

Working group : Lluis Fita, Yann Meurdesoif, Jan Polcher, Thomas Dubos, Josefine Ghattas

Branch created 18/01/2016 as a copy of ORCHIDEE/trunk revision 3115.

Planned work

1. XIOS2 and interpweight
   • Release of XIOS2
   • Grouping of current interpolation (read, aggregate_p, apply weights and bring back fractions) in ORCHIDEE into a new module (`interpweight')
2. Commit and test both steps on the trunk [done October 2016]
3. Open a new branch for the interpolation with XIOS2
4. Work on the interpolation (doing nothing for the calendar) with a new `interpweight_XIOS' but with XIOS2
5. Validation off-line
6. Implement it on the current version on the trunk
7. Coupling to DYNAMICO

interpweight

This module is used in the re-organization done in the section of the code responsible of the interpolation of the morphology datat (PFT, albedo, LAI, ...) from different files done by L. Fita in the mark of the HEAT project. Previous version of the interpolation of these files, proceed by three basic steps common for all the files:

1. Modification of values from file
2. Creation of land/sea mask using values from file
3. Interpolate values

This module serves to generalize all these steps and construct a common interface independently of the file. Each step is generalize and the different processes almost specific for each file in the previous form are introduced as different methodologies. Thus now it provides:

1. Modification of initial values: Removing that wrong values from the file
   • `interpweight_modifying_input[1/2/3/4]D': subroutines to modify 1D, 2D, 3D or 4D input data from file, using different methods: none
2. Masking input: Compute ’on fly’ the values of the land/sea mask taking the values in the file
   • `interpweight_masking_input[1/2/3/4]D': subroutines to compute the mask using data from input file using different methods: nomask, mbelow, mabove, msumrange, var
3. Area weights: Get coincident areas from input file to the target projection using `aggregate_p'
   • No changes on it
4. Interpolate: Use obtained areas and interpolate values at each grid point at the same format as it will be provided by XIOS
   • `interpweight_provide_fractions[1/2/3/4]D': Perform interpolation for fraction/category data
   • `interpweight_provide_interpolation[2/4]D': Perform interpolation for continuous data
   • `variableusetypes': Variable to provide the values along the additional dimension to perform the interpolation. (e.g.: in case the 13 pfts are non consecutive: 1,3,6,18,23,34,35,39,48,...)
5. A new variable is added which explains the ‘availability’ of data to perform the interpolation (0, 1). When it is negative it means that there was no data for that grid point

   availability = SUM(area_source)/Area_target_grid_point
   

Availability information for the vegetation map `maxvegetfrac' found in PFTmap.nc

A document with more detail and the working notes is available from here interpweight working notes​ A short presentation with the general aims of the module is also available here interpweight: Re-organization of ORCHIDEE’s interpolation​

This module will disappear once interpolation will be done throughout XIOS

Ciclad compilation

Compilation of the DYNAMICO-ORCDHIEE system in CICLAD is available here ORCHIDEE-DYNAMICO/DynOrCiclad