import netCDF4
import wrf


# 读取数据###############################################################开始
data_nc = netCDF4.Dataset(r"./data/WRF/wrfout_d01_2023-01-15_12%3A00%3A00")
# print(data_nc)
# print(data_nc.variables['Times'])
# print(wrf.getvar(data_nc, "Times"))
slp = wrf.getvar(data_nc, "slp", timeidx=wrf.ALL_TIMES)  # 读取计算海表面气压
# print(slp)
# 统计预报时间的次数
time_list = slp.coords['Time'].values.size
k = 1
for i in range(time_list):
    print(str(k)+':'+str(slp.coords['Time'].values[i]))
    k = k+1