def series_run(wy):
wind_redist_obj = wind_redistribution(res=500)
sw_ds = short_wave_ds(res=500)
end_month_list = [6, 7, 7, 7, 7, 7, 7, 7, 7, 7, 6, 4]
end_day_list = [30, 31, 31, 31, 31, 31, 31, 31, 31, 31, 30, 30]
i = wy - 2001
start_date = date(wy, 4, 1)
end_date = date(wy, end_month_list[i], end_day_list[i])
temp_date = start_date
while temp_date <= end_date:
swe = delta_swe_daily(temp_date, sw_ds, wind_redist_obj, res=500)
temp_date += timedelta(days=1)
swe = None
gc.collect()
def test():
for hour in range(12, 19):
wind_redist_obj = wind_redistribution(res=500)
sw_ds = short_wave_ds(res=500)
test_forcing = delta_swe_hourly(2006, 5, 24, hour, sw_ds, wind_redist_obj, res=500)
test_forcing.visual_net_rad()