python计算扰动位势高度（He）

什夜0.0

用python计算扰动位势高度（He），根据定义是从原始位势高度中减去同时的纬向平均位势高度（0°N–40°N，180°W–180°E），程序如下，但是得到的图像六月份还算正常，七八月份相差很大，请教各位大神这是什么原因呢？如有回答感激不尽！！


import matplotlib.pyplot as plt
import cartopy.crs as ccrs
import xarray as xr
import os
import numpy as np

def readat(year):
    pth = '/home/uranus/ERA-5/Daily/h500'
    da = xr.open_dataset(os.path.join(pth,'hgt.'+str(year)+'.nc'))
    dat = da['hgt'].loc[str(year)+'-05-01':str(year)+'-09-30',40:0,:]
    lon = dat['lon'].values
    lat = dat['lat'].values
    dat = dat.values
    dat = dat - np.moveaxis(np.tile(np.mean(dat,2),(len(lon),1,1)),0,2)

    return dat, lon, lat

if __name__=='__main__':
    for year in range(1979,2021):
       # print(year)
        dat, lon, lat = readat(year)
        if year==1979:
            datc = dat
        else:
            datc += dat
    datc = datc/(2020-1979+1)
    #print(datc.shape)
    da = xr.DataArray(datc,name='hgt',dims=('time','lat','lon'),
            coords={'time':np.arange(153),'lat':lat,'lon':lon})
    da.to_netcdf('climhe1.nc',encoding={'hgt':{'_FillValue':-9999.},
        'lon':{'_FillValue':-9999.},'lat':{'_FillValue':-9999.}})