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]) ax.set_ylim([30, 46]) ax.set_xlabel("Lon").set_fontsize(12) ax.set_ylabel("Lat").set_fontsize(12) ax.ticklabel_format(fontsize=10) ax.text(
print idate1 idate2 = timerequestors.Daily_req(year, month, day) print idate2 # Variable name VARLIST = ['P_l'] #'N3n','O2o'] read_adjusted = [True] #,False,False] # MASK of the domain TheMask = Mask( "/pico/home/usera07ogs/a07ogs00/OPA/V2C/etc/static-data/MED1672_cut/MASK/meshmask.nc" ) nav_lev = TheMask.zlevels layer = Layer(0, 200) #layer of the Float profile???? # new depth from 0 to 200 meters with 1 meter of resolution NewPres = np.linspace(0, 200, 201) dimnewpress = len(NewPres) # data interpolated on the vertical Z f = open(idate0 + "." + VARLIST[0] + "_arg_mis.dat", "w") # OUTPUT x il 3DVAR iniz = " \n" f.writelines(iniz) # LIST of profiles for the selected variable in VARLIST # in the interval idate1.time_interval in the mediterranean area Profilelist_1 = bio_float.FloatSelector(LOVFLOATVARS[VARLIST[0]], idate1.time_interval, OGS.med) TL = TimeList.fromfilenames(idate2.time_interval, INPUTDIR,
iFrame = 0 for k in indexes: t = TL.Timelist[k] inputfile = INPUTDIR + t.strftime("%Y%m%d") + \ "_d-OC_CNR-L4-CHL-MedOC4_SAM_7KM-MED-REP-v02.nc" print inputfile CHL = Sat.readfromfile(inputfile,'lchlm') SAT_3D[iFrame,:,:] = CHL iFrame +=1 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]) ax.set_ylim([30,46]) ax.set_xlabel('Lon').set_fontsize(12) ax.set_ylabel('Lat').set_fontsize(12) ax.ticklabel_format(fontsize=10) ax.text(-4,44.5,var + ' [mg /m^3]', \ horizontalalignment='left', \ verticalalignment='center', \