[ncl-talk] gsn_csm_pres_hgt_streamline problem

[Ioana Colfescu]

Hi,
I'm trying to plot the Walker cell and I m using  the
gsn_csm_pres_hgt_streamline function.
I read the previous posts about "errors" in plotting the streamlines and I
scaled the w but I still get a weird result - please see attached.
My problem is the W doesn't seem right and the direction of the streamlines
is not what it should be.
Below is my code - could you please point of what you think I'm doing wrong
on this one ? ( the plot is from a two plots version of the script but the
script below is the same).
Thanks.

begin
system("clear")
exp_name = (/"BTAL","BTALnEU","BTALnNA"/)
ens_nr = (/"1","2","3","4","5","6","7","8"/)

fw    = addfile("jjas_W__5-40N_30lev_BTAL.nc","r")
fu    = addfile("jjas_U__5-40N_30lev_BTAL.nc","r")
w = fw->w
u = fu->w

printVarSummary(w)
wavg1= dim_avg_n_Wrap(dim_avg_n(w,0),0)
uavg1= dim_avg_n_Wrap(dim_avg_n(u,0),0)
scale = abs(uavg1)/abs(wavg1_b)
wavg1 = scale*wavg1_b  ;

 wks=gsn_open_wks("x11","w_u_vertical_profile_"+lat1+"N")
 plot = new(1,graphic)
 res=True
 res at gsnDraw  = False                          ; don't draw
 res at gsnFrame = False
 res at cnLinesOn=False
 res at cnFillOn          = True                ; turn on color fill
 res at cnFillPalette     = "BlWhRe" ; choose color map
; res at cnLineLabelsOn    = False               ; turn off line labels
; res at lbLabelStride     = 2                   ; label every other box
 res at lbLabelBarOn         = False
plot(0)=gsn_csm_pres_hgt_streamline(wks,uavg1(0:5,{lat1},:),uavg1(0:5,{lat1},:),wavg1(0:5,{lat1},:),res)
 gsn_panel(wks,plot,(/1,1/),False)



*Ioana Colfescu *

*Research Scientist *

*National Centre for Atmospheric Science *

*Address: National Centre for Atmospheric Science, Fairbairn House, 71-75
Clarendon Road, Leeds, LS2 9PH
<https://maps.google.com/?q=71-75+Clarendon+Road,+Leeds,+LS2+9PH&entry=gmail&source=g>*
*University of Leeds *
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20190926/02035708/attachment-0001.html>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: w_u_vertical_profile_0N.pdf
Type: application/pdf
Size: 297610 bytes
Desc: not available
URL: <http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20190926/02035708/attachment-0001.pdf>