...
First, we define the parameter and time step for the cross section then call flexpart_filter() to extract the data. The result is a fieldset with units of "kg m**-3", which we scale explicitly to "ng m**-3" units for plotting (automatic units scaling does not work for cross sections).
| Code Block | ||
|---|---|---|
| ||
dIn="result_fwd/"
inFile=dIn & "conc_s001.grib"
#Define level, parameter and step
lev=-1 #all levels
par="mdc"
step=48
#Get fields for all levels for a given step
g=flexpart_filter(source: inFile,
param: par,
levType: "hl",
step: step)
#Scale into ng/m3 units
g=g*1E12
|
...