Fiche : dpu

Description

Dans constantes_soil.f90 sont définies les 'dpu' utilisés pour l'hydro 11 couches (hydrol.f90). Dans la version 1.9.4.1, on utilise 1 valeur de dpu par type de sol (dpu=vecteur de 3 élts). Dans la version LMD, on utilise le paramètre (scalaire) dpu_max.

Choix

On souhaite travailler avec une valeur de dpu identique quelque soit le type de sol. On retient la formulation de la version LMD. Doit-on unifier dpu_cste (utilisé pour choisnel) et dpu_max (utilisé pour le 11 couches) et n'en faire qu'un seul paramètre ? A moyen terme, on souhaitera travailler avec des valeurs spatialisées de dpu et donc lire une carte. Didier Solyga a initié ce travail de spatialisation de certains paramètres. Cette spatialisation ne sera pas considéré dans la version mergée 11-couches.

Mise en oeuvre

• dans constantes_soil.f90
  • suppression de dpu

    ! Total depth of soil reservoir (m)
    REAL(r_std),PARAMETER,DIMENSION(nstm) :: dpu = &
     & (/ 2.0_r_std, 2.0_r_std, 2.0_r_std /)
    ! dpu must be constant over the different soil types
    

• ajout de dpu_max

    REAL(r_std),SAVE :: dpu_max=2.0_r_std
  

• dans hydrol.f90
  • remplacement des variables dpu(jst) par dpu_max

Tests

none