if 'r5d' in var_stat and scens == 'historical':
of = p_stats.r5d(fn_hist_prmm, StComHis, EnComHis, model)
print "created outfile %s\n" % of
if 'r5d' in var_stat and scens != 'historical':
of = p_stats.r5d(fn_rcp_prmm, StYrsFut, EnYrsFut, model)
print "created outfile %s\n" % of
# simple daily precip intensity index
if 'sdii' in var_stat and scens == 'historical':
of = p_stats.sdii(fn_hist_prmm, StComHis, EnComHis, model)
print "created outfile %s\n" % of
if 'sdii' in var_stat and scens != 'historical':
of = p_stats.sdii(fn_rcp_prmm, StYrsFut, EnYrsFut, model)
print "created outfile %s\n" % of
# Pct of time precip exceeds ref pd 90th percentile (wet day values)
# calculate % of precip due to this too
if 'r90p' in var_stat:
ofpr90_ref = p_stats.R90ref(fn_hist_pr, StRefHis, EnRefHis)
print "created outfile %s\n" % ofpr90_ref
of = p_stats.R90P(fn_hist_pr, ofpr90_ref, StComHis, EnComHis)
print "created outfile %s\n" % of
of = p_stats.R90P(fn_rcp_pr, ofpr90_ref, StYrsFut, EnYrsFut)
print "created outfile %s\n" % of
of = p_stats.R90PTOT(fn_hist_pr, ofpr90_ref, StComHis, EnComHis)
print "created outfile %s\n" % of
of = p_stats.R90PTOT(fn_rcp_pr, ofpr90_ref, StYrsFut, EnYrsFut)
print "created outfile %s\n" % of
print '\n\n... done!\n'
print "created outfile %s\n" % (of)
#Max consec 5 day precip
if 'r5d' in vars:
of = Pstats.R5D(fn20pr,StRef,EnRef)
print "created outfile %s\n" % (of)
for j in range(len(StYrs)):
of = Pstats.R5D(fn21pr,StYrs[j],EnYrs[j])
print "created outfile %s\n" % (of)
#Pct of time precip exceeds ref pd 90th percentile (wet day values)
#calculate % of precip due to this too
if 'r90p' in vars:
ofpr90_ref = Pstats.R90ref(fn20pr,StRefh,EnRefh)
print "created outfile %s\n" % (ofpr90_ref)
of = Pstats.R90P(fn20pr,ofpr90_ref,StRef,EnRef)
print "created outfile %s\n" % (of)
for j in range(len(StYrs)):
of = Pstats.R90P(fn21pr,ofpr90_ref,StYrs[j],EnYrs[j])
print "created outfile %s\n" % (of)
of = Pstats.R90PTOT(fn20pr,ofpr90_ref,StRef,EnRef)
print "created outfile %s\n" % (of)
for j in range(len(StYrs)):
of = Pstats.R90PTOT(fn21pr,ofpr90_ref,StYrs[j],EnYrs[j])
print "created outfile %s\n" % (of)
#simple daily precip intensity index
if 'sdii' in vars:
of = Pstats.SDII(fn20pr,StRef,EnRef)
print "created outfile %s\n" % (of)
for j in range(len(StYrs)):