Code |
Description |
filter999 |
removes bad data points 999 from an array and replaces them with NaN |
filternan |
code to remove NaN from matrix M |
paddednum2str |
checks length of string is two characters long & if not
pads with a zero. Used to get days and months of 2 characters. |
scatter_plot |
function to plot scatter plots, one plot per page
and optionally print some basic statistics in the corner. |
brunt_p |
brunt_p calculates the brunt-vaisalla frequency for theta and thetae
using p as a vertical coordinate, theta, thetae in K, p (mb) |
brunt_z |
brunt_z calculates the brunt-vaisalla frequency for theta and thetae
using z as a vertical coordinate and assuming theta(sonde_no, data) |
hydro_pressure |
given z=height (m), t=temperature vectors (K), plus a ground pressure (mb) this function calculates the hydrostatic pressure (mb) |
hydro_height |
given pressure (mb), temp vectors (K), and scalar ground pressure
calculates the hydrostatic height |
rh_to_rhi |
converts rh to rh wrt ice, (%)
given rh (%), T (oC), and p (mb) |
rhi_to_rh |
converts rh wrt ice to rh, (%) given rhi (%), T (degC), and p (mb) |
rhi_to_rh_e |
converts rh wrt ice to rh (%) and e (mb) given rhi (%), t (oC), and p (mb) |
t_to_qsat |
generates saturated specific humidity (g/kg) wrt water and ice given T (degC), and p (mb) |
t_to_theta |
t_to_theta generates theta from T (degC) and p (mb) vectors |
thermo_es |
thermo_es calculates saturation vapour pressures (mb) wrt water and ice from T (degC) from the polynomial approximation of Lowe (see Sargent 1980) |
thermo_rh |
generates thermodynamic variables from T (degC) p(mb) rh(%) vectors; output [theta, thetae, q, qsat, qsati] = thermo_rh(t,p,rh) |
thermo_td |
generates thermodynamic variables from T (degC) p(mb) Td (degC) vectors; output [theta, thetae, q, qsat, qsati] = thermo_td(t,p,td) |
thermo2_td |
generates thermodynamic variables from T (degC) p(mb) Td (degC) vectors; output [theta, thetae, q, qsat, qsati] = thermo_td(t,p,td). Note: Same as thermo_td, except a fix for NaNs included. |
thermo_td_ecmwf |
generates thermodynamic variables from T (oC) p(mb) Td(oC); the modification for ec data is because their td at 2m field is
T_dewpoint if T >0 and T_frostpoint if T <0; output [theta, thetae, q, qsat, qsati], in K, K, g/kg, g/kg, g/kg |
theta_to_t |
generates T (K) from theta (K) and p(mb) vectors |
uv_to_vectormean |
takes vectors (u,v) and calculates the vector mean wind speed and wind direction |
uv_to_wswd |
takes vectors (u,v) and converts to wind speed and wind direction (deg from North) vectors |
wswd_to_uv |
takes wind speed and wind direction (deg from North) and converts to vectors (u,v) |