垂直速度斜剖面图

seanice

沿(33N,112E)做垂直速度斜剖面图
之前没画过斜剖面图，网上找到的模板非常有限，希望有大佬帮忙看看

load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"   
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl"   
;************************************************
begin

  f = addfile("omega.2021.nc","r")
  p1=1000
  p2=200

  it=200
  omega = (f->omega(it,{p1:p2},{10:30},{300:325}))
  var_zm= omega(:,{33},{112})

  ;var_zm= omega(:,{33},{112})


; f = addfile("uwnd.2021.nc","r")
;var = short2flt( f->uwnd )
; var = f->uwnd
;var_zm = dim_avg_n_Wrap(var,2) ;纬向平均
  ;Vadv =var_zm(:,:,{112})

;;**plot**************************************
; wks = gsn_open_wks ("eps", "plot-cross-section")
wks = gsn_open_wks ("png", "omexie01")
gsn_define_colormap(wks,"rainbow+gray")

res = True
res@gsnDraw = False
res@gsnFrame = False
res@gsnLeftString = ""
res@gsnRightString = ""

res@tmXBMode = "Explicit"
;res@tmXBValues = (/-180,-120,-60,0,60,120,180/)
res@tmXBValues              = (/50,150,250,350,450,550/)   ;指定经纬度
res@tmXBLabels              = (/"10N,55W","15N,50W","20N,45W","25N,40W","30N,35W","35N,30W"/)
;res@tmXBValues = (/0,60,120,180,240,300,360/)
;res@tmXBLabels = (/"180~W~o~E~W","120~W~o~E~W","60~W~o~E~W","0","60~W~o~E~E","120~W~o~E~E","180~W~o~E~E"/)
;res@tmXBLabels = (/"0","60~W~o~E~E","120~W~o~E~E","180~W~o~E~E","240~W~o~E~E","300~W~o~E~E","360~W~o~E~E"/)
res@tmXBLabelFontHeightF = 0.02
res@tmYLLabelFontHeightF = 0.02

res@tiYAxisString = "Level (hPa)"
;res@tmYBValues = (/20,200,400,600,800,1000/)
;res@cnFillOn = False


  ;res  = True                                      ; plot mods desired
; res@tiMainString     = ""
  ;res@vcRefMagnitudeF  = 50.0                       ; add a reference vector
  ;res@vcRefLengthF     = 0.2                      ; what the ref length is
  ;res@vcGlyphStyle     = "CurlyVector"             ; turn on curly vectors

  res                       = True     ; plot mods desired
  res@gsnMaximize           = True     ; Make plot fit the frame
  res@cnFillOn              = True     ; turn on color fill
  res@cnLinesOn             = False    ; turn of contour lines
  res@cnLevelSpacingF       = 0.05    ; contour spacing
  res@cnFillPalette         = "BlAqGrYeOrRe"
  ;res@lbOrientation         = "Vertical"
  res@gsnAddCyclic          = False    ; data already has cyclic poin

; res@cnLevelSelectionMode = "ManualLevels" ; manually set the contour levels with the following 3 resources
;res@cnMinLevelValF  = -10.0         ; set the minimum contour level
;res@cnMaxLevelValF  = 10.           ; set the maximum contour level
;res@cnLevelSpacingF = 0.00025          ; set the interval between contours

   ;base=gsn_csm_pres_hgt(wks,var_zm,res) ;1 月的纬向平均 u
   base=gsn_csm_pres_hgt(wks,var_zm,res)
   ;base = gsn_csm_contour(wks,var_zm,res)
;base = ColorNegDashZeroPosContour(base,"blue","black","red")
;小于 0 用蓝色虚线，大于 0 用红色实线，0 值用黑线

draw(base)
frame(wks)
end