- 积分
- 166
- 贡献
-
- 精华
- 在线时间
- 小时
- 注册时间
- 2021-9-23
- 最后登录
- 1970-1-1
|
楼主 |
发表于 2023-10-18 20:36:09
|
显示全部楼层
本帖最后由 爱吃三七花的LY 于 2023-10-18 20:40 编辑
以下是处理变量部分的程序
#%%sst场
f = xr.open_dataset('E:\DATA\sst.mnmean.V5.nc')
#转换经度坐标,将0-360,转换为-180~180
f.coords['lon'] = (f.coords['lon'] + 180) % 360 - 180
f = f.sortby(f.lon)
# 截取热带大西洋区域(20w-0,20S-20N)3,4,5月份海温数据MAM
f_data = f.sst.loc[f.time.dt.month.isin([3,4,5])].loc['1979-01-01':'2019-01-01',20:-20,-80:20]
lon = f_data.lon
lat = f_data.lat
sst=f_data
#%%降水
f1 = xr.open_dataset('E:\\DATA\\precip.mon.mean.nc')
#转换经度坐标,将0-360,转换为-180~180
f1.coords['lon'] = (f1.coords['lon'] + 180) % 360 - 180
f1 = f1.sortby(f1.lon)
f1 = f1.sortby(f1.lat, ascending=False)#纬度翻转
# 截取热带大西洋区域(20w-0,20S-20N)3,4,5月份降水数据MAM
f1_data = f1.precip.loc[f1.time.dt.month.isin([3,4,5])].loc['1979-01-01':'2019-01-01',20:-20,-80:20]
lon1 = f1_data.lon
lat1 = f1_data.lat
precip=f1_data
图片分别是两个变量的属性,发现这两个时间范围和长度都一样,所以很疑惑是哪里出错了。
|
|