0

您好，我收到此错误：ValueError: setting an array element with a sequence。当我运行此代码时：

import metpy.calc as mpcalc

from metpy.units import units

import cartopy.crs as ccrs
import cartopy.feature as cfeat
import os

from netCDF4 import Dataset as netcdf_dataset


fname = os.path.join(config["repo_data_dir"],
                     'netcdf', "nam_218_20170625_1200_031.nc"
                     )

dataset = netcdf_dataset(fname)

lat = dataset.variables['gridlat_0'][:]
lon = dataset.variables['gridlon_0'][:]
sfct=units.degC*(np.mean(dataset.variables['TMP_P0_L103_GLC0'][:,:,:],0)-273.15)
dp=units.degC*(dataset.variables['DPT_P0_L103_GLC0'][:]-273.15)
sfcp=units.hPa*(dataset.variables['PRES_P0_L103_GLC0'][:]/100.0)

lcl_pressure=np.empty(((428),(614)))
lcl_temperature=np.empty(((428),(614)))
lclht=np.empty(((428),(614)))
for i in range(428):
    for j in range(614):
        lcl_pressure[i,j], lcl_temperature[i,j] = mpcalc.lcl(sfcp[i,j],sfct[i,j],dp[i,j])
        lclht[i,j]=mpcalc.pressure_to_height_std(lcl_pressure[i,j]).to('feet')

有问题的行是这样的：

lcl_pressure[i,j], lcl_temperature[i,j] = mpcalc.lcl(sfcp[i,j],sfct[i,j],dp[i,j])

因为我不确定这是如何导致它的，因为有问题的行适用于它所说的不同代码

lcl_pressure[i,j], lcl_temperature[i,j] = mpcalc.lcl(p[0],t[0],dp[0])

因为第二个版本中的 p、t 和 dp 也是数组

4

0 回答 0