<div dir="ltr"><div class="gmail_default" style="font-size:small">Matt,</div><div class="gmail_default" style="font-size:small"><br></div><div class="gmail_default" style="font-size:small">NCL is plotting this data correctly as far as I can tell.</div><div class="gmail_default" style="font-size:small"><br></div><div class="gmail_default" style="font-size:small">When I&#39;m debugging a WRF plot, I usually like to plot the HGT variable, because then I can usually tell if the map projection is correct.  I also like to plot the lat/lon locations (from XLAT and XLONG on the file) as filled dots on the plot, so I can see where the actual plot area is.</div><div class="gmail_default" style="font-size:small"><br></div><div class="gmail_default" style="font-size:small">I took your script and modified it to do both a native filled contour/map (calling wrf_map_resources) and a non native filled contour/map (using XLAT/XLONG read off the file).</div><div class="gmail_default" style="font-size:small"><br></div><div class="gmail_default" style="font-size:small">Both plots indicate that this data is indeed over the eastern part of the U.S. See attached modified scripts and PNGs.</div><div class="gmail_default" style="font-size:small"><br></div><div class="gmail_default" style="font-size:small">--Mary</div><div class="gmail_default" style="font-size:small"><br></div></div><div class="gmail_extra"><br><div class="gmail_quote">On Tue, Feb 7, 2017 at 4:41 PM, Matt Masarik <span dir="ltr">&lt;<a href="mailto:mattmasarik@boisestate.edu" target="_blank">mattmasarik@boisestate.edu</a>&gt;</span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr">ps, I&#39;m working on Yellowstone using the following modules;<div><br></div><div>ncarenv/1.0</div><div>ncarbinlibs/1.1</div><div>intel/12.1.5</div><div>ncarcompilers/1.0</div><div>netcdf/4.3.0</div><div>ncl/6.3.0</div><div><br></div></div><div class="HOEnZb"><div class="h5"><div class="gmail_extra"><br><div class="gmail_quote">On Tue, Feb 7, 2017 at 4:34 PM, Matt Masarik <span dir="ltr">&lt;<a href="mailto:mattmasarik@boisestate.edu" target="_blank">mattmasarik@boisestate.edu</a>&gt;</span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr">I tried plotting some WRF data today with an NCL script that has worked correctly previously, and found that the location of the domain is incorrect.  The domain should be centered over / containing Idaho (verified with ncview), but the output from the ncl script shows the domain located over Florida, containing the SE US (attached).  <div><br></div><div>Here&#39;s the (editted) script:</div><div>------------------------------<wbr>------------------------------<wbr>------------------------------<wbr>----------------</div><div><pre style="color:rgb(0,0,0);word-wrap:break-word;white-space:pre-wrap">;-----------------------------<wbr>------------------------------<wbr>-----------
; wrf_gsn_6.ncl
;-----------------------------<wbr>------------------------------<wbr>-----------
; Concepts illustrated:
;   - Using gsn_csm scripts to plot WRF-ARW data
;   - Overlaying filled contours, and vectors on a map
;   - Setting the correct WRF map projection using wrf_map_resources
;   - Setting lots of resources to customize a WRF plot
;-----------------------------<wbr>------------------------------<wbr>-----------
; This script is meant to show the difference between plotting WRF
; data using wrf_xxx scripts, and using gsn_csm_xxx scripts.
;-----------------------------<wbr>------------------------------<wbr>-----------
; You can safely ignore these warning messages:
;
; warning:start_lat is not a valid resource in wrf_gsn_contour at this time
; warning:start_lon is not a valid resource in wrf_gsn_contour at this time
; warning:end_lat is not a valid resource in wrf_gsn_contour at this time
; warning:end_lon is not a valid resource in wrf_gsn_contour at this time
; warning:mpNestTime is not a valid resource in map at this time
;-----------------------------<wbr>------------------------------<wbr>-----------
load &quot;$NCARG_ROOT/lib/ncarg/nclscri<wbr>pts/csm/shea_util.ncl&quot;


begin

;---Parameters
  input_file         = &quot;wrfout_d01_2015-02-05_00:00:<a href="http://00.nc" target="_blank">0<wbr>0.nc</a>&quot;
  input_dir          = &quot;/glade/u/home/mmasarik/WRF-NC<wbr>L/input&quot;
  input_path         = input_dir + &quot;/&quot; +  input_file

  plot_name          = &quot;streamlines_station-points&quot;



;---Open WRF output file
  a                  = addfile(input_path,&quot;r&quot;)


;--Create Workstation
  wks                = gsn_open_wks(&quot;png&quot;,&quot;wrf_map&quot;)


;--Start resources
  res                                = True
  res@gsnFrame                       = False
  res@gsnDraw                        = False
  res@gsnLeftString                  = &quot;&quot;
  res@gsnRightString                 = &quot;&quot;
  res@tfDoNDCOverlay                 = True         ; necesssary


;---Map plot
  map_res                             = res
  map_res                             = wrf_map_resources(a,map_res)
  map_res@gsnFrame                    = False
  map_res@gsnDraw                     = False
  map_res@gsnStringFontHeightF        = 0.015
  map_res@tiMainJust                  = &quot;CenterLeft&quot;
  map_res@tiMainPosition              = &quot;Left&quot;
  map_res@tiMainOffsetYF              = -0.01


;---Set map resources based on projection on WRF output file
  map_res@mpFillOn                    = False
  map_res@mpOutlineOn                 = True
  map_res@mpDataBaseVersion           = &quot;MediumRes&quot;
  map_res@mpOutlineDrawOrder          = &quot;PostDraw&quot;
  map_res@mpUSStateLineColor          = &quot;Black&quot;
  map_res@mpPerimLineColor            = &quot;Black&quot;
  map_res@mpNationalLineColor         = &quot;Black&quot;
  map_res@mpLimbLineColor             = &quot;Black&quot;
  map_res@mpGridLineColor             = &quot;Black&quot;
  map_res@mpGeophysicalLineColor      = &quot;Black&quot;
  map_res@mpUSStateLineThickness<wbr>F     = 3.0
  map_res@mpNationalLineThicknes<wbr>sF    = 3.0
  map_res@mpGeophysicalLineThick<wbr>nessF = 3.0

  map = gsn_csm_map(wks,map_res)



;---Overlay plots on map and draw.
  draw(map)
  frame(wks)


end</pre><pre style="color:rgb(0,0,0);word-wrap:break-word;white-space:pre-wrap"><span style="font-family:arial,sans-serif;color:rgb(34,34,34)">------------------------------<wbr>------------------------------<wbr>------------------------------<wbr>----------------</span></pre></div><div><br></div><div>The only warnings I get are the one&#39;s mentioned above in the script.  What am I doing wrong here?</div><div><br></div><div>Thanks!</div><span class="m_-306333007365841813HOEnZb"><font color="#888888"><div><br></div><div>matt</div><div><br></div><div><br></div></font></span></div>
</blockquote></div><br></div>
</div></div><br>______________________________<wbr>_________________<br>
ncl-talk mailing list<br>
<a href="mailto:ncl-talk@ucar.edu">ncl-talk@ucar.edu</a><br>
List instructions, subscriber options, unsubscribe:<br>
<a href="http://mailman.ucar.edu/mailman/listinfo/ncl-talk" rel="noreferrer" target="_blank">http://mailman.ucar.edu/<wbr>mailman/listinfo/ncl-talk</a><br>
<br></blockquote></div><br></div>