...
| Code Block |
|---|
t_fc96 = read("/homepath/graphicsto/cgkhome/metview/training/Analysis Viewsday_1/analysis views/t_fc96.grib")
t_fc24 = read("/homepath/graphicsto/cgkhome/metview/training/Analysis Viewsday_1/analysis views/t_fc24.grib")
diff = t_fc96 - t_fc24
pos = mcont(
legend : "on",
contour_level_selection_type : "level_list",
contour_max_level : 10,
contour_min_level : 0.1,
contour_level_list : [0.5,1,2,4,10],
contour_shade : "on",
contour_shade_method : "area_fill",
contour_shade_max_level_colour : "red",
contour_shade_min_level_colour : "orange_yellow",
contour_shade_colour_direction : "clockwise"
)
neg = mcont(
legend : "on",
contour_level_selection_type : "level_list",
contour_max_level : -0.5,
contour_min_level : -10,
contour_shade_max_level : -0.1,
contour_shade_min_level : -10,
contour_level_list : [-10,-4,-2,-1,-0.5],
contour_shade : "on",
contour_shade_method : "area_fill",
contour_shade_max_level_colour : "greenish_yellow",
contour_shade_min_level_colour : "blue",
contour_shade_colour_direction : "clockwise"
)
xs_europe = mxsectview(
line : [55,-6,43,16]
)
plot(xs_europe,diff,neg,pos) |
...