Re: WRF vertical velocity vectors on vertical cross section

Jonthan,

You are calling vinth2p with only 4 parameters, and it needs nine. I
think the error message you are getting
is not telling you that you need 9 vectors, but that you need 9 input
variables (or parameters).

Please see the documentation which has some examples:

http://www.ncl.ucar.edu/Document/Functions/Built-in/vinth2p.shtml

--Mary

> 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.
>
> <vert_xsec_with_w.rtf>_______________________________________________
> ncl-talk mailing list
> List instructions, subscriber options, unsubscribe:
> http://mailman.ucar.edu/mailman/listinfo/ncl-talk

_______________________________________________
ncl-talk mailing list
List instructions, subscriber options, unsubscribe:
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
Received on Wed Nov 18 12:06:33 2009