...
Next, we define the contouring. The "feflux" fields are automatically scaled into "ng have the units of "kg m**-2 " for contouring (see here for details) s**-1" but with the current value range it is better to use "g ng m**-2 s**-1" units would better fit for contouring. To achieve it we simply multiply the "tcmdfeflux" fieldset with 10001012:
| Code Block | ||
|---|---|---|
| ||
g=g*1E12 |
The contour definition itself goes like this:
...
| Code Block | ||||||
|---|---|---|---|---|---|---|
| ||||||
#The contour levels
cont_list=[1,10,50,100,150,200,250,500,750,1000,2000,4000]
#Define contour shading
flux_shade = mcont(
legend : "on",
contour : "off",
contour_level_selection_type : "level_list",
contour_level_list : cont_list,
contour_label : "off",
contour_shade : "on",
contour_shade_method : "area_fill",
contour_shade_max_level_colour : "red",
contour_shade_min_level_colour : "RGB(0.14,0.37,0.86)",
contour_shade_colour_direction : "clockwise",
contour_method : "linear"
) |
...