Merge-Hydro Branch

The Hydrology branch has been included into the svn trunk of ORCHIDEE. See information here.

This version will merge the version of ORCHIDEE developed at LMD with the 11-layer hydrological scheme into the latest version of the trunk ( Based on the differences between the LMD version and those from which it originates (either 1_9 or 1_9_4_1) (see the attachments), several notes have been written in order to prepare the merge. Here below are the links to these notes.

The hydrology group is : isabelle.gouttevin, jan.polcher, aducharne, aurelien.campoy, matthieu.guimberteau, gerhard.krinner,, catherine.ottle, pierre.brender,, frederique.cheruy, marie-alice.foujols, patricia.cadule.
Everybody with a svn login can(since 12 april 2012) download this branch but only people in the group can commit changes.

How to download ORCHIDEE Hydrology branch for use in offline configuration

  1. Get modipsl
    > svn co modipsl
  1. Get the ORCHIDEE TOOLS-box. You will be prompted for username and password. Use personal svn login (firstname.lastname) or old sechiba cvs login.
    > cd modipsl/util
    > svn co svn://
  2. Get the hydrology version. svn_login can be your personal svn login or sechiba
    > TOOLS/recup_my_ORCHIDEE svn_login branches/Hydrology

How to download ORCHIDEE Hydrology branch for use in coupled LMDZOR configuration

  1. Get modipsl
    > svn co modipsl
  1. Modifiy in modipsl/util/mod.def to get ORCHIDEE Hydrology sources instead of default.
    In modipsl/util/mod.def change following line :
    #-C- LMDZOR_v4  ORCHIDEE                    orchidee_1_9_5   2  .   modeles
    #-C- LMDZOR_v4  branches/Hydrology/ORCHIDEE HEAD             15 .  modeles
  1. Extract configuration LMDZOR_v4 and compile, for more information see wiki LMDZOR_v4
    cd modipsl/util
    ./model LMDZOR_v4       => you will be asked for passwd for sechiba
    cd ../config/LMDZOR_v4
  1. Some changes are needed to run this version using Choisnel 2 layers. No default set up is done yet for using CWRR.
    • Add in PARAM/orchidee.def :
    • Add in orchidee.card, in section [InitialStateFiles] following line (and add ", \" at line above):
      (${R_INIT}/SRF/${config_UserChoices_TagName}/, .)
    • NB! For the moment this version bug with default set up using RIVER_ROUTING=y. To have it runnig, change in orchidee.def to RIVER_ROUTING=n.

Notes to prepare the merge

Note on CWRR

Note on evapnu

Note on dpu

Note on flag couple

Note on logz0

Note on LAI

Note on Soil Map

Note on Interpolations

Martial notes for the merge

Testing the Merge-Hydro version

Matthieu Guimberteau: Off-line tests in the Amazon with bug corrections

Aurélien Campoy (Feb 2012) - test peformed off-line using 50 years of SAFRAN met forcing for the SIRTA site, without routing, without STOMATE : compares reasonably well with the CWRR version of ORCHIDEE before the merge

Jan Polcher (March 2012) - Off-line simulation over Europe using the WATCH forcing data set : simulation performed, to be analyzed

Nicolas Vuichard (March 2012) - off-line tests at FLUXNET sites using STOMATE : ongoing. See details

Fabienne Maignan (Feb-Jul 2012) - global off-line simulation

  • forcing: ERA-Interim regular grid 0.72° 3-hourly
  • computer: LSCE/obelix*
  • 16 procs
    • segmentation fault SIGSEGV here:
       Projection arrays for Slope map                      : 
       nbvmax =           16
       We will work with      1000105  points of the fine grid
       Aggregate_2d : Slope map                     
       aggregate_2D nbvmax =           16 max used =           16
      The problem is in slowproc_slope after the call to aggregate_vec:
          DO ib = 1, nbpt
            !- Reinfiltration coefficient due to the slope: Calculation with parameteres maxlope_ro 
            slopecoef = zero
      !@BUG SOURCE      DO WHILE ( sub_area(ib,idi) > zero ) !@END BUG SOURCE 
             DO WHILE ( sub_area(ib,idi) > zero .AND. idi .LE. nbvmax) 
               ip = sub_index(ib,idi,1)
               jp = sub_index(ib,idi,2)
               slopecoef = slopecoef + MIN(slopemap(ip,jp)/slope_noreinf, un) * sub_area(ib,idi)
               idi = idi +1
      idi may become greater than nbvmax, which causes the segmentation fault. The problem is resolved by the correction proposed in the above code.
    • stop in hydrol_soil_infilt:
      Error in the calculation of infilt tot -4.176114610371373E-008
       k, ji, jst, mc   1.17784270216726        2.16374183793683             2394
                 1  0.277241865204730
      FATAL ERROR FROM ROUTINE hydrol_soil_infilt
       --> We will STOP after hydrol_soil_infilt.
      Fatal error from IOIPSL. STOP in ipslerr with code
      This is caused by slightly negative values in the ERA-Interim precipitations, as already mentionned by Nicolas Vuichard. Proposed (and validated) modification in intersurf_main_2d and intersurf_main_1d:
          WHERE(zprecip_rain(:) .LT. 0.)
          IF (check_INPUTS) THEN
  • test option CHECK_CWRR: OK (no stop, much longer ~*3)
  • Comparison Choisnel against 196: on-going

Frédérique Cheruy - global simulation coupled to LMDz : to be done

Last modified 6 years ago Last modified on 03/26/15 10:42:47

Attachments (8)

Download all attachments as: .zip