[ncl-talk] Can't overlay maps on a NAM (native grid)

[David Brown]

I am attaching a modified script that works on your data (except you need
to figure out what you want to do about the titles). Most of your script
was okay except for a few crucial points.

Main points to remember:
All the resource warning messages come from setting map resources for a
plot that does not include a map -- the warnings are not fatal, these
resources will simply be ignored. So don't copy the resources from the
contour map to the plain contour. Instead start over with a new resource
variable. You had
zres = res     ; copies all the attributes of "res" to "zres"
Instead
zres = True

Since you are doing a native plot (tfDoNDCOverlay = True) you also need to
set this for the overlay plot (contour_z):

zres at tfDoNDCOverlay = True

Finally whenever you call the "draw" routine, you also need to call
"frame", so:

draw(contour_rh)
frame(wks)

Hope this helps.
 -dave




On Wed, Nov 23, 2016 at 12:58 PM, Barry Lynn <barry.h.lynn at gmail.com> wrote:

> Dave: Please see below.
>
> But first: Here are the three components needed to test the program.
>
> 1) The file, 2) the ncl program, and 3 the get_ij.f, which has to be
> wrapped.
>
> Dave: You wrote:  "The problem is that both your individual plots are
> drawn using a map projection. You should draw the contour_z field without
> the map
> projection using "gsn_csm_contour". Then you will be able to do the
> overlay without error."
>
> If I remove the "map" from the second gsn_csm_contour command, I just get
> those warnings again and no overlay.
>
> Thank you,​
>  namanl_218_20130208_1800_000.grb
> <https://drive.google.com/file/d/0B3E3GjUeKpBdQWhSVFpISWM2SkE/view?usp=drive_web>
> ​
>
> On Wed, Nov 23, 2016 at 9:26 PM, Barry Lynn <barry.h.lynn at gmail.com>
> wrote:
>
>> Hello:
>>
>> I am trying to overlay the height field from the NAM on the relative
>> humidity.
>>
>> I succeeded to zoom and project the NAM fields correctly, and I succeeded
>> to make two separate plots, both of which were projected correctly.
>>
>> I consulted these two pages.
>>
>> http://www.ncl.ucar.edu/Applications/Scripts/lcnative_2.ncl
>>
>> and:
>>
>> http://www.ncl.ucar.edu/Applications/Scripts/conOncon_2.ncl
>>
>> I can create test.00001, and test.00002 using the attached ncl code.
>>
>> However, I cannot create the overlay.
>>
>> I try to do this using the last two lines of the code:
>>
>>   overlay(contour_rh,contour_z)
>>   draw(contour_rh)
>>
>>
>> But, I get this error.
>>
>> fatal:NhlAddOverlay: plot class mapPlotClass cannot be overlay plot member
>>
>> If I try to change gsn_csm_contour_map to gsn_csm_contour to then invoke
>> the overlay procedure, I get:
>>
>>
>>
>> warning:mpProjection is not a valid resource in test_contour at this time
>> warning:mpLimitMode is not a valid resource in test_contour at this time
>> warning:mpLeftCornerLonF is not a valid resource in test_contour at this
>> time
>> warning:mpRightCornerLonF is not a valid resource in test_contour at this
>> time
>> warning:mpLeftCornerLatF is not a valid resource in test_contour at this
>> time
>> warning:mpRightCornerLatF is not a valid resource in test_contour at this
>> time
>> warning:mpLambertMeridianF is not a valid resource in test_contour at
>> this time
>> warning:mpLambertParallel1F is not a valid resource in test_contour at
>> this time
>> warning:mpLambertParallel2F is not a valid resource in test_contour at
>> this time
>> warning:mpOutlineDrawOrder is not a valid resource in test_contour at
>> this time
>> warning:mpFillOn is not a valid resource in test_contour at this time
>> warning:mpOutlineBoundarySets is not a valid resource in test_contour at
>> this time
>> fatal:NhlAddOverlay: plot class mapPlotClass cannot be overlay plot member
>>
>> and the projection is lost.
>>
>> Your suggestions are appreciated.
>>
>> Barry Lynn
>>
>>
>> --
>> 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
>>
>
>
>
> --
> 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/20161123/0c6ad40a/attachment.html 
-------------- next part --------------
A non-text attachment was scrubbed...
Name: test.000003.png
Type: image/png
Size: 249338 bytes
Desc: not available
Url : http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20161123/0c6ad40a/attachment.png 
-------------- next part --------------
A non-text attachment was scrubbed...
Name: plot_native_mod.ncl
Type: application/octet-stream
Size: 6705 bytes
Desc: not available
Url : http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20161123/0c6ad40a/attachment.obj