想用ncl作图，用fnl格式为grib2的资料画整层的水汽通量图，有脚本，画不出图，求指点

liutong

load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/contributed.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/shea_util.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/skewt_func.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/wind_rose.ncl"
;*********************************************************
begin
plot=new(2,graphic)                              
grb_file = addfile("D:/NCL/cygwin/ziliao/fnl_20150629_18_00.grib2","r")
print
;500hgt*********************************************************************
  tt       = grb_file->TMP_3_ISBL_10({850},:,:)
  rhprs    = grb_file->R_H_3_ISBL_10({850},:,:)
  u        = grb_file->U_GRD_3_ISBL_10({850},:,:)
  v        = grb_file->V_GRD_3_ISBL_10({850},:,:)
  tc       = tt-273.16
  lev      = tt@lv_ISBL3
  td=tc-((14.55+0.114*tc)*(1-0.01*rhprs)+(((2.5+0.007*tc)*(1-0.01*rhprs))^3)+(15.9+0.37*tc)*((1-0.01*rhprs)^14))
  vapr=6.112*exp((17.67*td)/(td+243.5))
  ee=vapr*1.001+(lev-100)/900*0.0034
  mixr=0.62137*(ee/(lev-ee))*1000/9.8
  qu=u*mixr
  qv=v*mixr

  mm=sqrt(qu^2+qv^2)
   do i=0,359
    do j=0,180
         if (mm(j,i).lt.15)then
              mm(j,i)=0
          end if
    end do
   end do

  copy_VarCoords(tt,qu)
  copy_VarCoords(tt,qv)
  copy_VarCoords(rhprs,mm)
  copy_VarMeta(tt,qu)
  copy_VarMeta(tt,qv)
  copy_VarMeta(rhprs,mm)

  wks =gsn_open_wks("pdf","fnl_vapour_flux")
  gsn_define_colormap(wks,"BlAqGrYeOrReVi200")   
  vcres = True
  vcres@gsnAddCyclic  = False
  vcres@gsnDraw      = False                    ; don't draw yet
  vcres@gsnFrame     = False                        ; don't advance frame yet

  vcres@mpDataSetName         = "Earth..3"   ; This new database contains  
  vcres@mpDataBaseVersion     = "MediumRes"  ; Medium resolution database
  vcres@mpOutlineOn           = True         ; Turn on map outlines
  vcres@mpProjection          = "CylindricalEquidistant"
  vcres@tiMainString          =""
  vcres@mpGeophysicalLineThicknessF = 0.35      
; double the thickness of geophysical boundaries
  vcres@mpNationalLineThicknessF    = 0.5     
; double the thickness of national boundaries
  vcres@pmTickMarkDisplayMode = "Always"
  vcres@mpMinLatF         =  -25                         ;鏈

andrewsoong

啥错？

liutong

andrewsoong 发表于 2016-4-25 20:59
啥错？

脚本没发完
vcres@mpMaxLatF         =  50                        ;鏈

liutong

andrewsoong 发表于 2016-4-25 20:59
啥错？

fatal:["Execute.c":6332]:variable (TMP_3_ISBL_10) is not in file (grb_file)
fatal:["Execute.c":8578]:Execute: Error occurred at or near line 12 in file /vaper_flux_wind (1).ncl

andrewsoong

liutong 发表于 2016-4-25 21:02
fatal:["Execute.c":6332]:variable (TMP_3_ISBL_10) is not in file (grb_file)
fatal:["Execute.c":85 ...

这个变量不在，你是不是写错变量名了

liutong

andrewsoong 发表于 2016-4-25 21:02
这个变量不在，你是不是写错变量名了

检查了好多遍检查不出来，我把脚本发给你，能帮忙给看一下吗？

andrewsoong

liutong 发表于 2016-4-25 21:05
检查了好多遍检查不出来，我把脚本发给你，能帮忙给看一下吗？

你去这个QQ群里讨论吧176212282

liutong

好的，谢谢哈

好多鱼

很好的帖子，留在以后用了。

Lucia

请问楼主 脚本怎么改好了的