[ncl-talk] gc_latlon

Dennis Shea shea at ucar.edu
Mon Mar 18 09:16:01 MDT 2019


As noted by Soma:  a nice function

*https://www.ncl.ucar.edu/Document/Functions/Contributed/geolocation_circle.shtml*
<https://www.ncl.ucar.edu/Document/Functions/Contributed/geolocation_circle.shtml>

*Available in version 6.5.0
<https://www.ncl.ucar.edu/prev_releases.shtml#6.5.0> and later.*


*https://www.ncl.ucar.edu/Applications/polyg.shtml
<https://www.ncl.ucar.edu/Applications/polyg.shtml>*

*Examples 24 & 25*

*Y*ou can specify 500km as your 'circle radius'.
Use in conjunction with* gc_inout
<https://www.ncl.ucar.edu/Document/Functions/Built-in/gc_inout.shtml> t*o
mask data inside a great circle.
========

The following was created prior to the introduction of *geolocation_circle*.
It offers an alternative approach:


*https://www.ncl.ucar.edu/Applications/mask.shtml
<https://www.ncl.ucar.edu/Applications/mask.shtml>*

*Example:  Katrina_circle.ncl*

On Mon, Mar 18, 2019 at 1:49 AM Md. Jalal Uddin <dmjalal90 at gmail.com> wrote:

> I want to show the rainfall only within 500 km. Please look at the picture
> (within the circle).
>
> [image: jan.png]
>
> On Mon, Mar 18, 2019 at 3:35 PM Soma Roy <somaroy892 at gmail.com> wrote:
>
>> Try to use geolocation circle and based on that cut out the inner values
>> by using mask option as available in ncl pages.
>>
>> Hope this idea help you...
>>
>> Cheers,
>> Soma
>>
>> On Mon, Mar 18, 2019, 12:59 Md. Jalal Uddin <dmjalal90 at gmail.com wrote:
>>
>>> Hi all,
>>> I wanted to extract rainfall within 500 km from certain lat and lon. I
>>> am not getting the right one. any suggestions?
>>>
>>> f  = addfile("jan500.nc", "r")
>>>  prc =f->pre(:,{75:100},{0:30})
>>>  lon = f->longitude
>>>  lat =f->latitude
>>>  prc_reorder = prc(time|:,latitude|:,longitude|:)
>>>
>>>  prc_reorder&latitude = lat
>>>  prc_reorder&longitude = lon
>>>
>>> lat1 =(/5.5/)
>>> lon1 = (/87.0/)
>>>
>>> do i = 0,dimsizes(lat)-1
>>>       do j = 0,dimsizes(lon)-1
>>>       if gc_latlon(lat1,lon1,i,j,2,4) .lt.500 then
>>>           prc_reorder(:,i,j) =  prc_reorder(:,i,j)
>>>               end if
>>>       end do
>>>  end do
>>> ;;;;
>>> Jalal
>>> --
>>> *Md. Jalal Uddin*
>>> MSc in Applied Meteorology (English Language)
>>> Nanjing University of Information, Science and Technology, China
>>> B.Sc. in Disaster Management (Hons.)
>>> Patuakhali Science and Technology University, Bangladesh.
>>> Cell: +8613260859092, +8801792052662
>>> Web: www.dmjalal90.weebly.com
>>> Facebook: jalal.hossen.39
>>> LinkedIn: https://bd.linkedin.com/in/md-jalal-uddin-80a026b0
>>> Twitter: dmjalal90
>>> Skype: dmjalal90
>>> _______________________________________________
>>> ncl-talk mailing list
>>> ncl-talk at ucar.edu
>>> List instructions, subscriber options, unsubscribe:
>>> http://mailman.ucar.edu/mailman/listinfo/ncl-talk
>>>
>>
>
> --
> *Md. Jalal Uddin*
> MSc in Applied Meteorology (English Language)
> Nanjing University of Information, Science and Technology, China
> B.Sc. in Disaster Management (Hons.)
> Patuakhali Science and Technology University, Bangladesh.
> Cell: +8613260859092, +8801792052662
> Web: www.dmjalal90.weebly.com
> Facebook: jalal.hossen.39
> LinkedIn: https://bd.linkedin.com/in/md-jalal-uddin-80a026b0
> Twitter: dmjalal90
> Skype: dmjalal90
> _______________________________________________
> 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/20190318/6a3081a0/attachment-0001.html>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: jan.png
Type: image/png
Size: 307126 bytes
Desc: not available
URL: <http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20190318/6a3081a0/attachment-0001.png>


More information about the ncl-talk mailing list