of = t_stats.cd18(fn_hist_tas, StComHis, EnComHis, model)
print "created outfile %s\n" % of
if 'cd18' in var_stat and scens != 'historical':
of = t_stats.cd18(fn_rcp_tas, StYrsFut, EnYrsFut, model)
print "created outfile %s\n" % of
# 90th percentile Tmax - one value per year
if 'tx90' in var_stat:
of = t_stats.TX90(fn_hist_tx, StComHis, EnComHis)
print "created outfile %s\n" % of
of = t_stats.TX90(fn_rcp_tx, StYrsFut, EnYrsFut)
print "created outfile %s\n" % of
# Generalized version: Pct of time T doesn't exceeds ref pd Nth percentile
if 'tx90p' in var_stat:
oftx90p_ref = t_stats.Tref(fn_hist_tx, 90, 'tasmax', StRefHis, EnRefHis)
print "created outfile %s\n" % oftx90p_ref
of = t_stats.TP(fn_hist_tx, 90, 'tasmax', oftx90p_ref, StComHis, EnComHis)
of = t_stats.TP(fn_rcp_tx, 90, 'tasmax', oftx90p_ref, StYrsFut, EnYrsFut)
print "created outfile %s\n" % of
# Generalized version: Pct of time T doesn't exceeds ref pd Nth percentile
if 'tx10p' in var_stat:
oftx90p_ref = t_stats.Tref(fn_hist_tx, 10, 'tasmax', StRefHis, EnRefHis)
print "created outfile %s\n" % oftx90p_ref
of = t_stats.TP(fn_hist_tx, 10, 'tasmax', oftx90p_ref, StComHis, EnComHis)
of = t_stats.TP(fn_rcp_tx, 10, 'tasmax', oftx90p_ref, StYrsFut, EnYrsFut)
print "created outfile %s\n" % of
# Generalized version: Pct of time T doesn't exceeds ref pd Nth percentile
if 'tn90p' in var_stat:
print "created outfile %s\n" % (of)
for j in range(len(StYrs)):
of = Tstats.TNN(fn21tn,StYrs[j],EnYrs[j])
print "created outfile %s\n" % (of)
#90th percentile Tmax - one value per year
if 'tx90' in vars:
of = Tstats.TX90(fn20tx,StRef,EnRef)
print "created outfile %s\n" % (of)
for j in range(len(StYrs)):
of = Tstats.TX90(fn21tx,StYrs[j],EnYrs[j])
print "created outfile %s\n" % (of)
#Generalized version: Pct of time T doesn't exceeds ref pd Nth percentile
if 'tx90p' in vars:
oftx90p_ref = Tstats.Tref(fn20tx,90,'tasmax',StRefh,EnRefh)
print "created outfile %s\n" % (oftx90p_ref)
of = Tstats.TP(fn20tx,90,'tasmax',oftx90p_ref,StRef,EnRef)
for j in range(len(StYrs)):
of = Tstats.TP(fn21tx,90,'tasmax',oftx90p_ref,StYrs[j],EnYrs[j])
print "created outfile %s\n" % (of)
#Generalized version: Pct of time T doesn't exceeds ref pd Nth percentile
if 'tx10p' in vars:
oftx90p_ref = Tstats.Tref(fn20tx,10,'tasmax',StRefh,EnRefh)
print "created outfile %s\n" % (oftx90p_ref)
of = Tstats.TP(fn20tx,10,'tasmax',oftx90p_ref,StRef,EnRef)
for j in range(len(StYrs)):
of = Tstats.TP(fn21tx,10,'tasmax',oftx90p_ref,StYrs[j],EnYrs[j])
print "created outfile %s\n" % (of)
#Generalized version: Pct of time T doesn't exceeds ref pd Nth percentile