def validate(directory, HUC8, NWISgage, start, end, temp = True,
snow = True, hydrology = True, warmup = 1, upcomids = []):
run_dates = (datetime.datetime(start, 1, 1),
datetime.datetime(end, 1, 1))
process_dates = (datetime.datetime(start + warmup, 1, 1),
datetime.datetime(end, 1, 1))
f = '{}/{}/calibrations/{}/hspfmodel'.format(directory, HUC8, NWISgage)
folder = '{}/{}/validations/{}'.format(directory, HUC8, NWISgage)
with open(f, 'rb') as m: hspfmodel = pickle.load(m)
# figure out the external targets needed
targets = ['reach_outvolume', 'groundwater', 'water_state',
'evaporation', 'runoff', 'snow_state', 'snowpack',
'snowfall']
# build the input wdmfile
hspfmodel.build_wdminfile()
# create the UCI file and the output WDM file
hspfmodel.build_uci(targets, run_dates[0], run_dates[1], temp = temp,
snow = snow, hydrology = hydrology)
hspfmodel.messagepath = None
# run it
hspfmodel.run(verbose = True)
# add the simulation parameters to the calibrator
postprocessor = Postprocessor(hspfmodel, process_dates)
# calculate the errors and make a calibration report
postprocessor.get_hspexp_parameters(verbose = False)
postprocessor.get_calibration(verbose = False, vverbose = False)
# plot everything and save to file
if not os.path.isdir('{}/{}/validations'.format(directory, HUC8)):
os.mkdir('{}/{}/validations'.format(directory, HUC8))
if not os.path.isdir(folder): os.mkdir(folder)
output = folder + '/calibration_report.csv'
postprocessor.calculate_errors(verbose = False, output = output)
output = folder + '/snow'
postprocessor.plot_snow(output = output, show = False)
output = folder + '/snowcalibration'
postprocessor.plot_snowcalibration(output = output, show = False)
output = folder + '/monthlyhydrograph'
postprocessor.plot_hydrograph(output = output, show = False)
output = folder + '/dayofyear'
postprocessor.plot_dayofyear(output = output, show = False)
output = folder + '/averagewaterbudget'
postprocessor.plot_waterbudget(output = output, show = False)
output = folder + '/waterbudgets'
postprocessor.plot_allwaterbudgets(output = output, show = False)
output = folder + '/runoff'
postprocessor.plot_runoff(output = output, show = False)
output = folder + '/storms'
postprocessor.plot_storms(output = output, show = False)
output = folder + '/calibration'
postprocessor.plot_calibration(output = output, show = False)
postprocessor.close()
def validate(directory, HUC8, NWISgage, start, end, temp=True, snow=True, hydrology=True, warmup=1, upcomids=[]):
run_dates = (datetime.datetime(start, 1, 1), datetime.datetime(end, 1, 1))
process_dates = (datetime.datetime(start + warmup, 1, 1), datetime.datetime(end, 1, 1))
f = "{}/{}/calibrations/{}/hspfmodel".format(directory, HUC8, NWISgage)
folder = "{}/{}/validations/{}".format(directory, HUC8, NWISgage)
with open(f, "rb") as m:
hspfmodel = pickle.load(m)
# figure out the external targets needed
targets = [
"reach_outvolume",
"groundwater",
"water_state",
"evaporation",
"runoff",
"snow_state",
"snowpack",
"snowfall",
]
# build the input wdmfile
hspfmodel.build_wdminfile()
# create the UCI file and the output WDM file
hspfmodel.build_uci(targets, run_dates[0], run_dates[1], temp=temp, snow=snow, hydrology=hydrology)
hspfmodel.messagepath = None
# run it
hspfmodel.run(verbose=True)
# add the simulation parameters to the calibrator
postprocessor = Postprocessor(hspfmodel, process_dates)
# calculate the errors and make a calibration report
postprocessor.get_hspexp_parameters(verbose=False)
postprocessor.get_calibration(verbose=False, vverbose=False)
# plot everything and save to file
if not os.path.isdir("{}/{}/validations".format(directory, HUC8)):
os.mkdir("{}/{}/validations".format(directory, HUC8))
if not os.path.isdir(folder):
os.mkdir(folder)
output = folder + "/calibration_report.csv"
postprocessor.calculate_errors(verbose=False, output=output)
output = folder + "/snow"
postprocessor.plot_snow(output=output, show=False)
output = folder + "/snowcalibration"
postprocessor.plot_snowcalibration(output=output, show=False)
output = folder + "/monthlyhydrograph"
postprocessor.plot_hydrograph(output=output, show=False)
output = folder + "/dayofyear"
postprocessor.plot_dayofyear(output=output, show=False)
output = folder + "/averagewaterbudget"
postprocessor.plot_waterbudget(output=output, show=False)
output = folder + "/waterbudgets"
postprocessor.plot_allwaterbudgets(output=output, show=False)
output = folder + "/runoff"
postprocessor.plot_runoff(output=output, show=False)
output = folder + "/storms"
postprocessor.plot_storms(output=output, show=False)
output = folder + "/calibration"
postprocessor.plot_calibration(output=output, show=False)
postprocessor.close()
