Page History

For preparations and running the simulation needed for this tutorial click here ...

...

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" that , which we need to scale explicitly convert to "ng m**-3" units for plotting since the (automatic units scaling only works for map based plotsdoes not work for cross sections).

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

...