爱气象,爱气象家园! 

气象家园

 找回密码
 立即注册

QQ登录

只需一步,快速开始

搜索
查看: 6951|回复: 7

叠加风矢量图报错

[复制链接]
发表于 2022-9-21 20:15:35 | 显示全部楼层 |阅读模式

登录后查看更多精彩内容~

您需要 登录 才可以下载或查看,没有帐号?立即注册

x
这是怎么回事呢: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









密码修改失败请联系微信:mofangbao
 楼主| 发表于 2022-9-27 20:35:48 | 显示全部楼层
解决啦,是经纬度的格点范围不一样,我直接手动设置了他们的范围,让他们一样,就可以出图了
密码修改失败请联系微信:mofangbao
回复 支持 1 反对 0

使用道具 举报

发表于 2023-4-13 16:09:12 | 显示全部楼层
Augleo 发表于 2022-9-27 20:35
解决啦,是经纬度的格点范围不一样,我直接手动设置了他们的范围,让他们一样,就可以出图了

您好,请问怎样设置格点范围呢?
密码修改失败请联系微信:mofangbao
回复 支持 反对

使用道具 举报

发表于 2023-4-13 17:21:45 | 显示全部楼层
zyj1210 发表于 2023-4-13 16:09
您好,请问怎样设置格点范围呢?

  res@mpMinLatF
  res@mpMaxLatF
  res@mpMinLonF
  res@mpMaxLonF
密码修改失败请联系微信:mofangbao
回复 支持 反对

使用道具 举报

发表于 2023-4-25 11:25:38 | 显示全部楼层
KasuganoSora 发表于 2023-4-13 17:21
res@mpMinLatF
  res@mpMaxLatF
  res@mpMinLonF

好的,谢谢!我发现我报错的原因是之前读的是“U”这个变量,改为读取“U10”这个变量就正确啦
密码修改失败请联系微信:mofangbao
回复 支持 反对

使用道具 举报

发表于 2024-5-12 20:39:47 | 显示全部楼层
KasuganoSora 发表于 2023-4-13 17:21
res@mpMinLatF
  res@mpMaxLatF
  res@mpMinLonF

楼主您好,我也出现了同样的报错信息,画图时加上格点范围也依旧报错,请问楼主还知道其他解决办法吗
密码修改失败请联系微信:mofangbao
回复 支持 反对

使用道具 举报

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

这个应该是画图UV变量大小不一致,问题应该在数据上
密码修改失败请联系微信:mofangbao
回复 支持 反对

使用道具 举报

发表于 2024-5-13 17:48:45 | 显示全部楼层
KasuganoSora 发表于 2024-5-13 17:06
这个应该是画图UV变量大小不一致,问题应该在数据上

好嘟,谢谢楼主答疑解惑,我再检查一下;我记得我好像是改了一下参考箭头长度,就没有报错了
密码修改失败请联系微信:mofangbao
回复 支持 反对

使用道具 举报

您需要登录后才可以回帖 登录 | 立即注册

本版积分规则

Copyright ©2011-2014 bbs.06climate.com All Rights Reserved.  Powered by Discuz! (京ICP-10201084)

本站信息均由会员发表,不代表气象家园立场,禁止在本站发表与国家法律相抵触言论

快速回复 返回顶部 返回列表