<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="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/<wbr>nclscripts/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-<wbr>NCL/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@<wbr>mpUSStateLineThicknessF     = 3.0
  map_res@<wbr>mpNationalLineThicknessF    = 3.0
  map_res@<wbr>mpGeophysicalLineThicknessF = 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="HOEnZb"><font color="#888888"><div><br></div><div>matt</div><div><br></div><div><br></div></font></span></div>
</blockquote></div><br></div>