Changeset 343 for branches/ORCHIDEE_EXT/ORCHIDEE/src_parameters/pft_parameters.f90

Timestamp:

2011-07-26T11:04:04+02:00 (13 years ago)

Author:

didier.solyga

Message:

Delete all default initializations of arrays by zero. Externalized llaimax and height_presc because they are read even if NOT(impveg). Activate getin when the corresponding flag is activated

File:

• branches/ORCHIDEE_EXT/ORCHIDEE/src_parameters/pft_parameters.f90 (modified) (2 diffs)

branches/ORCHIDEE_EXT/ORCHIDEE/src_parameters/pft_parameters.f90

@@ -512 +512 @@
 
    !
-   ! 1. Initialisation !! not all the parameters are initialized 
-   !
-
-   !-
-   ! Vegetation structure 
-   !-
-   !
-   ! 1 .Sechiba
-   !
-   veget_ori_fixed_test_1(:) = zero
-   llaimax(:) = zero
-   llaimin(:) = zero
-   height_presc(:) = zero
-   rveg_pft(:) = zero
-   !
-   ! 2 .Stomate
-   !
-   leaf_tab(:) = zero_int
-   sla(:) = zero   
-   !-
-   ! Evapotranspiration -  sechiba
-   !-
-   rstruct_const(:) = zero
-   kzero(:) = zero
-   !-
-   ! Water - sechiba
-   !-
-   wmax_veg(:) = zero
-   humcste(:) = zero
-   throughfall_by_pft(:) = zero
-   !-
-   ! Albedo - sechiba
-   !-
-   snowa_ini(:) = zero
-   snowa_dec(:) = zero
-   alb_leaf_vis(:) = zero
-   alb_leaf_nir(:) = zero
-   !-
-   ! Soil - vegetation
-   !-
-   pref_soil_veg_sand(:) = zero_int
-   pref_soil_veg_loan(:) = zero_int
-   pref_soil_veg_clay(:) = zero_int
-   !-
-   ! Photosynthesis
-   !-
-   !
-   ! 1 .CO2
-   !
-   gsslope(:) = zero
-   gsoffset(:) = zero
-   vcmax_fix(:) = zero
-   vjmax_fix(:) = zero
-   co2_tmin_fix(:) = zero
-   co2_topt_fix(:) = zero
-   co2_tmax_fix(:) = zero
-   ext_coeff(:) = zero
-   !
-   ! 2 .Stomate
-   !
-   vcmax_opt(:) = zero
-   vjmax_opt(:) = zero
-   tphoto_min_a(:) = zero
-   tphoto_min_b(:) = zero
-   tphoto_min_c(:) = zero
-   tphoto_opt_a(:) = zero
-   tphoto_opt_b(:) = zero
-   tphoto_opt_c(:) = zero
-   tphoto_max_a(:) = zero
-   tphoto_max_b(:) = zero
-   tphoto_max_c(:) = zero
-   !-
-   ! Respiration - stomate
-   !-
-   maint_resp_slope_c(:) = zero
-   maint_resp_slope_b(:) = zero
-   maint_resp_slope_a(:) = zero
-   cm_zero_leaf(:) = zero
-   cm_zero_sapabove(:) = zero
-   cm_zero_sapbelow(:) = zero
-   cm_zero_heartabove(:) = zero
-   cm_zero_heartbelow(:) = zero
-   cm_zero_root(:) = zero
-   cm_zero_fruit(:) = zero
-   cm_zero_carbres(:) = zero
-   !-
-   ! Fire - stomate
-   !-
-   ! 
-   flam(:) = zero
-   resist(:) = zero
-   !-
-   ! Flux - LUC
-   !-
-   coeff_lcchange_1(:) = zero
-   coeff_lcchange_10(:) = zero
-   coeff_lcchange_100(:) = zero
-   !-
-   ! Phenology
-   !-
-   !
-   ! 1 .Stomate
-   !
-   lai_max(:) = zero
-   pheno_type(:) = zero_int
-   !
-   ! 2. Leaf Onset
-   !
-   pheno_gdd_crit_c(:) = zero
-   pheno_gdd_crit_b(:) = zero
-   pheno_gdd_crit_a(:) = zero
-   ngd_crit(:) = zero
-   ncdgdd_temp(:) = zero
-   hum_frac(:) = zero
-   lowgpp_time(:) = zero
-   hum_min_time(:) = zero
-   tau_sap(:) = zero
-   tau_fruit(:) = zero
-   ecureuil(:) = zero
-   alloc_min(:) = zero
-   alloc_max(:) = zero
-   demi_alloc(:) = zero  
-  !>> DS new for merge in the trunk 15/06/2011
-   leaflife_tab(:) = zero
-   !
-   ! 3. Senescence
-   !
-   leaffall(:) = zero
-   leafagecrit(:) = zero
-   senescence_hum(:) = zero
-   nosenescence_hum(:) = zero
-   max_turnover_time(:) = zero   
-   min_turnover_time(:) = zero 
-   min_leaf_age_for_senescence(:) = zero 
-   senescence_temp_c(:) = zero 
-   senescence_temp_b(:) = zero 
-   senescence_temp_a(:) = zero 
-   !-
-   ! DGVM
-   !-
-   residence_time(:) = zero
-   tmin_crit(:) = zero
-   tcm_crit(:) = zero
-   !-
-   !  Internal_parameters
-   !-
-   lai_initmin(:) = zero
-   bm_sapl(:,:) = zero
-   migrate(:) = zero
-   maxdia(:) = zero
-   cn_sapl(:) = zero
-   leaf_timecst(:) = zero  
-
-
-   !
-   ! 2. Correspondance between the PFTs values and thes MTCs values 
+   ! 1. Correspondance between the PFTs values and thes MTCs values 
    !
  
    DO j= 1, nvm
-
 
       PFT_name(j) = MTC_name(pft_to_mtc(j))
@@ -1099 +943 @@
   IF(first_call) THEN
 
-     ! No calling to getin for veget_ori_fixed_test_1, llaimax and height_presc
-     ! use of setvar in slowproc.f90
-
      !-
      ! Vegetation structure
      !-
+     CALL getin_p('SECHIBA_LAI',llaimax)
      CALL getin_p('LLAIMIN',llaimin)
+     CALL getin_p('SLOWPROC_HEIGHT', height_presc)
      CALL getin('TYPE_OF_LAI',type_of_lai)
      CALL getin_p('IS_TREE',is_tree)