[ncl-talk] A simple plot question (sort of)

[Mary Haley]

Hi Barry,

If  plots are not displaying properly after they are overlaid on another,
usually the first thing I suspect is that the two plots are not in the same
data space.

But first, is this data rectilinear? The reason I ask is that I see you are
reading lat/lon off the file, but they are not used in the plotting.  If
you're not sure what I mean by rectilinear, please see the description at
the top of this page:

http://www.ncl.ucar.edu/Applications/plot_data_on_map.shtml

Then, scroll down to the first example, which explains how to tell if your
data is rectilinear or not.

Your filled contours over the map are roughly in the range 10 to 60 lat,
and -150 to -55 lon.  You didn't include the plot of the 700 mb contours,
but I would look at the X/Y axis of the plot of this data to see what range
it's in.  Also, if you look at your data with:

printVarSummary(pett_fgen)
printVarSummary(z_plane)

this should provide information about coordinate arrays and their ranges,
*if you have rectilinear data* This is another way to double-check that the
two datasets are in the same data space.

If you continue to have problems, can you email ncl-talk with the output
results from the two printVarSummary statements, and also include the plot
of z_plane?

Thanks,

--Mary



On Tue, Nov 22, 2016 at 10:33 AM, Barry Lynn <barry.h.lynn at gmail.com> wrote:

> Hello:
>
> I have moved away from using wrf interfaces for plotting and am simply
> trying to plot relative humidity and height at 700 mb from the NAM data
> (yesterday, Guido explained untangled my question about getting the
> variables I needed from the files themselves).
>
> If I comment out overlay and draw(contour_fgen, below, and uncomment
> draw(contour_z), I get a plot without any country/state line references,
> but I get perfect 700 mb contours.
>
> However, if I run the program as is, I get only a single height line
> superimposed on the relative humidity map (as shown).
>
> Clearly, something I have done is acting at cross purposes, but I various
> tweaks I've made (adding/removing code) have not solved the problem.
>
> This may be a simple fix, but it would be nice to know what it is.
>
> Thank you.
>
>   overlay(contour_fgen,contour_z)
>   draw(contour_fgen)   ; this will draw both the line and filled contours
> over a map
> ; draw(contour_z)   ; this will draw both the line and filled contours
> over a map
>
>
> --
> Barry H. Lynn, Ph.D
> Senior Lecturer,
> The Institute of the Earth Science,
> The Hebrew University of Jerusalem,
> Givat Ram, Jerusalem 91904, Israel
> Tel: 972 547 231 170
> Fax: (972)-25662581
>
> C.E.O, Weather It Is, LTD
> Weather and Climate Focus
> http://weather-it-is.com
> Jerusalem, Israel
> Local: 02 930 9525
> Cell: 054 7 231 170
> Int-IS: x972 2 930 9525
> US 914 432 3108
>
> _______________________________________________
> ncl-talk mailing list
> ncl-talk at ucar.edu
> List instructions, subscriber options, unsubscribe:
> http://mailman.ucar.edu/mailman/listinfo/ncl-talk
>
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20161128/f560093a/attachment.html