Comments (5)
I think this because the IPSL-CM5A-LR model has a grid that includes the poles but has an even number of points in the latitude direction. The underlying software assumes that girds with an even number of points do not include the poles (and conversely grids with an odd number of latitude points must include the poles).
Regridding your data to a suitable grid may be an easy option since you are using iris.
from windspharm.
Thank you!
I tried regridding using: u.regrid(v, iris.analysis.Linear())
according to the example
rotated_air_temp = global_air_temp.regrid(rotated_psl, iris.analysis.Linear()) in http://scitools.org.uk/iris/docs/latest/userguide/interpolation_and_regridding.html,
but I am not sure this is correct. The code seems to take ages or to be stuck doing this...
from windspharm.
It may take a long time to regrid if you have a lot of data. You may want to look into caching the interpolator.
from windspharm.
Thank you. I tried interpolating..
delta_latitude = 180/96.0
sample_points = [('longitude', u.coord('longitude').points),('latitude', np.linspace(90 - 0.5 *
delta_latitude,-90 + 0.5 * delta_latitude,96))]
result = u.interpolate(sample_points, iris.analysis.Linear())
print(result.summary(shorten=True))
which gives me as result
eastward_wind / (m s-1) (time: 13140; latitude: 96; longitude: 96)
northward_wind / (m s-1) (time: 13140; latitude: 96; longitude: 96)
but now I get the following error message when I want to compute w = VectorWind(u,v) (for u = result):
raise ValueError('u and v cannot contain missing values')
ValueError: u and v cannot contain missing values
How should I proceed?
from windspharm.
I've taken a look at the data on JASMIN, it looks like the horizontal interpolation step is successful (and necessary), but the problem is that this data set really does contain missing values. Try plotting the first level at the first time, you'll see what I mean. Windspharm cannot handle data with missing values.
The missing values are there because the data are already interpolated onto pressure levels, and the lower pressure levels intersect the surface in regions of high topography. Some model output chooses to extrapolate below ground (e.g., ERA) but in this case that has not been done.
If I leave out the lowest few pressure levels (the ones that intersect the surface) you can successfully create the VectorWind
instance. If you really want the vorticity on the lower levels too, you'll either need to interpolate to fill in the gaps yourself or use a finite difference method to compute the vorticity.
from windspharm.
Related Issues (20)
- Conda Install Packages Not Found HOT 3
- Enhancement - including APIs for adding filters HOT 1
- Curious issue creating WindVector from 2D arrays HOT 2
- divergence is not right? HOT 4
- Calculating shear and stretch deform using windspharm HOT 4
- take divergence and vorticity as input HOT 1
- Issue attempting to import Windspharm in script HOT 2
- Importing JRA-55 files with windspharm
- Problem with computing vorticity of NWP model output
- ValueError: Cube 'U velocity' must contain a single 1D x coordinate. HOT 1
- Typo in windspharm.xarray docs
- windspharm future
- prep_data occur TypeError: unhashable type: 'list' HOT 1
- Possible issue with windspharm package's calculations of irrotational winds
- installation on Apple silicon HOT 2
- windspharm.iris.VectorWind fails with Numpy v1.24 HOT 1
- use windspharm on a Mercator grid? HOT 1
- numpy.float64 object cannot be interpreted as an index HOT 1
- Class instance fails if u and v are just lat,long arrays with no other (eg time) dimension HOT 3
- Memory conservation: legfunc kwarg HOT 2
Recommend Projects
-
React
A declarative, efficient, and flexible JavaScript library for building user interfaces.
-
Vue.js
🖖 Vue.js is a progressive, incrementally-adoptable JavaScript framework for building UI on the web.
-
Typescript
TypeScript is a superset of JavaScript that compiles to clean JavaScript output.
-
TensorFlow
An Open Source Machine Learning Framework for Everyone
-
Django
The Web framework for perfectionists with deadlines.
-
Laravel
A PHP framework for web artisans
-
D3
Bring data to life with SVG, Canvas and HTML. 📊📈🎉
-
Recommend Topics
-
javascript
JavaScript (JS) is a lightweight interpreted programming language with first-class functions.
-
web
Some thing interesting about web. New door for the world.
-
server
A server is a program made to process requests and deliver data to clients.
-
Machine learning
Machine learning is a way of modeling and interpreting data that allows a piece of software to respond intelligently.
-
Visualization
Some thing interesting about visualization, use data art
-
Game
Some thing interesting about game, make everyone happy.
Recommend Org
-
Facebook
We are working to build community through open source technology. NB: members must have two-factor auth.
-
Microsoft
Open source projects and samples from Microsoft.
-
Google
Google ❤️ Open Source for everyone.
-
Alibaba
Alibaba Open Source for everyone
-
D3
Data-Driven Documents codes.
-
Tencent
China tencent open source team.
from windspharm.