# [ncl-talk] Trying to draw circles with different radius in Km around a station

Dennis Shea shea at ucar.edu
Thu Feb 8 07:45:59 MST 2018

```An alternative:

nggog can be used. Sample code segments follow. For small radii Npts=100 is
fine but for large radii and nice circles, Npts=360 might be better.

;=========================================================================
; On a sphere with radius 6371 km, one degree is approximately 110.57 km.
;---
; Of course, the earth is an oblate spheroid:
;      http://www.longitudestore.com/how-big-is-one-gps-degree.html
;=========================================================================
; https://www.ncl.ucar.edu/Document/Functions/Built-in/nggcog.shtml
;=========================================================================

;---One station

slat     = 0.0                  ; station lat
slon     = 0.0                  ;         lon
km2deg   = 1/110.57             ; conversion factor: deg ==> km
; one degree is approximately 110.57 km.

Npts = 100                      ; # of points at each radii (arbitrary)

print("")

print(Rlat(nr,:)+"   "+ Rlon(nr,:))
print("========")
end do

+++++++++++++++++++++++++++++++++++++++++++++++++++++
;---Multiple stations

slat     = (/ 0.0, 50.0     /)  ; station lat
slon     = (/ 0.0, 30.0     /)  ;         lon
km2deg   = 1/110.57                ; conversion factor: deg ==> km
; one degree is
approximately 110.57 km.

Nsta     = dimsizes(slat)            ; # stations
Npts     = 100                            ; # of points at each radii
(arbitrary; eg 360)

Rlon = new((/Nsta, Nrad, Npts/), "float")

do ns=0,Nsta-1
end do
end do

On Thu, Feb 8, 2018 at 2:38 AM, Karin Meier-Fleischer <
meier-fleischer at dkrz.de> wrote:

> Hi Soumik,
>
> you can use polymarkers to plot the circles around lat/lon locations. Have
> a look at example newcolor_4.ncl at http://ncl.ucar.edu/
> Applications/Scripts/newcolor_4.ncl
>
> I've attached an example script using random data for marker size and
> lat/lon locations overlayed on filled contour map.
>
> Hope this helps.
>
> -Karin
>
> Am 08.02.18 um 00:55 schrieb Soumik Basu:
>
> Hello NCL-users,
>
> I am trying to draw circles around a station with different radius. For
> example 50 Km, 100 Km, 150 Km etc.
>
> So, I am wondering if there is any way to do it in NCL.
>
> I looked at the "nggcog" function but from the description what I
> understand is that it can draw circles along a great circle only not
> according to a given distance.
>
> Thank you for any help.
>
> Thanks,
> Soumik
>
> --
> --
> “We’ve all got both light and dark inside us. What matters is the part we
> choose to act on. That’s who we really are.” – J.K. Rowling
>
> ************************************************************
> **************************************
> Dr. Soumik Basu
> Post Doctoral Research Faculty
> International Arctic Research Center, UAF, Fairbanks, AK, USA
> PhD in Atmospheric Sciences
> M.Sc. in Atmospheric Sciences
> Email: suvro05 at gmail.com
> website: http://soumikbasu.weebly.com/
>
> ************************************************************
> ***************************************
>
>
> _______________________________________________
> ncl-talk mailing listncl-talk at ucar.edu
> List instructions, subscriber options, unsubscribe:http://mailman.ucar.edu/mailman/listinfo/ncl-talk
>
>
>
> _______________________________________________
> 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/20180208/9fefbb99/attachment.html>
```