changed NOAA QC flags into numerical values (0/1)

da/ct/obs.py

@@ -91,6 +91,7 @@ class CtObservations(Observation):
         flags       = ncf.variables['NOAA_QC_flags'].get_value().take(subselect,axis=0)
         flags       = [s.tostring().lower() for s in flags]
         flags       = map(strip,flags)
+        flags       = [int(f == '...') for f in flags]
         dummy       = ncf.close()
 
         msg         = "Successfully read data from obs file (%s)"%self.ObsFilename ; logging.debug(msg)
@@ -225,7 +226,7 @@ class CtObservations(Observation):
                 obs.may_reject      = SiteInfo[obs.code]['may_reject']
             else:
                 msg= "Observation NOT found (%s), please check sites.rc file  (%s)  !!!" % (obs.code, self.SitesFile,) ; logging.warning(msg)
-                obs.flag            = 1
+                obs.flag            = 99
 
             # Add SiteInfo dictionary to the Observation object for future use