回答此问题可获得 20 贡献值,回答如果被采纳可获得 50 分。
<p>我正在使用NetCDF格式的一些大气模型输出。我想计算850百帕风的温度平流。我试着使用metpy函数
<a href="https://unidata.github.io/MetPy/latest/api/generated/metpy.calc.advection.html#metpy.calc.advection" rel="nofollow noreferrer">https://unidata.github.io/MetPy/latest/api/generated/metpy.calc.advection.html#metpy.calc.advection</a>
同样的。我附加了我的代码的一部分。请帮助解决此错误。所有必需的库都将更新为其最新版本</p>
<p><code>import xarray as xr</code>
<code>import numpy as np</code>
<code>import matplotlib.pyplot as plt</code>
<code>from metpy.units import units</code>
<code>import cartopy.crs as ccrs</code>
<code>import cartopy.feature as cfeature</code>
<code>import metpy.calc as mpcalc</code></p>
<p><code>temp = xr.open_dataset('ts.nc')</code>
<code>u10 = xr.open_dataset('u10.nc')</code>
<code>v10 = xr.open_dataset('v10.nc')</code></p>
<p><code>lat = temp['lat']</code>
<code>lon = temp['lon']</code>
<code>lonn, latt = np.meshgrid(lon,lat)</code></p>
<p><code>dx, dy = mpcalc.lat_lon_grid_deltas(lon, lat, )</code></p>
<p><code>time = np.array(temp.time[0:9])</code></p>
<p><code>T1 = temp['ts'][0:9,:,:]</code></p>
<p><code>u10_1 = u10['u10'][0:9,:,:]</code><br/>
<code>v10_1 = v10['v10'][0:9,:,:]</code></p>
<p><code>T11 = xr.DataArray(data = T1.values, dims = ['time', 'lat', 'lon'], coords=dict(lat=lat, lon = lon, time = time))</code>
<code>u10_11 = xr.DataArray(data = u10_1.values, dims = ['time', 'lat', 'lon'], coords=dict(lat=lat, lon = lon, time = time))</code>
<code>v10_11 = xr.DataArray(data = v10_1.values, dims = ['time', 'lat', 'lon'], coords=dict(lat=lat, lon = lon, time = time))</code></p>
<p><code>advec = mpcalc.advection(T1, [u10_1, v10_1], (dx, dy))</code></p>
<p>[守则][1]</p>
<p><strong>错误:属性错误:crs属性不可用。</strong></p>
<p>如果需要更多的参考资料,请告诉我</p>