叠加风矢量图报错

Augleo

这是怎么回事呢：fatal:Unable to Destroy (Bad PID#-4)
fatal:VectorFieldInitialize: dimensions of vfUDataArray and vfVDataArray are inconsistent
fatal:Unable to initialize layer-Can't Create
fatal:Unable to access object with id:-4
fatal:PID #-4 can't be found in NhlSetValues
fatal:PID #-4 can't be found in NhlSetValues
fatal:Unable to access object with id:-4
warning:Unable to add DataItem "(null)" to DataList "vcVectorFieldData"
fatal:CompileDataList:DataList has no valid members

代码如下：；取文件
so2_files = systemfunc("ls "+outdir+"so2*")
f1 = addfiles(so2_files,"r")
V_files = systemfunc("ls "+outdir2+"v2*")
f2 = addfiles(V_files,"r")
U_files = systemfunc("ls "+outdir3+"u2*")
f3 = addfiles(U_files,"r")
f = addfile("/home/login/chenxsh/kc/data.wrfchem/wrfout_d01_2020-01-01_01:00:00.nc","r")
；读取变量
so2 = f->so2
u1=f->U
v1=f->V
lat1 = f->XLAT
lon1 = f->XLONG
SO2_00=new(dimsizes(so2),typeof(so2))
V_00=new(dimsizes(v1),typeof(v1))
U_00=new(dimsizes(u1),typeof(u1))
do k = 1,70,3；循环读取相应变量
        SO2_1=f1[k]->SO2
        SO2_2=f1[k+1]->SO2
        SO2_3=f1[k+2]->SO2
        SO2_00=(SO2_1+SO2_2+SO2_3)/3
        SO2_00=SO2_00#1000
        V_1=f2[k]->V(0,0,:,:)
        V_2=f2[k+1]->V(0,0,:,:)
        V_3=f2[k+2]->V(0,0,:,:)
        V_00=(V_1+V_2+V_3)/3
        U_1=f3[k]->U(0,0,:,:)
        U_2=f3[k+1]->U(0,0,:,:)
        U_3=f3[k+2]->U(0,0,:,:)
        U_00=(U_1+U_2+U_3)/3

vdef=-9999
  nx = 162
  ny = 119;
SO2S = new((/1,1,ny,nx/),float,vdef)
  SO2S2 = new((/1,ny,nx/),float,vdef)
VS = new((/1,1,ny,nx/),float,vdef)
VS2 = new((/1,ny,nx/),float,vdef)
US = new((/1,1,ny,nx/),float,vdef)
US2 = new((/1,ny,nx/),float,vdef)
  lat = new((/ny,nx/),float,vdef)
  lon = new((/ny,nx/),float,vdef)
  lat(:,:) = lat1(0,:,:)
  lon(:,:) = lon1(0,:,:)
  lat!0="lat"
  lat!1="lon"
  lon!0="lat"
  lon!1="lon"
  lon@units="degree_east"
  lat@units="degree_north"
; printVarSummary(lat1)
  SO2S(0,0,:,:) = SO2_00(0,0,:,:)
  SO2S2(0,:,:) = SO2S(0,0,:,:)
  SO2S2@lat2d = lat
  SO2S2@lon2d = lon
  SO2S2!0 = "lev"
  SO2S2!1 = "lat"
  SO2S2!2 = "lon"
type = "png"
  type@wkWidth = 1000.*2.
  type@wkHeight = 1000.*2.*0.85
; dir_fig = "fig.wrfout_d01_010100/"
  ;wks = gsn_open_wks(type,"dir_fig")
SO2S2@long_name = "so2 concentration(ug/m3) at "+k
;fname = "so2 at 20200101 1:00-3:00"
wks = gsn_open_wks(type,k+"wind")
gsn_define_colormap(wks,"WhiteBlueGreenYellowRed")

这里省略了一些res
plot =gsn_csm_contour_map(wks,SO2S2(0,:,:),res1)
draw(plot)
vplot =gsn_csm_vector(wks,U_00,V_00,resv)
overlay(plot,vplot)
frame(wks)
end do


end

Augleo

解决啦，是经纬度的格点范围不一样，我直接手动设置了他们的范围，让他们一样，就可以出图了

zyj1210

Augleo 发表于 2022-9-27 20:35
解决啦，是经纬度的格点范围不一样，我直接手动设置了他们的范围，让他们一样，就可以出图了

您好，请问怎样设置格点范围呢？

KasuganoSora

zyj1210 发表于 2023-4-13 16:09
您好，请问怎样设置格点范围呢？

  res@mpMinLatF
  res@mpMaxLatF
  res@mpMinLonF
  res@mpMaxLonF

zyj1210

KasuganoSora 发表于 2023-4-13 17:21
res@mpMinLatF
  res@mpMaxLatF
  res@mpMinLonF

好的，谢谢！我发现我报错的原因是之前读的是“U”这个变量，改为读取“U10”这个变量就正确啦

Jyy1108

KasuganoSora 发表于 2023-4-13 17:21
res@mpMinLatF
  res@mpMaxLatF
  res@mpMinLonF

楼主您好，我也出现了同样的报错信息，画图时加上格点范围也依旧报错，请问楼主还知道其他解决办法吗

KasuganoSora

Jyy1108 发表于 2024-5-12 20:39
楼主您好，我也出现了同样的报错信息，画图时加上格点范围也依旧报错，请问楼主还知道其他解决办法吗

这个应该是画图UV变量大小不一致，问题应该在数据上

Jyy1108

KasuganoSora 发表于 2024-5-13 17:06
这个应该是画图UV变量大小不一致，问题应该在数据上

好嘟，谢谢楼主答疑解惑，我再检查一下；我记得我好像是改了一下参考箭头长度，就没有报错了