<div dir="ltr"><div>To whom it may concern,</div><div><br></div><div> I am trying to plot a pressure/height vs. latitude plot of zonal winds with WV wind vectors overlaid on top. I've been trying to follow the example shown in '<a href="https://www.ncl.ucar.edu/Applications/Scripts/vector_5.ncl">https://www.ncl.ucar.edu/Applications/Scripts/vector_5.ncl</a>' and use the 'gsn_csm_pres_hgt_vector' function, but every time the script gets to the point where it calls that function it just hangs and does nothing. I can successfully run 'gsn_csm_pres_hgt' and produce a plot with all of my variables so I know that they are formatted properly. I am using NCL version 6.6.2 so I don't think it is a version issue. I've pasted my plotting code below. Any advice would be greatly appreciated.</div><div><br></div><div>Plotting Script (the script hangs on the first plot(0) call):</div><div><i><br>  wks = gsn_open_wks("png", "/Users/cmaloney/Other_NCL_scripts/images/zonal_avgs/diff_plots/" + runa + "_" + runDef + "seasonal_wind_anomalies")<br>  plot = new(2,graphic)<br>  <br>  <br>  res = True<br>  res@gsnDraw         = False<br>  res@gsnFrame                = False <br>  res@cnFillOn         = True                <br>  res@cnLinesOn        = False              <br>    res@trYMinF          = 1.0<br>    res@cnFillPalette     = "BlueDarkRed18"    <br>    res@lbTitleString    = "Delta (m/sec)"<br>    res@lbTitlePosition  = "Bottom"  <br>    res@tiYAxisString     = "Pressure (hPa)"<br>    res@tmYLMode = "Explicit"<br>    res@tmYLValues = (/ 850, 500, 300, 200, 150, 100, 70, 50, 30, 10, 1 /)<br>    res@tmYLLabels = (/ "850", "500", "300", "200", "150", "100", "70", "50", "30", "10", "1" /)<br>    res@cnLevelSelectionMode = "ManualLevels"<br>    res@cnMaxLevelValF       = 2<br>    res@cnMinLevelValF       = -2<br>    res@cnLevelSpacingF      = 0.1<br>    res@pmLabelBarOrthogonalPosF = 0.025<br>    res@vcMapDirection = False<br>    res@vcGlyphStyle    = "CurlyVector"      ; turn on curly vectors<br>  <br> <br>    plot(0) = gsn_csm_pres_hgt_vector(wks, Zonal_DJF_U_anomaly, Zonal_DJF_vstar_anomaly, Zonal_DJF_wstar_anomaly, res)<br>    <br>    plot(1) = gsn_csm_pres_hgt_vector(wks, Zonal_JJA_U_anomaly, Zonal_JJA_vstar_anomaly, Zonal_JJA_wstar_anomaly, res)<br><br>  resP                  = True<br>  resP@gsnFrame         = False                  <br>  resP@gsnPanelLabelBar = True                  <br>  resP@lbLabelFontHeightF = 0.015 <br>  resP@gsnPanelBottom   = 0.2 <br>  resP@gsnPanelTop      = 0<br>  resP@gsnPanelYWhiteSpacePercent = 5.<br>  resP@gsnPanelXWhiteSpacePercent = 5.<br>  resP@gsnPanelMainString = "Wind Anomalies"    <br>  gsn_panel(wks,plot,(/1,2/),resP)<br> </i></div><div><i>    draw(plot)<br></i></div><div><i>    frame(wks)</i></div><div><br></div><div><br></div><div><br></div><div>Thanks for your help,</div><div>- Chris</div><div><br></div><div>-- <br><div dir="ltr" class="gmail_signature" data-smartmail="gmail_signature"><div dir="ltr"><div><div dir="ltr"><div><div><div><div><div>Christopher Maloney                      <br></div></div></div></div></div><div><div><div><div><div><div>Professional Research Associate<br></div><div>CIRES at the University of Colorado/NOAA Chemical Sciences Lab</div><div><br></div></div></div></div></div></div></div></div></div></div></div></div>