<div dir="ltr">Hi Mary,<div><br></div><div>thanks for working on this and sending editted scripts. I used the scripts and got correct output from each. User error was probably my issue...</div><div><br></div><div>I also appreciate the tip on plotting the XLAT/XLONG as dots, haven't done that before..</div><div><br></div><div>Thanks again!</div><div>Cheers,</div><div><br></div><div>Matt</div><div><br></div></div><div class="gmail_extra"><br><div class="gmail_quote">On Mon, Feb 13, 2017 at 9:45 AM, Mary Haley <span dir="ltr"><<a href="mailto:haley@ucar.edu" target="_blank">haley@ucar.edu</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><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'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><span class="HOEnZb"><font color="#888888"><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></font></span></div><div class="gmail_extra"><br><div class="gmail_quote"><div><div class="h5">On Tue, Feb 7, 2017 at 4:41 PM, Matt Masarik <span dir="ltr"><<a href="mailto:mattmasarik@boisestate.edu" target="_blank">mattmasarik@boisestate.edu</a>></span> wrote:<br></div></div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div><div class="h5"><div dir="ltr">ps, I'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="m_8792142955542432848HOEnZb"><div class="m_8792142955542432848h5"><div class="gmail_extra"><br><div class="gmail_quote">On Tue, Feb 7, 2017 at 4:34 PM, Matt Masarik <span dir="ltr"><<a href="mailto:mattmasarik@boisestate.edu" target="_blank">mattmasarik@boisestate.edu</a>></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'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 "$NCARG_ROOT/lib/ncarg/nclscri<wbr>pts/csm/shea_util.ncl"
begin
;---Parameters
input_file = "wrfout_d01_2015-02-05_00:00:<a href="http://00.nc" target="_blank">0<wbr>0.nc</a>"
input_dir = "/glade/u/home/mmasarik/WRF-NC<wbr>L/input"
input_path = input_dir + "/" + input_file
plot_name = "streamlines_station-points"
;---Open WRF output file
a = addfile(input_path,"r")
;--Create Workstation
wks = gsn_open_wks("png","wrf_map")
;--Start resources
res = True
res@gsnFrame = False
res@gsnDraw = False
res@gsnLeftString = ""
res@gsnRightString = ""
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 = "CenterLeft"
map_res@tiMainPosition = "Left"
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 = "MediumRes"
map_res@mpOutlineDrawOrder = "PostDraw"
map_res@mpUSStateLineColor = "Black"
map_res@mpPerimLineColor = "Black"
map_res@mpNationalLineColor = "Black"
map_res@mpLimbLineColor = "Black"
map_res@mpGridLineColor = "Black"
map_res@mpGeophysicalLineColor = "Black"
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's mentioned above in the script. What am I doing wrong here?</div><div><br></div><div>Thanks!</div><span class="m_8792142955542432848m_-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></div></div><span class="">______________________________<wbr>_________________<br>
ncl-talk mailing list<br>
<a href="mailto:ncl-talk@ucar.edu" target="_blank">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/mailma<wbr>n/listinfo/ncl-talk</a><br>
<br></span></blockquote></div><br></div>
</blockquote></div><br></div>