jkrasting / momlevel

Sea Level Routines for the Modular Ocean Model (MOM)
Other
6 stars 6 forks source link

Return potential vorticity in different units #21

Closed jkrasting closed 1 year ago

jkrasting commented 2 years ago

An option to return potential vorticity in units of $10^{14} cm^{-1} sec^{-1}$ would be useful. This conversion is done by:

pvort = np.abs(((pvort/100)*1.e14))

It would allow one to generate plots similar to this one from Talley and McCarthy:

image

http://sam.ucsd.edu/papers/WOCE_potvor_abs.html

The colormap and levels for this figure are:

colors = ["#ffffff",
          "#66ff00",
          "#66ff00",
          "#66ff00",
          "#66ff00",
          "#66ff00",
          "#99ff00",
          "#99ff00",
          "#99ff00",
          "#99ff00",
          "#99ff00",
          "#99ff00",
          "#99ff00",
          "#ccff00",
          "#ccff00",
          "#ccff00",
          "#ccff00",
          "#ffccff",
          "#ffccff",
          "#ff99ff",
          "#ff99ff",
          "#ff4cff",
          "#ff4cff",
          "#ff4cff",
          "#ff4cff",
          "#ff4cff",
          "#ff00ff",
          "#ff00ff",
          "#ff00ff",
          "#ff00ff",
         ]

levels = [0,0.0001,1,2,3,4,5,6,7,8,9,10,15,20,30,40,50,60,70,80,90,100,120,140,160,180,200,300,400,2000]
jkrasting commented 1 year ago

Addressed by e57d263. Closing