RH的计算

转自https://zhuanlan.zhihu.com/p/713234717

Temp = xr.open_dataset(r'E:\Sencond_droughandphenelogy\ERA5_temp_d2m2000-2020.nc')['t2m'][:,:720,:]
Demp = xr.open_dataset(r'E:\Sencond_droughandphenelogy\ERA5_temp_d2m2000-2020.nc')['d2m'][:,:720,:]
Zeros = np.zeros((252,720,1440))

for i in range(Demp.shape[0]):
    Td = np.array(Demp[i,:,:])
    Ta = np.array(Temp[i,:,:])
    RH1 = (17.269*(Td - 273.3))/Td
    RH2 = (17.269*(Ta - 273.3))/Ta

    Zeros[i,:,:] = np.exp(RH1 - RH2)