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

Fri Nov 25 04:41:25 MST 2016

Hi!

The resources

   gsnLeftStringParallelPosF
   gsnLeftStringOrthogonalPosF

   gsnRightStringParallelPosF
   gsnRightStringOrthogonalPosF

set the horizontal and vertical position for the "long_name" and  
"unit" respectively.

Setting them (Parallel~ 0.5 and Orthogonal~ -1.5) before plotting the  
Relative Humidity plot does the job.

Giannis Tegoulias

Quoting Barry Lynn <barry.h.lynn at gmail.com>:

> Hi:
>
> I haven't been able to move the title Relative Humidity to the label bar
> below. Nor, can I move it to the right.  Part of the problem is that I
> don't know how to reference them, so setting their position with title
> position information doesn't help.
>
> The name Relative Humidity, like the name Geopotential Height, comes from
> the attribute list.
>
> Usually, titles (data attributes) are automatically placed correctly on the
> page.
>
> I have consulted pages like:
>
> https://www.ncl.ucar.edu/Applications/title.shtml
>
> and I have read about variable attributes,
>
> https://www.ncl.ucar.edu/Document/Manuals/Ref_Manual/NclVariables.shtml
>
> *Moreover: the page that deals with Overlays shows an example where the
> variable names are space correctly, but never referenced (or their
> positions defined).*
>
> http://www.ncl.ucar.edu/Applications/Scripts/overlay_1.ncl
>
> I have made other searches using various key words, but this situation
> seems to be unusual (as resulting from the use of overlay).
>
> I realize that this may be a trivial thing, but it is not unimportant.
>
> Thank you for solving this problem.
>
> Barry Lynn
>
> On Thu, Nov 24, 2016 at 3:13 AM, David Brown <dbrown at ucar.edu> wrote:
>
>> 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
>>>
>>>
>>
>
>
> --
> 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