Re: Nesting Domain

Thanks a lot, Adam and Mary ! It works perfectly now.

I just have a quick question that I expect have State boundary on the graph.

As shown on the graph, you can see the boundary of the salt lake in Utah.
However, the upper right corner, there should have boundary between Utah and
Wyoming but it does not show up. I am not sure why. Could you please help to
answer?

Thank you!
Guan

--------------------------------------------------
From: "Adam Phillips" <asphilli_at_cgd.ucar.edu>
Sent: Wednesday, December 03, 2008 7:29 PM
To: "gs6r" <gs6r_at_virginia.edu>
Cc: "Mary Haley" <haley_at_ucar.edu>; <ncl-talk_at_ucar.edu>
Subject: 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
Best Regards,

Guan Song
Graduate Student
Reply email: gs6r_at_virginia.edu
New Clark 272
Dept. of Environmental Sciences
University of Virginia
Received on Sat Jun 07 2008 - 17:21:32 MDT