ncl画wrf输出风场的涡度图

yusw

请问大家是怎么用wrfout风场数据画涡度图的，在官网上的脚本必须要求数据是经纬度一维格式的，wrf输出的数据经纬度是二维的，大神们是怎么处理的？

1649518749

;   Example script to produce Vorticity plots from WRF ARW model data
;   Novemner 2008

load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/wrf/WRFUserARW.ncl"

begin
;
; The WRF ARW input file.  
; This needs to have a ".nc" appended, so just do it.
  a = addfile("../wrfout_d01_2000-01-24_12:00:00.nc","r")

; We generate plots, but what kind do we prefer?； type = "x11"
; type = "pdf"
type = "ps"
; type = "ncgm"
  wks = gsn_open_wks(type,"plt_pv")

; Set some basic resources
  res = True
  res@MainTitle                   = "REAL-TIME WRF"

  pltres = True
  mpres = True


;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;

; What times and how many time steps are in the data set?
  times  = wrf_user_getvar(a,"times",-1) ; get times in the file
  ntimes = dimsizes(times)         ; number of times in the file

;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;

  do it = 0,ntimes-1,2             ; TIME LOOP

    print("Working on time: " + times(it) )
    res@TimeLabel = times(it)   ; Set Valid time to use on plots


  ; Get the data
    pvo   = wrf_user_getvar(a,"pvo",it)
    avo   = wrf_user_getvar(a,"avo",it)
    p     = wrf_user_getvar(a,"pressure",it)


  ; Interpolate to pressure
    pv_plane = wrf_user_intrp3d(pvo,p,"h",300.,0,False)
    av_plane = wrf_user_intrp3d(avo,p,"h",500.,0,False)


    ; Plotting options
      opts = res                        
      opts@cnFillOn = True  
      opts@gsnSpreadColorEnd = -3  ; End third from the last color in color map
      opts@ContourParameters = (/ 0., 100., 10./)
      contour_a = wrf_contour(a,wks,av_plane,opts)
      opts@ContourParameters = (/ 0., 5., 1./)
      contour = wrf_contour(a,wks,pv_plane,opts)
      delete(opts)

    ; MAKE PLOTS                                       
      plot = wrf_map_overlays(a,wks,(/contour_a/),pltres,mpres)
      plot = wrf_map_overlays(a,wks,(/contour/),pltres,mpres)


  end do        ; END OF TIME LOOP

;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;

end

zongqing123

很棒，学习了

sdhjdfjsb

这个是绝对涡度吗？

supersupergirl

请问相对涡度应该怎么计算呀

yusw

1649518749 发表于 2016-8-15 17:30
;   Example script to produce Vorticity plots from WRF ARW model data
;   Novemner 2008

查了这个函数，是绝对涡度哦，请问有画过相对涡度的吗？

1649518749

supersupergirl 发表于 2019-8-11 11:30
请问相对涡度应该怎么计算呀

http://www.ncl.ucar.edu/Document ... -in/uv2vr_cfd.shtml

1649518749

yusw 发表于 2019-8-30 10:46
查了这个函数，是绝对涡度哦，请问有画过相对涡度的吗？

http://www.ncl.ucar.edu/Document ... -in/uv2vr_cfd.shtml

supersupergirl

1649518749 发表于 2019-8-30 17:00
http://www.ncl.ucar.edu/Document/Functions/Built-in/uv2vr_cfd.shtml

好的 感谢！

yusw

1649518749 发表于 2019-8-30 17:03
http://www.ncl.ucar.edu/Document/Functions/Built-in/uv2vr_cfd.shtml

这个函数要求计算的uv风场，是等经纬格点的数据，wrfout输出的为等网格距的数据，直接用这个函数算不了