- 积分
- 67296
- 贡献
-
- 精华
- 在线时间
- 小时
- 注册时间
- 2011-7-23
- 最后登录
- 1970-1-1
|
发表于 2021-7-28 18:20:05
|
显示全部楼层
pint不算是数据格式,pint本身就是单位(千克、米、秒……)管理包。
metpy的优势在于在以前的单纯数据计算时,加入了单位,不再需要人为地解决开尔文vs摄氏度、克vs千克、米vs千米。你只要告诉函数你的数据是什么单位的就行了,函数内部会归一化。
mpcalc.divergence( u , v , dx=ddx , dy=ddy ,xdim=-1 ,ydim=-2)
其中u和v就不说了。dx是位置参数关键字,ddx是你算出来的差分距离(我数学不好,你应该能懂我的意思,比如1°分辨率的是110km),纬向同。
ats = ds.latitude.data
lons = ds.longitude.data
dx, dy = mpcalc.lat_lon_grid_deltas(lons, lats)
得到的就是这里的ddx,ddy。参照楼上“完整”代码。
PS. 我现在越来越觉得xarray比netCDF4好用了,当然也有些特例。 |
|