...
First, we define the parameter and time step for the cross section then call mvl_flexpart_read_hlfilter() 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.
| 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=mvl_flexpart_read_hl(inFile,par,lev,step,1filter(source: inFile, param: par, levType: "hl", step: step) #Scale into ng/m3 units g=g*1E12 |
...