Changeset 428 for branches


Ignore:
Timestamp:
2011-08-31T11:26:58+02:00 (13 years ago)
Author:
didier.solyga
Message:

Create a new subroutine in pft_parameters for the parameters read whatever flag is activated

Location:
branches/ORCHIDEE_EXT/ORCHIDEE
Files:
2 edited

Legend:

Unmodified
Added
Removed

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

    r425 r428  
    10251025!= 
    10261026! 
    1027  SUBROUTINE getin_sechiba_pft_parameters 
    1028  
     1027 SUBROUTINE getin_pft_parameters 
     1028    
    10291029   IMPLICIT NONE 
    10301030   
    10311031   LOGICAL, SAVE ::  first_call = .TRUE. 
    10321032 
    1033   IF(first_call) THEN 
     1033   IF (first_call) THEN 
    10341034 
    10351035     !- 
    … …  
    11871187     !Config  Units =  
    11881188     CALL getin_p('EXT_COEFF',ext_coeff) 
    1189      !- 
    1190      ! Evapotranspiration -  sechiba 
    1191      !- 
    1192      ! 
    1193      !Config  Key  = RSTRUCT_CONST 
    1194      !Config  Desc = Structural resistance : one for each vegetation type 
    1195      !Config  if  = OK_SECHIBA 
    1196      !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 
    1197      !Config  Help = 
    1198      !Config  Units =  
    1199      CALL getin_p('RSTRUCT_CONST',rstruct_const) 
    1200      ! 
    1201      !Config  Key  = KZERO 
    1202      !Config  Desc =  A vegetation dependent constant used in the calculation of the surface resistance. 
    1203      !Config  if  = OK_SECHIBA 
    1204      !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  
    1205      !Config  Help = 
    1206      !Config  Units =  
    1207      CALL getin_p('KZERO',kzero) 
    1208      ! 
    1209      ! Ajouts Nathalie - le 28 Mars 2006 - sur conseils Fred Hourdin 
    1210      ! 
    1211      !Config  Key  = RVEG_PFT 
    1212      !Config  Desc = Artificial parameter to increase or decrease canopy resistance. 
    1213      !Config  if  = OK_SECHIBA 
    1214      !Config  Def  = 1. 
    1215      !Config  Help = This parameter is set by PFT. 
    1216      !Config  Units =  
    1217      CALL getin_p('RVEG_PFT', rveg_pft)     
     1189 
    12181190     !- 
    12191191     ! Water-hydrology - sechiba 
    12201192     !- 
    1221      ! 
    1222      !Config  Key  = WMAX_VEG 
    1223      !Config  Desc = Maximum field capacity for each of the vegetations (Temporary): max quantity of water 
    1224      !Config  if  = OK_SECHIBA 
    1225      !Config  Def  =  150., 150., 150., 150., 150., 150., 150.,150., 150., 150., 150., 150., 150. 
    1226      !Config  Help = 
    1227      !Config  Units = [Kg/M3] 
    1228      CALL getin_p('WMAX_VEG',wmax_veg) 
    12291193     ! 
    12301194     !Config  Key  = HYDROL_HUMCSTE 
    … …  
    12371201     !Config  Units =  
    12381202     CALL getin_p('HYDROL_HUMCSTE', humcste) 
     1203 
     1204     !- 
     1205     ! Soil - vegetation 
     1206     !- 
     1207     ! 
     1208     !Config  Key  = PREF_SOIL_VEG_SAND 
     1209     !Config  Desc = Table which contains the correlation between the soil types and vegetation type 
     1210     !Config  if  = OK_SECHIBA 
     1211     !Config  Def  = 1, 3, 3, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2 
     1212     !Config  Help = first layer of the soil 
     1213     !Config  Units =  
     1214     CALL getin_p('PREF_SOIL_VEG_SAND',pref_soil_veg_sand) 
     1215     ! 
     1216     !Config  Key  = PREF_SOIL_VEG_LOAN 
     1217     !Config  Desc = Table which contains the correlation between the soil types and vegetation type 
     1218     !Config  if  = OK_SECHIBA 
     1219     !Config  Def  = 2, 2, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3 
     1220     !Config  Help = second layer of the soil 
     1221     !Config  Units =  
     1222     CALL getin_p('PREF_SOIL_VEG_LOAN',pref_soil_veg_loan)  
     1223      ! 
     1224     !Config  Key  = PREF_SOIL_VEG_CLAY 
     1225     !Config  Desc = Table which contains the correlation between the soil types and vegetation type 
     1226     !Config  if  = OK_SECHIBA 
     1227     !Config  Def  = 3, 1, 1, 1, 1, 1 ,1 ,1 ,1 ,1 ,1 ,1, 1  
     1228     !Config  Help = third layer of the soil  
     1229     !Config  Units =         
     1230     CALL getin_p('PREF_SOIL_VEG_CLAY',pref_soil_veg_clay) 
     1231 
     1232     first_call = .FALSE. 
     1233 
     1234   ENDIF 
     1235 
     1236 END SUBROUTINE getin_pft_parameters 
     1237! 
     1238!= 
     1239! 
     1240 SUBROUTINE getin_sechiba_pft_parameters 
     1241 
     1242   IMPLICIT NONE 
     1243   
     1244   LOGICAL, SAVE ::  first_call = .TRUE. 
     1245 
     1246  IF (first_call) THEN 
     1247     !- 
     1248     ! Evapotranspiration -  sechiba 
     1249     !- 
     1250     ! 
     1251     !Config  Key  = RSTRUCT_CONST 
     1252     !Config  Desc = Structural resistance : one for each vegetation type 
     1253     !Config  if  = OK_SECHIBA 
     1254     !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 
     1255     !Config  Help = 
     1256     !Config  Units =  
     1257     CALL getin_p('RSTRUCT_CONST',rstruct_const) 
     1258     ! 
     1259     !Config  Key  = KZERO 
     1260     !Config  Desc =  A vegetation dependent constant used in the calculation of the surface resistance. 
     1261     !Config  if  = OK_SECHIBA 
     1262     !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  
     1263     !Config  Help = 
     1264     !Config  Units =  
     1265     CALL getin_p('KZERO',kzero) 
     1266     ! 
     1267     ! Ajouts Nathalie - le 28 Mars 2006 - sur conseils Fred Hourdin 
     1268     ! 
     1269     !Config  Key  = RVEG_PFT 
     1270     !Config  Desc = Artificial parameter to increase or decrease canopy resistance. 
     1271     !Config  if  = OK_SECHIBA 
     1272     !Config  Def  = 1. 
     1273     !Config  Help = This parameter is set by PFT. 
     1274     !Config  Units =  
     1275     CALL getin_p('RVEG_PFT', rveg_pft)     
     1276 
     1277     !- 
     1278     ! Water-hydrology - sechiba 
     1279     !- 
     1280     ! 
     1281     !Config  Key  = WMAX_VEG 
     1282     !Config  Desc = Maximum field capacity for each of the vegetations (Temporary): max quantity of water 
     1283     !Config  if  = OK_SECHIBA 
     1284     !Config  Def  =  150., 150., 150., 150., 150., 150., 150.,150., 150., 150., 150., 150., 150. 
     1285     !Config  Help = 
     1286     !Config  Units = [Kg/M3] 
     1287     CALL getin_p('WMAX_VEG',wmax_veg) 
    12391288     ! 
    12401289     !Config  Key  = PERCENT_THROUGHFALL_PFT 
    … …  
    12461295     !Config  Units = [Kg/M3] 
    12471296     CALL getin_p('PERCENT_TROUGHFALL_PFT',throughfall_by_pft) 
     1297 
    12481298     !- 
    12491299     ! Albedo - sechiba 
    … …  
    12811331     !Config  Units =  
    12821332     CALL getin_p('ALB_LEAF_NIR',alb_leaf_nir) 
    1283      !- 
    1284      ! Soil - vegetation 
    1285      !- 
    1286      ! 
    1287      !Config  Key  = PREF_SOIL_VEG_SAND 
    1288      !Config  Desc = Table which contains the correlation between the soil types and vegetation type 
    1289      !Config  if  = OK_SECHIBA 
    1290      !Config  Def  = 1, 3, 3, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2 
    1291      !Config  Help = first layer of the soil 
    1292      !Config  Units =  
    1293      CALL getin_p('PREF_SOIL_VEG_SAND',pref_soil_veg_sand) 
    1294      ! 
    1295      !Config  Key  = PREF_SOIL_VEG_LOAN 
    1296      !Config  Desc = Table which contains the correlation between the soil types and vegetation type 
    1297      !Config  if  = OK_SECHIBA 
    1298      !Config  Def  = 2, 2, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3 
    1299      !Config  Help = second layer of the soil 
    1300      !Config  Units =  
    1301      CALL getin_p('PREF_SOIL_VEG_LOAN',pref_soil_veg_loan)  
    1302       ! 
    1303      !Config  Key  = PREF_SOIL_VEG_CLAY 
    1304      !Config  Desc = Table which contains the correlation between the soil types and vegetation type 
    1305      !Config  if  = OK_SECHIBA 
    1306      !Config  Def  = 3, 1, 1, 1, 1, 1 ,1 ,1 ,1 ,1 ,1 ,1, 1  
    1307      !Config  Help = third layer of the soil  
    1308      !Config  Units =         
    1309      CALL getin_p('PREF_SOIL_VEG_CLAY',pref_soil_veg_clay) 
    13101333 
    13111334     first_call = .FALSE. 
    … …  
    13231346  LOGICAL, SAVE ::  first_call = .TRUE. 
    13241347 
    1325   IF(first_call) THEN 
     1348  IF (first_call) THEN 
    13261349 
    13271350     !- 
    … …  
    13441367     !Config  Units = [m**2/gC] 
    13451368     CALL getin_p('SLA',sla) 
     1369 
    13461370     !- 
    13471371     ! Photosynthesis 
    … …  
    14351459     !Config  Units = NONE 
    14361460     CALL getin_p('TPHOTO_MAX_C',tphoto_max_c) 
     1461 
    14371462     !- 
    14381463     ! Respiration - stomate 
    … …  
    15261551     !Config  Units = [g/g/day]  
    15271552     CALL getin_p('CM_ZERO_CARBRES',cm_zero_carbres) 
     1553 
    15281554     !- 
    15291555     ! Fire - stomate 
    … …  
    15451571     !Config  Units = 
    15461572     CALL getin_p('RESIST',resist) 
     1573 
    15471574     !- 
    15481575     ! Flux - LUC 
    … …  
    15721599     !Config  Units = 
    15731600     CALL getin_p('COEFF_LCCHANGE_100',coeff_lcchange_100) 
     1601 
    15741602     !- 
    15751603     ! Phenology 
    … …  
    15991627     !Config  Units = NONE 
    16001628     CALL getin_p('PHENO_TYPE',pheno_type) 
     1629 
    16011630     !- 
    16021631     ! Phenology : Leaf Onset 
    … …  
    18171846     !Config  Units = NONE 
    18181847     CALL getin_p('SENESCENCE_TEMP_A',senescence_temp_a) 
     1848 
    18191849     !- 
    18201850     ! DGVM 

  • branches/ORCHIDEE_EXT/ORCHIDEE/src_sechiba/intersurf.f90

    r426 r428  
    29472947 
    29482948    ! 6. Read the parameters in the run.def file  according the flags 
     2949 
     2950    IF (impose_param ) THEN  
     2951       CALL getin_pft_parameters 
     2952    ENDIF 
    29492953 
    29502954    IF ( control_flags%ok_sechiba ) THEN 
Note: See TracChangeset for help on using the changeset viewer.