def Outputs(inputpath, outputpath):
print 'Generating outputs'
# Load pickeled model results
p = open('_pickle/scenario', 'r')
scenario = pickle.load(p)
p.close()
p = open('_pickle/filenames', 'r')
filenames = pickle.load(p)
p.close()
p = open('_pickle/reefmap', 'r')
reefmap = pickle.load(p)
p.close()
p = open('_pickle/params', 'r')
params = pickle.load(p)
p.close()
p = open('_pickle/mcavs', 'r')
mcavs = pickle.load(p)
p.close()
p = open('_pickle/mcCIs', 'r')
mcCIs = pickle.load(p)
p.close()
p = open('_pickle/mcavsSR', 'r')
mcavsSR = pickle.load(p)
p.close()
p = open('_pickle/mcCIsSR', 'r')
mcCIsSR = pickle.load(p)
p.close()
p = open('_pickle/cellavs', 'r')
cellavs = pickle.load(p)
p.close()
p = open('_pickle/forcing_nut', 'r')
forcing_nut = pickle.load(p)
p.close()
p = open('_pickle/forcing_sed', 'r')
forcing_sed = pickle.load(p)
p.close()
p = open('_pickle/forcing_hurr', 'r')
forcing_hurr = pickle.load(p)
p.close()
p = open('_pickle/forcing_cm', 'r')
forcing_cm = pickle.load(p)
p.close()
p = open('_pickle/forcing_f', 'r')
forcing_f = pickle.load(p)
p.close()
p = open('_pickle/forcing_df', 'r')
forcing_df = pickle.load(p)
p.close()
p = open('_pickle/hlog_cat', 'r')
hlog_cat = pickle.load(p)
p.close()
p = open('_pickle/hlog_sr', 'r')
hlog_sr = pickle.load(p)
p.close()
p = open('_pickle/cmlog_sr', 'r')
cmlog_sr = pickle.load(p)
p.close()
p = open('_pickle/defectruns', 'r')
defectruns = pickle.load(p)
p.close()
p = open('_pickle/timing', 'r')
timing = pickle.load(p)
p.close()
p = open('_pickle/coralmaps', 'r')
coralmaps = pickle.load(p)
p.close()
# Pathnames for output files
pathmeans = outputpath + scenario + '_means'
pathmeansSR = outputpath + scenario + '_meansSR'
pathspatial = outputpath + scenario + '_spatial'
pathlog = outputpath + scenario + '_log.txt'
pathlogF = outputpath + scenario + '_logF.txt'
pathparamsout = outputpath + scenario + '_params.txt'
pathoptionsout = outputpath + scenario + '_options.txt'
pathcoralmap = outputpath + scenario + '_coralmap'
output = Output()
runf = params.runs - len(defectruns)
# Write output files
output.meansCSV(pathmeans, params, mcavs, runf)
output.meansSRCSV(pathmeansSR, params, reefmap, mcavsSR, runf)
output.spatialTXT(pathspatial, variables, params, reefmap, cellavs, 100, 1)
output.meansCDF(pathmeans, reefmap, params, variables, runf, mcavs)
output.meansSRCDF(pathmeansSR, reefmap, params, variables, runf, mcavsSR)
output.writeLog(pathlog, version, scenario, params, reefmap,
filenames['paramfile'], filenames['ivfile'], timing,
defectruns)
output.writeLogF(pathlogF, params, forcing_f, forcing_nut, forcing_sed,
forcing_hurr, forcing_cm, forcing_df,
filenames['paramfile'], filenames['ivfile'],
filenames['ffile'], filenames['nfile'],
filenames['sfile'], filenames['hfile'],
filenames['cmfile'], filenames['dffile'], hlog_cat,
hlog_sr, cmlog_sr, runf)
output.coralmapTXT(pathcoralmap, reefmap, coralmaps, 100)
#output.writeParams(pathparamsout,filenames['paramfile'])
#output.writeParams(pathoptionsout,inputpath + '/CORSET_options.txt')
# Create and save graphics
display = Display()
display.plotEndSR(reefmap, params, output)
display.plotDynamics(reefmap, params, output)
display.plotDynamicsSR(reefmap, params, output)
display.plotCatchSR(reefmap, params, output)
display.savePlot(outputpath, scenario)
