Hello,
I want to put WRF vertical velocity (w) vectors on vertical cross
section plots.
I have attached the script. I attempt to use the the outline from
wrfCrosssection3.ncl. I add the following additional lines:
u = wrf_user_getvar(a, "ua",it) ; u averaged to mass points ;
u = u * 1.94386 ; kts
u@units = "kts"
w = wrf_user_getvar(a, "wa",it) ; w averaged to mass points ;
w = w * 1.94386 ; kts
w@units = "kts"
in the loop I add:
wp = vinth2p(w,p,1,opts)
up = vinth2p(u,p,1,opts)
but it says nine vectors are needed
Then to actually print the vectors I add these lines:
; Plotting options for w wind vectors
opts_w = res
opts_w@vcRefAnnoOrthogonalPosF = -0.165 ; move
reference vec up
opts_w@pmLabelBarOrthogonalPosF = -0.05 ; move
color label bar up
opts_w@vcRefMagnitudeF = 0.04 ; define
vector ref mag
opts_w@vcRefLengthF = 0.045 ; define
length of vec ref
I add up and wp to this line
wrf_overlay(wks,(/contour_o3, up, wp/),True)
----------------------------------------------------------------
This message was sent using IMP, the Internet Messaging Program.
_______________________________________________
ncl-talk mailing list
List instructions, subscriber options, unsubscribe:
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
This archive was generated by hypermail 2.1.8 : Mon Nov 23 2009 - 12:24:02 MST