--------------- 1. MAJOR CHANGE --------------- ISBA-CC spinup procedure, DIF coupling, and code optimization (see sections 3) Add parameterization on soil organic carbon impact on isba hydro and thermo parameters (CSOC option see section 4) Add the possibility to read a permafrost map and change some parameter according to this map (see section 5) Add isba water budget and new outputs using new LWATER_BUDGET key in NAM_DIAG_ISBAn (see section 7) Add (or change) some isba DIF outputs using new LSURF_MISC_DIF key in NAM_DIAG_ISBAn (see section 8) Averaged gird-box soil liquid and ice water content diag in kg/m2 and m3/m3 --------------- 2. MINOR CHANGE --------------- SOM abreviation (soil organic matter) expressed in % changed to SOC (soil organic carbon) expressed in kg/m² EVAPCOR snow correction not taken into account with ISBA-DIF in hydro_soildif.F90 ZTAUICE (:) = MAX(PTSTEP,PTAUICE(:)) for ISBA-FR in isba_fluxes.F90 = Caracteristic time for ice in force-restore if Ptstep > 3300s Waiting for MEB : Active Topmodel soil moisture adjustement if ISBA-DIF in isba_sgh_update.F90 Add LDIAG_FA_NOCOMPACT in NAM_IO_OFFLINE : no compaction for diag file (default = false) If true, fa diag files not compacted. -------------------------- 3. MAJOR CHANGE in ISBA-CC or A-gs -------------------------- ISBA-DIF coupling Top surface = 0.1 m => rapid carbon pool Sub surface = 1.0 m => slow carbon pool Code optimization Remove all power operation with real-type and replace by EXP(LOG) (more rapid). Attention, PQDGAMM, PQDGMES and PQDAMAX now expressed in Log of Q10 (constant in time) Add XCC_NITRO and XBIOMASST_LIM in MODD_CO2V_PAR and remove ZCC_NITRO and ZBIOMASST_LIM from nitro_decline.F90 Move ISBA-CC fluxes diag from diag_misc to diag_evap and new ISBA-CC outputs : GPP_ISBA, GPPC_ISBA, GPP_PATCH, GPPC_PATCH = gross primary production over tile nature (averaged, cumulated and/or per patch) R_AUTO_ISBA, RC_AUTO_ISBA, R_AUTO_PATCH, RC_AUTO_PATCH = autotrophic respiration over tile nature (averaged, cumulated and/or per patch) R_ECO_ISBA, RC_ECO_ISBA, R_ECO_PATCH, RC_ECO_PATCH = ecosystem respiration over tile nature (averaged, cumulated and/or per patch) IF PROVAR_TO_DIAG : IF CPHOTO=='NIT'.OR.CPHOTO=='NCB' BIOM??_ISBA = grid-average Biomass_reservoir for each?? reservoirs (kgDM/m2) IF CRESPSL=='CNT' LIT??_??_ISBA = grid-average Litter_pool for each ?? reservoirs (gC/m2) SCARB??_ISBA = grid-average Soil_carbon_pool for each ?? reservoirs (gC/m2) LIGSTR??_ISBA = grid-average Ratio_Lignin/Carbon_in_structural_litter (-) Minor: Add or change some dimension in comments New Soil and wood carbon spin up procedure (see C. Delire and/or spinup_max.F90 for more details): In new NAM_SPINUP_CARB: LOGICAL :: LSPINUPCARBS ! T: do the soil carb spinup, F: no (default) LOGICAL :: LSPINUPCARBW ! T: do the wood carb spinup, F: no (default) REAL :: XSPINMAXS ! max number of times where CARBON_SOIL subroutine is called for each timestep in simulation during acceleration procedure REAL :: XSPINMAXW ! max number of times where the wood carbon is accelerated INTEGER :: NNBYEARSPINS ! nbr years needed to reaches soil equilibrium INTEGER :: NNBYEARSPINW ! nbr years needed to reaches wood equilibrium Exemple : &NAM_SPINUP_CARB LSPINUPCARBS = T , => soil carb spinup activated LSPINUPCARBW = T , => wood carb spinup activated XSPINMAXS = 200. , => CARBON_SOIL called 200 in one time step NNBYEARSPINS = 100 , => during 100 years (equivalent 20 000 years soil carb spinup) XSPINMAXW = 10. , => wood carbon called 10 in one time step NNBYEARSPINW = 15 => during 15 years (equivalent 150 years soil carb spinup) / To be removed : main_carb_spinup.F90 and main_wood_spinup.F90 ----------------------------------- 4. MAJOR CHANGE due new CSOC option ----------------------------------- Only with ISBA-DIF ISBA parameterizations for soil thermal and hydraulic properties are modified to accommodate both mineral and organic carbon soils according to observations from Boelter (1969). Details in isba_soc_parameters.F90 Activated if CSOC=SGH in NAM_SGH_ISBAn and if SOC data given in PGD fields via NAM_ISBA In NAM_ISBA : CHARACTER(LEN=28) :: YSOC_TOP ! file name for organic carbon expressed in kg/m2 CHARACTER(LEN=28) :: YSOC_SUB ! file name for organic carbon expressed in kg/m2 CHARACTER(LEN=6) :: YSOCFILETYPE ! organic carbon data file type REAL :: XUNIF_SOC_TOP ! uniform value of organic carbon (kg/m2) REAL :: XUNIF_SOC_SUB ! uniform value of organic carbon (kg/m2) LOGICAL :: LIMP_SOC ! Imposed maps of organic carbon from same resolution pgd file Can be used with soc_top.dir (.hdr) and soc_sub.dir (.hdr) at 1km resolution, ex : &NAM_ISBA YSOC_TOP = 'soc_top' , YSOC_SUB = 'soc_sub' , YSOCFILETYPE = 'DIRECT' ------------------------------------- 5. MAJOR CHANGE due to permafrost map ------------------------------------- If a permafrost map is given in pgd, soil depth put 12m if permafrost area fraction is > to 0.25 (permafrost_depth.F90) In addition and waiting for new vegtypes with explicit representation of tundra, over the same area (fraction > 0.25), the jackson root extinction parameter is fixed to 0.914 (observed value for tundra) with ISBA-DIF In NAM_ISBA : CHARACTER(LEN=28) :: YPERM ! file name for permafrost distribution expressed in fraction CHARACTER(LEN=6) :: YPERMFILETYPE ! permafrost distribution data file type REAL :: XUNIF_PERM ! uniform value of permafrost distribution expressed in fraction LOGICAL :: LIMP_PERM ! Imposed maps of permafrost distribution from same resolution pgd file Can be used with perm_glo_10km.dir (.hdr) at 10km resolution (from data at 0.5° resolution), ex: &NAM_ISBA YPERM = 'perm_glo_10km' , YPERMFILETYPE = 'DIRECT' -------------- 6. NEW OUTPUTS -------------- MISC : ------ ALT_ISBA or ALT_PATCH = active layer thickness over permafrost (m) if DIF FLT_ISBA or FLT_PATCH = frozen layer thickness over non-permafrost (m) if DIF LAI_ISBA = grid-box averaged leaf area index (m2/m2) PGD (diag_pgd_isba) : --------------------- XFRACSOC = Soil organic carbon fraction FLUXES (diag_seb_isba) : ------------------------ LESL_ISBA, LESLC_ISBA, LESL_PATCH, LESLC_PATCH = liquid water evaporation over snow (ES, Crocus) for tile nature (averaged, cumulated and/or per patch) IRRIG_ISBA, IRRIG_PATCH, IRRIGC_ISBA, IRRIGC_PATCH = irrigation rate (averaged, cumulated and/or per patch) if LAGRIP (agricultural option) --------------- PROVAR_TO_DIAG -------------- WG_??_ISBA = averaged soil liquid water content (m3/m3) (with ?? = each layer) SOILM_??_ISBA = averaged soil moisture (kg/m2) WGI_??_ISBA = averaged soil ice water content (m3/m3) SOILI_??_ISBA = averaged soil ice content (kg/m2) --------------------------------------- 7. NEW OUTPUTS FOR LWATER_BUDGET = True --------------------------------------- 2 new routines : isba_budget.F90 and isba_budget_init.F90 RAINF_ISBA, RAINFC_ISBA = input rainfall rate (averaged or cumulated) SNOWF_ISBA, SNOWFC_ISBA = input snowfall rate (averaged or cumulated) DWG_ISBA, DWGC_ISBA, DWG_PATCH, DWGC_PATCH = change in liquid soil moisture (averaged or cumulated) DWGI_ISBA, DWGIC_ISBA, DWGI_PATCH, DWGIC_PATCH = change in solid soil moisture (averaged or cumulated) DWR_ISBA, DWRC_ISBA, DWR_PATCH, DWRC_PATCH = change in water on canopy (averaged or cumulated) DSWE_ISBA, DSWEC_ISBA, DSWE_PATCH, DSWEC_PATCH = change in snow water equivalent (averaged or cumulated) WATBUD_ISBA, WATBUDC_ISBA, WATBUD_PATCH, WATBUDC_PATCH = isba water budget as residue (averaged or cumulated) Note that WATBUD is computed in isba_budget.F90 as : DWG + DWGI + DWR + DSWE = input fluxes - output fluxes in kg/m2/s input fluxes = RAINF + SNOWF + IFLOOD + IRRIG_FLUX output fluxes = EVAP + DRAIN + RUNOFF + PFLOOD + ICEFLUX --------------------------------------- 8. NEW OUTPUTS FOR LSURF_MISC_DIF = True ---------------------------------------- TSWI_R_ISBA =total soil wetness index over the root zone (-) WGTOT_R_ISBA =total water content (liquid+solid) over the root zone (kg/m2) WGI_R_ISBA =total ice content (solid) over the root zone (kg/m2) TSWI_S_ISBA =total soil wetness index over the surface (-) WG_S_ISBA =liquid water content over the surface (m3/m3) WGI_S_ISBA =ice content over the surface (m3/m3) TSWI_D2_ISBA =total soil wetness index over comparable FR-DG2 reservoir (-) WG_D2_ISBA =liquid water content over comparable FR-DG2 reservoir (m3/m3) WGI_D2_ISBA =ice content over comparable FR-DG2 reservoir (m3/m3) TSWI_D3_ISBA =total soil wetness index over comparable FR-DG3 reservoir (-) WG_D3_ISBA =liquid water content over comparable FR-DG3 reservoir (m3/m3) WGI_D3_ISBA =ice content over comparable FR-DG3 reservoir (m3/m3)