CHANGELOG for Biome-BGCMuSo v6.2 2021-10-08 by Dóra HIDY -------------------------------- Changes from 6.1 to 6.2-b1 -------------------------------- - NEW DEVELOPMENT: calculation of potential evaporation (using stomatal conductance without soil moisture limitation) - BUG FIXING in cutdown_to_litter.c (litter fluxes of belowground cut-down material) - BUG FIXING in dayphen.c (incorrect value assignment of the variable offday) - BUG FIXING in fertilizing_init.c: free fertilizing date variables - DEVELOPMENT of tipping.c: tipping diffusion limitation in case of saturation - also percolation - DEVELOPMENT of Richards method: implementation of finer spatial resolution - calculation of actual value of hydrological diffusivity and conductivity -------------------------------- Changes from 6.2-b1 to 6.2-b2 -------------------------------- - BUG FIXING in richards.c: Richards finer spatial resolution - BUG FIXING in fertilizing_init: data reading - DEVELOPMENT of thinning: thinning also affects the root biomass of the trees (to the same extent as aboveground biomass) new output variables: root biomass to THN (carbon fluxes: [1187-1198], nitrogen fluxes: [2061-2072]) - BUG FIXING in potential evaporation calculation in canopy_et.c (using potential stom.cond. instead of actual) - MODIFICATION: instead of the calculation of hydraulic conductivity and diffusivity and the starting/end point of the simulation (finer temporal resolution), only averaged value is calculated (eliminated output variables for start/end hyrd.conduct/diffus [2559-2598], instead: averaged hyrd.conduct/diffus [2559-2578] -------------------------------- Changes from 6.2-b2 to 6.2-b5 -------------------------------- - DEVELOPMENT of actual evaporation calculation - presence of groundwater has also effect on pontential infiltration new output variable: soilw_from_GW0 [237] - DEVELOPMENT in groundwater.c: special case for GW on the ground - DEVELOPMENT of optional input reading: using year information in onday/offday/mortality/conductance/groundwater - DEVELOPMENT of Richards method: BUG FIXING of finer spatial resolution -------------------------------- Changes from 6.2-b5 to 6.2-b6 -------------------------------- - BUG FIXING: small bug with relatively large impact in daily_allocation.c - C allocation calculation is corrected in N-limited situations - DEVELOPMENT of Richards: BUG FIXING of finer spatial resolution - MODIFICATION: ceasing output variable: GDD previous day - GDDpre [77] - BUG FIXING in output_map_init.c: variables [2045-2071] were overwritten with [2405-2432]. - DEVELOPMENT of groundwater simulation recharging of soil moisture in groundwater (net loss for the system) if groundwater table is in actual layer (above lower boundary): GWlayer = actual layer, lower layers are charged soil layers below the groundwater table are saturated - net water gain from soil system new output variables for recharging fluxes and pools: GW_recharge [227-236], groundwater_snk [129], GWlayer [2797], GWD [2798], CapillFringe [2799], GWeff [2931-2940] thickness of capillary zone is calculated by the model based on soil type (and literautre data) CF is regulated based on the ratio of actual potential evaporation and theroetically maximum of potential evaporation two methods for GW-calculations method1: increasing of VWCfc, water from GW (net gain) method2: water from GW (net gain) depending on the position of capillary zone - BUG FIXING in phenology.c: first day after EMERGNECE period - BUG FIXING in multilayer_hydrolprocess.c: managing of POND WATER from soil and special bottom layer -------------------------------- Changes from 6.2-b6 to 6.2-b7 -------------------------------- - BUG FIXING in soilstress_calculation.c (calculation of m_SWCstressLENGTH paramter) - DEVELOPMENT: standing dead biomass turns into litter pool on the day of planting new output parameter: cumulative C loss due to planting (standing dead biomass turns into litter pool) cum_Closs_PLT [3117] - DEVELOPMENT in senescence.c: control step to avoid negative SMSI value -------------------------------- Changes from 6.2-b7 to 6.2-b8 -------------------------------- - BUG FIXING in planting.c: initialization of fruitC_HRV and vegC_HRV variables on the day of planting -------------------------------- Changes from 6.2-b8 to 6.2-b9 -------------------------------- - MODIFICATION of groundwater.c: no dependence of actual value of capillary fringe on soil evaporation - MODIFICATION in output_map_init.c: completing of the name of mineralized N variables - NH4-N/NO3-N instead of NH4 and NO3 -------------------------------- Changes from 6.2-b9 to 6.2-b10 -------------------------------- - DEVELOPMENT: introduction of extra parameters for soil calculations: optional extraSOIparam.txt file (similar to WPM or GW file) parameters of diffusion calculation in tipping.c parameters of temperature and pH function of nitrification in multilayer_sminn.c parameters of temperature and pH function of decomposition in decomp.c parameters of surface residue/mulch calculation in mulch.c parameters of groundwater calculations in groundwater.c - DEVELOPMENT: residue/mulch simulation is developed: estimation of aboveground litrC and cwdC - based on lit/cwd fluxes new output variables cwdc_above and litC_above [564-565] - MODIFICATION: new pH function of nitrification - MODIFICATION: modified T response function for decomposition and nitrification calculation - BUG FIXING in multilayer_rootdepth: in case of germ_layer=0 -------------------------------- Changes from 6.2-b10 to 6.2-b11 -------------------------------- - DEVELOPMENT of Richards method: optional finer spatial resolution (SHCM_flag=0) or optional normal spatial resolution (SHCM_flag=3) calculation of actual value of hydrological diffusivity and conductivity - DEVELOPMENT: option to REVERT to original response functions of nitrification and decomposition using extraSOIparam.txt - DEVELOPMENT: retranslocation of N during senescence of leaves, stems and fruits - result: varying C:N ratio of plant pools (dead+live) new output variables: SNSC_to_retrans [2433], leafSNSCgenprog_to_retrans [2434], leaf/froot/fruit/softstem C:N [3118-3121] - DEVELOPMENT of GW: GW above surface (negative GW-data) - > turns into pond water new output variable: GW_to_pondw [238] -------------------------------- Changes from 6.2-b11 to 6.2-b12 -------------------------------- - BUG FIXING in output_map_init.c: ts_nitrif variables [2589-2594] - MODIFICATION: of SOI parameter settings: BD can be set independent from critical VWC values - DEVELOPMENT: extension of log file with information about using of extraSOIparameters.txt and decompostion/nitrification Tresponse function - BUG FIXING: in groundwater.c: updating of increased VWC_fc in case of groundwater - BUG FIXING: in groundwater.c: using of GW_flag - BUG FIXING: in multilayer_rootdepth.c: upper limit for rooting depth (MRD from EPC file) - DEVELOPMENT: code optimalization - maximum rootzone is calculated only on first day instead of every day n_maxrootlayers is calculated in firstday.c instead of multilayer_rootdepth.c - MODIFICATION: rooting zone variables (VWC_RZ, PSI_RZ) are calculated not only vegetation period, but on every simulation day - MODIFICATION: default value of minimum temperature for decomposition and nitrification is changed from -10 to -3 Celsius -------------------------------- Changes from 6.2-b11 to 6.2-b12 -------------------------------- - DEVELOPMENT: new output variables for testing GSI-simulation - MODIFICATION: SHCM_flag (in INI file) - 0: tipping with diffusion, 1: tipping without diffusion, 2: Richards with normal spatial discretization, 3. Richards with finer spatial discretization -------------------------------- Changes from 6.2-b12 to 6.2-b13 -------------------------------- - MODIFICATION: default value of Tmin of decomposition and nitrification - change default value from -10 to -3 Celsius - MODIFICATION: calculation or rootzone value - instead of actual, maximum rooting depth values - BUG FIXING: calculation of assim_SScoeff in photosynthesis.c - DEVELOPMENT: heatsum growing indexes are available as output variable new output variables: tmin_index, vpd_index, dayl_index, gsi_indexAVG, heatsum_index [25-29] -------------------------------- Changes from 6.2-b13 to 6.2-b14 -------------------------------- - BUG FIXING: calculation of assim_SScoeff in photosynthesis.c - returning to b12 logic - DEVELOPMENT: heatsum is available as output variable [30] -------------------------------- Changes from 6.2-b14 to 6.2-b15 -------------------------------- - BUG FIXING: calculation of number of maximum rooting layers (problem occured only in case of crop rotation) - WARNING: in case of crop rotation every EPC file must contain information about germination and emergence phenophase (problem occured only if germination and emergence phenophase wa set to 0) - DEVELOPMENT: new output variables - soil organic carbon content in % layer by layer [3122-3126] -------------------------------- Changes from 6.2-b15 to 6.2-b16 -------------------------------- - BUG FIXING: in pond water filling - DEVELOPMENT: new output variable prcp_to_pondw [239]: water flux from precipitation to pond water) iterstep_Richards [130]: Number of iteration step using Richards-method - DEVELOPMENT: if RCN number is equal to 0 - no runoff - DEVELOPMENT (MAJOR CHANGE): option to use extended SOI file if extraSOIparameters.txt doesn't exist (see User's Guide) -------------------------------- Changes from 6.2-b16 to 6.2-b18 -------------------------------- - BUG FIXING: in pond water filling and runoff calculation - DEVELOPMENT: reorganized SOI file (new, logical order of SOI parameters) - DEVELOPMENT: new SOI parameter - critical WFPS value for denitrification -------------------------------- Changes from 6.2-b18 to 6.219 -------------------------------- - BUG FIXING: in pond water evaporation - DEVELOPMENT: new output variables to summarize management related input data FRZamount [3500] IRGamount [3501] plantDensity [3502] -------------------------------- Changes from 6.2-b18 to 6.2-b19 -------------------------------- - BUG FIXING: in bare soil evaporation - BUG FIXING: in bare precision control of coarse woody debris pool -------------------------------- Changes from 6.2-b19 to 6.2-b20 -------------------------------- - BUG FIXING: in actual evaporation calculation (dsr calculation) - BUG FIXING: in mortality calculations (m_livestemn_to_cwdn) - BUG FIXING: in precision control calculation (resetting cwdc below a critical value) -------------------------------- Changes from 6.2-b20 to 6.2-b21 -------------------------------- - DEVELOPMENT: new output file - .econout for crop managment information - BUG FIXING: in mortality calculations (m_livestemn_to_cwdn) - MODIFICATION: ceasing of i2o output struct (FRZamount,IRGamount,plantDensity) - it has become redundant due to econout file - MODIFICATION: in conditional management section of INI file (new: first befVWCratio_condIRG, second: befSMSIratio_condIRG) - BUG FIXING: in decomposition calculations (upper limit for potential litter loss) - MODIFICATION: changing the sequence of hydrological routines - first evaporation, after percolation - DEVELOPMENT: new output variables Amount of conditional irrigation [131] CN ratio error [566] - MODIFICATION: RZ variables refer to actual (not to maximum) rooting zone - DEVELOMENT:new output file for economic variables (.econout)- only in case of harvesting! - BUG FIXING: in litter decomposition problem: k value can be greater than 1 -> decomposition of given pool can larger than pool -> negative pool fixation: limitation of decomposition flux - BUG FIXING: in harvest (also leaf transfer pool is eliminated in case of harvest) - MODIFICATION: order of rountines in multilayer_hydrolprocess.c - now: evaporation before percolation - MODIFICATION: RZ variables refer to actual (not to maximum) rooting zone - DEVELOPMENT: new output data - CNratioERR [566] - MODIFICATION: order of lines in conditional management section of INI file (conditional IRRIGATION) condIRG_flag=1: condIRG based on SWC ratio, condIRG_flag=2: condIRG based on SMSI - DEVELOPMENT: control of conditional management - in case of SMSI-method, if SWC after condIRG is less than SWC before condIRG (due to the false settings) - > no conditional irrigation - MODIFICATION: unit of fruitC_HRV and vegC_HRV variables (558-559; from kgC/m2 to kgDM/t) - DEVELPMENT: new output variable - Cumulative amount of conditional irrigation [condIRGsrc;131] - DEVELOPMENT: soilEvapCUM1 can't be negative - control in potEVAP_to_actEVAP.c - BUG FIXING: control of CWD in precision_control.c - avoiding negative CWD pools in case of evergreen forests - BUG FIXING: in calculation of soilEvapCUM2 (in case of evaporation limitation) - BUG FIXING: in sprop_init.c - initialization of Capillary Fringe (resetting only if DATA GAP) -------------------------------- Changes from 6.2-b21 to 6.2-b22 -------------------------------- - MODIFICATION: in unit of fruitC_HRV and vegC_HRV (tC/ha to kgC/m2) - BUG FIXING: in calculation of primary product and secondary product in econout file - BUG FIXING: in calculation infiltration limit in pondANDrunoff.c - DEVELOPMENT: differentiation of Hortonian (prcp_to_runoff) and Dunnian (soilw_to_runoff) new output variables: pondw_to_runoff [240] - DEVELOPMENT: improving of Richard's method new output variables: infilt_to_pondw [241], infilt_to_soilw [242] -------------------------------- Changes from 6.2-b22 to 6.2-b23 -------------------------------- - DEVELOPMENT: new output variable for the end of wintering - BUG FIXING: in richards.c no extra limitation of soil evaporation and transpiration pond water calculation outside the internal cycle - MODIFICATION: soil hydrological calculation method flag in INI (0: tipping, 1: Richards, 2: tipping without diffusion) -------------------------------- Changes from 6.2-b23 to 6.2-b24 -------------------------------- - BUG FIXING: calculation of harvest index (DM of stem was missing from the equation) - BUG FIXING: in richards.c limitation of soil evaporation (only in top soil layer) pond water limitation - correction of typing error - WARNING HANDLING: elimination of variables set but not used -------------------------------- Changes from 6.2-b24 to 6.2-b25 -------------------------------- - BUG FIXING: balance error elimination in case of flux limitation - DEVELOPMENT: new output variables rootlength [2579], Rooting depth of the phenphase1-phenphase7 [2580-2586] - BUG FIXING: balance error elimination in planting.c (resetting of storage variables)