[ncl-talk] Radar data plotting, no MAP_PROJ attribute
Ryan Connelly
rconne01 at gmail.com
Wed Jul 20 15:55:55 MDT 2016
Got it. Thanks!
On Wed, Jul 20, 2016 at 5:50 PM, David Brown <dbrown at ucar.edu> wrote:
> Set mpOutlineBoundarySets to "USStates" or "GeophysicalAndUSStates".
> If you want counties use "AllBoundaries"
> -dave
>
> On Wed, Jul 20, 2016 at 3:23 PM, Ryan Connelly <rconne01 at gmail.com> wrote:
> > Alex,
> >
> > Thanks! That and a few other changes got it working. My last problem
> now
> > is that I can't get state outlines to plot, even though I think I have
> the
> > necessary resources set:
> >
> > radres at mpUSStateLineColor = "Black"
> > radres at mpUSStateLineThicknessF = 2.0
> >
> > I get coastal outlines, but nothing else.
> >
> > Thanks,
> > Ryan
> >
> > On Wed, Jul 20, 2016 at 3:44 PM, Alexander Schaefer <aschae11 at uncc.edu>
> > wrote:
> >>
> >> Ryan,
> >>
> >> You might try looking here:
> >>
> >> http://www.ncl.ucar.edu/Document/Graphics/overlays_on_map.shtml
> >>
> >> Check out point 2 talking about lat2d and lon2d.
> >>
> >> It might be as simple a fix as
> >>
> >> rd_plane at lat2d = rlat
> >> rd_plane at lon2d = rlon
> >>
> >> It should associate lat/lon point wise like you are looking for.
> >>
> >> Cheers,
> >> -Alex
> >>
> >>
> >>
> >>
> >> On Jul 20, 2016, at 3:31 PM, Ryan Connelly <rconne01 at gmail.com> wrote:
> >>
> >> Hi Mary,
> >>
> >> Thanks for getting back to me. Sorry about the delay on my end.
> >>
> >> printVarSummary of the rlat (or rlon) array shows
> >>
> >> Variable: rlat
> >> Type: float
> >> Total Size: 131044 bytes
> >> 32761 values
> >> Number of Dimensions: 2
> >> Dimensions and sizes: [y0 | 181] x [x0 | 181]
> >> Coordinates:
> >> y0: [-120..119.94]
> >> x0: [-120..119.94]
> >> Number Of Attributes: 3
> >> _FillValue : 41.94297
> >> standard_name : latitude
> >> units : degrees_north
> >>
> >>
> >> and for rd_plane shows
> >>
> >>
> >> Variable: rd_plane
> >> Type: float
> >> Total Size: 131044 bytes
> >> 32761 values
> >> Number of Dimensions: 2
> >> Dimensions and sizes: [y0 | 181] x [x0 | 181]
> >> Coordinates:
> >> y0: [-120..119.94]
> >> x0: [-120..119.94]
> >> Number Of Attributes: 12
> >> z0 : 2
> >> time : 1424628563
> >> grid_mapping : grid_mapping_0
> >> coordinates : lon0 lat0
> >> units : dBZ
> >> long_name : radar_reflectivity
> >> standard_name : REF
> >> _FillValue : 9.96921e+36
> >> valid_max : 25.3757
> >> valid_min : -9.187451
> >>
> >> printMinMax shows:
> >>
> >> (0) latitude: min=39.779 max=41.943
> >> (0) longitude: min=-74.313 max=-71.4156
> >>
> >> The values of rlat and rlon are the correct latitude and longitude
> values,
> >> but the trouble seems to be in matching the i,j'th point in rd_plane
> with
> >> the lat,lon contained in the i'th point of rlon and j'th point of rlat,
> so
> >> that the code can plot radar reflectivity in the correct place.
> >>
> >> Thanks,
> >> Ryan
> >>
> >> On Sun, Jul 10, 2016 at 11:53 PM, Mary Haley <haley at ucar.edu> wrote:
> >>>
> >>> Hi Ryan,
> >>>
> >>> You can't use wrf_xxxx functions to plot anything but WRF-ARW data.
> >>> These functions assume that you are using WRF files with names like
> >>> "wrfout_d01_2008-09-28_00:00:00" which have a bunch of global
> attributes in
> >>> them that define the native map projection of the WRF data. This is
> why you
> >>> were getting an error about MAP_PROJ being undefined. This is one of
> the
> >>> required global attributes.
> >>>
> >>> If you are trying to plot the data over a map, then you will need to
> >>> examine your file and your data variable (rd_plane) to determine if it
> is
> >>> also on some native map projection, or whether you need to provide
> lat/lon
> >>> coordinates for plotting. I then suggest using the gsn_csm_contour or
> >>> gsn_csm_contour_map functions to generate just contours, or contours
> over a
> >>> map.
> >>>
> >>> Please see this page for some examples on how to plot data over a map
> >>> using NCL:
> >>>
> >>> http://www.ncl.ucar.edu/Applications/plot_data_on_map.shtml
> >>>
> >>> You might also want to look at this page:
> >>>
> >>> http://www.ncl.ucar.edu/Applications/radar.shtml
> >>>
> >>>
> >>> --Mary
> >>>
> >>>
> >>> On Fri, Jul 8, 2016 at 4:35 PM, Ryan Connelly <rconne01 at gmail.com>
> wrote:
> >>>>
> >>>> Hi,
> >>>>
> >>>> I have Level II radar data in netcdf format (using radx and the steps
> >>>> outlined here:
> >>>>
> https://publish.illinois.edu/snesbitt/resources/radar-data-processing-with-radx/
> ).
> >>>>
> >>>> I want to overlay it on a base map, in the same way I would with WRF
> >>>> data. Here's a highly condensed version of what I have:
> >>>>
> >>>> rf = addfile(diri+fili,"r")
> >>>> rd = rf->REF(:,:,:,:) ; float REF(time, z0, y0,
> x0) ;
> >>>> RHO, ZDR
> >>>>
> >>>> rd_plane = rd(0,0,:,:)
> >>>>
> >>>> contour = gsn_csm_contour(wks,rd_plane,radres)
> >>>> plot = wrf_map_overlays(rf,wks,(/contour/),radres,mpres)
> >>>>
> >>>> I was just assuming the netcdf had location data packaged into it,
> but I
> >>>> get the error:
> >>>>
> >>>> (0) wrf_map: Error: no MAP_PROJ attribute in input file
> >>>> fatal:overlay: missing value as input, can't continue
> >>>>
> >>>> I tried adding in a wrf file I had laying around and using it for map
> >>>> plotting, and while that got be a basemap, by radar data was
> obviously not
> >>>> matched geographically (it showed the KOKX RDA over New York City!).
> >>>>
> >>>> Thanks,
> >>>> Ryan
> >>>>
> >>>> --
> >>>> Ryan Connelly
> >>>> M.S. Student in Atmospheric Sciences, Stony Brook University
> >>>> B.S. in Meteorology with Minors in Mathematics and GIS, Valparaiso
> >>>> University
> >>>> rconne01 at gmail.com
> >>>> ryan.connelly at stonybrook.edu
> >>>>
> >>>> _______________________________________________
> >>>> ncl-talk mailing list
> >>>> ncl-talk at ucar.edu
> >>>> List instructions, subscriber options, unsubscribe:
> >>>> http://mailman.ucar.edu/mailman/listinfo/ncl-talk
> >>>>
> >>>
> >>
> >>
> >>
> >> --
> >> Ryan Connelly
> >> M.S. Student in Atmospheric Sciences, Stony Brook University
> >> B.S. in Meteorology with Minors in Mathematics and GIS, Valparaiso
> >> University
> >> rconne01 at gmail.com
> >> ryan.connelly at stonybrook.edu
> >> _______________________________________________
> >> ncl-talk mailing list
> >> ncl-talk at ucar.edu
> >> List instructions, subscriber options, unsubscribe:
> >> http://mailman.ucar.edu/mailman/listinfo/ncl-talk
> >>
> >>
> >
> >
> >
> > --
> > Ryan Connelly
> > M.S. Student in Atmospheric Sciences, Stony Brook University
> > B.S. in Meteorology with Minors in Mathematics and GIS, Valparaiso
> > University
> > rconne01 at gmail.com
> > ryan.connelly at stonybrook.edu
> >
> > _______________________________________________
> > ncl-talk mailing list
> > ncl-talk at ucar.edu
> > List instructions, subscriber options, unsubscribe:
> > http://mailman.ucar.edu/mailman/listinfo/ncl-talk
> >
>
--
Ryan Connelly
M.S. Student in Atmospheric Sciences, Stony Brook University
B.S. in Meteorology with Minors in Mathematics and GIS, Valparaiso
University
rconne01 at gmail.com
ryan.connelly at stonybrook.edu
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20160720/c4d4d78a/attachment.html
More information about the ncl-talk
mailing list