Re: add topography in the lon-height cross section plot

Hi Erika,
Here is a short script that masks a 3D (lev,lat,lon) array by using a
surface pressure array (lat,lon) and the where function:

; untested: Make sure the following code does what you think it is doing.
; uwnd = 3D array dimensioned lev x lat x lon, lev is in units of hPa
; ps = 2D array dimensioned lat x lon

lev = ywnd&lev
do gg = 0,dimsizes(lev)-1
      print("Now working on "+lev(gg)+"hPa")
      uwnd(gg,:,:) = where(ps.ge.lev(gg),uwnd@_FillValue, uwnd(gg,:,:))
end do

Adam







On Tue, Jul 15, 2014 at 9:29 AM, Erika Folova <e.folova@gmail.com> wrote:

> Thanks Adam and Tabish,
>
> Unfortunately I am not using WRF output, so i cant used this function
>
> http://www2.mmm.ucar.edu/wrf/OnLineTutorial/Graphics/NCL/Examples/CROSS_SECTION/wrf_CrossSection1.htm
>
> Adam, would you please demonstrate it with a small script how to add the
> topography/ topographic line in the hovmueller diagram using pressure
> data?? I still have no idea how do execute it. Any easy way for instance to
> mask with the topography data?
>
>
> On Tue, Jul 15, 2014 at 5:20 PM, Adam Phillips <asphilli@ucar.edu> wrote:
>
>> Hi Erika,
>> You could use gsn_add_polygon to add a filled polygon at the bottom of
>> your Hovmoller plot. The question would be how to get the topographic line
>> that you will fill beneath. (Use the surface pressure field to create one?)
>>
>> Another graphical alternative is to again use a surface pressure field
>> and mask your flux data when it falls beneath the surface pressure.
>>
>> Hope that all helps! If not, or if you have further questions, please
>> respond to the ncl-talk email list only.
>> Adam
>>
>>
>> On Sun, Jul 13, 2014 at 1:55 AM, Erika Folova <e.folova@gmail.com> wrote:
>>
>>> Dear NCL,
>>>
>>> Anyone can suggest me how can I add the topography in the lon-height
>>> cross section plot as i attached in
>>> this email. Is there any trick to do that since the data that i used
>>> have a complete data till 1000 hPa.
>>> Thank you
>>>
>>> Erika,
>>>
>>> _______________________________________________
>>> ncl-talk mailing list
>>> List instructions, subscriber options, unsubscribe:
>>> http://mailman.ucar.edu/mailman/listinfo/ncl-talk
>>>
>>>
>>
>>
>> --
>> Adam Phillips
>> Associate Scientist, Climate and Global Dynamics Division, NCAR
>> www.cgd.ucar.edu/staff/asphilli/ 303-497-1726
>>
>> <http://www.cgd.ucar.edu/staff/asphilli>
>>
>
>


--
Adam Phillips
Associate Scientist, Climate and Global Dynamics Division, NCAR
www.cgd.ucar.edu/staff/asphilli/ 303-497-1726

<http://www.cgd.ucar.edu/staff/asphilli>

Received on Tue Jul 15 06:11:57 2014