读取/home/mw/input/fnl6384/fnl/fnl_20210101_00_00.grib2文件,计算2021年1月1日0时,北半球(含赤道)地表以上 2 米(👀 注意:不是 2 米以上)温度最高值和最低值之差(保留整数)。
提示:注意是北半球!可以用类似于 NumPy 中的切片方法,截取出北半球部分,然后再进行计算。具体的截取方法,请查看 xarray 的用户手册(Getting Started (xarray.dev)xarray 的用户手册(,从中寻找。
import xarray as xr
ds = xr.open_dataset("/home/mw/input/fnl6384/fnl/fnl_20210101_00_00.grib2",engine="pynio") # 读取数据
tem=ds['TMP_P0_L103_GLL0'][0,][0:91,:] #第二个方括号索引是第0维:地表2米数据;第三个方括号索引是经纬度坐标(赤道及以北)。
t= tem.max()-tem.min() #最大值-最小值
a1=int(t)
a1