Index: wflow-py/wflow/wflow_adapt.py
===================================================================
diff -u -r80d9c22b6f90a1d35c03f205993929f444f638e8 -raf44f8f40198557b3bd5c324342c1372dfa52225
--- wflow-py/wflow/wflow_adapt.py (.../wflow_adapt.py) (revision 80d9c22b6f90a1d35c03f205993929f444f638e8)
+++ wflow-py/wflow/wflow_adapt.py (.../wflow_adapt.py) (revision af44f8f40198557b3bd5c324342c1372dfa52225)
@@ -308,24 +308,28 @@
Converts a .tss file to a PI-xml file
"""
-
+ missval = "-999.0"
+
#try:
tss,header = pcrut.readtss(tssfile)
+
#except:
# logger.error("Tss file not found or corrupt: ", tssfile)
# return
# Add dummpy first timesteps
if len(tss.shape) > 1:
dumm = tss[0,:].copy()
- dumm[:] = -999
+ dumm[:] = -999.0
tss = pcrut.numpy.vstack((dumm,tss))
else:
dumm = tss.copy()
- dumm[:] = -999
+ dumm[:] = -999.0
tss = pcrut.numpy.vstack((dumm,tss))
-
+ # replace NaN with missing values
+ tss[pcrut.numpy.isnan(tss)] = missval
+
trange = timedelta(seconds=timestep * (tss.shape[0]))
extraday = timedelta(seconds=timestep)
@@ -334,12 +338,10 @@
#Edate = Sdate + trange - extraday
Sdate = Sdate + extraday
Edate = Sdate + trange - extraday - extraday
-
-
+
Sdatestr = Sdate.strftime('%Y-%m-%d')
Stimestr = Sdate.strftime('%H:%M:%S')
-
-
+
Edatestr = Edate.strftime('%Y-%m-%d')
Etimestr = Edate.strftime('%H:%M:%S')
ofile = open(xmlfile,'w')
@@ -358,7 +360,7 @@
ofile.write("\n")
ofile.write("\n")
ofile.write("\n")
- ofile.write("-999.0\n")
+ ofile.write(""+str(missval)+"\n")
ofile.write("" + header[count-1] + "\n")
ofile.write("\n")
# add data here