Index: doc/_build/html/plots/calibplot.py
===================================================================
diff -u -r504c0ea7310f40dd1dc984a080e2e6fa68682491 -r2f17058524a2631c08051465fe6f54f0a8d098c4
--- doc/_build/html/plots/calibplot.py (.../calibplot.py) (revision 504c0ea7310f40dd1dc984a080e2e6fa68682491)
+++ doc/_build/html/plots/calibplot.py (.../calibplot.py) (revision 2f17058524a2631c08051465fe6f54f0a8d098c4)
@@ -1,47 +1,47 @@
-# -*- coding: utf-8 -*-
-"""
-Created on Mon Jun 18 14:42:47 2012
-
-@author: schelle
-"""
-
-
-import numpy
-import matplotlib as mpl
-#mpl.use('Agg')
-import matplotlib.pyplot as plt
-import datetime
-import pandas
-import pcrut
-from stats import *
-
-Case = "rhine_sbm"
-#Case = "rhineNew"
-warmup = 500 # timesteps to skip in warmup phase for calculations
-cooldown = 50
-
-location = 7
-obs,head=pcrut.readtss(Case + "/calib_new.tss")
-obs = obs[1:3601,:]
-
-
-for location in range(0,14):
- sim,hd=pcrut.readtss(Case + "/newcalib/run.tss")
- #shift one day!!!!
- pers = numpy.size(obs,axis=0)
-
- a = get_nash_sutcliffe(obs[warmup:len(obs)-cooldown,location],sim[warmup:len(sim)-cooldown,location],NoData=numpy.nan)
-
- #trange = pandas.DatetimeIndex(datetime.datetime(1985,1,1),periods=pers,offset=pandas.DateOffset())
- ts = pandas.Series(obs[:,location],index=pandas.date_range('1/1/1985',periods=pers))
- tssim = pandas.Series(sim[:,location],index=pandas.date_range('1/1/1985',periods=pers))
-
- plt.figure(location)
- plt.autoscale(enable=True)
- ts.plot(label='Observed',color='blue')
- plt.autoscale(enable=False)
- tssim.plot(label='Simulated',color='black')
- plt.title(head[location] + ": NS = " + str(a[0]))
- plt.legend()
-
-
+# -*- coding: utf-8 -*-
+"""
+Created on Mon Jun 18 14:42:47 2012
+
+@author: schelle
+"""
+
+
+import numpy
+import matplotlib as mpl
+#mpl.use('Agg')
+import matplotlib.pyplot as plt
+import datetime
+import pandas
+import pcrut
+from stats import *
+
+Case = "rhine_sbm"
+#Case = "rhineNew"
+warmup = 500 # timesteps to skip in warmup phase for calculations
+cooldown = 50
+
+location = 7
+obs,head=pcrut.readtss(Case + "/calib_new.tss")
+obs = obs[1:3601,:]
+
+
+for location in range(0,14):
+ sim,hd=pcrut.readtss(Case + "/newcalib/run.tss")
+ #shift one day!!!!
+ pers = numpy.size(obs,axis=0)
+
+ a = get_nash_sutcliffe(obs[warmup:len(obs)-cooldown,location],sim[warmup:len(sim)-cooldown,location],NoData=numpy.nan)
+
+ #trange = pandas.DatetimeIndex(datetime.datetime(1985,1,1),periods=pers,offset=pandas.DateOffset())
+ ts = pandas.Series(obs[:,location],index=pandas.date_range('1/1/1985',periods=pers))
+ tssim = pandas.Series(sim[:,location],index=pandas.date_range('1/1/1985',periods=pers))
+
+ plt.figure(location)
+ plt.autoscale(enable=True)
+ ts.plot(label='Observed',color='blue')
+ plt.autoscale(enable=False)
+ tssim.plot(label='Simulated',color='black')
+ plt.title(head[location] + ": NS = " + str(a[0]))
+ plt.legend()
+
+