Changeset 713

Timestamp:

05/04/16 13:50:22 (8 years ago)

Author:

llmlod

Message:

Check that empty orbit are not save in netcdf

File:

• trunk/src/file_Tbcorr_daily.py (modified) (2 diffs)

trunk/src/file_Tbcorr_daily.py

@@ -13 +13 @@
 import os
 from pylab import date2num
-from numpy import arange, size, hstack, dstack, dtype, argsort, mean, nanmean, isnan
+from numpy import arange, size, hstack, dstack, dtype, argsort, mean, nanmean, isnan, floor
 from dateutil.rrule import rrule, DAILY
 from useCorrTb_daily import applycorr2orb, SatName, areaname
@@ -215 +215 @@
         for cc, ich in enumerate(ch):
             # Index to remove orbites completly 
-            iorb = (mean(tbcorr[cc,:,:].squeeze(),axis=0) != fvalue)
+            iorb = (floor(mean(tbcorr[cc,:,:].squeeze(),axis=0)) > fvalue)
             nsat = satvec[iorb].size
             # Output file name