= 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 [[BR]]
   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 [[BR]]
   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 [[BR]]
   lrwxrwxrwx 1 p86ipsl dsm        22 Oct 20  2007 soils_param_usda.nc -> soils_param_usdatop.nc

== Programmation ==
{{{
    !Config  Key  = 'PREF_SOIL_VEG' 
    !Config  Desc = The soil tile number for each vegetation
    !Config  Def  = 1, 2, 2, 2, 2, 2, 2, 2, 2, 3, 3, 3, 3
    !Config  Help = Gives the number of the soil tile on which we will
    !Config         put each vegetation. This allows to divide the hydrological column
    pref_soil_veg = (/ 1, 2, 2, 2, 2, 2, 2, 2, 2, 3, 3, 3, 3 /)
    CALL getin_p('PREF_SOIL_VEG', pref_soil_veg)

[...]
       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(ji,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
}}}
Le pref_soil_veg perd donc une dimension par rapport à son ancienne définition :
{{{
    ! Select the preferences for the distribution of the 13 PFTs on the 3 soil types.
    !
    vegsoil_dist='EQUI'
    !
    SELECT CASE(vegsoil_dist)
       !
       ! A reasonable choice
       !
    CASE('MAXR')
       pref_soil_veg(:,1)  = (/ 1, 3, 3, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2 /)
       pref_soil_veg(:,2)  = (/ 2, 2, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3 /)
       pref_soil_veg(:,3)  = (/ 3, 1, 1, 1, 1, 1 ,1 ,1 ,1 ,1 ,1 ,1, 1 /)
       !
       ! Current default : equidistribution.
       !
    CASE('EQUI')
       !
       pref_soil_veg(:,:) = zero
       !
    CASE DEFAULT
       !
       WRITE(numout,*) 'The vegetation soil type preference you have chosen does not exist'
       WRITE(numout,*) 'You chose :', vegsoil_dist
       STOP 'slowproc_init'
       !
    END SELECT
}}}
== tests ==