Commit 426c61b2 authored by Nils Brüggemann's avatar Nils Brüggemann

notebooks/examp_atm_converting_omega_to_w.ipynb: Added notebook to convert w[Pa/s] to w[m/s].

parent f7a179eb
......@@ -104,7 +104,7 @@
"cell_type": "markdown",
"metadata": {},
"source": [
"Load pressure (time dependent) and geometric height (time independent):"
"Load pressure, density and vert. vel. [Pa/s]:"
]
},
{
......@@ -152,6 +152,13 @@
"wap, it_ave = pyic.time_average(IcD_atm3d, 'wap', t1=t1, t2=t2, iz='all')"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"Convert wap from [Pa/s] to wa [m/s]:"
]
},
{
"cell_type": "code",
"execution_count": 6,
......@@ -162,6 +169,13 @@
"wa = -wap/(rho*grav)"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"Interoplate to regular grid."
]
},
{
"cell_type": "code",
"execution_count": 7,
......@@ -171,6 +185,13 @@
"lon, lat, wahi = pyic.interp_to_rectgrid(wa, fpath_ckdtree_atm, coordinates='clat clon')"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"Plot the data at model level:"
]
},
{
"cell_type": "code",
"execution_count": 8,
......@@ -217,7 +238,7 @@
"source": [
"## Interpolate to pressure level\n",
"\n",
"Define pressure levels for which to calculate the geopotential height:"
"Define pressure levels on which we want to interpolate wa, and speciy 500hPa level:"
]
},
{
......@@ -230,6 +251,13 @@
"ip500 = np.argmin((IcD_atm3d.plevc-500e2)**2)"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"Interpolate from model levels to pressure levels:"
]
},
{
"cell_type": "code",
"execution_count": 10,
......@@ -240,6 +268,13 @@
"wavi = wa[ind_lev,icall]*fac+wa[ind_lev+1,icall]*(1.-fac)"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"Interpolate horizontally to rectangular grid:"
]
},
{
"cell_type": "code",
"execution_count": 11,
......@@ -249,6 +284,13 @@
"lon, lat, wavihi = pyic.interp_to_rectgrid(wavi, fpath_ckdtree_atm, coordinates='clat clon')"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"Plot the data:"
]
},
{
"cell_type": "code",
"execution_count": 13,
......
Markdown is supported
0%
or
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment