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
of = t_stats.HWDI(fn_hist_tx, oftxnorm_ref, StComHis, EnComHis)
print "created outfile %s\n" % of
of = t_stats.HWDI(fn_rcp_tx, oftxnorm_ref, StYrsFut, EnYrsFut)
print "created outfile %s\n" % of
###########################################################################
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:
oftxnorm_ref = Tstats.TXnorm(fn20tx,StRefh,EnRefh)
print "created outfile %s\n" % (oftxnorm_ref)
of = Tstats.HWDI(fn20tx,oftxnorm_ref,StRef,EnRef)
print "created outfile %s\n" % (of)
for j in range(len(StYrs)):
of = Tstats.HWDI(fn21tx,oftxnorm_ref,StYrs[j],EnYrs[j])
print "created outfile %s\n" % (of)