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)
# MY_YEAR = TimeInterval('20130101','20140101',"%Y%m%d")
req_label = Timestart[0:4] #'Ave:2013'
req = requestors.Generic_req(TI)
indexes, weights = TL.select(req)
nFrames = len(indexes)
SAT_3D = np.zeros((nFrames, jpj, jpi), np.float32)
for iFrame, k in enumerate(indexes):
t = TL.Timelist[k]
inputfile = INPUTDIR + t.strftime("%Y%m%d") + "_d-OC_CNR-L4-CHL-MedOC4_SAM_7KM-MED-REP-v02.nc"
CHL = Sat.readfromfile(inputfile, "lchlm")
SAT_3D[iFrame, :, :] = CHL
Sat2d = Sat.averager(SAT_3D)
masknan = TheMask.mask_at_level(0)
Sat2d[~masknan] = np.NaN
var = "SATchl"
layer = Layer(0, 10)
fig, ax = mapplot(
{"varname": var, "clim": [0, 0.4], "layer": layer, "data": Sat2d, "date": "annual"},
fig=None,
ax=None,
mask=TheMask,
coastline_lon=clon,
coastline_lat=clat,
)
ax.set_xlim([-5, 36])
