|
Erosion new variables
|
Functions/Subroutines | |
| subroutine | erosion (minErosWindSped) |
| Main entry point to erosion submodel. | |
| subroutine | calcAeroRoug (ridgSpacParaWind, ridgHght, randRoug, aeroAnemFlg, resiLAI, resiSAI, resiAvgHght, cropLAI, cropSAI, cropHght, cropRowSpac, cropFurrFlg, aeroRougRidg, aeroRougRand, aeroRougSurf, aeroRougCnpy, aeroRougStatHght, bioDragCoef) |
| Calculates aerodynamic roughness parameter, which is aeroRougSurf if there is no standing biomass or aeroRougCnpy if standing biomass is present. | |
| subroutine | calcSurfThrsFricVelc (soilFracDiamLt84, soilLayrAggDen, soilCrstFrac, soilLayrRock, soilLoosCovFrac, bioFracFlatCovr, aeroRougSurf, soilSurfWatrCont, soilLayrWiltPt, fricVelc, soilFracDiamLt84ic, asoilLayrRock, soilMassAvalDelt, thrsFricVelc, thrsFricVelcTrap, fricVelcMod, soilFracDiamLt84mn, soilMassAvalAggSurf, soilMassAvalAggSurfmx, thrsFricVelcBareSoil, thrsFricVelcSurfCov, thrsFricVelcSurfWet, thrsFricVelcSurfAgg, surfFracNoEmit) |
| Calculate threshold soil surface friction velocity as a function of ag size dist., aerodynamic roughness, crust, rock, & flat biomass cover,and soil surface wetness. | |
| real function | calcFricVelc (anemHght, aeroRougStatHght, avgWindSped, aeroRougCnpy, bioDragCoef) |
| Calculates friction velocity, given station anemometer height, surface roughness, wind speed, and aerodynamic roughness. | |
| subroutine | updGridPts (threFricVelUpdFlg, avgWindSped, wind_dir, ntstep, intstep, fricThrsVelcRatio) |
| to update aeroRougSurf at each grid point; To update soil friction velocity on each grid point and modify it for barriers and hills; To initialize en thresh. and cp thresh. fr. velocites on grid; To calculate max ratios of friction velocity to threshold friction velocity | |
| subroutine | calcSoilLossGain (time, flg) |
| To calc loss/dep of saltation/creep, susp. and PM-10 at cells To call sumOutput to calc. qo, qsso, q10o for each cell To calc. deposition in the boundary cells of sim. region To update the threshold friction velocity as the loose material depletes upwind and increases downwind. | |
| subroutine | setupGrid |
| to calculate grid size and spacing for EROSION. grid size assumes outer points are outside field boundary a distance ix/2 to calculate number of grid points for EROSION. A max 'interior' square grid of 29X29 is assigned-no barriers A max 'interior' rectangular grid of 59X59 is assigned barriers to assign subregion index no. to each grid point. | |
| subroutine | updSurf (fricVelc, thrsFricVelc, fricVelcMod, soilFracDiamLt84ic, asoilLayrRock, soilMassAvalAggSurfmx, soilMassAvalAggSurf, soilFracDiamLt84mn, soilFracDiamLt84) |
| calc. minimum erodible fraction (soilFracDiamLt84mn) needed to stop emission on aggregated portion of surface at current threshold friction velocity that occurs when soilMassAvalAggSurf is removed. calc. potential mobile soil reservoir (soilMassAvalAggSurfmx) for a smooth surface 0.8 m/s friction velocity and soilFracDiamLt84ic that armors clod surface. calc. available mobile reservoir (soilMassAvalAggSurf) of soilFracDiamLt84 for current surface based on wus-thrsFricVelc ratio for | |
| subroutine | calcEmit (ounit, ws, hhr) |
| To calc the emissions for each time step of the input wind speed The emissions for EPA are the suspension component with units kg m-2 s-1. To write out a file in the format: 12 blank col, yr, mo, day, hr, soucename, emissionrate. | |
| subroutine | calcWind (wind_dir, prev_dir) |
| Calc. wind angle on the sim. region Calc. sweep sequence for update of grid cells calc. ridge spacing parallel the wind. | |
| subroutine | calcPM10 (soilLayrAggSta, soilCrstStab, soilLayrClay, soilLayrSand, awzypt, coefAbraAgg, coefAbraCrst, soilPM10FracAbraSusp, soilPM10FracEmitSusp, soilPM10FracSaltBrkSusp) |
| Calculates abrasion coefficients and PM10 fractions in sources of suspended soil. | |
| subroutine | initSubr |
| Input subregion values of variables from other submodels to the grid points of the erosion submodel which erosion changes Initialize output grid array Calc. soil fraction of 4 dia. from asd, & rr shelter angles. | |
| subroutine | calcOpenFldFricVelc |
| Calculates the fraction of open field friction velocity in shelter for 8 cardinal wind directions 0, 45, ...315 at all interior nodes. Assigns the 8 fractons calculated at each node to as 3-d array (W0br(i,j,k)) for all nodes inside sim. region. | |
| real function | calcDist (x1, y1, x2, y2) |
| Calculates distance between two points given xy coordinates. | |
| real function | calcAngToBarr (x1, y1, x2, y2) |
| Calculates angle of vector from starting point to destination point. | |
| real function | adjFricVelcForBarr (distToBarr, barrHght, effPor, barrWid) |
| Calculates reduction in friction velocity due to upwind barriers (ranges: porosity 0 to 0.9, distance: -5*barrHght to 50*barrHght) | |
| subroutine | zeroGrid |
| Sets the grid output arrays to zero. | |
| subroutine | sumOutput (flg, fricVelc, thrsFricVelc, thrsFricVelcTrap, soilFracDiaLt10, soilFracDiamLt84, soilFracDiaLt200, soilCrstThck, soilCrstFrac, soilLoosCovFrac, soilLoosMass, ridgHght, orgRidgSpac, ridgSpacParaWind, randRoug, soilLayrClay, soilLayrSand, sfvfs, soilLayrRock, brsai, bioMassHght, bioFracFlatCovr, time, coefAbraAgg, coefAbraCrst, soilPM10FracAbraSusp, soilPM10FracEmitSusp, soilPM10FracSaltBrkSusp, lx, qi, qssi, q10i, i, j, imax, jmax, soilMassAvalAggSurf, soilMassAvalDelt, soilFracDiamLt84mn, soilFracDiamLt84ic, soilFracDiaLt10ic, asoilLayrRock, soilMassAvalAggSurfmx, qo, qsso, q10o) |
| Calculates the saltation/creep, suspension, and PM-10 discharge from each control volume. | |
| real function | calcSoilMassFrac (aggDistGeoMean, aggDistGeoStd, aggDistLoLmt, aggDistUpLmt, aggDistInpDiam) |
| Calculates soil mass fraction (soilMassFrac) < diameter (aggDistInpDiam) given modified lognormal distribution parameters. | |
| subroutine | calcZeroPlanDisp (ridgHght, cropRowSpac, cropFurrFlg, aeroAnemFlg, resiLAI, resiSAI, resiAvgHght, cropLAI, cropSAI, cropHght, dispHghtWeatStat, zeroPlanDisp) |
| Set field zero plane displacement equal to Anemometer zero plane displacement when anemometer at the field site. | |
| subroutine | anemometer_init |
| Provides initial default values to weather station variables. | |
| real function adjFricVelcForBarr | ( | real | distToBarr, |
| real | barrHght, | ||
| real | effPor, | ||
| real | barrWid | ||
| ) |
Calculates reduction in friction velocity due to upwind barriers (ranges: porosity 0 to 0.9, distance: -5*barrHght to 50*barrHght)
| distToBarr | - distance from barrier |
| barrHght | - barrier height |
| effPor | - effective porosity |
| barrWid | - barrier width |
| subroutine anemometer_init | ( | ) |
Provides initial default values to weather station variables.
| subroutine calcAeroRoug | ( | real | ridgSpacParaWind, |
| real | ridgHght, | ||
| real | randRoug, | ||
| integer | aeroAnemFlg, | ||
| real | resiLAI, | ||
| real | resiSAI, | ||
| real | resiAvgHght, | ||
| real | cropLAI, | ||
| real | cropSAI, | ||
| real | cropHght, | ||
| real | cropRowSpac, | ||
| integer | cropFurrFlg, | ||
| real | aeroRougRidg, | ||
| real | aeroRougRand, | ||
| real | aeroRougSurf, | ||
| real | aeroRougCnpy, | ||
| real | aeroRougStatHght, | ||
| real | bioDragCoef | ||
| ) |
Calculates aerodynamic roughness parameter, which is aeroRougSurf if there is no standing biomass or aeroRougCnpy if standing biomass is present.
| ridgSpacParaWind | - row/dike spacing parallel the wind (mm) |
| ridgHght | - ridge height (mm) |
| aeroAnemFlg | - flag=0 - anemometer at station flag=1 - anemometer at field |
| randRoug | - random roughness (mm) |
| resiLAI | - residue leaf area index (total)(m2/m2) |
| resiSAI | - residue stem area index (total)(m2/m2) |
| resiAvgHght | - composite average residue height (m) |
| cropLAI | - crop leaf area index (m2/m2) |
| cropSAI | - crop stem area index (m2/m2) |
| cropHght | - crop height (m) |
| cropRowSpac | - crop row spacing (m) |
| cropFurrFlg | - flag=0 - crop planted in furrow bottom flag=1 - crop planted on ridge top |
| aeroRougSurf | - aerodynamic roughness of surface below canopy (mm) |
| aeroRougRidg | - aerodynamic roughness of ridge |
| aeroRougRand | - aerodynamic roughness of random roughness |
| aeroRougCnpy | - aerodynamic roughness length of canopy (mm) |
| aeroRougStatHght | - aerodynamic roughness at anemom. site (mm) |
| bioDragCoef | - biomass drag coefficient |
uses biodrag()
| real function calcAngToBarr | ( | real | x1, |
| real | y1, | ||
| real | x2, | ||
| real | y2 | ||
| ) |
Calculates angle of vector from starting point to destination point.
| x1 | - x cooridinate of destination point (barrier) |
| y1 | - y cooridinate of destination point (barrier) |
| x2 | - x cooridinate of starting point (target cell) |
| y2 | - y cooridinate of starting point (target cell) |
| real function calcDist | ( | real | x1, |
| real | y1, | ||
| real | x2, | ||
| real | y2 | ||
| ) |
Calculates distance between two points given xy coordinates.
| x1 | - x coordinate of first point |
| y1 | - y coordinate of first point |
| x2 | - x coordinate of second point |
| y2 | - y coordinate of second point |
| subroutine calcEmit | ( | integer | ounit, |
| real | ws, | ||
| real | hhr | ||
| ) |
To calc the emissions for each time step of the input wind speed The emissions for EPA are the suspension component with units kg m-2 s-1. To write out a file in the format: 12 blank col, yr, mo, day, hr, soucename, emissionrate.
| real function calcFricVelc | ( | real | anemHght, |
| real | aeroRougStatHght, | ||
| real | avgWindSped, | ||
| real | aeroRougCnpy, | ||
| real | bioDragCoef | ||
| ) |
Calculates friction velocity, given station anemometer height, surface roughness, wind speed, and aerodynamic roughness.
| anemHght | - parameter, anemometer height of input wind speed (m). |
| aeroRougStatHght | - parameter, surface aerodynamic roughness at input wind speed location (mm). |
| avgWindSped | - input wind speed driving EROSION submodel (m/s). |
| aeroRougCnpy | - subregion aerodynamic roughness (mm). |
| bioDragCoef | - biomass drag coefficient  |
Friction velocity at the weather stations (E-70) 
Translate friction velocity at station to friction velocity at the field (E-71)
If biomass is present translate friction velocity at station to friction velocity at the field (E-72) 
Adjust for above to below canopy when standing biomass is present (E-73) 
| subroutine calcOpenFldFricVelc | ( | ) |
Calculates the fraction of open field friction velocity in shelter for 8 cardinal wind directions 0, 45, ...315 at all interior nodes. Assigns the 8 fractons calculated at each node to as 3-d array (W0br(i,j,k)) for all nodes inside sim. region.
| subroutine calcPM10 | ( | real | soilLayrAggSta, |
| real | soilCrstStab, | ||
| real | soilLayrClay, | ||
| real | soilLayrSand, | ||
| real | awzypt, | ||
| real | coefAbraAgg, | ||
| real | coefAbraCrst, | ||
| real | soilPM10FracAbraSusp, | ||
| real | soilPM10FracEmitSusp, | ||
| real | soilPM10FracSaltBrkSusp | ||
| ) |
Calculates abrasion coefficients and PM10 fractions in sources of suspended soil.
| soilLayrAggSta | = aggreg. stability [Ln(J/Kg)] |
| soilCrstStab | = crust stabitlity [Ln(J/Kg)] |
| soilLayrClay | = soil surface fraction clay |
| soilLayrSand | = soil surface faction sand |
| awzypt | = annual average precipitation (mm) |
| coefAbraAgg | = coefficent of abrasion of aggregates (1/m) |
| coefAbraCrst | = coefficient of abrasion of crust (1/m) |
| soilPM10FracAbraSusp | = fraction pm10 in abraded supension size soil |
| soilPM10FracEmitSusp | = fraction pm10 in emitted suspension size soil |
| soilPM10FracSaltBrkSusp | = fraction pm10 in breakage from saltion size soil |
| subroutine calcSoilLossGain | ( | real | time, |
| integer | flg | ||
| ) |
To calc loss/dep of saltation/creep, susp. and PM-10 at cells To call sumOutput to calc. qo, qsso, q10o for each cell To calc. deposition in the boundary cells of sim. region To update the threshold friction velocity as the loose material depletes upwind and increases downwind.
| time | - time interval (seconds) |
| flg | - surface update flag (1=on, 0=off) |
| real function calcSoilMassFrac | ( | real | aggDistGeoMean, |
| real | aggDistGeoStd, | ||
| real | aggDistLoLmt, | ||
| real | aggDistUpLmt, | ||
| real | aggDistInpDiam | ||
| ) |
Calculates soil mass fraction (soilMassFrac) < diameter (aggDistInpDiam) given modified lognormal distribution parameters.
| aggDistGeoMean | - aggregate distribution geometric mean diameter (mm). |
| aggDistGeoStd | - aggregate distribution geometric standard deviation. |
| aggDistLoLmt | - aggregate distribution lower limit (mm). |
| aggDistUpLmt | - aggregate distribution upper limit (mm). |
| aggDistInpDiam | - soil diameter in distribution (mm) |
| soilMassFrac | - soil mass fraction < aggDistInpDiam |
| subroutine calcSurfThrsFricVelc | ( | real | soilFracDiamLt84, |
| real | soilLayrAggDen, | ||
| real | soilCrstFrac, | ||
| real | soilLayrRock, | ||
| real | soilLoosCovFrac, | ||
| real | bioFracFlatCovr, | ||
| real | aeroRougSurf, | ||
| real | soilSurfWatrCont, | ||
| real | soilLayrWiltPt, | ||
| real | fricVelc, | ||
| real | soilFracDiamLt84ic, | ||
| real | asoilLayrRock, | ||
| real | soilMassAvalDelt, | ||
| real | thrsFricVelc, | ||
| real | thrsFricVelcTrap, | ||
| real | fricVelcMod, | ||
| real | soilFracDiamLt84mn, | ||
| real | soilMassAvalAggSurf, | ||
| real | soilMassAvalAggSurfmx, | ||
| real | thrsFricVelcBareSoil, | ||
| real | thrsFricVelcSurfCov, | ||
| real | thrsFricVelcSurfWet, | ||
| real | thrsFricVelcSurfAgg, | ||
| real | surfFracNoEmit | ||
| ) |
Calculate threshold soil surface friction velocity as a function of ag size dist., aerodynamic roughness, crust, rock, & flat biomass cover,and soil surface wetness.
| sfd84 | - soil mass fraction in surface layer < 0.84 mm |
| soilLayrAggDen | - aggregate density (Mg/m^3) |
| soilCrstFrac | - fraction of crust cover. |
| soilLayrRock | - updated surface vol. rock > 2.0 mm (m^3/m^3). |
| soilLoosCovFrac | - soil fraction loose material cover on crust (m^3/m^3) |
| bioFracFlatCovr | - biomass fraction of flat cover (m^2/m^2) |
| aeroRougSurf | - aerodynamic roughness length of surface below canopy(mm) |
| soilSurfWatrCont | - soil water content on mass basis (at surface) (kg/kg). |
| soilLayrWiltPt | - soil water content on mass basis, at -1.5 MPa (kg/kg) |
| fricVelc | - Soil surface friction velocity (m/s) |
| soilFracDiamLt84ic | - surface soil fraction <0.84 mm initial condition |
| asoilLayrRock | - initial surface soil volume roc fraction |
| soilMassAvalDelt | - mobile soil mass change from erosion of aggregated |
| surface | (kg/m^2) (+)= gain, (-)= loss from surface. |
| thrsFricVelc | - friction velocity theshold for en (m/s) |
| thrsFricVelcTrap | - friction velocity threshold of tp and trans. cap.(m/s) |
| soilFracDiamLt84mn | - surface soil fraction <0.84 mm where wust= fricVelc of ag.sfc. |
| soilMassAvalAggSurf | - mobile soil reservoir of initial aggregated sfc.(kg/m^2) |
| soilMassAvalAggSurfmx | - max mobile soil reservoir of aggregateed sfc.(kg/m^2) |
| thrsFricVelcBareSoil | - bare soil threshold friction velocity |
| thrsFricVelcSurfCov | - surface cover addition to bare soil threshold friction velocity |
| thrsFricVelcSurfWet | - surface wetness addition to bare soil threshold friction velocity |
| thrsFricVelcSurfAgg | - aggregate density addition to bare soil threshold friction velocity |
| surfFracNoEmit | - fraction bare surface that does not emit |
| fricVelcMod | - calculate threshold frictional velocity of bare, smooth surface |
| subroutine calcWind | ( | real | wind_dir, |
| real | prev_dir | ||
| ) |
Calc. wind angle on the sim. region Calc. sweep sequence for update of grid cells calc. ridge spacing parallel the wind.
| wind_dir | - direction of the wind in degrees from north |
| prev_dir | - previously computed direction of the wind |
| subroutine calcZeroPlanDisp | ( | real | ridgHght, |
| real | cropRowSpac, | ||
| integer | cropFurrFlg, | ||
| integer | aeroAnemFlg, | ||
| real | resiLAI, | ||
| real | resiSAI, | ||
| real | resiAvgHght, | ||
| real | cropLAI, | ||
| real | cropSAI, | ||
| real | cropHght, | ||
| real | dispHghtWeatStat, | ||
| real | zeroPlanDisp | ||
| ) |
Set field zero plane displacement equal to Anemometer zero plane displacement when anemometer at the field site.
| subroutine erosion | ( | real | minErosWindSped | ) |
Main entry point to erosion submodel.
| minErosWindSped | Minimum erosive wind speed (m/s) to evaluate for erosion loss |
| subroutine initSubr | ( | ) |
Input subregion values of variables from other submodels to the grid points of the erosion submodel which erosion changes Initialize output grid array Calc. soil fraction of 4 dia. from asd, & rr shelter angles.
| subroutine setupGrid | ( | ) |
to calculate grid size and spacing for EROSION. grid size assumes outer points are outside field boundary a distance ix/2 to calculate number of grid points for EROSION. A max 'interior' square grid of 29X29 is assigned-no barriers A max 'interior' rectangular grid of 59X59 is assigned barriers to assign subregion index no. to each grid point.
| subroutine sumOutput | ( | integer | flg, |
| real | fricVelc, | ||
| real | thrsFricVelc, | ||
| real | thrsFricVelcTrap, | ||
| real | soilFracDiaLt10, | ||
| real | soilFracDiamLt84, | ||
| real | soilFracDiaLt200, | ||
| real | soilCrstThck, | ||
| real | soilCrstFrac, | ||
| real | soilLoosCovFrac, | ||
| real | soilLoosMass, | ||
| real | ridgHght, | ||
| real | orgRidgSpac, | ||
| real | ridgSpacParaWind, | ||
| real | randRoug, | ||
| real | soilLayrClay, | ||
| real | soilLayrSand, | ||
| real | sfvfs, | ||
| real | soilLayrRock, | ||
| real | brsai, | ||
| real | bioMassHght, | ||
| real | bioFracFlatCovr, | ||
| real | time, | ||
| real | coefAbraAgg, | ||
| real | coefAbraCrst, | ||
| real | soilPM10FracAbraSusp, | ||
| real | soilPM10FracEmitSusp, | ||
| real | soilPM10FracSaltBrkSusp, | ||
| real | lx, | ||
| real | qi, | ||
| real | qssi, | ||
| real | q10i, | ||
| integer | i, | ||
| integer | j, | ||
| integer | imax, | ||
| integer | jmax, | ||
| real | soilMassAvalAggSurf, | ||
| real | soilMassAvalDelt, | ||
| real | soilFracDiamLt84mn, | ||
| real | soilFracDiamLt84ic, | ||
| real | soilFracDiaLt10ic, | ||
| real | asoilLayrRock, | ||
| real | soilMassAvalAggSurfmx, | ||
| real | qo, | ||
| real | qsso, | ||
| real | q10o | ||
| ) |
Calculates the saltation/creep, suspension, and PM-10 discharge from each control volume.
| fricVelc | = friction velocity (m/s) |
| thrsFricVelc | = friction velocity threhold at emission (m/s) |
| thrsFricVelcTrap | = friction velocity threshold at transport cap.(m/s) |
| soilFracDiaLt10 | = soil fractions less than 0.1 |
| soilFracDiamLt84 | = soil fractions less than 0.84 mm |
| soilFracDiaLt200 | = soil fractions less than 2.0 mm |
| soilFracDiamLt84ic | = soil surface fraction less than 0.84 initially |
| soilFracDiaLt10ic | = soil surface fraction less than 0.10 initially |
| asoilLayrRock | = soil surface volume rock at start of event |
| soilCrstFrac | = soil fraction area crusted |
| soilLoosCovFrac | = soil fraction area of loose soil (only on crust) |
| sarrc | = soil angle r. roughness weibull c parm |
| ridgHght | = soil ridge height (mm) |
| ridgSpacParaWind | = soil ridge spacing parallel wind direction (mm) |
| sxrg | = soil ridge spacing (mm) |
| soilLayrAggSta | = soil aggregate stability (ln(J/kg)) |
| soilCrstStab | = soil crust stability (ln(j/kg)) |
| soilLayrClay | = soil fraction clay by mass |
| soilLayrSand | = soil fraction sand by mass |
| sfvfs | = soil fraction very fine sand by mass (0.05-0.1 mm) |
| soilLayrRock | = soil fraction rock >2.0 mm by volume |
| brsai | = biomass stem area index |
| bioMassHght | = biomass height (m) |
| bioFracFlatCovr | = biomass fraction flat cover |
| coefAbraAgg | = coeficient of aggregate abrasion (1/m) |
| coefAbraCrst | = coefficient of crust abrasion (1/m) |
| soilPM10FracAbraSusp | = soil fraction of pm10 in abraded suspension size |
| soilPM10FracEmitSusp | = soil fraction of pm10 in emitted suspension size |
| soilPM10FracSaltBrkSusp | = soil fraction of pm10 in saltion/creep breakage |
| lx | = node spacing in x-direction (m) |
| qi,qssi,q10i | = input to C.V. of saltion, creep, pm-10 (kg/m*2) |
| qo,qsso,q10o | = output from C.V. of saltion, creep, pm-10 (kg/m*s) |
| szc | = tmp variable for prior crust thickness (mm) |
| sz | = tmp variable for roughness (mm) |
| crlos | = factor that decreases loose cover with roughness |
| soilMassAvalDelt | = change in loose mass on aggregated sfc. (kg/m^2) |
| dmcld | = change in clod mass on aggregated sfc. (kg/m^2) |
| szv | = change in height based on volume change (mm) |
| sacd | = specific area of clods per unit volume (mm^2/mm^3) |
| flg | = surface update flag, 1=update surface, 0=do not update surface |
| soilCrstThck | = consolidated crust thickness |
| soilLoosMass | = Amount of loose erodible material on crusted fraction of soil surface (kg/m^2) |
| orgRidgSpac | = Ridge spacing at beginning of day |
| randRoug | = soil random roughness (mm) |
| time | = time interval (seconds) |
| i,j | = grid cell indices |
| imax | = no. grid intervals in x-direction |
| jmax | = no. grid intervals in y-direction. |
| soilMassAvalAggSurf | = mobile soil reservoir of initial aggregated sfc.(kg/m^2) |
| soilFracDiamLt84mn | = soil fraction <0.84 mm where wust= fricVelc of ag.sfc. |
| soilMassAvalAggSurfmx | = max mobile soil reservoir of aggregateed sfc.(kg/m^2) |
| subroutine updGridPts | ( | integer | threFricVelUpdFlg, |
| real | avgWindSped, | ||
| real | wind_dir, | ||
| integer | ntstep, | ||
| integer | intstep, | ||
| real | fricThrsVelcRatio | ||
| ) |
to update aeroRougSurf at each grid point; To update soil friction velocity on each grid point and modify it for barriers and hills; To initialize en thresh. and cp thresh. fr. velocites on grid; To calculate max ratios of friction velocity to threshold friction velocity
| intstep | - current index of ntstep thru time |
| ntstep | - max. no. of time steps in day |
| fricThrsVelcRatio | - max ratio of friction velocity to thresh. friction vel. |
| wind_dir | - direction of wind (degrees from north) |
| threFricVelUpdFlg | - dummy |
| avgWindSped | - dummy |
| subroutine updSurf | ( | real | fricVelc, |
| real | thrsFricVelc, | ||
| real | fricVelcMod, | ||
| real | soilFracDiamLt84ic, | ||
| real | asoilLayrRock, | ||
| real | soilMassAvalAggSurfmx, | ||
| real | soilMassAvalAggSurf, | ||
| real | soilFracDiamLt84mn, | ||
| real | soilFracDiamLt84 | ||
| ) |
calc. minimum erodible fraction (soilFracDiamLt84mn) needed to stop emission on aggregated portion of surface at current threshold friction velocity that occurs when soilMassAvalAggSurf is removed. calc. potential mobile soil reservoir (soilMassAvalAggSurfmx) for a smooth surface 0.8 m/s friction velocity and soilFracDiamLt84ic that armors clod surface. calc. available mobile reservoir (soilMassAvalAggSurf) of soilFracDiamLt84 for current surface based on wus-thrsFricVelc ratio for
| fricVelc | - friction velocity (m/s) |
| thrsFricVelc | - friction velocity theshold for en (m/s) |
| fricVelcMod | - threhold friction velocity = fricVelc minus flat biomass and wetness effects (m/s) |
| soilFracDiamLt84ic | - surface soil fraction <0.84 mm initial condition |
| asoilLayrRock | - surface soil volume rock (m^3/m^3) |
| soilMassAvalAggSurfmx | - max mobile soil reservoir of aggregateed sfc.(kg/m^2) |
| soilMassAvalAggSurf | - potential mobile soil reservoir of aggregated sfc.(kg/m^2) |
| soilFracDiamLt84mn | - surface soil fraction <0.84 mm where thrsFricVelc= fricVelc of ag.sfc. |
| soilFracDiamLt84 | - soil mass fraction in surface layer < 0.84 mm |
| subroutine zeroGrid | ( | ) |
Sets the grid output arrays to zero.
1.8.0
