CODEWITHSUNDEEP
pythonmatplotlib
Fir Nor
Fir Nor

Reputation: 193

Cross section with orography

First of all, I've looked at these links:

Contouring non-uniform 2d data in python/matplotlib above terrain

Cross-section plotting with terrain as background

Matplotlib axis labels: how to find out where they will be located?

but unsuccessfully solve my problem.

As you can see from this image (ignore the wrong title, it's actually zonal winds):

zonal winds with hybrid level

where the black line is the orography and the shaded contour is zonal wind values, the winds value at near sea level, at for example 101.5 E, is under the terrain/orography instead of on top.

At this point, the data used to produce that plot are:

  1. zonal winds data with 814x70 shape. (814-longitude, 70-height)
  2. orography with 814x1 shape (basically a list of terrain height)

I wonder if I can manipulate something, to plot the data following the terrain. Some suggest grid interpolation or interpolation in general, but I've been scratching my head on how to get it works.

Upvotes: 4

Views: 2586

Answers (1)

Shakka
Shakka

Reputation: 51

I've had this problem too - it is quite tricky to get the terrain to follow orography. I tried manually adding in orography using the surface altitude variable I had in my model output but it was too computationally intensive (at least with the number of model levels and the resolution I have).

I don't know what format your files are in, or what kind of data you're using but have you tried looking at the module iris? It includes a method of adding orography into profiles when it loads files (HybridHeightFactory). Link here: http://scitools.org.uk/iris/

Upvotes: 0

Related Questions