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
###########################################################################
# Monthly total precip
if 'ptot' in var_stat and scens == 'historical':
of = p_stats.ptot(fn_hist_prmm, StComHis, EnComHis, model)
print "created outfile %s\n" % of
if 'ptot' in var_stat and scens != 'historical':
of = p_stats.ptot(fn_rcp_prmm, StYrsFut, EnYrsFut, model)
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)
#cooling degree days
if 'cd18' in vars:
of = Tstats.CD18(fn20tg,StRef,EnRef)
print "created outfile %s\n" % (of)
for j in range(len(StYrs)):
of = Tstats.CD18(fn21tg,StYrs[j],EnYrs[j])
print "created outfile %s\n" % (of)