<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>----------------------------------------------------------------------------------------------------------</div><div><pre style="color:rgb(0,0,0);word-wrap:break-word;white-space:pre-wrap">;----------------------------------------------------------------------
; wrf_gsn_6.ncl
;----------------------------------------------------------------------
; 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
;----------------------------------------------------------------------
; This script is meant to show the difference between plotting WRF
; data using wrf_xxx scripts, and using gsn_csm_xxx scripts.
;----------------------------------------------------------------------
; 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
;----------------------------------------------------------------------
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/shea_util.ncl"
begin
;---Parameters
input_file = "wrfout_d01_2015-02-05_00:00:<a href="http://00.nc">00.nc</a>"
input_dir = "/glade/u/home/mmasarik/WRF-NCL/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@mpUSStateLineThicknessF = 3.0
map_res@mpNationalLineThicknessF = 3.0
map_res@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)">----------------------------------------------------------------------------------------------------------</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><div><br></div><div>matt</div><div><br></div><div><br></div></div>