@@ -40,7 +40,7 @@ def hsu(rainfall, cleaning_threshold, tilt, pm2_5, pm10,
4040 Concentration of airborne particulate matter (PM) with
4141 aerodynamicdiameter less than 10 microns. [g/m^3]
4242
43- depo_veloc : dict, default {'2_5': 0.4 , '10': 0.09 }
43+ depo_veloc : dict, default {'2_5': 0.0009 , '10': 0.004 }
4444 Deposition or settling velocity of particulates. [m/s]
4545
4646 rain_accum_period : Timedelta, default 1 hour
@@ -69,15 +69,15 @@ def hsu(rainfall, cleaning_threshold, tilt, pm2_5, pm10,
6969
7070 # never use mutable input arguments
7171 if depo_veloc is None :
72- depo_veloc = {'2_5' : 0.004 , '10' : 0.0009 }
72+ depo_veloc = {'2_5' : 0.0009 , '10' : 0.004 }
7373
7474 # accumulate rainfall into periods for comparison with threshold
7575 accum_rain = rainfall .rolling (rain_accum_period , closed = 'right' ).sum ()
7676 # cleaning is True for intervals with rainfall greater than threshold
7777 cleaning_times = accum_rain .index [accum_rain >= cleaning_threshold ]
7878
7979 horiz_mass_rate = pm2_5 * depo_veloc ['2_5' ]\
80- + np .maximum (pm10 - pm2_5 , 0. ) * depo_veloc ['10' ]
80+ + np .maximum (pm10 - pm2_5 , 0. ) * depo_veloc ['10' ] * 3600
8181 tilted_mass_rate = horiz_mass_rate * cosd (tilt ) # assuming no rain
8282
8383 # tms -> tilt_mass_rate
0 commit comments