TI = TimeInterval(Timestart,Time__end,"%Y%m%d")
TLCheck = TimeList.fromfilenames(TI, CHECKDIR,"*.nc",prefix='',dateformat='%Y%m%d')
#IonamesFile = '../../../postproc/IOnames_sat.cci.xml'
#IOname = IOnames.IOnames(IonamesFile)
suffix = os.path.basename(TLCheck.filelist[0])[8:]
WEEK_reqs=TLCheck.getWeeklyList(2)
jpi = Sat.masks.KD490mesh.jpi
jpj = Sat.masks.KD490mesh.jpj
for req in WEEK_reqs:
outfile = req.string + suffix
outpathfile = WEEKLYDIR + outfile
conditionToSkip = (os.path.exists(outpathfile)) and (not reset)
if conditionToSkip: continue
print outfile
ii, w = TLCheck.select(req)
nFiles = len(ii)
M = np.zeros((nFiles,jpj,jpi),np.float32)
for iFrame, j in enumerate(ii):
inputfile = TLCheck.filelist[j]
Kext = Sat.readfromfile(inputfile,'KD490')
M[iFrame,:,:] = Kext
Kext_OUT = Sat.averager(M)
Sat.dump_KD490_nativefile(outpathfile, Kext_OUT)
#CHL_IN[581:,164:] = Sat.fillValue # BLACK SEA
cloudsLandTIME = CHL_IN == Sat.fillValue
cloudlandsCLIM = DAILY_REF_MEAN == Sat.fillValue
CHL_OUT = CHL_IN.copy()
CHL_OUT[cloudsLandTIME] = Sat.fillValue
CHL_OUT[cloudlandsCLIM] = Sat.fillValue
counter_refNAN = (~cloudsLandTIME & cloudlandsCLIM).sum(axis=None)
outOfRange = np.abs(CHL_IN - DAILY_REF_MEAN) > DAILY_REF_STD *2.0
outOfRange[cloudsLandTIME | cloudlandsCLIM ] = False
counter_elim = outOfRange.sum(axis = None)
CHL_OUT[outOfRange] = Sat.fillValue
print filename
print 'Rejection: after check', counter_elim, ' values'
print 'rejected for NAN in Climatology', counter_refNAN, ' values'
Sat.dump_KD490_nativefile(outfile, CHL_OUT)
