Re: Nesting Domain

Hi Guan,

I just ran your script, and I see the problem. You are not explicitly
setting the contour levels to be drawn for each plot. Thus, NCL is
determining what contours to draw and what colors to shade for each of
the 3 plots. So your nested grids were all being drawn with different
contours/colors, and that was making it look like the interior domains
were had faded colors.

Try adding these two lines to your cnres resource list:
cnres_at_cnLevelSelectionMode = "ExplicitLevels"
cnres_at_cnLevels = ispan(1000,3400,200)

That will fix the problem..

My thanks to Mary for pointing me on the right track with debugging this.
Adam

gs6r wrote:
> Hey
>
> Thanks for your suggestions.
>
> Adam, I have followed your instruction to modify the code. Apparently,
> there is no map variable to draw, so I change it to draw(plot1) only.
>
> datatondc(plot1,xbox,ybox,x_out,y_out)
> gsn_polyline_ndc(wks,x_out,y_out,lnres)
> draw(map)
> frame(wks)
> end
>
> The results does not show complete graph at all, which shows that half
> of topographic map is out of the paper range. I do not understand.
> Anyway, it does show a graph with half topographic data.
>
> Mary, yes, as in the script, I only use the first domain ' terrain
> height and overlay the two rectangular box indicting the nesting domains.
>
> However, I have come out a new script with many modifications, the same
> thing occurs.
>
> The new code is attached.
>
> Thanks a lot,
> Guan
>
>
> --------------------------------------------------
> From: "Adam Phillips" <asphilli_at_cgd.ucar.edu>
> Sent: Wednesday, December 03, 2008 12:25 PM
> To: "Mary Haley" <haley_at_ucar.edu>
> Cc: "gs6r" <gs6r_at_virginia.edu>; <ncl-talk_at_ucar.edu>
> Subject: Re: Nesting Domain
>
>> Hi Guan,
>>
>> Looking at your plot, I believe you are pointing out that when you
>> draw a square using gsn_polyline_ndc (to denote one of your nested
>> grids), that the hgtd01 contour field appears lighter within the
>> square than outside the square.
>>
>> My first guess is that for whatever reason the overlay line:
>> map=wrf_map_overlays(find01,wks,plot1,cnres,mapres)
>> is causing the problem.
>>
>> To test, I would replace this:
>> map=wrf_map_overlays(find01,wks,plot1,cnres,mapres)
>> datatondc(map,xbox,ybox,x_out,y_out)
>> gsn_polyline_ndc(wks,x_out,y_out,lnres)
>> gsres = True
>> ........<snip>.....
>> draw(map)
>> frame(wks)
>> end
>>
>> with this:
>> datatondc(plot1,xbox,ybox,x_out,y_out)
>> gsn_polyline_ndc(wks,x_out,y_out,lnres)
>> draw(map)
>> frame(wks)
>> end
>>
>> Let us know if that changes the plot.
>> Adam
>>
>> Mary Haley wrote:
>>> Guan,
>>>
>>> I'm a little confused by this plot. I thought there should be two
>>> separate contour fields, but in your script, I think I only see one,
>>> and that's the one for hgtd01.
>>>
>>> My suggestion was going to be to make sure that both of your contour
>>> plots are represented by the same contour levels, but since I don't
>>> see a second plot, I'm not sure if this is the problem.
>>>
>>> --Mary
>>>
>>> On Dec 3, 2008, at 9:23 AM, gs6r wrote:
>>>
>>>> <Slideshow.gif> Play slideshow
>>>> <http://g.msn.com/5meen_us/171?path=/photomail/{58d07d8b-6741-43b3-86aa-824cafd8c3e0}&image=CEF991426176A654!896&imagehi=CEF991426176A654!894&CID=-3532632718406408620>
>>>> <SaveImages.gif> Save all photos Want to save all these photos at
>>>> once? Learn how
>>>> <http://help.live.com/help.aspx?mkt=en-us&project=WL_Mailv2&querytype=keyword&query=segami_lla_evas>
>>>>
>>>>
>>>> Online pictures are available for 30 days
>>>>
>>>> *Hey *
>>>> ** *I am trying to make a topographic map with my parent domain and
>>>> two rectangular boxes indicating two nesting domains. *
>>>> ** *However, one I overlay the two boxes, the color of inner domains
>>>> becomes lighter, which shows the discontinuous of the topographic
>>>> data. That does not make any senses. So I am wondering how can I fix
>>>> the problem.*
>>>> ** *I attached the graph and the code. Could you please help me?*
>>>> *Thanks a lot,*
>>>> *Guan*
>>>> <189.jpg>
>>>> <http://jmboig.blu.livefilestore.com/y1pIerLshmLfAy7w1YpbT_u0n70zvSuUU2NOXqOXwzWUdXOYft6mg2BFqVjaAP8Z6EVqZK_lsHuZE4/topo.jpg?download>
>>>>
>>>>
>>>> *Best Regards,*
>>>> *Guan Song
>>>> Graduate Student
>>>> Reply email: gs6r_at_virginia.edu <mailto:gs6r_at_virginia.edu>
>>>> New Clark 272
>>>> Dept. of Environmental Sciences
>>>> University of Virginia*
>>>> Get Windows Live Mail to create your own photo e-mails
>>>> <http://g.msn.com/5meen_us/175>
>>>>
>>>> <Slideshow.gif><SaveImages.gif><189.jpg><domain.ncl>_______________________________________________
>>>>
>>>> ncl-talk mailing list
>>>> List instructions, subscriber options, unsubscribe:
>>>> http://mailman.ucar.edu/mailman/listinfo/ncl-talk
>>>
>>>
>>> ------------------------------------------------------------------------
>>>
>>> _______________________________________________
>>> ncl-talk mailing list
>>> List instructions, subscriber options, unsubscribe:
>>> http://mailman.ucar.edu/mailman/listinfo/ncl-talk
>>
>> --
>> --------------------------------------------------------------
>> Adam Phillips asphilli_at_ucar.edu
>> National Center for Atmospheric Research tel: (303) 497-1726
>> ESSL/CGD/CAS fax: (303) 497-1333
>> P.O. Box 3000 Boulder, CO 80307-3000 http://www.cgd.ucar.edu/cas/asphilli
> Best Regards,
>
> Guan Song
> Graduate Student
> Reply email: gs6r_at_virginia.edu
> New Clark 272
> Dept. of Environmental Sciences
> University of Virginia
>
>
> ------------------------------------------------------------------------
>
> _______________________________________________
> ncl-talk mailing list
> List instructions, subscriber options, unsubscribe:
> http://mailman.ucar.edu/mailman/listinfo/ncl-talk

-- 
--------------------------------------------------------------
Adam Phillips			             asphilli_at_ucar.edu
National Center for Atmospheric Research   tel: (303) 497-1726
ESSL/CGD/CAS                               fax: (303) 497-1333
P.O. Box 3000				
Boulder, CO 80307-3000	  http://www.cgd.ucar.edu/cas/asphilli