[ncl-talk] wrf domain location
Matt Masarik
mattmasarik at boisestate.edu
Tue Feb 7 16:41:12 MST 2017
ps, I'm working on Yellowstone using the following modules;
ncarenv/1.0
ncarbinlibs/1.1
intel/12.1.5
ncarcompilers/1.0
netcdf/4.3.0
ncl/6.3.0
On Tue, Feb 7, 2017 at 4:34 PM, Matt Masarik <mattmasarik at boisestate.edu>
wrote:
> 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).
>
> Here's the (editted) script:
> ------------------------------------------------------------
> ----------------------------------------------
>
> ;----------------------------------------------------------------------
> ; 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:00.nc"
> 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 at gsnFrame = False
> res at gsnDraw = False
> res at gsnLeftString = ""
> res at gsnRightString = ""
> res at tfDoNDCOverlay = True ; necesssary
>
>
> ;---Map plot
> map_res = res
> map_res = wrf_map_resources(a,map_res)
> map_res at gsnFrame = False
> map_res at gsnDraw = False
> map_res at gsnStringFontHeightF = 0.015
> map_res at tiMainJust = "CenterLeft"
> map_res at tiMainPosition = "Left"
> map_res at tiMainOffsetYF = -0.01
>
>
> ;---Set map resources based on projection on WRF output file
> map_res at mpFillOn = False
> map_res at mpOutlineOn = True
> map_res at mpDataBaseVersion = "MediumRes"
> map_res at mpOutlineDrawOrder = "PostDraw"
> map_res at mpUSStateLineColor = "Black"
> map_res at mpPerimLineColor = "Black"
> map_res at mpNationalLineColor = "Black"
> map_res at mpLimbLineColor = "Black"
> map_res at mpGridLineColor = "Black"
> map_res at mpGeophysicalLineColor = "Black"
> map_res at mpUSStateLineThicknessF = 3.0
> map_res at mpNationalLineThicknessF = 3.0
> map_res at mpGeophysicalLineThicknessF = 3.0
>
> map = gsn_csm_map(wks,map_res)
>
>
>
> ;---Overlay plots on map and draw.
> draw(map)
> frame(wks)
>
>
> end
>
> ----------------------------------------------------------------------------------------------------------
>
>
> The only warnings I get are the one's mentioned above in the script. What
> am I doing wrong here?
>
> Thanks!
>
> matt
>
>
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20170207/86556b20/attachment.html
More information about the ncl-talk
mailing list