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,