wiki:DevelopmentActivities/MergeHydro/Soilmap_note

Version 5 (modified by mmaipsl, 10 years ago) (diff)

--

Introduce new soil map and soil classification

Description

Parameters for soil type distribution.

L'ancien nstm devient ncsm.

Default soil texture distribution in the following order : COARSE, MEDIUM, FINE ! For FAO Classification

INTEGER(i_std),PARAMETER
nscm_fao=3
! For USDA Classification
INTEGER(i_std),PARAMETER
nscm_usda=12
INTEGER(i_std), SAVE
nscm=nscm_fao

Affected modules

constantes_soil, slowproc, hydrol, hydrolc ?

standard Zobler soils_param.nc : Transfered from the ISLSCP CD-rom on 19990916 at 165208.686

  • soilcolor : soil color types from the Henderson-Sellers and Wilson dataset.
  • soiltext : Soil texture from Zobler 86

soils_param_zobler.nc

Ce qui est déjà dans les dépôts : /dmnfs/cont003/p86ipsl/IGCM/INIT/SRF/OL2:

  • -rw-r--r-- 2 p86ipsl dsm 1037820 Nov 19 2003 soils_param.nc
    lrwxrwxrwx 1 p86ipsl dsm 14 Oct 20 2007 soils_param_zobler.nc -> soils_param.nc
  • -rw-r--r-- 2 p86ipsl dsm 112000952 Feb 9 2006 soils_param_faotop.nc
    lrwxrwxrwx 1 p86ipsl dsm 21 Oct 20 2007 soils_param_fao.nc -> soils_param_faotop.nc
  • -rw-r--r-- 2 p86ipsl dsm 149334740 Sep 21 2006 soils_param_fao2.nc
  • -rw-r--r-- 2 p86ipsl dsm 112000956 Feb 9 2006 soils_param_usdatop.nc
    lrwxrwxrwx 1 p86ipsl dsm 22 Oct 20 2007 soils_param_usda.nc -> soils_param_usdatop.nc

Programmation

       soilclass=val_exp
       CALL setvar (soilclass, val_exp, 'SOIL_FRACTIONS', soilclass_default)
       njsc(:) = 0
       soiltile(:,:) = zero
       DO ji = 1, kjpindex
          njsc(ji) = MAXLOC(soilclass(ji,:),1)
          soiltile(:,1) = SUM(frac_nobio(ji,:))
       ENDDO
       DO jv = 1, nvm
          jst = pref_soil_veg(jv)
          DO ji = 1, kjpindex
             soiltile(ji,jst) = soiltile(ji,jst) + veget(ji,jv)
          ENDDO
       ENDDO

tests