求职，ncl绘制垂直廓线问题--gsn_csm_pres_hgt_vector函数已解决

litg04

请教大家一个问题
如图所示，我用16年7月的fnl资料画了同时刻两种图，一张图中是850hPaV风场和温度，另一张是117E经度垂直剖面图，箭头是V-W风场，阴影为温度剖面图用的函数是gsn_csm_pres_hgt_vector。
问题是，在117E位置上，第一张图上30~40N之间吹得是南风，而在剖面图上相同纬度上确是北风，找了好久也没发现问题出在哪里
还请大家不吝赐教啊{:eb511:}


下面是画第二张图的程序：
files = systemfunc ("ls data/fnl/fnl*.grib2")
   f    = addfiles (files,"r")              
  ListSetType (f,"cat")
  t=new((/20,26,181,360/),"float")
  v=new((/20,26,181,360/),"float")
  omega=new((/20,21,181,360/),"float")
  do i=0,19
read_t=f->TMP_P0_L100_GLL0({100000:1000},:,:)
  read_v=f->VGRD_P0_L100_GLL0({100000:1000},:,:)
read_omega=f->VVEL_P0_L100_GLL0({100000:1000},:,:)
t(i,:,:,:)=read_t
v(i,:,:,:)=read_v
omega(i,:,:,:)=read_omega
  end do
  printVarSummary(v)
  
air=t-273.15
air!0="time"
air!1="lev"
air!2="lat"
air!3="lon"
air&lev=t&lv_ISBL0
air&lat=t&lat_0
air&lon=t&lon_0

lev_omega=omega&lv_ISBL6

scale_s=fabs(avg(v(:,{lev_omega},:,{117.5}))/avg(omega(:,:,:,{117.5})*9.8))
omega_s=-omega*scale_s
omega_s!1="level"
omega_s!2="lat"
omega_s!3="lon"
omega_s&level=omega&lv_ISBL6
omega_s&lat=omega&lat_0
omega_s&lon=omega&lon_0

r es=True
res@gsnMaximize=True
res@gsnAddCyclic         = False
  res@gsnDraw      = False                        ; don't draw yet
  res@gsnFrame     = False                        ; don't advance frame yet
res@gsnLeftString =""
  res@gsnRightString =""
  res@tiMainString =""
res@cnLineLabelsOn  = False                ; turn off line labels
res@lbLabelBarOn = False
  res@cnFillOn        = True                 ; turn on color fill
  res@cnInfoLabelOn = False
  res@cnLevelSelectionMode = "ExplicitLevels"
  res@cnLevels= (/-35,-30,-25,-20,-15,-10,-5,-3,-2,-1,0,1,2,3,5,10,15,20,25,30,35/)
  res@cnFillPalette  = "BlueWhiteOrangeRed"
  res@lbLabelStride   = 2                    ; every other color

  res@vcLineArrowThicknessF  = 10.0
  res@vcRefAnnoOrthogonalPosF  = -1.0            ;ref的位置
  res@vcRefAnnoString1On      = True
  res@vcRefAnnoString2On      = False
  res@vcRefMagnitudeF         = 10                 ; make vectors larger
  res@vcRefLengthF            = 0.04               ; ref vec length
  res@vcGlyphStyle            = "CurlyVector"     ; turn on curly vectors
  res@vcMinDistanceF          =0.01                ;箭头疏密
res@vcFillArrowWidthF       =10000

  plot=gsn_csm_pres_hgt_vector(wks,air(it,{lev_omega},{50:10},{117.5}),\
     v(it,{lev_omega},{50:10},{117.5}),omega_s(it,{lev_omega},{50:10},{117.5}),res)

lovjiangyp

膜拜膜拜膜拜

litg04

问题已经解决，\(^o^)/~

自导自演

请问楼主怎么解决的

yeah...

litg04 发表于 2016-11-17 16:48
问题已经解决，\(^o^)/~

请问是如何解决的啊 我也遇到了同样的问题 不同的是我的w是相反的

寰殇丶天使

楼主你几年前的帖子安装WRF:安装WRF遇到问题
这一堆个:REGISTRY WARNING: o3rad: o3rad_b is not a variable or number; ignoring it......
怎么解决的啊

litg04

yeah... 发表于 2016-12-3 13:21
请问是如何解决的啊 我也遇到了同样的问题 不同的是我的w是相反的

是个低级错误，路径的问题。程序本身没有错误

litg04

寰殇丶天使 发表于 2016-12-9 15:57
楼主你几年前的帖子安装WRF:安装WRF遇到问题
这一堆个:REGISTRY WARNING: o3rad: o3rad_b is not a variab ...

记不大清了，貌似是编译器配置的问题吧，用的intel的编译器，编译器配置好了后全部用intel的编译器重新编译下试试看，记不大请了

abc258886

看看嫩不能用上！！！！！！！！！！！

吴阿好大人

请问楼主怎么解决的