public void CalcEo_AtmosphericPotential(MetData Met, CanopyData Canopy) { //Get Eo and assign it to the public Eo field for this object. //private void soilwat2_priestly_taylor() // { double albedo; //! albedo taking into account plant material double eeq; //! equilibrium evaporation rate (mm) double wt_ave_temp; //! weighted mean temperature for the day (oC) //* ******* calculate potential evaporation from soil surface (eos) ****** // ! find equilibrium evap rate as a // ! function of radiation, albedo, and temp. albedo = cons.max_albedo - (cons.max_albedo - salb) * (1.0 - Canopy.cover_green_sum); // ! wt_ave_temp is mean temp, weighted towards max. wt_ave_temp = (0.60 * Met.maxt) + (0.40 * Met.mint); eeq = Met.radn * 23.8846 * (0.000204 - 0.000183 * albedo) * (wt_ave_temp + 29.0); //! find potential evapotranspiration (eo) from equilibrium evap rate Eo = eeq * soilwat2_eeq_fac(Met); // } }
public void CalcCoverForRunoff(CanopyData Canopy, SurfaceCoverData SurfaceCover) { //private void soilwat2_cover_surface_runoff() // { //This does NOT calculate runoff. It calculates an effective cover that is used for runoff. //In the process event this is called before the soilwat2_runoff. //*+ Purpose //* calculate the effective runoff cover //*+ Assumptions //* Assumes that if canopy height is negative it is missing. //*+ Mission Statement //* Calculate the Effective Runoff surface Cover double canopyfact; //! canopy factor (0-1) int crop; //! crop number double effective_crop_cover; //! effective crop cover (0-1) double cover_surface_crop; //! efective total cover (0-1) //! cover cn response from perfect - ML & dms 7-7-95 //! nb. perfect assumed crop canopy was 1/2 effect of mulch //! This allows the taller canopies to have less effect on runoff //! and the cover close to ground to have full effect (jngh) //! weight effectiveness of crop canopies //! 0 (no effect) to 1 (full effect) cover_surface_crop = 0.0; for (crop = 0; crop < Canopy.NumberOfCrops; crop++) { if (Canopy.canopy_height[crop] >= 0.0) { bool bDidInterpolate; canopyfact = MathUtilities.LinearInterpReal(Canopy.canopy_height[crop], cons.canopy_fact_height, cons.canopy_fact, out bDidInterpolate); } else { canopyfact = cons.canopy_fact_default; } effective_crop_cover = Canopy.cover_tot[crop] * canopyfact; cover_surface_crop = add_cover(cover_surface_crop, effective_crop_cover); } //! add cover known to affect runoff //! ie residue with canopy shading residue cover_surface_runoff = add_cover(cover_surface_crop, SurfaceCover.surfaceom_cover); }
public void CalcCoverForRunoff(CanopyData Canopy, SurfaceCoverData SurfaceCover) { //private void soilwat2_cover_surface_runoff() // { //This does NOT calculate runoff. It calculates an effective cover that is used for runoff. //In the process event this is called before the soilwat2_runoff. //*+ Purpose //* calculate the effective runoff cover //*+ Assumptions //* Assumes that if canopy height is negative it is missing. //*+ Mission Statement //* Calculate the Effective Runoff surface Cover double canopyfact; //! canopy factor (0-1) int crop; //! crop number double effective_crop_cover; //! effective crop cover (0-1) double cover_surface_crop; //! efective total cover (0-1) //! cover cn response from perfect - ML & dms 7-7-95 //! nb. perfect assumed crop canopy was 1/2 effect of mulch //! This allows the taller canopies to have less effect on runoff //! and the cover close to ground to have full effect (jngh) //! weight effectiveness of crop canopies //! 0 (no effect) to 1 (full effect) cover_surface_crop = 0.0; for (crop = 0; crop < Canopy.NumberOfCrops; crop++) { if (Canopy.canopy_height[crop] >= 0.0) { bool bDidInterpolate; canopyfact = MathUtilities.LinearInterpReal(Canopy.canopy_height[crop], cons.canopy_fact_height, cons.canopy_fact, out bDidInterpolate); } else { canopyfact = cons.canopy_fact_default; } effective_crop_cover = Canopy.cover_tot[crop] * canopyfact; cover_surface_crop = add_cover(cover_surface_crop, effective_crop_cover); } //! add cover known to affect runoff //! ie residue with canopy shading residue cover_surface_runoff = add_cover(cover_surface_crop, SurfaceCover.surfaceom_cover); }
private void OnSimulationCommencing(object sender, EventArgs e) { SaveModuleConstants(); //daily inputs met = new MetData(); irrig = new IrrigData(); canopy = new CanopyData(); surfaceCover = new SurfaceCoverData(); //optional daily inputs runon = 0.0; interception = 0.0; residueinterception = 0.0; if (Soil.Thickness != null) { try { SoilObject = new SoilWaterSoil(constants, Soil); //constructor can throw an Exception surfaceFactory = new SurfaceFactory(); surface = surfaceFactory.GetSurface(SoilObject, Clock); //constructor can throw an Exception (Evap class constructor too) //optional inputs (array) inflow_lat = null; } catch (Exception Ex) { throw new ApsimXException(this, Ex.Message); //catch any constructor Exceptions and rethrow as an ApsimXException. } } else { throw new ApsimXException(this, "SoilWater module has detected that the Soil has no layers."); } }
public void CalcEos_EoReducedDueToShading(double Eo, CanopyData Canopy, SurfaceCoverData SurfaceCover) { //private void soilwat2_pot_soil_evaporation() // { //eos -> ! (output) potential soil evap after modification for crop cover & residue_w double eos_canopy_fract; //! fraction of potential soil evaporation limited by crop canopy (mm) double eos_residue_fract; //! fraction of potential soil evaporation limited by crop residue (mm) //! 1. get potential soil water evaporation //!---------------------------------------+ //! reduce Eo to that under plant CANOPY <DMS June 95> //!---------------------------------------+ //! Based on Adams, Arkin & Ritchie (1976) Soil Sci. Soc. Am. J. 40:436- //! Reduction in potential soil evaporation under a canopy is determined //! the "% shade" (ie cover) of the crop canopy - this should include th //! green & dead canopy ie. the total canopy cover (but NOT near/on-grou //! residues). From fig. 5 & eqn 2. <dms June 95> //! Default value for c%canopy_eos_coef = 1.7 //! ...minimum reduction (at cover =0.0) is 1.0 //! ...maximum reduction (at cover =1.0) is 0.183. eos_canopy_fract = Math.Exp(-1 * cons.canopy_eos_coef * Canopy.cover_tot_sum); // !-----------------------------------------------+ // ! reduce Eo under canopy to that under mulch <DMS June 95> // !-----------------------------------------------+ // !1a. adjust potential soil evaporation to account for // ! the effects of surface residue (Adams et al, 1975) // ! as used in Perfect // ! BUT taking into account that residue can be a mix of // ! residues from various crop types <dms june 95> if (SurfaceCover.surfaceom_cover >= 1.0) { //! We test for 100% to avoid log function failure. //! The algorithm applied here approaches 0 as cover approaches //! 100% and so we use zero in this case. eos_residue_fract = 0.0; } else { //! Calculate coefficient of residue_wt effect on reducing first //! stage soil evaporation rate //! estimate 1st stage soil evap reduction power of //! mixed residues from the area of mixed residues. //! [DM. Silburn unpublished data, June 95 ] //! <temporary value - will reproduce Adams et al 75 effect> //! c%A_to_evap_fact = 0.00022 / 0.0005 = 0.44 eos_residue_fract = Math.Pow((1.0 - SurfaceCover.surfaceom_cover), cons.A_to_evap_fact); } //! Reduce potential soil evap under canopy to that under residue (mulch) Eos = Eo * eos_canopy_fract * eos_residue_fract; }
public void CalcEo_AtmosphericPotential(MetData Met, CanopyData Canopy) { //Get Eo and assign it to the public Eo field for this object. //private void soilwat2_priestly_taylor() // { double albedo; //! albedo taking into account plant material double eeq; //! equilibrium evaporation rate (mm) double wt_ave_temp; //! weighted mean temperature for the day (oC) //* ******* calculate potential evaporation from soil surface (eos) ****** // ! find equilibrium evap rate as a // ! function of radiation, albedo, and temp. albedo = cons.max_albedo - (cons.max_albedo - salb) * (1.0 - Canopy.cover_green_sum); // ! wt_ave_temp is mean temp, weighted towards max. wt_ave_temp = (0.60 * Met.maxt) + (0.40 * Met.mint); eeq = Met.radn * 23.8846 * (0.000204 - 0.000183 * albedo) * (wt_ave_temp + 29.0); //! find potential evapotranspiration (eo) from equilibrium evap rate Eo = eeq * soilwat2_eeq_fac(Met); // } }
public void CalcEos_EoReducedDueToShading(CanopyData Canopy, SurfaceCoverData SurfaceCover) { //private void soilwat2_pot_soil_evaporation() // { //eos -> ! (output) potential soil evap after modification for crop cover & residue_w double eos_canopy_fract; //! fraction of potential soil evaporation limited by crop canopy (mm) double eos_residue_fract; //! fraction of potential soil evaporation limited by crop residue (mm) //! 1. get potential soil water evaporation //!---------------------------------------+ //! reduce Eo to that under plant CANOPY <DMS June 95> //!---------------------------------------+ //! Based on Adams, Arkin & Ritchie (1976) Soil Sci. Soc. Am. J. 40:436- //! Reduction in potential soil evaporation under a canopy is determined //! the "% shade" (ie cover) of the crop canopy - this should include th //! green & dead canopy ie. the total canopy cover (but NOT near/on-grou //! residues). From fig. 5 & eqn 2. <dms June 95> //! Default value for c%canopy_eos_coef = 1.7 //! ...minimum reduction (at cover =0.0) is 1.0 //! ...maximum reduction (at cover =1.0) is 0.183. eos_canopy_fract = Math.Exp(-1 * cons.canopy_eos_coef * Canopy.cover_tot_sum); // !-----------------------------------------------+ // ! reduce Eo under canopy to that under mulch <DMS June 95> // !-----------------------------------------------+ // !1a. adjust potential soil evaporation to account for // ! the effects of surface residue (Adams et al, 1975) // ! as used in Perfect // ! BUT taking into account that residue can be a mix of // ! residues from various crop types <dms june 95> if (SurfaceCover.surfaceom_cover >= 1.0) { //! We test for 100% to avoid log function failure. //! The algorithm applied here approaches 0 as cover approaches //! 100% and so we use zero in this case. eos_residue_fract = 0.0; } else { //! Calculate coefficient of residue_wt effect on reducing first //! stage soil evaporation rate //! estimate 1st stage soil evap reduction power of //! mixed residues from the area of mixed residues. //! [DM. Silburn unpublished data, June 95 ] //! <temporary value - will reproduce Adams et al 75 effect> //! c%A_to_evap_fact = 0.00022 / 0.0005 = 0.44 eos_residue_fract = Math.Pow((1.0 - SurfaceCover.surfaceom_cover), cons.A_to_evap_fact); } //! Reduce potential soil evap under canopy to that under residue (mulch) Eos = Eo * eos_canopy_fract * eos_residue_fract; }