Changeset 350


Ignore:
Timestamp:
2011-07-29T14:04:49+02:00 (13 years ago)
Author:
didier.solyga
Message:

Add labels for the externalized pft parameters. Update obsolete aspects of F77. Replace the last 0.0 by zero in stomate

Location:
branches/ORCHIDEE_EXT/ORCHIDEE
Files:
8 edited

Legend:

Unmodified
Added
Removed

  • branches/ORCHIDEE_EXT/ORCHIDEE/src_parameters/pft_parameters.f90

    r343 r350  
    411411       
    412412      ! 3.2 Reading of the conrrespondance table in the .def file 
     413      ! 
     414      !Config  Key  = PFT_TO_MTC 
     415      !Config  Desc = correspondance array linking a PFT to MTC 
     416      !Config  if  = ANYTIME 
     417      !Config  Def  = 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13 
     418      !Config  Help = 
     419      !Config  Units = NONE 
    413420      CALL getin_p('PFT_TO_MTC',pft_to_mtc) 
    414421       
    … …  
    457464 
    458465      ! 5.1 Read the name of the PFTs given by the user 
     466      ! 
     467      !Config  Key = PFT_NAME 
     468      !Config  Desc = Name of a PFT 
     469      !Config  if = ANYTIME 
     470      !Config  Def = bared ground, tropical broad-leaved evergreen, tropical broad-leaved raingreen, temperate needleleaf evergreen,temperate broad-leaved evergreen 
     471      !              temperate broad-leaved summergreen, boreal needleleaf evergreen, boreal broad-leaved summergreen, boreal needleleaf summergreen, 
     472      !              C3 grass, C4 grass, C3 agriculture, C4 agriculture     
     473      !Config  Help = the user can name the new PFTs he/she introducing for new species 
     474      !Config  Units = NONE 
    459475      CALL getin('PFT_NAME',pft_name) 
    460476 
    … …  
    946962     ! Vegetation structure 
    947963     !- 
     964     ! 
     965     !Config  Key  = SECHIBA_LAI 
     966     !Config  Desc = laimax for maximum lai see also type of lai interpolation 
     967     !Config  if  =  IMPOSE_VEG 
     968     !Config  Def  = 0., 8., 8., 4., 4.5, 4.5, 4., 4.5, 4., 2., 2., 2., 2. 
     969     !Config  Help = Values of lai used for interpolation of the lai map 
     970     !Config  Units =  
    948971     CALL getin_p('SECHIBA_LAI',llaimax) 
     972     ! 
     973     !Config  Key  = LLAIMIN 
     974     !Config  Desc = laimin for minimum lai see also type of lai interpolation 
     975     !Config  if  = OK_SECHIBA 
     976     !Config  Def  = 0., 8., 0., 4., 4.5, 0., 4., 0., 0., 0., 0., 0., 0. 
     977     !Config  Help = 
     978     !Config  Units =  
    949979     CALL getin_p('LLAIMIN',llaimin) 
     980     ! 
     981     !Config  Key  = SLOWPROC_HEIGHT 
     982     !Config  Desc = prescribed height of vegetation : one for each vegetation type 
     983     !Config  if  = OK_SECHIBA 
     984     !Config  Def  = 0.,30.,30.,20.,20.,20.,15.,15.,15.,.5,.6,1.,1. 
     985     !Config  Help = 
     986     !Config  Units = Meters (m) ? 
    950987     CALL getin_p('SLOWPROC_HEIGHT', height_presc) 
     988     ! 
     989     !Config  Key  = TYPE_OF_LAI 
     990     !Config  Desc = Type of behaviour of the LAI evolution algorithm for each vegetation type : mean or interp 
     991     !Config  if  = OK_SECHIBA 
     992     !Config  Def  = inter','inter','inter','inter','inter','inter','inter','inter','inter','inter','inter','inter','inter' 
     993     !Config  Help = 
     994     !Config  Units = NONE 
    951995     CALL getin('TYPE_OF_LAI',type_of_lai) 
     996     ! 
     997     !Config  Key  = IS_TREE 
     998     !Config  Desc = Is the vegetation type a tree ? 
     999     !Config  if  = OK_SECHIBA 
     1000     !Config  Def  = n, y, y, y, y, y, y, y, y, n, n, n, n 
     1001     !Config  Help = 
     1002     !Config  Units = NONE 
    9521003     CALL getin_p('IS_TREE',is_tree) 
     1004     ! 
     1005     !Config  Key  = NATURAL 
     1006     !Config  Desc = natural?  
     1007     !Config  if  = OK_SECHIBA 
     1008     !Config  Def  = y, y, y, y, y, y, y, y, y, y, y, n, n  
     1009     !Config  Help = 
     1010     !Config  Units = NONE 
    9531011     CALL getin_p('NATURAL',natural) 
    9541012 
    9551013     !>> DS new for merge in the trunk   ! 15/06/2011  
    9561014     ! Add for writing history files in stomate_lpj.f90 'treeFracPrimDec' and 'treeFracPrimEver' 
     1015     ! 
     1016     !Config  Key  = IS_DECIDUOUS 
     1017     !Config  Desc = is PFT deciduous ? 
     1018     !Config  if  = OK_SECHIBA 
     1019     !Config  Def  = n, n, y, n, n, y, n, y, y, n, n, n, n 
     1020     !Config  Help = 
     1021     !Config  Units = NONE 
    9571022     CALL getin('IS_DECIDUOUS',is_deciduous) 
    958      CALL getin('IS_EVERGREEN',is_evergreen)   
     1023     ! 
     1024     !Config  Key  = IS_EVERGREEN 
     1025     !Config  Desc = is PFT evergreen ? 
     1026     !Config  if  = OK_SECHIBA 
     1027     !Config  Def  = n, y, n, y, y, n, y, n, n, n, n, n, n 
     1028     !Config  Help = 
     1029     !Config  Units = NONE  
     1030     CALL getin('IS_EVERGREEN',is_evergreen) 
     1031     ! 
     1032     !Config  Key  = IS_C3 
     1033     !Config  Desc = is PFT C3 ? 
     1034     !Config  if  = OK_SECHIBA 
     1035     !Config  Def  = n, n, n, n, n, n, n, n, n, n, y, n, y, n 
     1036     !Config  Help = 
     1037     !Config  Units = NONE    
    9591038     CALL getin_p('IS_C3',is_c3)    
    9601039 
    … …  
    9621041     ! Photosynthesis 
    9631042     !- 
     1043     ! 
     1044     !Config  Key  = IS_C4 
     1045     !Config  Desc = flag for C4 vegetation types 
     1046     !Config  if  = OK_SECHIBA 
     1047     !Config  Def  = n, n, n, n, n, n, n, n, n, n, n, y, n, y 
     1048     !Config  Help = 
     1049     !Config  Units = NONE 
    9641050     CALL getin_p('IS_C4',is_c4) 
     1051     ! 
     1052     !Config  Key  = GSSLOPE 
     1053     !Config  Desc = Slope of the gs/A relation (Ball & al. 
     1054     !Config  if  = OK_SECHIBA AND OK_CO2 
     1055     !Config  Def  = 0., 9., 9., 9., 9., 9., 9., 9., 9., 9., 3., 9., 3. 
     1056     !Config  Help = 
     1057     !Config  Units =  
    9651058     CALL getin_p('GSSLOPE',gsslope) 
     1059     ! 
     1060     !Config  Key  = GSOFFSET 
     1061     !Config  Desc = intercept of the gs/A relation (Ball & al.) 
     1062     !Config  if  = OK_SECHIBA 
     1063     !Config  Def  = 0.0,  0.01, 0.01, 0.01, 0.01, 0.01, 0.01, 0.01, 0.01, 0.01, 0.03, 0.01, 0.03 
     1064     !Config  Help = 
     1065     !Config  Units =  
    9661066     CALL getin_p('GSOFFSET',gsoffset) 
     1067     ! 
     1068     !Config  Key  = VCMAX_FIX 
     1069     !Config  Desc = values used for vcmax when STOMATE is not activated 
     1070     !Config  if  =  OK_SECHIBA 
     1071     !Config  Def  =  0., 40., 50., 30., 35., 40.,30., 40., 35., 60., 60., 70., 70. 
     1072     !Config  Help = 
     1073     !Config  Units =  
    9671074     CALL getin_p('VCMAX_FIX',vcmax_fix) 
     1075     ! 
     1076     !Config  Key  = VJMAX_FIX 
     1077     !Config  Desc = values used for vjmax when STOMATE is not activated 
     1078     !Config  if  = OK_SECHIBA 
     1079     !Config  Def  =  0., 80., 100., 60., 70., 80.,  60., 80., 70., 120., 120., 140., 140. 
     1080     !Config  Help = 
     1081     !Config  Units =  
    9681082     CALL getin_p('VJMAX_FIX',vjmax_fix) 
     1083     ! 
     1084     !Config  Key  = CO2_TMIN_FIX 
     1085     !Config  Desc = values used for photosynthesis tmin when STOMATE is not activated 
     1086     !Config  if  = OK_SECHIBA  
     1087     !Config  Def  = 0.,  2.,  2., -4., -3., -2., -4., -4., -4., -5.,  6., -5.,  6. 
     1088     !Config  Help = 
     1089     !Config  Units =  
    9691090     CALL getin_p('CO2_TMIN_FIX',co2_tmin_fix) 
     1091     ! 
     1092     !Config  Key  = CO2_TOPT_FIX  
     1093     !Config  Desc = values used for photosynthesis topt when STOMATE is not activated  
     1094     !Config  if  = OK_SECHIBA 
     1095     !Config  Def  =  0., 27.5, 27.5, 17.5, 25., 20.,17.5, 17.5, 17.5, 20.,  32.5, 20.,  32.5 
     1096     !Config  Help = 
     1097     !Config  Units =  
    9701098     CALL getin_p('CO2_TOPT_FIX',co2_topt_fix) 
     1099     ! 
     1100     !Config  Key  = CO2_TMAX_FIX 
     1101     !Config  Desc = values used for photosynthesis tmax when STOMATE is not activated  
     1102     !Config  if  = OK_SECHIBA 
     1103     !Config  Def  = 0., 55., 55., 38., 48., 38.,38., 38., 38., 45., 55., 45., 55. 
     1104     !Config  Help = 
     1105     !Config  Units =  
    9711106     CALL getin_p('CO2_TMAX_FIX',co2_tmax_fix) 
     1107     ! 
     1108     !Config  Key  = EXT_COEFF 
     1109     !Config  Desc = extinction coefficient of the Monsi&Seaki relationship (1953) 
     1110     !Config  if  = OK_SECHIBA OR OK_STOMATE 
     1111     !Config  Def  = .5, .5, .5, .5, .5, .5, .5, .5, .5, .5, .5, .5, .5 
     1112     !Config  Help = 
     1113     !Config  Units =  
    9721114     CALL getin_p('EXT_COEFF',ext_coeff) 
    9731115     !- 
    9741116     ! Evapotranspiration -  sechiba 
    9751117     !- 
     1118     ! 
     1119     !Config  Key  = RSTRUCT_CONST 
     1120     !Config  Desc = Structural resistance : one for each vegetation type 
     1121     !Config  if  = OK_SECHIBA 
     1122     !Config  Def  = 0.0, 25.0, 25.0, 25.0, 25.0, 25.0, 25.0, 25.0, 25.0,  2.5,  2.0,  2.0,  2.0 
     1123     !Config  Help = 
     1124     !Config  Units =  
    9761125     CALL getin_p('RSTRUCT_CONST',rstruct_const) 
     1126     ! 
     1127     !Config  Key  = KZERO 
     1128     !Config  Desc =  A vegetation dependent constant used in the calculation of the surface resistance. 
     1129     !Config  if  = OK_SECHIBA 
     1130     !Config  Def  = 0.0, 12.E-5, 12.E-5, 12.e-5, 12.e-5, 25.e-5, 12.e-5,25.e-5, 25.e-5, 30.e-5, 30.e-5, 30.e-5, 30.e-5  
     1131     !Config  Help = 
     1132     !Config  Units =  
    9771133     CALL getin_p('KZERO',kzero) 
    9781134     ! 
    … …  
    9811137     !Config Key  = RVEG_PFT 
    9821138     !Config Desc = Artificial parameter to increase or decrease canopy resistance. 
     1139     !Config  if  = OK_SECHIBA 
    9831140     !Config Def  = 1. 
    9841141     !Config Help = This parameter is set by PFT. 
     1142     !Config  Units =  
    9851143     CALL getin_p('RVEG_PFT', rveg_pft)     
    9861144     !- 
    9871145     ! Water-hydrology - sechiba 
    9881146     !- 
     1147     ! 
     1148     !Config  Key  = WMAX_VEG 
     1149     !Config  Desc = Maximum field capacity for each of the vegetations (Temporary): max quantity of water 
     1150     !Config  if  = OK_SECHIBA 
     1151     !Config  Def  =  150., 150., 150., 150., 150., 150., 150.,150., 150., 150., 150., 150., 150. 
     1152     !Config  Help = 
     1153     !Config  Units = Kg/M3 
    9891154     CALL getin_p('WMAX_VEG',wmax_veg) 
     1155     ! 
    9901156     !Config Key  = HYDROL_HUMCSTE 
    9911157     !Config Desc = Root profile 
    9921158     !Config Def  = 5., .8, .8, 1., .8, .8, 1., 1., .8, 4., 4., 4., 4. 
     1159     !Config  if  = OK_SECHIBA 
    9931160     !Config Help = Default values were defined for 2 meters soil depth. 
    9941161     !Config        For 4 meters soil depth, you may use those ones : 
    9951162     !Config        5., .4, .4, 1., .8, .8, 1., 1., .8, 4., 1., 4., 1. 
     1163     !Config  Units =  
    9961164     CALL getin_p('HYDROL_HUMCSTE', humcste) 
     1165     ! 
    9971166     !Config  Key  = PERCENT_THROUGHFALL_PFT 
    9981167     !Config  Desc = Percent by PFT of precip that is not intercepted by the canopy 
     1168     !Config  if  = OK_SECHIBA OR OK_CWRR 
    9991169     !Config  Def  = 30. 30. 30. 30. 30. 30. 30. 30. 30. 30. 30. 30. 30. 
    10001170     !Config  Help = During one rainfall event, PERCENT_THROUGHFALL_PFT% of the incident rainfall 
    10011171     !Config         will get directly to the ground without being intercepted, for each PFT. 
     1172     !Config  Units = Kg/M3 
    10021173     CALL getin_p('PERCENT_TROUGHFALL_PFT',throughfall_by_pft) 
    10031174     !- 
    10041175     ! Albedo - sechiba 
    10051176     !- 
     1177     ! 
     1178     !Config  Key  = SNOWA_INI 
     1179     !Config  Desc = Initial snow albedo value for each vegetation type as it will be used in condveg_snow 
     1180     !Config  if  = OK_SECHIBA 
     1181     !Config  Def  = 0.35, 0.,   0.,   0.14, 0.14,0.14, 0.14, 0.14, 0.14, 0.18,0.18, 0.18, 0.18 
     1182     !Config  Help = Values are from the Thesis of S. Chalita (1992) 
     1183     !Config  Units =  
    10061184     CALL getin_p('SNOWA_INI',snowa_ini) 
     1185     ! 
     1186     !Config  Key  = SNOWA_DEC 
     1187     !Config  Desc = Decay rate of snow albedo value for each vegetation type as it will be used in condveg_snow 
     1188     !Config  if  =  OK_SECHIBA 
     1189     !Config  Def  = 0.45, 0.,  0., 0.06, 0.06, 0.11, 0.06, 0.11, 0.11, 0.52,0.52, 0.52, 0.52 
     1190     !Config  Help = Values are from the Thesis of S. Chalita (1992) 
     1191     !Config  Units =  
    10071192     CALL getin_p('SNOWA_DEC',snowa_dec) 
     1193     ! 
     1194     !Config  Key  = ALB_LEAF_VIS 
     1195     !Config  Desc = leaf albedo of vegetation type, visible albedo 
     1196     !Config  if  = OK_SECHIBA 
     1197     !Config  Def  = .00, .04, .06, .06, .06,.06, .06, .06, .06, .10, .10, .10, .10 
     1198     !Config  Help = 
     1199     !Config  Units =  
    10081200     CALL getin_p('ALB_LEAF_VIS',alb_leaf_vis) 
     1201     ! 
     1202     !Config  Key  = ALB_LEAF_NIR 
     1203     !Config  Desc = leaf albedo of vegetation type, near infrared albedo 
     1204     !Config  if  = OK_SECHIBA 
     1205     !Config  Def  =.00, .20, .22, .22, .22,.22, .22, .22, .22, .30,.30, .30, .30  
     1206     !Config  Help = 
     1207     !Config  Units =  
    10091208     CALL getin_p('ALB_LEAF_NIR',alb_leaf_nir) 
    10101209     !- 
    10111210     ! Soil - vegetation 
    10121211     !- 
     1212     ! 
     1213     !Config  Key  = PREF_SOIL_VEG_SAND 
     1214     !Config  Desc = Table which contains the correlation between the soil types and vegetation type 
     1215     !Config  if  = OK_SECHIBA 
     1216     !Config  Def  = 1, 3, 3, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2 
     1217     !Config  Help = first layer of the soil 
     1218     !Config  Units =  
    10131219     CALL getin_p('PREF_SOIL_VEG_SAND',pref_soil_veg_sand) 
    1014      CALL getin_p('PREF_SOIL_VEG_LOAN',pref_soil_veg_loan)          
     1220     ! 
     1221     !Config  Key  = PREF_SOIL_VEG_LOAN 
     1222     !Config  Desc = Table which contains the correlation between the soil types and vegetation type 
     1223     !Config  if  = OK_SECHIBA 
     1224     !Config  Def  = 2, 2, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3 
     1225     !Config  Help = second layer of the soil 
     1226     !Config  Units =  
     1227     CALL getin_p('PREF_SOIL_VEG_LOAN',pref_soil_veg_loan)  
     1228      ! 
     1229     !Config  Key  = PREF_SOIL_VEG_CLAY 
     1230     !Config  Desc = Table which contains the correlation between the soil types and vegetation type 
     1231     !Config  if  = OK_SECHIBA 
     1232     !Config  Def  = 3, 1, 1, 1, 1, 1 ,1 ,1 ,1 ,1 ,1 ,1, 1  
     1233     !Config  Help = third layer of the soil  
     1234     !Config  Units =         
    10151235     CALL getin_p('PREF_SOIL_VEG_CLAY',pref_soil_veg_clay) 
    10161236 
    … …  
    10311251  IF(first_call) THEN 
    10321252 
    1033       !- 
    1034       ! Vegetation structure 
    1035       !- 
    1036       CALL getin_p('LEAF_TAB',leaf_tab) 
    1037       CALL getin_p('SLA',sla) 
    1038       !- 
    1039       ! Photosynthesis 
    1040       !- 
    1041       CALL getin_p('VCMAX_OPT',vcmax_opt) 
    1042       CALL getin_p('VJMAX_OPT',vjmax_opt) 
    1043       CALL getin_p('TPHOTO_MIN_A',tphoto_min_a) 
    1044       CALL getin_p('TPHOTO_MIN_B',tphoto_min_b) 
    1045       CALL getin_p('TPHOTO_MIN_C',tphoto_min_c) 
    1046       CALL getin_p('TPHOTO_OPT_A',tphoto_opt_a) 
    1047       CALL getin_p('TPHOTO_OPT_B',tphoto_opt_b) 
    1048       CALL getin_p('TPHOTO_OPT_C',tphoto_opt_c) 
    1049       CALL getin_p('TPHOTO_MAX_A',tphoto_max_a) 
    1050       CALL getin_p('TPHOTO_MAX_B',tphoto_max_b) 
    1051       CALL getin_p('TPHOTO_MAX_C',tphoto_max_c) 
    1052       !- 
    1053       ! Respiration - stomate 
    1054       !- 
    1055       CALL getin_p('MAINT_RESP_SLOPE_C',maint_resp_slope_c)  
    1056       CALL getin_p('MAINT_RESP_SLOPE_B',maint_resp_slope_b) 
    1057       CALL getin_p('MAINT_RESP_SLOPE_A',maint_resp_slope_a) 
    1058       CALL getin_p('CM_ZERO_LEAF',cm_zero_leaf) 
    1059       CALL getin_p('CM_ZERO_SAPABOVE',cm_zero_sapabove) 
    1060       CALL getin_p('CM_ZERO_SAPBELOW',cm_zero_sapbelow) 
    1061       CALL getin_p('CM_ZERO_HEARTABOVE',cm_zero_heartabove) 
    1062       CALL getin_p('CM_ZERO_HEARTBELOW',cm_zero_heartbelow) 
    1063       CALL getin_p('CM_ZERO_ROOT',cm_zero_root) 
    1064       CALL getin_p('CM_ZERO_FRUIT',cm_zero_fruit) 
    1065       CALL getin_p('CM_ZERO_CARBRES',cm_zero_carbres) 
    1066       !- 
    1067       ! Fire - stomate 
    1068       !- 
    1069       CALL getin_p('FLAM',flam) 
    1070       CALL getin_p('RESIST',resist) 
    1071       !- 
    1072       ! Flux - LUC 
    1073       !- 
    1074       CALL getin_p('COEFF_LCCHANGE_1',coeff_lcchange_1) 
    1075       CALL getin_p('COEFF_LCCHANGE_10',coeff_lcchange_10) 
    1076       CALL getin_p('COEFF_LCCHANGE_100',coeff_lcchange_100) 
    1077       !- 
    1078       ! Phenology 
    1079       !- 
    1080       CALL getin_p('LAI_MAX',lai_max) 
    1081       CALL getin('PHENO_MODEL',pheno_model) 
    1082       CALL getin_p('PHENO_TYPE',pheno_type) 
    1083       !- 
    1084       ! Phenology : Leaf Onset 
    1085       !- 
    1086       CALL getin_p('PHENO_GDD_CRIT_C',pheno_gdd_crit_c) 
    1087       CALL getin_p('PHENO_GDD_CRIT_B',pheno_gdd_crit_b) 
    1088       CALL getin_p('PHENO_GDD_CRIT_A',pheno_gdd_crit_a) 
    1089       CALL getin_p('NGD_CRIT',ngd_crit) 
    1090       CALL getin_p('NCDGDD_TEMP', ncdgdd_temp) 
    1091       CALL getin_p('HUM_FRAC', hum_frac) 
    1092       CALL getin_p('LOWGPP_TIME', lowgpp_time) 
    1093       CALL getin_p('HUM_MIN_TIME', hum_min_time) 
    1094       CALL getin_p('TAU_SAP',tau_sap) 
    1095       CALL getin_p('TAU_FRUIT',tau_fruit) 
    1096       CALL getin_p('ECUREUIL',ecureuil) 
    1097       CALL getin_p('ALLOC_MIN',alloc_min) 
    1098       CALL getin_p('ALLOC_MAX',alloc_max) 
    1099       CALL getin_p('DEMI_ALLOC',demi_alloc) 
    1100  
    1101       !>> DS new for merge in the trunk 
    1102       ! 15/06/2011 : add leaflife_mtc for the new formalism used for calculate sla 
    1103       CALL getin_p('LEAFLIFE_TAB',leaflife_tab) 
    1104  
    1105       !- 
    1106       ! Phenology : Senescence 
    1107       !- 
    1108       CALL getin_p('LEAFFALL',leaffall) 
    1109       CALL getin_p('LEAFAGECRIT',leafagecrit)   
    1110       CALL getin('SENESCENCE_TYPE', senescence_type)  
    1111       CALL getin_p('SENESCENCE_HUM', senescence_hum) 
    1112       CALL getin_p('NOSENESCENCE_HUM', nosenescence_hum)  
    1113       CALL getin_p('MAX_TURNOVER_TIME',max_turnover_time) 
    1114       CALL getin_p('MIN_TURNOVER_TIME',min_turnover_time) 
    1115       CALL getin_p('MIN_LEAF_AGE_FOR_SENESCENCE', min_leaf_age_for_senescence) 
    1116       CALL getin_p('SENESCENCE_TEMP_C',senescence_temp_c) 
    1117       CALL getin_p('SENESCENCE_TEMP_B',senescence_temp_b) 
    1118       CALL getin_p('SENESCENCE_TEMP_A',senescence_temp_a) 
    1119       !- 
    1120       ! DGVM 
    1121       !- 
    1122       CALL getin_p('RESIDENCE_TIME',residence_time) 
    1123       CALL getin_p('TMIN_CRIT',tmin_crit) 
    1124       CALL getin_p('TCM_CRIT',tcm_crit) 
     1253     !- 
     1254     ! Vegetation structure 
     1255     !- 
     1256     ! 
     1257     !Config  Key  = LEAF_TAB 
     1258     !Config  Desc =  leaf type : 1=broad leaved tree, 2=needle leaved tree, 3=grass 4=bare ground 
     1259     !Config  if  = OK_STOMATE 
     1260     !Config  Def  =  4, 1, 1, 2, 1, 1, 2, 1, 2, 3, 3, 3, 3  
     1261     !Config  Help =  
     1262     !Config  Units = NONE 
     1263     CALL getin_p('LEAF_TAB',leaf_tab) 
     1264     ! 
     1265     !Config  Key  = SLA 
     1266     !Config  Desc =  specif leaf area  
     1267     !Config  if  = OK_STOMATE 
     1268     !Config  Def  = 1.5E-2, 1.53E-2,  2.6E-2, 9.26E-3,  2E-2,  2.6E-2,  9.26E-3, 2.6E-2,  1.9E-2,  2.6E-2, 2.6E-2, 2.6E-2, 2.6E-2 
     1269     !Config  Help = 
     1270     !Config  Units = m**2/gC 
     1271     CALL getin_p('SLA',sla) 
     1272     !- 
     1273     ! Photosynthesis 
     1274     !- 
     1275     ! 
     1276     !Config  Key  = VCMAX_OPT 
     1277     !Config  Desc =  Maximum rate of carboxylation 
     1278     !Config  if  = OK_STOMATE 
     1279     !Config  Def  =  undef, 65., 65., 35., 45., 55., 35., 45., 35., 70., 70.,70., 70. 
     1280     !Config  Help = 
     1281     !Config  Units = % 
     1282     CALL getin_p('VCMAX_OPT',vcmax_opt) 
     1283     ! 
     1284     !Config  Key  = VJMAX_OPT 
     1285     !Config  Desc = Maximum rate of RUbp regeneration 
     1286     !Config  if  =  OK_STOMATE 
     1287     !Config  Def  =  undef, 130., 130., 70., 80., 110., 70., 90., 70., 160.,160.,200., 200. 
     1288     !Config  Help = 
     1289     !Config  Units = % 
     1290     CALL getin_p('VJMAX_OPT',vjmax_opt) 
     1291     ! 
     1292     !Config  Key  = TPHOTO_MIN_A 
     1293     !Config  Desc =  minimum photosynthesis temperature, constant a of ax^2+bx+c (deg C), tabulated 
     1294     !Config  if  = OK_STOMATE 
     1295     !Config  Def  =  undef,  0., 0., 0., 0., 0., 0.,  0., 0.,  0.0025, 0., 0., 0.  
     1296     !Config  Help = 
     1297     !Config  Units = NONE 
     1298     CALL getin_p('TPHOTO_MIN_A',tphoto_min_a) 
     1299     ! 
     1300     !Config  Key  = TPHOTO_MIN_B 
     1301     !Config  Desc = minimum photosynthesis temperature, constant b of ax^2+bx+c (deg C), tabulated 
     1302     !Config  if  = OK_STOMATE 
     1303     !Config  Def  =  undef,  0.,  0., 0., 0., 0., 0., 0., 0., 0.1, 0.,0.,0. 
     1304     !Config  Help = 
     1305     !Config  Units = NONE 
     1306     CALL getin_p('TPHOTO_MIN_B',tphoto_min_b) 
     1307     ! 
     1308     !Config  Key  = TPHOTO_MIN_C 
     1309     !Config  Desc = minimum photosynthesis temperature, constant c of ax^2+bx+c (deg C), tabulated 
     1310     !Config  if  = OK_STOMATE 
     1311     !Config  Def  =  undef,  2., 2., -4., -3.,-2.,-4., -4.,-4.,-3.25, 13.,-5.,13. 
     1312     !Config  Help = 
     1313     !Config  Units = NONE 
     1314     CALL getin_p('TPHOTO_MIN_C',tphoto_min_c) 
     1315     ! 
     1316     !Config  Key  = TPHOTO_OPT_A 
     1317     !Config  Desc = optimum photosynthesis temperature, constant a of ax^2+bx+c (deg C), tabulated 
     1318     !Config  if  = OK_STOMATE 
     1319     !Config  Def  =  undef, 0., 0., 0., 0.,0.,0.,0.,0.,0.0025,0.,0.,0.  
     1320     !Config  Help =  
     1321     !Config  Units = NONE 
     1322     CALL getin_p('TPHOTO_OPT_A',tphoto_opt_a) 
     1323     ! 
     1324     !Config  Key  = TPHOTO_OPT_B 
     1325     !Config  Desc = optimum photosynthesis temperature, constant b of ax^2+bx+c (deg C), tabulated 
     1326     !Config  if  = OK_STOMATE 
     1327     !Config  Def  =  undef, 0.,0.,0.,0.,0.,0., 0.,0.,0.25,0.,0.,0.   
     1328     !Config  Help = 
     1329     !Config  Units = NONE 
     1330     CALL getin_p('TPHOTO_OPT_B',tphoto_opt_b) 
     1331     ! 
     1332     !Config  Key  = TPHOTO_OPT_C 
     1333     !Config  Desc = optimum photosynthesis temperature, constant c of ax^2+bx+c (deg C), tabulated 
     1334     !Config  if  = OK_STOMATE 
     1335     !Config  Def  = undef, 37., 37., 25., 32., 26., 25., 25., 25., 27.25, 36., 30., 36. 
     1336     !Config  Help = 
     1337     !Config  Units = NONE 
     1338     CALL getin_p('TPHOTO_OPT_C',tphoto_opt_c) 
     1339     ! 
     1340     !Config  Key  = TPHOTO_MAX_A 
     1341     !Config  Desc = maximum photosynthesis temperature, constant a of ax^2+bx+c (deg C), tabulated 
     1342     !Config  if  = OK_STOMATE 
     1343     !Config  Def  = undef,  0., 0., 0., 0., 0., 0., 0., 0., 0.00375, 0., 0., 0. 
     1344     !Config  Help = 
     1345     !Config  Units = NONE 
     1346     CALL getin_p('TPHOTO_MAX_A',tphoto_max_a) 
     1347     ! 
     1348     !Config  Key  = TPHOTO_MAX_B 
     1349     !Config  Desc = maximum photosynthesis temperature, constant b of ax^2+bx+c (deg C), tabulated 
     1350     !Config  if  = OK_STOMATE 
     1351     !Config  Def  = undef, 0., 0., 0., 0., 0., 0., 0., 0.,0.35, 0., 0., 0.    
     1352     !Config  Help = 
     1353     !Config  Units = NONE 
     1354     CALL getin_p('TPHOTO_MAX_B',tphoto_max_b) 
     1355     ! 
     1356     !Config  Key  = TPHOTO_MAX_C 
     1357     !Config  Desc = maximum photosynthesis temperature, constant c of ax^2+bx+c (deg C), tabulated  
     1358     !Config  if  = OK_STOMATE 
     1359     !Config  Def  = undef, 55., 55.,38., 48.,38.,38., 38., 38., 41.125, 55., 45., 55.   
     1360     !Config  Help = 
     1361     !Config  Units = NONE 
     1362     CALL getin_p('TPHOTO_MAX_C',tphoto_max_c) 
     1363     !- 
     1364     ! Respiration - stomate 
     1365     !- 
     1366     ! 
     1367     !Config  Key  = MAINT_RESP_SLOPE_C 
     1368     !Config  Desc = slope of maintenance respiration coefficient (1/K), constant c of aT^2+bT+c , tabulated 
     1369     !Config  if  = OK_STOMATE 
     1370     !Config  Def  = undef, .12, .12,.16,.16,.16,.16,.16,.16,.16,.12,.16,.12  
     1371     !Config  Help = 
     1372     !Config  Units = NONE 
     1373     CALL getin_p('MAINT_RESP_SLOPE_C',maint_resp_slope_c)  
     1374     ! 
     1375     !Config  Key  = MAINT_RESP_SLOPE_B 
     1376     !Config  Desc = slope of maintenance respiration coefficient (1/K), constant b of aT^2+bT+c , tabulated 
     1377     !Config  if  = OK_STOMATE 
     1378     !Config  Def  = undef,.0,.0,.0,.0,.0,.0,.0,.0, -.00133,.0, -.00133,.0  
     1379     !Config  Help = 
     1380     !Config  Units = NONE 
     1381     CALL getin_p('MAINT_RESP_SLOPE_B',maint_resp_slope_b) 
     1382     ! 
     1383     !Config  Key  = MAINT_RESP_SLOPE_A 
     1384     !Config  Desc = slope of maintenance respiration coefficient (1/K), constant a of aT^2+bT+c , tabulated 
     1385     !Config  if  = OK_STOMATE 
     1386     !Config  Def  = undef,.0,.0, .0,.0,.0,.0,.0,.0,.0,.0,.0,.0     
     1387     !Config  Help = 
     1388     !Config  Units = NONE 
     1389     CALL getin_p('MAINT_RESP_SLOPE_A',maint_resp_slope_a) 
     1390     ! 
     1391     !Config  Key  = CM_ZERO_LEAF 
     1392     !Config  Desc =  maintenance respiration coefficient at 0 deg C, for leaves, tabulated 
     1393     !Config  if  = OK_STOMATE 
     1394     !Config  Def  = undef,  2.35E-3, 2.62E-3, 1.01E-3, 2.35E-3, 2.62E-3, 1.01E-3,2.62E-3, 2.05E-3, 2.62E-3, 2.62E-3, 2.62E-3, 2.62E-3 
     1395     !Config  Help = 
     1396     !Config  Units = g/g/day 
     1397     CALL getin_p('CM_ZERO_LEAF',cm_zero_leaf) 
     1398     ! 
     1399     !Config  Key  = CM_ZERO_SAPABOVE 
     1400     !Config  Desc = maintenance respiration coefficient at 0 deg C,for sapwood above, tabulated 
     1401     !Config  if  = OK_STOMATE 
     1402     !Config  Def  = undef,  1.19E-4, 1.19E-4, 1.19E-4, 1.19E-4, 1.19E-4, 1.19E-4,1.19E-4, 1.19E-4, 1.19E-4, 1.19E-4, 1.19E-4, 1.19E-4 
     1403     !Config  Help = 
     1404     !Config  Units = g/g/day 
     1405     CALL getin_p('CM_ZERO_SAPABOVE',cm_zero_sapabove) 
     1406     ! 
     1407     !Config  Key  = CM_ZERO_SAPBELOW 
     1408     !Config  Desc = maintenance respiration coefficient (g/g/day) at 0 deg C, for sapwood below, tabulated 
     1409     !Config  if  = OK_STOMATE 
     1410     !Config  Def  = undef,   1.19E-4, 1.19E-4, 1.19E-4, 1.19E-4, 1.19E-4, 1.19E-4,1.19E-4, 1.19E-4, 1.19E-4, 1.19E-4, 1.19E-4, 1.19E-4  
     1411     !Config  Help = 
     1412     !Config  Units = g/g/day 
     1413     CALL getin_p('CM_ZERO_SAPBELOW',cm_zero_sapbelow) 
     1414     ! 
     1415     !Config  Key  = CM_ZERO_HEARTABOVE 
     1416     !Config  Desc = maintenance respiration coefficient at 0 deg C, for heartwood above, tabulated 
     1417     !Config  if  = OK_STOMATE  
     1418     !Config  Def  =  undef, 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.  
     1419     !Config  Help = 
     1420     !Config  Units = g/g/day 
     1421     CALL getin_p('CM_ZERO_HEARTABOVE',cm_zero_heartabove) 
     1422     ! 
     1423     !Config  Key  = CM_ZERO_HEARTBELOW 
     1424     !Config  Desc = maintenance respiration coefficient at 0 deg C,for heartwood below, tabulated 
     1425     !Config  if  = OK_STOMATE  
     1426     !Config  Def  = undef, 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.  
     1427     !Config  Help = 
     1428     !Config  Units = g/g/day  
     1429     CALL getin_p('CM_ZERO_HEARTBELOW',cm_zero_heartbelow) 
     1430     ! 
     1431     !Config  Key  = CM_ZERO_ROOT 
     1432     !Config  Desc = maintenance respiration coefficient at 0 deg C, for roots, tabulated 
     1433     !Config  if  = OK_STOMATE 
     1434     !Config  Def  = undef,1.67E-3, 1.67E-3, 1.67E-3, 1.67E-3, 1.67E-3, 1.67E-3,1.67E-3, 1.67E-3, 1.67E-3, 1.67E-3, 1.67E-3, 1.67E-3 
     1435     !Config  Help = 
     1436     !Config  Units = g/g/day  
     1437     CALL getin_p('CM_ZERO_ROOT',cm_zero_root) 
     1438     ! 
     1439     !Config  Key  = CM_ZERO_FRUIT 
     1440     !Config  Desc = maintenance respiration coefficient at 0 deg C, for fruits, tabulated 
     1441     !Config  if  = OK_STOMATE 
     1442     !Config  Def  = undef, 1.19E-4, 1.19E-4, 1.19E-4, 1.19E-4, 1.19E-4, 1.19E-4,1.19E-4, 1.19E-4, 1.19E-4, 1.19E-4, 1.19E-4, 1.19E-4     
     1443     !Config  Help = 
     1444     !Config  Units = g/g/day  
     1445     CALL getin_p('CM_ZERO_FRUIT',cm_zero_fruit) 
     1446     ! 
     1447     !Config  Key  = CM_ZERO_CARBRES 
     1448     !Config  Desc = maintenance respiration coefficient at 0 deg C, for carbohydrate reserve, tabulated 
     1449     !Config  if  = OK_STOMATE 
     1450     !Config  Def  = undef, 1.19E-4, 1.19E-4, 1.19E-4, 1.19E-4, 1.19E-4, 1.19E-4,1.19E-4, 1.19E-4, 1.19E-4, 1.19E-4, 1.19E-4, 1.19E-4 
     1451     !Config  Help = 
     1452     !Config  Units = g/g/day  
     1453     CALL getin_p('CM_ZERO_CARBRES',cm_zero_carbres) 
     1454     !- 
     1455     ! Fire - stomate 
     1456     !- 
     1457     ! 
     1458     !Config  Key  = FLAM 
     1459     !Config  Desc = flamability: critical fraction of water holding capacity 
     1460     !Config  if  = OK_STOMATE 
     1461     !Config  Def  = undef,.15,.25,.25,.25,.25,.25,.25,.25,.25,.25,.35,.35 
     1462     !Config  Help = 
     1463     !Config  Units = 
     1464     CALL getin_p('FLAM',flam) 
     1465     ! 
     1466     !Config  Key  = RESIST 
     1467     !Config  Desc = fire resistance 
     1468     !Config  if  = OK_STOMATE 
     1469     !Config  Def  = undef, .95,.90,.12,.50,.12,.12,.12,.12,.0,.0,.0,.0  
     1470     !Config  Help = 
     1471     !Config  Units = 
     1472     CALL getin_p('RESIST',resist) 
     1473     !- 
     1474     ! Flux - LUC 
     1475     !- 
     1476     ! 
     1477     !Config  Key  = COEFF_LCCHANGE_1 
     1478     !Config  Desc = Coeff of biomass export for the year 
     1479     !Config  if  = OK_STOMATE 
     1480     !Config  Def  = undef,0.597,0.597,0.597,0.597,0.597,0.597, 0.597, 0.597, 0.597, 0.597, 0.597, 0.597  
     1481     !Config  Help = 
     1482     !Config  Units = 
     1483     CALL getin_p('COEFF_LCCHANGE_1',coeff_lcchange_1) 
     1484     ! 
     1485     !Config  Key  = COEFF_LCCHANGE_10 
     1486     !Config  Desc = Coeff of biomass export for the decade 
     1487     !Config  if  = OK_STOMATE 
     1488     !Config  Def  = undef,0.403,0.403,0.299,0.299,0.299,0.299,0.299,0.299,0.299,0.403,0.299,0.403 
     1489     !Config  Help = 
     1490     !Config  Units = 
     1491     CALL getin_p('COEFF_LCCHANGE_10',coeff_lcchange_10) 
     1492     ! 
     1493     !Config  Key  = COEFF_LCCHANGE_100 
     1494     !Config  Desc = Coeff of biomass export for the century 
     1495     !Config  if  = OK_STOMATE 
     1496     !Config  Def  = undef, 0.,0.,0.104,0.104,0.104,0.104,0.104,0.104,0.104, 0.,0.104,0. 
     1497     !Config  Help = 
     1498     !Config  Units = 
     1499     CALL getin_p('COEFF_LCCHANGE_100',coeff_lcchange_100) 
     1500     !- 
     1501     ! Phenology 
     1502     !- 
     1503     ! 
     1504     !Config  Key  = LAI_MAX 
     1505     !Config  Desc = maximum LAI, PFT-specific 
     1506     !Config  if  = OK_STOMATE 
     1507     !Config  Def  = undef, 7., 7., 5., 5., 5.,4.5, 4.5, 3.0, 2.5, 2.5, 5.,5.  
     1508     !Config  Help = 
     1509     !Config  Units = 
     1510     CALL getin_p('LAI_MAX',lai_max) 
     1511     ! 
     1512     !Config  Key  = PHENO_MODEL 
     1513     !Config  Desc = which phenology model is used? (tabulated)  
     1514     !Config  if  = OK_STOMATE 
     1515     !Config  Def  =  'none  ', 'none  ', 'moi   ', 'none  ','none  ','ncdgdd','none  ','ncdgdd','ngd   ','moigdd','moigdd','moigdd','moigdd' 
     1516     !Config  Help = 
     1517     !Config  Units = NONE 
     1518     CALL getin('PHENO_MODEL',pheno_model) 
     1519     ! 
     1520     !Config  Key  = PHENO_TYPE 
     1521     !Config  Desc = type of phenology, 0=bare ground 1=evergreen,  2=summergreen,  3=raingreen,  4=perennial 
     1522     !Config  if  = OK_STOMATE 
     1523     !Config  Def  =  0, 1, 3, 1, 1, 2, 1, 2, 2, 4, 4, 2, 3 
     1524     !Config  Help = 
     1525     !Config  Units = NONE 
     1526     CALL getin_p('PHENO_TYPE',pheno_type) 
     1527     !- 
     1528     ! Phenology : Leaf Onset 
     1529     !- 
     1530     ! 
     1531     !Config  Key  = PHENO_GDD_CRIT_C 
     1532     !Config  Desc = critical gdd, tabulated (C), constant c of aT^2+bT+c 
     1533     !Config  if  = OK_STOMATE 
     1534     !Config  Def  = undef, undef, undef, undef, undef, undef, undef, undef, undef, 270., 400., 125., 400. 
     1535     !Config  Help = 
     1536     !Config  Units = NONE 
     1537     CALL getin_p('PHENO_GDD_CRIT_C',pheno_gdd_crit_c) 
     1538     ! 
     1539     !Config  Key  = PHENO_GDD_CRIT_B 
     1540     !Config  Desc = critical gdd, tabulated (C), constant b of aT^2+bT+c 
     1541     !Config  if  = OK_STOMATE 
     1542     !Config  Def  = undef, undef, undef, undef, undef, undef, undef,undef, undef, 6.25, 0., 0., 0. 
     1543     !Config  Help = 
     1544     !Config  Units = NONE 
     1545     CALL getin_p('PHENO_GDD_CRIT_B',pheno_gdd_crit_b) 
     1546     ! 
     1547     !Config  Key  = PHENO_GDD_CRIT_A 
     1548     !Config  Desc = critical gdd, tabulated (C), constant a of aT^2+bT+c 
     1549     !Config  if  = OK_STOMATE 
     1550     !Config  Def  = undef, undef, undef, undef, undef, undef, undef, undef, undef, 0.03125,  0., 0., 0. 
     1551     !Config  Help = 
     1552     !Config  Units = NONE 
     1553     CALL getin_p('PHENO_GDD_CRIT_A',pheno_gdd_crit_a) 
     1554     ! 
     1555     !Config  Key  = NGD_CRIT 
     1556     !Config  Desc =  critical ngd, tabulated. Threshold -5 degrees 
     1557     !Config  if  = OK_STOMATE 
     1558     !Config  Def  = undef, undef, undef, undef, undef, undef, undef, 0., undef, undef, undef, undef, undef 
     1559     !Config  Help = 
     1560     !Config  Units = 
     1561     CALL getin_p('NGD_CRIT',ngd_crit) 
     1562     ! 
     1563     !Config  Key  = NCDGDD_TEMP 
     1564     !Config  Desc = critical temperature for the ncd vs. gdd function in phenology 
     1565     !Config  if  = OK_STOMATE 
     1566     !Config  Def  = undef, undef, undef, undef, undef, 5., undef, 0., undef, undef, undef, undef, undef 
     1567     !Config  Help = 
     1568     !Config  Units = celsius degrees (C)  ? 
     1569     CALL getin_p('NCDGDD_TEMP', ncdgdd_temp) 
     1570     ! 
     1571     !Config  Key  = HUM_FRAC 
     1572     !Config  Desc = critical humidity (relative to min/max) for phenology 
     1573     !Config  if  = OK_STOMATE 
     1574     !Config  Def  = undef, undef, .5, undef, undef, undef, undef, undef,  undef, .5, .5, .5,.5      
     1575     !Config  Help = 
     1576     !Config  Units = 
     1577     CALL getin_p('HUM_FRAC', hum_frac) 
     1578     ! 
     1579     !Config  Key  = LOWGPP_TIME 
     1580     !Config  Desc = minimum duration of dormance (d) for phenology  
     1581     !Config  if  = OK_STOMATE 
     1582     !Config  Def  =  undef, undef, 30., undef, undef, 30., undef, 30., 30., 30., 30., 30., 30.   
     1583     !Config  Help = 
     1584     !Config  Units = days 
     1585     CALL getin_p('LOWGPP_TIME', lowgpp_time) 
     1586     ! 
     1587     !Config  Key  = HUM_MIN_TIME 
     1588     !Config  Desc = minimum time elapsed since moisture minimum 
     1589     !Config  if  = OK_STOMATE 
     1590     !Config  Def  =  undef, undef, 50., undef, undef, undef, undef, undef, undef, 35., 35., 75., 75. 
     1591     !Config  Help = 
     1592     !Config  Units = days (d) 
     1593     CALL getin_p('HUM_MIN_TIME', hum_min_time) 
     1594     ! 
     1595     !Config  Key  = TAU_SAP 
     1596     !Config  Desc = sapwood -> heartwood conversion time 
     1597     !Config  if  = OK_STOMATE 
     1598     !Config  Def  = undef, 730., 730., 730., 730., 730., 730., 730., 730., undef, undef, undef, undef 
     1599     !Config  Help = 
     1600     !Config  Units = days (d) 
     1601     CALL getin_p('TAU_SAP',tau_sap) 
     1602     ! 
     1603     !Config  Key  = TAU_FRUIT 
     1604     !Config  Desc =  fruit lifetime 
     1605     !Config  if  = OK_STOMATE 
     1606     !Config  Def  = undef, 90., 90., 90., 90., 90., 90., 90., 90., undef, undef, undef, undef 
     1607     !Config  Help = 
     1608     !Config  Units = days (d) 
     1609     CALL getin_p('TAU_FRUIT',tau_fruit) 
     1610     ! 
     1611     !Config  Key  = ECUREUIL 
     1612     !Config  Desc = fraction of primary leaf and root allocation put into reserve 
     1613     !Config  if  = OK_STOMATE 
     1614     !Config  Def  = undef, .0, 1.,.0,.0, 1., .0,1., 1., 1., 1., 1., 1. 
     1615     !Config  Help = 
     1616     !Config  Units = NONE 
     1617     CALL getin_p('ECUREUIL',ecureuil) 
     1618     ! 
     1619     !Config  Key  = ALLOC_MIN 
     1620     !Config  Desc = allocation above/below = f(age) - 30/01/04 NV/JO/PF 
     1621     !Config  if  = OK_STOMATE 
     1622     !Config  Def  = undef,  0.2, 0.2, 0.2, 0.2, 0.2, 0.2, 0.2, 0.2, undef, undef, undef, undef  
     1623     !Config  Help = 
     1624     !Config  Units = 
     1625     CALL getin_p('ALLOC_MIN',alloc_min) 
     1626     ! 
     1627     !Config  Key  = ALLOC_MAX 
     1628     !Config  Desc =  
     1629     !Config  if  = OK_STOMATE 
     1630     !Config  Def  = undef, 0.8, 0.8, 0.8, 0.8, 0.8, 0.8, 0.8, 0.8, undef, undef, undef, undef 
     1631     !Config  Help = 
     1632     !Config  Units = 
     1633     CALL getin_p('ALLOC_MAX',alloc_max) 
     1634     ! 
     1635     !Config  Key  = DEMI_ALLOC  
     1636     !Config  Desc =  
     1637     !Config  if  = OK_STOMATE 
     1638     !Config  Def  = undef, 5., 5., 5., 5., 5., 5., 5., 5., undef, undef, undef, undef 
     1639     !Config  Help = 
     1640     !Config  Units = 
     1641     CALL getin_p('DEMI_ALLOC',demi_alloc) 
     1642 
     1643     !>> DS new for merge in the trunk 
     1644     ! 15/06/2011 : add leaflife_mtc for the new formalism used for calculate sla 
     1645     ! 
     1646     !Config  Key  = LEAFLIFE_TAB 
     1647     !Config  Desc =  
     1648     !Config  if  = OK_STOMATE 
     1649     !Config  Def  = undef, .5, 2.,.33, 1., 2.,.33, 2., 2., 2., 2., 2., 2.  
     1650     !Config  Help = 
     1651     !Config  Units = 
     1652     CALL getin_p('LEAFLIFE_TAB',leaflife_tab) 
     1653 
     1654     !- 
     1655     ! Phenology : Senescence 
     1656     !- 
     1657     ! 
     1658     !Config  Key  = LEAFFALL 
     1659     !Config  Desc = length of death of leaves, tabulated  
     1660     !Config  if  = OK_STOMATE 
     1661     !Config  Def  = undef,     undef,     10.,   undef,   undef,     10.,   undef,10.,     10.,     10.,     10.,     10.,     10.  
     1662     !Config  Help = 
     1663     !Config  Units = days (d) 
     1664     CALL getin_p('LEAFFALL',leaffall) 
     1665     ! 
     1666     !Config  Key  = LEAFAGECRIT 
     1667     !Config  Desc = critical leaf age, tabulated (d) 
     1668     !Config  if  = OK_STOMATE 
     1669     !Config  Def  = undef,     730.,    180.,    910.,    730.,    180.,    910.,180.,    180.,    120.,    120.,    90.,    90.   
     1670     !Config  Help = 
     1671     !Config  Units = days (d) 
     1672     CALL getin_p('LEAFAGECRIT',leafagecrit)  
     1673     ! 
     1674     !Config  Key  = SENESCENCE_TYPE 
     1675     !Config  Desc = type of senescence, tabulated 
     1676     !Config  if  = OK_STOMATE 
     1677     !Config  Def  =  'none  ',  'none  ',   'dry   ',  'none  ',  'none  ', 'cold  ',  'none  ',   'cold  ',  'cold  ',  'mixed ','mixed ',  'mixed ',   'mixed '  
     1678     !Config  Help = 
     1679     !Config  Units = NONE 
     1680     CALL getin('SENESCENCE_TYPE', senescence_type)  
     1681     ! 
     1682     !Config  Key  = SENESCENCE_HUM 
     1683     !Config  Desc = critical relative moisture availability for senescence 
     1684     !Config  if  = OK_STOMATE 
     1685     !Config  Def  = undef, undef, .3, undef, undef, undef, undef, undef, undef, .2, .2, .3, .2  
     1686     !Config  Help = 
     1687     !Config  Units =  
     1688     CALL getin_p('SENESCENCE_HUM', senescence_hum) 
     1689     ! 
     1690     !Config  Key  = NOSENESCENCE_HUM 
     1691     !Config  Desc = relative moisture availability above which there is no humidity-related senescence 
     1692     !Config  if  = OK_STOMATE 
     1693     !Config  Def  = undef, undef, .8, undef, undef, undef, undef, undef, undef, .3, .3, .3, .3  
     1694     !Config  Help = 
     1695     !Config  Units = 
     1696     CALL getin_p('NOSENESCENCE_HUM', nosenescence_hum)  
     1697     ! 
     1698     !Config  Key  = MAX_TURNOVER_TIME 
     1699     !Config  Desc = maximum turnover time for grasse 
     1700     !Config  if  = OK_STOMATE 
     1701     !Config  Def  = undef,  undef, undef, undef, undef, undef, undef, undef, undef,  80.,  80., 80., 80.  
     1702     !Config  Help = 
     1703     !Config  Units = days (d) 
     1704     CALL getin_p('MAX_TURNOVER_TIME',max_turnover_time) 
     1705     ! 
     1706     !Config  Key  = MIN_TURNOVER_TIME 
     1707     !Config  Desc = minimum turnover time for grasse  
     1708     !Config  if  = OK_STOMATE 
     1709     !Config  Def  = undef, undef, undef, undef, undef, undef, undef, undef, undef, 10., 10., 10., 10.  
     1710     !Config  Help = 
     1711     !Config  Units = days (d) 
     1712     CALL getin_p('MIN_TURNOVER_TIME',min_turnover_time) 
     1713     ! 
     1714     !Config  Key  = MIN_LEAF_AGE_FOR_SENESCENCE 
     1715     !Config  Desc =  minimum leaf age to allow senescence g 
     1716     !Config  if  = OK_STOMATE 
     1717     !Config  Def  = undef, undef, 90., undef, undef, 90., undef, 60., 60., 30., 30., 30., 30. 
     1718     !Config  Help = 
     1719     !Config  Units = days ? 
     1720     CALL getin_p('MIN_LEAF_AGE_FOR_SENESCENCE', min_leaf_age_for_senescence) 
     1721     ! 
     1722     !Config  Key  = SENESCENCE_TEMP_C 
     1723     !Config  Desc = critical temperature for senescence (C), constant c of aT^2+bT+c, tabulated 
     1724     !Config  if  = OK_STOMATE 
     1725     !Config  Def  = undef, undef, undef, undef, undef, 12., undef, 7., 2., -1.375, 5., 5., 10. 
     1726     !Config  Help = 
     1727     !Config  Units = NONE 
     1728     CALL getin_p('SENESCENCE_TEMP_C',senescence_temp_c) 
     1729     ! 
     1730     !Config  Key  = SENESCENCE_TEMP_B 
     1731     !Config  Desc = critical temperature for senescence (C), constant b of aT^2+bT+c ,tabulated 
     1732     !Config  if  = OK_STOMATE  
     1733     !Config  Def  =  undef, undef, undef, undef, undef, 0., undef, 0., 0., .1, 0., 0., 0. 
     1734     !Config  Help = 
     1735     !Config  Units = NONE 
     1736     CALL getin_p('SENESCENCE_TEMP_B',senescence_temp_b) 
     1737     ! 
     1738     !Config  Key  = SENESCENCE_TEMP_A 
     1739     !Config  Desc = critical temperature for senescence (C), constant a of aT^2+bT+c , tabulated 
     1740     !Config  if  = OK_STOMATE 
     1741     !Config  Def  = undef, undef, undef, undef, undef, 0., undef, 0., 0.,.00375, 0., 0., 0.  
     1742     !Config  Help = 
     1743     !Config  Units = NONE 
     1744     CALL getin_p('SENESCENCE_TEMP_A',senescence_temp_a) 
     1745     !- 
     1746     ! DGVM 
     1747     !- 
     1748     ! 
     1749     !Config  Key  = RESIDENCE_TIME 
     1750     !Config  Desc =  residence time of trees 
     1751     !Config  if  = OK_DGVM AND .NOT. LPJ_GAP_CONST_MORT 
     1752     !Config  Def  =  undef, 30.0, 30.0, 40.0, 40.0, 40.0, 80.0, 80.0, 80.0, 0.0, 0.0, 0.0, 0.0  
     1753     !Config  Help = 
     1754     !Config  Units = years (y) 
     1755     CALL getin_p('RESIDENCE_TIME',residence_time) 
     1756     ! 
     1757     !Config  Key  = TMIN_CRIT 
     1758     !Config  Desc = critical tmin, tabulated (C) 
     1759     !Config  if  = OK_STOMATE 
     1760     !Config  Def  = undef,  0.0, 0.0, -30.0, -14.0, -30.0, -45.0, -45.0, undef, undef, undef, undef, undef 
     1761     !Config  Help = 
     1762     !Config  Units = Celsius degrees 
     1763     CALL getin_p('TMIN_CRIT',tmin_crit) 
     1764     ! 
     1765     !Config  Key  = TCM_CRIT 
     1766     !Config  Desc = critical tcm, tabulated (C) 
     1767     !Config  if  = OK_STOMATE 
     1768     !Config  Def  =  undef, undef, undef, 5.0, 15.5, 15.5, -8.0, -8.0, -8.0, undef, undef,  undef, undef 
     1769     !Config  Help = 
     1770     !Config  Units = Celsius degrees 
     1771     CALL getin_p('TCM_CRIT',tcm_crit) 
    11251772       
    11261773     first_call = .FALSE. 

  • branches/ORCHIDEE_EXT/ORCHIDEE/src_stomate/lpj_gap.f90

    r257 r350  
    138138             ELSEWHERE 
    139139 
    140                 vigour(:) = 0.0 
     140                vigour(:) = zero 
    141141 
    142142             ENDWHERE 

  • branches/ORCHIDEE_EXT/ORCHIDEE/src_stomate/stomate_lcchange.f90

    r257 r350  
    297297          flux10(i,m)     =  flux10(i,m-1) 
    298298 
    299           IF (prod10(i,m) .LT. 1.0) prod10(i,m) = 0.0 
     299          IF (prod10(i,m) .LT. 1.0) prod10(i,m) = zero 
    300300!MM => quid de prod10_total ???  
    301301       ENDDO 
    … …  
    314314          flux100(i,m)      =  flux100(i,m-1) 
    315315 
    316           IF (prod100(i,m).LT.1.0) prod100(i,m) = 0.0 
     316          IF (prod100(i,m).LT.1.0) prod100(i,m) = zero 
    317317 
    318318       ENDDO 

  • branches/ORCHIDEE_EXT/ORCHIDEE/src_stomate/stomate_litter.f90

    r257 r350  
    126126    REAL(r_std), DIMENSION(npts)                                       :: zdiff_min 
    127127    ! for messages 
    128     CHARACTER*10, DIMENSION(nlitt)                                    :: litter_str 
    129     CHARACTER*22, DIMENSION(nparts)                                   :: part_str 
    130     CHARACTER*7, DIMENSION(ncarb)                                     :: carbon_str 
    131     CHARACTER*5, DIMENSION(nlevs)                                     :: level_str 
     128    CHARACTER(LEN=10), DIMENSION(nlitt)                                    :: litter_str 
     129    CHARACTER(LEN=22), DIMENSION(nparts)                                   :: part_str 
     130    CHARACTER(LEN=7), DIMENSION(ncarb)                                     :: carbon_str 
     131    CHARACTER(LEN=5), DIMENSION(nlevs)                                     :: level_str 
    132132    ! Indices 
    133133    INTEGER(i_std)                                                    :: i,j,k,l,m 
    … …  
    445445    ! 4.2.2 integrate over the nbdl levels 
    446446 
    447     tsoil_decomp(:) = 0.0 
     447    tsoil_decomp(:) = zero 
    448448 
    449449    DO l = 1, nbdl 
    … …  
    477477    ! 5.2.2 integrate over the nbdl levels 
    478478 
    479     soilhum_decomp(:) = 0.0 
     479    soilhum_decomp(:) = zero 
    480480 
    481481    DO l = 1, nbdl 

  • branches/ORCHIDEE_EXT/ORCHIDEE/src_stomate/stomate_phenology.f90

    r304 r350  
    331331             ! 
    332332 
    333              when_growthinit(i,j) = 0.0 
     333             when_growthinit(i,j) = zero 
    334334 
    335335          ENDIF    ! start of the growing season 

  • branches/ORCHIDEE_EXT/ORCHIDEE/src_stomate/stomate_resp.f90

    r257 r350  
    121121       ! 1.3.2 integrate over the nbdl levels 
    122122 
    123        t_root(:,j) = 0.0 
     123       t_root(:,j) = zero 
    124124 
    125125       DO l = 1, nbdl 

  • branches/ORCHIDEE_EXT/ORCHIDEE/src_stomate/stomate_soilcarbon.f90

    r257 r350  
    8080    REAL(r_std), DIMENSION(npts,ncarb,ncarb)                           :: flux 
    8181    ! for messages 
    82     CHARACTER*7, DIMENSION(ncarb)                                     :: carbon_str 
     82    CHARACTER(LEN=7), DIMENSION(ncarb)                                     :: carbon_str 
    8383    ! Indices 
    8484    INTEGER(i_std)                                                    :: k,kk,m 

  • branches/ORCHIDEE_EXT/ORCHIDEE/src_stomate/stomate_turnover.f90

    r257 r350  
    166166    ! 
    167167 
    168     leaf_meanage(:,:) = 0.0 
     168    leaf_meanage(:,:) = zero 
    169169 
    170170    DO m = 1, nleafages 
    … …  
    348348 
    349349       ! initialize leaf mass change in age class 
    350        delta_lm(:,:) = 0.0 
     350       delta_lm(:,:) = zero 
    351351 
    352352       IF ( tree(j) ) THEN 
    … …  
    451451       DO m = 1, nleafages 
    452452 
    453           WHERE ( biomass(:,j,ileaf) .GT. 0.0 ) 
     453          WHERE ( biomass(:,j,ileaf) .GT. zero ) 
    454454             leaf_frac(:,j,m) = ( leaf_frac(:,j,m)*lm_old(:) + delta_lm(:,m) ) / biomass(:,j,ileaf) 
    455455          ELSEWHERE 
Note: See TracChangeset for help on using the changeset viewer.