mjucker/aostools

Helper functions for postprocessing and analysis of netCDF data. Also includes I/O routines to seamlessly work with pv_atmos.

aostools's Issues

Time mean of d(theta_bar)/dp?

Is there any reason for the following time mean of d(theta_bar)/dp?

aostools/climate.py

Line 482 in 5b88663

dthdp = np.nanmean(dthdp,axis=0)[np.newaxis,:]

I did not find any justification for the time mean in the literature [1,2].

[1] Andrews, D. G., J. R. Holton, and C. B. Leovy, 1987: Middle Atmosphere Dynamics. Academic Press.
[2] Edmon, H. J., B. J. Hoskins, and M. E. McIntyre, 1980: Eliassen-Palm Cross Sections for the Troposphere. J. Atmos. Sci., 37, 2600–2616.

eof() should accept multi-dimensional arrays

As is, eof() only accepts arrays of shape [time x space]. If we want a 2d EOF (for instance), the array would have shape [time x space1 x space2], and the function won't work. There is a simple solution with numpy.reshape, but this should be included in the function rather than up to the user before calling the function.

constant p0

Hi Martin, I have been playing around with the "ComputeMeridionalPVGrad" function and I think the call to the .constants module may be slightly off. The call says "import p0_Pa as p0" but I think it should just be "p0". Thanks!

Difference between uref and u

Hi,
Thanks for developing this awesome tool.

I would like to ask,
Please, what is the difference between uref and u when using the ComputeWaveActivityFlux function?
Since there are no tutorials on it yet, please, could you clarify this?
I have monthly u and v wind datasets. Please, see my code below.

data = xr.open_dataset(data1)
lon=data['lon']
lat=data['lat']
u = data['uwnd']
v = data['vwnd']
pres = data['level']
tref = data['air'] + 273.15 ##must be in Kelvin
hgt = data['hgt']

WAF = cli.ComputeWaveActivityFlux(phi_or_u = u,phiref_or_v = v,uref=u,vref=v,lat='infer',lon='infer',pres='infer',tref=tref,is_anomaly=False,qg=False,use_windspharm=False)

Thank you