Page History

...

First, we define the parameter and time step for the cross section then call mvl_flexpart_read_hl() to extract the data. The result is a fieldset with units of "kg m**-3" that we need to explicitly convert to "ng m**-3" units for plotting since the automatic units scaling only works for map based plots.

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=mvl_flexpart_read_hl(inFile,par,lev,step,1)

#Scale into ng/m3 units
g=g*1E12

...