Changeset 7432 for branches/ORCHIDEE_2_2/ORCHIDEE/src_parameters/constantes.f90

Timestamp:

2022-01-17T19:18:33+01:00 (2 years ago)

Author:

agnes.ducharne

Message:

Changes to make the IMPOSE_SOILT mode functional. This mode is not anymore dependent on IMPOSE_VEG, so we can impose soil
properties even when we read a vegetation map. A new output variable is added to export "ksref" before being vertically modified.

With IMPOSE_SOILT, we can either impose a certain USDA texture (via SOIL_FRACTIONS) and this propagates to all the soil parameters, but we can as well impose specific parameters in isolation (with a default txture as Loam). This committ was tested with various configs, and it also works if we use a restart file, and we can even change the soil texture map in this case.

File:

• branches/ORCHIDEE_2_2/ORCHIDEE/src_parameters/constantes.f90 (modified) (3 diffs)

branches/ORCHIDEE_2_2/ORCHIDEE/src_parameters/constantes.f90

@@ -214 +214 @@
     CALL getin_p('IMPOSE_VEG', impveg)
 
-    IF (impveg) THEN
+    !IF (impveg) THEN
        !Config Key   = IMPOSE_SOILT
        !Config Desc  = Should the soil type be prescribed ?
@@ -225 +225 @@
        !Config Units = [FLAG]
        CALL getin_p('IMPOSE_SOILT', impsoilt)     
-    ENDIF
+    !ENDIF
 
     !Config Key   = LAI_MAP
@@ -1004 +1004 @@
     !-
     !
-    !Config Key   = CLAYFRACTION_DEFAULT
-    !Config Desc  = default fraction of clay
-    !Config If    = OK_SECHIBA 
-    !Config Def   = 0.2 
-    !Config Help  = 
-    !Config Units = [-]   
-    CALL getin_p('CLAYFRACTION_DEFAULT',clayfraction_default)
-    !
-    !Config Key   = SANDFRACTION_DEFAULT
-    !Config Desc  = default fraction of sand
-    !Config If    = OK_SECHIBA 
-    !Config Def   = 0.4 
-    !Config Help  = 
-    !Config Units = [-]   
-    CALL getin_p('SANDFRACTION_DEFAULT',sandfraction_default)
-    !
-    !Config Key   = SILTFRACTION_DEFAULT
-    !Config Desc  = default fraction of silt
-    !Config If    = OK_SECHIBA 
-    !Config Def   = 0.4 
-    !Config Help  = 
-    !Config Units = [-]   
-    CALL getin_p('SILTFRACTION_DEFAULT',siltfraction_default)
-
-
-    IF ( ABS(clayfraction_default+sandfraction_default+siltfraction_default-1) > min_sechiba) THEN
-       WRITE(numout,*) 'Incoherence found. clayfraction_default=', clayfraction_default, ' sandfraction_default=',&
-            sandfraction_default,' siltfraction_default=',siltfraction_default
-       CALL ipslerr_p(3,"config_sechiba_parameters",&
-            "Inconsistecy between CLAYFRACTION_DEFAULT, SANDFRACTION_DEFAULT and SILTFRACTION_DEFAULT set in run.def",&
-            "The sum should be equal 1 but this is not the case.","Modify run.def and restart the model")
-    END IF
-    !
     !Config Key   = MIN_VEGFRAC 
     !Config Desc  = Minimal fraction of mesh a vegetation type can occupy