print "created outfile %s\n" % oftx90p_ref
of = t_stats.TP(fn_hist_tn, 90, 'tasmin', oftx90p_ref, StComHis, EnComHis)
of = t_stats.TP(fn_rcp_tn, 90, 'tasmin', oftx90p_ref, StYrsFut, EnYrsFut)
print "created outfile %s\n" % of
# Generalized version: Pct of time T doesn't exceeds ref pd Nth percentile
if 'tn10p' in var_stat:
oftx90p_ref = t_stats.Tref(fn_hist_tn, 10, 'tasmin', StRefHis, EnRefHis)
print "created outfile %s\n" % oftx90p_ref
of = t_stats.TP(fn_hist_tn, 10, 'tasmin', oftx90p_ref, StComHis, EnComHis)
of = t_stats.TP(fn_rcp_tn, 10, 'tasmin', oftx90p_ref, StYrsFut, EnYrsFut)
print "created outfile %s\n" % of
# Frost days
if 'fd' in var_stat:
of = t_stats.FD(fn_hist_tn, StComHis, EnComHis)
print "created outfile %s\n" % of
of = t_stats.FD(fn_rcp_tn, StYrsFut, EnYrsFut)
print "created outfile %s\n" % of
# Thermal growing season length
if 'gsl' in var_stat:
of = t_stats.GSL(fn_hist_tn, fn_hist_tx, lsm, StComHis, EnComHis)
print "created outfile %s\n" % of
of = t_stats.GSL(fn_rcp_tn, fn_rcp_tx, lsm, StYrsFut, EnYrsFut)
print "created outfile %s\n" % of
# Heat wave duration index wrt mean of reference_period
if 'hwdi' in var_stat:
oftxnorm_ref = t_stats.TXnorm(fn_hist_tx, StRefHis, EnRefHis)
print "created outfile %s\n" % oftxnorm_ref
for j in range(len(StYrs)):
of = Tstats.TP(fn21tn,90,'tasmin',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 'tn10p' in vars:
oftx90p_ref = Tstats.Tref(fn20tn,10,'tasmin',StRefh,EnRefh)
print "created outfile %s\n" % (oftx90p_ref)
of = Tstats.TP(fn20tn,10,'tasmin',oftx90p_ref,StRef,EnRef)
for j in range(len(StYrs)):
of = Tstats.TP(fn21tn,10,'tasmin',oftx90p_ref,StYrs[j],EnYrs[j])
print "created outfile %s\n" % (of)
#Frost days
if 'fd' in vars:
of = Tstats.FD(fn20tn,StRef,EnRef)
print "created outfile %s\n" % (of)
for j in range(len(StYrs)):
of = Tstats.FD(fn21tn,StYrs[j],EnYrs[j])
print "created outfile %s\n" % (of)
#Thermal growing season length
if 'gsl' in vars:
of = Tstats.GSL(fn20tn,fn20tx,lsm,StRef,EnRef)
print "created outfile %s\n" % (of)
for j in range(len(StYrs)):
of = Tstats.GSL(fn21tn,fn21tx,lsm,StYrs[j],EnYrs[j])
print "created outfile %s\n" % (of)
#Heat wave duration index wrt mean of reference_period
if 'hwdi' in vars: