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

Timestamp:

2022-07-20T11:30:43+02:00 (2 years ago)

Author:

josefine.ghattas

Message:

Integrated new irrigation scheme developed by Pedro Arboleda. See ticket #857
This corresponds to revsion 7708 of version pedro.arboleda/ORCHIDEE. Following differences were made but were not made on the pedro.arboleda/ORCHIDEE :

• argumet place in call to routing_wrapper_intialize changed order
• lines with only change in space were not taken
• some indentation changed
• set irrigation output as enalbled false if not do_irrigation

File:

• branches/ORCHIDEE_2_2/ORCHIDEE/src_parameters/constantes.f90 (modified) (1 diff)

branches/ORCHIDEE_2_2/ORCHIDEE/src_parameters/constantes.f90

@@ -879 +879 @@
 
 
+
+    !Crop irrigation
+    !Config Key   = IRRIG_DOSMAX
+    !Config Desc  = The maximum irrigation water injected per hour (kg.m^{-2}/hour)
+    !Config If    = DO_IRRIGATION AND old_irrig_scheme = FALSE
+    !Config Def   = 1.
+    !Config Help  =
+    !Config Units = [kg.m^{-2}/hour]
+    CALL getin_p('IRRIG_DOSMAX',irrig_dosmax)
+
+    !Config Key   = CUM_NROOT_THR
+    !Config Desc  = Cumulated nroot threshoold to define root zone, and calculate water deficit for irrigation
+    !Config If    = DO_IRRIGATION AND old_irrig_scheme = FALSE
+    !Config Def   = 0.90
+    !Config Help  =
+    !Config Units = []
+    CALL getin_p('CUM_NROOT_THR',cum_nroot_thr)
+
+    !Crop irrigation
+    !Config Key   = IRRIGATED_SOILTILE
+    !Config Desc  = Do we introduce a new soil tile for irrigated croplands?
+    !Config If    = DO_IRRIGATION AND old_irrig_scheme = FALSE
+    !Config Def   = n
+    !Config Help  =
+    !Config Units = [FLAG]
+    CALL getin_p('IRRIGATED_SOILTILE',irrigated_soiltile)
+
+    !
+    !Config Key   = OLD_IRRIG_SCHEME
+    !Config Desc  = Do we run with the old irrigation scheme?
+    !Config If    = DO_IRRIGATION
+    !Config Def   = n
+    !Config Help  =
+    !Config Units = [FLAG]
+    CALL getin_p('OLD_IRRIG_SCHEME',old_irrig_scheme)
+
+    !   To run with new irrigation scheme
+    !Config Key   = IRRIG_ST
+    !Config Desc  = Which is the soil tile with irrigation flux
+    !Config If    = DO_IRRIGATION and old_irrig_scheme = FALSE and IRRIGATED_SOILTILE = FALSE
+    !Config Def   = 3
+    !Config Help  =
+    !Config Units = []
+    CALL getin_p('IRRIG_ST',irrig_st)
+
+    !  To run with NEW irrigation scheme
+    !Config Key   = AVAIL_RESERVE
+    !Config Desc  = Max. fraction of routing reservoir volume that can be used for irrigation
+    !Config If    = DO_IRRIGATION and old_irrig_scheme = FALSE
+    !Config Def   = 0.9,0.0,0.9
+    !Config Help  =
+    !Config Units = []
+    CALL getin_p('AVAIL_RESERVE',avail_reserve)
+
+    !  To run with NEW irrigation scheme
+    !Config Key   = BETA_IRRIG
+    !Config Desc  = Threshold multiplier of Target SM to calculate root deficit
+    !Config If    = DO_IRRIGATION and old_irrig_scheme = FALSE
+    !Config Def   = 1.
+    !Config Help  =
+    !Config Units = []
+    CALL getin_p('BETA_IRRIG',beta_irrig)
+    !
+    !   To run with new irrigation scheme
+    !Config Key   = LAI_IRRIG_MIN
+    !Config Desc  = Min. LAI (Leaf Area Index) to trigger irrigation
+    !Config If    = DO_IRRIGATION and old_irrig_scheme = FALSE and IRRIGATED_SOILTILE = FALSE
+    !Config Def   = 0.1
+    !Config Help  =
+    !Config Units = [m2/m2]
+    CALL getin_p('LAI_IRRIG_MIN',lai_irrig_min)
+
+    !   To run with new irrigation scheme
+    !Config Key   = irrig_map_dynamic
+    !Config Desc  = Do we use a dynamic irrig map?
+    !Config If    = DO_IRRIGATION = T
+    !Config Def   = n
+    !Config Help  =
+    !Config Units = [FLAG]
+    CALL getin_p('IRRIG_MAP_DYNAMIC_FLAG',irrig_map_dynamic_flag)
+
+    !   To run with new irrigation scheme
+    !Config Key   = select_source_irrig
+    !Config Desc  = Do we use the new priorization scheme, based on maps of equipped area with surface water?
+    !Config If    = DO_IRRIGATION = T
+    !Config Def   = n
+    !Config Help  =
+    !Config Units = [FLAG]
+    CALL getin_p('SELECT_SOURCE_IRRIG',select_source_irrig)
+
+    !   Do we reinfiltrate all/part of runoff in crop soil tile?
+    !Config Key   = Reinfiltr_IrrigField
+    !Config Desc  = DO_IRRIGATION
+    !Config If    = DO_IRRIGATION = T
+    !Config Def   = n
+    !Config Help  =
+    !Config Units = [FLAG]
+    CALL getin_p('REINFILTR_IRRIGFIELD',Reinfiltr_IrrigField)
+
+    !   Externalized reinf_slope for reinfiltration in irrigated cropland
+    !Config Key   = reinf_slope_cropParam
+    !Config Desc  = DO_IRRIGATION
+    !Config If    = DO_IRRIGATION = T, REINFILTR_IRRIGFIELD = T
+    !Config Def   = 0.8
+    !Config Help  =
+    !Config Units = [FRACTION 0-1]
+    CALL getin_p('REINF_SLOPE_CROPPARAM',reinf_slope_cropParam)
+
+    !   Externalized for available volume to adduction
+    !Config Key   = a_stream_adduction  
+    !Config Desc  = DO_IRRIGATION
+    !Config If    = DO_IRRIGATION = T
+    !Config Def   = zero
+    !Config Help  =
+    !Config Units = [FRACTION 0-1]
+    CALL getin_p('A_STREAM_ADDUCTION',a_stream_adduction)
+
+
     !Surface resistance
     !