<div dir="ltr">Hello,<div><br></div><div>After some good help from some list members, I believe that I am setting up the resources for plotting WRF data without the WRF functions.   However, I still get no streamlines in my plot, or at least not in the correct place.</div><div><br></div><div>I am using NCL 6.1.2.  My WRF domain is close to CONUS.  The domain is LambertConformal and I am setting the mapping properties stored in the Global section of the NetCDF file to the appropriate variables for the LC mapping.   To draw the streamlines, I am basing my code off of the example here:</div><div><br></div><div><a href="https://www.ncl.ucar.edu/Applications/Scripts/gsn_stream_4.ncl" target="_blank">https://www.ncl.ucar.edu/Applications/Scripts/gsn_stream_4.ncl</a><br></div><div><br></div><div>I read in the data, massage it, then try and plot it.  When I plot the speed, u1, or v1 arrays with contour:</div><div><br></div><div><div>        contour = gsn_csm_contour_map(wks,speed,res)</div><div><br></div></div><div> I get what I would want, a contour of that variable with the bounds of the plot that of my domain.  When I try the same thing to create a streamline plot:</div><div><br></div><div>       plot = gsn_streamline_scalar_map(wks,u1,v1,speed,res)<br></div><div><br></div><div>I get nothing.  When I remove setting the Corners of the domain by commenting out the following:</div><div><br></div><div><div> res@mpRightCornerLatF                 = lat2d(0,0)</div><div> res@mpRightCornerLonF                = lon2d(0,0)</div><div> res@mpLeftCornerLatF                   = lat2d(nlat-1,nlon-1) </div><div> res@mpLeftCornerLonF                  = lon2d(nlat-1,nlon-1)</div><div> res@mpLimitMode                           = &quot;Corners&quot;</div></div><div><br></div><div>I get streamlines, but they are between 0E and 225E and 0N to 90N.  For the CONUS domain I am using the corners are about (138W,21N) and (298W,47N).  Something is being calculated, but not what I was expecting.</div><div><br></div><div>Any ideas are appreciated.</div><div><br></div><div>Thanks,</div><div>Craig</div><div><br></div><div><div>load &quot;$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl&quot;</div><div>load &quot;$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl&quot;</div><div>load &quot;$NCARG_ROOT/lib/ncarg/nclscripts/wrf/WRFUserARW.ncl&quot;</div><div><br></div><div>begin</div><div>;---Open WRF output file.</div><div>  dir      = &quot;/scratch2/portfolios/BMC/winds/Craig.Tierney/theia/wrf/windswrf/trunk/conusplus-wrf-6day/201308010000/nowfp-nocu/wrf/&quot;</div><div>  filename = &quot;wrfout_d01_2013-08-04_00:00:00&quot;</div><div>  a = addfile(dir + filename + &quot;.nc&quot;,&quot;r&quot;)</div><div><br></div><div>;---Read terrain height and lat/lon off file.</div><div>  it        = 0     ; first time step</div><div>  du1       = wrf_user_getvar(a,&quot;U10&quot;,it)    ; Terrain elevation</div><div>  dv1       = wrf_user_getvar(a,&quot;V10&quot;,it)    ; Terrain elevation</div><div>  dlon2d=wrf_user_getvar(a,&quot;XLONG&quot;,it)</div><div>  dlat2d=wrf_user_getvar(a,&quot;XLAT&quot;,it)</div><div><br></div><div>; Down sample hack</div><div> d=dimsizes(du1)</div><div> nx=d(0)</div><div> ny=d(1)</div><div> s=8</div><div> u1=du1(0:nx-1:s,0:ny-1:s)</div><div> v1=dv1(0:nx-1:s,0:ny-1:s)</div><div> lat2d=dlat2d(0:nx-1:s,0:ny-1:s)</div><div> lon2d=dlon2d(0:nx-1:s,0:ny-1:s)</div><div><br></div><div> u1@lon2d=lon2d</div><div> u1@lat2d=lat2d</div><div> v1@lon2d=lon2d</div><div> v1@lat2d=lat2d</div><div><br></div><div>  speed=sqrt(u1*u1+v1*v1)</div><div>  speed@lat2d=lat2d</div><div>  speed@lon2d=lon2d</div><div><br></div><div>  wks = gsn_open_wks(&quot;png&quot;,&quot;wrf_gsn&quot;)</div><div><br></div><div>;---Set some basic plot options</div><div>  res               = True</div><div><br></div><div>  res@gsnMaximize   = True   ; maximize plot in frame</div><div><br></div><div>  res@tiMainString  = filename</div><div><br></div><div>  ;res@mpProjection  = &quot;CylindricalEquidistant&quot;   </div><div>  res@mpProjection  = &quot;LambertConformal&quot; </div><div><br></div><div>      res@mpLambertParallel1F = a@TRUELAT1         ; two parallels</div><div>      res@mpLambertParallel2F = a@TRUELAT2</div><div>      res@mpLambertMeridianF  = a@STAND_LON        ; central meridian</div><div>      res@mpCenterLonF        = a@CEN_LON</div><div>      res@mpCenterLatF       = a@CEN_LAT</div><div><br></div><div>    d=dimsizes(u1)</div><div>    nlat=d(0)</div><div>    nlon=d(1)</div><div><br></div><div>;;;; Comment this section out, and we get streamlines, but in the wrong place</div><div>;;;; This is why they don&#39;t show up in the zoomed in plot</div><div>;    res@mpRightCornerLatF                 = lat2d(0,0)</div><div>;    res@mpRightCornerLonF                = lon2d(0,0)</div><div>;    res@mpLeftCornerLatF                   = lat2d(nlat-1,nlon-1)</div><div>;    res@mpLeftCornerLonF                  = lon2d(nlat-1,nlon-1)</div><div>;    res@mpLimitMode                           = &quot;Corners&quot;</div><div><br></div><div>;---Additional resources desired</div><div>  res@pmTickMarkDisplayMode = &quot;Always&quot;   ; nicer tickmarks</div><div><br></div><div>  res@mpDataBaseVersion     = &quot;MediumRes&quot;       ; better and more map outlines</div><div>  res@mpDataSetName         = &quot;Earth..4&quot;</div><div>  res@mpOutlineBoundarySets = &quot;AllBoundaries&quot;</div><div>  res@mpOutlineOn           = True</div><div><br></div><div>  res@lbOrientation         = &quot;Vertical&quot;</div><div>  res@tiMainOffsetYF        = -0.03           ; Move the title down</div><div><br></div><div>;; This is where we decide what to plot</div><div>  ptype=&quot;contour&quot;</div><div>;  ptype=&quot;streamline&quot;</div><div><br></div><div>  if (ptype .eq. &quot;contour&quot;) then</div><div><span class="" style="white-space:pre">        </span>print(&quot;Plotting: &quot;+ptype)</div><div><span class="" style="white-space:pre">        </span>res@cnFillOn      = True</div><div><span class="" style="white-space:pre">        </span>res@cnFillPalette = &quot;OceanLakeLandSnow&quot;</div><div><span class="" style="white-space:pre">        </span>res@cnLinesOn     = False</div><div><span class="" style="white-space:pre">        </span>;contour = gsn_csm_contour_map(wks,speed,res)</div><div><span class="" style="white-space:pre">        </span>;contour = gsn_csm_contour_map(wks,u1,res)</div><div><span class="" style="white-space:pre">        </span>contour = gsn_csm_contour_map(wks,v1,res)</div><div>  else</div><div><span class="" style="white-space:pre">        </span>if (ptype .eq. &quot;streamline&quot;) then</div><div><span class="" style="white-space:pre">                </span>print(&quot;Plotting: &quot;+ptype)</div><div><span class="" style="white-space:pre">                </span>res@pmTickMarkDisplayMode    = &quot;Always&quot;</div><div>  <span class="" style="white-space:pre">                </span>res@pmLabelBarDisplayMode    = &quot;Always&quot;</div><div><span class="" style="white-space:pre">                </span>res@pmLabelBarOrthogonalPosF = -0.02</div><div><span class="" style="white-space:pre">                </span>res@pmLabelBarWidthF         = 0.1</div><div><span class="" style="white-space:pre">                </span>res@lbPerimOn                = False</div><div><br></div><div><span class="" style="white-space:pre">                </span>plot = gsn_streamline_scalar_map(wks,u1,v1,speed,res)</div><div><span class="" style="white-space:pre">        </span>else</div><div><span class="" style="white-space:pre">                </span>print(&quot;Unable to plot either contour or streamline: &quot;+ptype)</div><div>        <span class="" style="white-space:pre">        </span>exit</div><div><span class="" style="white-space:pre">        </span>end if</div><div>  end if</div><div><br></div><div>end</div></div><div><br></div></div>