Пример #1
0
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,