[ncl-talk] IBTrACS: storm tracking problem

[Adam Phillips]

Hi Jalal,
As you switched to using gsn_add_polyline/gsn_add_polymarker, my previous
advice will not work. The gs resources do not allow one to specify multiple
colors for multiple entries. Looking at your script, it looks to me like
you are only drawing one storm path, but your graphic shows multiple.
stormLat/stormLon look like single dimensioned arrays. (Maybe I am
misreading your code though.)

In any case, assuming you have multiple tracks, you need to loop over each
track, setting each track a different color along the way:
; ntracks is the number of tracks to be plotted

colors = (/"red","blue","green","yellow"/)   ; and so on and so on
d1 = new(ntracks,graphic)
d2 = d1
do gg = 0,ntracks-1   ; loop over each track
    resp at gsLineColor = colors(gg)
    amres at gsMarkerColor = colors(gg)
    d1(gg) = gsn_add_polyline(wks,map,stormLon,stormLat,resp)
    d2(gg) = gsn_add_polymarker(wks,map,stormLon,stormLat,amres)
end do
draw(map)
frame(wks)

For the colors, you can alternatively go old school and read in a colormap
directly and just refer to color indices:
gsn_define_colormap(wks,"srip_reanalysis") ; only has 19 colors
; (no need to set a colors array)
do gg = 0,ntracks-1   ; loop over each track
    resp at gsLineColor = gg+2
    amres at gsMarkerColor = gg+2

Anyway, the above should at the very least be able to get you started. Good
luck!
Adam


On Mon, Jan 25, 2021 at 6:52 PM Md. Jalal Uddin <dmjalal90 at gmail.com> wrote:

> Sorry, I forgot to add the updated script. Please see attached.
>
> On Tue, 26 Jan 2021 at 07:59, Adam Phillips <asphilli at ucar.edu> wrote:
>
>> Hi Jalal,
>> You did not send an updated version of your script, so I can only address
>> what you are doing by looking at the script attached to your initial
>> ncl-talk email. Assuming that you are drawing all the tracks with a single
>> gsn_csm_xy plot call, you will need to set the color of each marker using
>> the xyMarkerColors resource.
>>
>> http://www.ncl.ucar.edu/Document/Graphics/Resources/xy.shtml#xyMarkerColors
>> An example of doing this can be seen in XY example 37 here:
>> http://www.ncl.ucar.edu/Applications/xy.shtml#ex37
>> Hope that helps!
>> Adam
>>
>>
>> On Mon, Jan 18, 2021 at 6:52 AM Md. Jalal Uddin via ncl-talk <
>> ncl-talk at mailman.ucar.edu> wrote:
>>
>>> Hi all,
>>>
>>> I have solved the previous issue. However, I am still facing a problem.
>>> All typhoons are getting one line (attached figure), but every typhoon
>>> should have a different line.
>>>
>>> Any suggestions would be appreciated.
>>>
>>> On Sun, 17 Jan 2021 at 12:28, Md. Jalal Uddin <dmjalal90 at gmail.com>
>>> wrote:
>>>
>>>> Hi all,
>>>>
>>>> I am getting a map for only one typhoon instead of all (attached
>>>> figure). I followed the ibtracs_1.ncl example (
>>>> https://www.ncl.ucar.edu/Applications/Scripts/ibtracs_1.ncl). I need
>>>> to modify the following code.
>>>>
>>>> ID      = 2
>>>>  year  = 2019
>>>>
>>>> if( ID.ge.2 ) then
>>>>     stormInd = ind( ( year.eq.stormYear ).and.( ID.eq.stormID ) )
>>>>   else
>>>>     stormInd = ind( stormID.gt.ID )
>>>>   end if
>>>>
>>>>   stormLat = Latitude(stormInd)
>>>>   stormLon = Longitude(stormInd)
>>>>
>>>> The full code and data have been added in the attachment.
>>>>
>>>> Best regards,
>>>> Jalal
>>>> --
>>>> *Md. Jalal Uddin*
>>>> Founder and Director of Research Society (
>>>> https://researchsociety20.org/founder-and-director/)
>>>> PhD candidate - Atmospheric Physics and Atmospheric Environment (NUIST)
>>>> M.Sc. in Applied Meteorology (NUIST)
>>>> B.Sc. in Disaster Management (PSTU)
>>>> Email: dmjalal90 at nuist.edu.cn
>>>> Personal website: https://dmjalal90.weebly.com/
>>>> Research gate profile: https://www.researchgate.net/profile/Md_Uddin125
>>>>
>>>> YouTube channel: http://www.youtube.com/c/MdJalalUddindmjalal90
>>>> Cell: +8613260859092, +8801792052662
>>>>
>>>
>>>
>>> --
>>> *Md. Jalal Uddin*
>>> Founder and Director of Research Society (
>>> https://researchsociety20.org/founder-and-director/)
>>> PhD candidate - Atmospheric Physics and Atmospheric Environment (NUIST)
>>> M.Sc. in Applied Meteorology (NUIST)
>>> B.Sc. in Disaster Management (PSTU)
>>> Email: dmjalal90 at nuist.edu.cn
>>> Personal website: https://dmjalal90.weebly.com/
>>> Research gate profile: https://www.researchgate.net/profile/Md_Uddin125
>>>
>>> YouTube channel: http://www.youtube.com/c/MdJalalUddindmjalal90
>>> Cell: +8613260859092, +8801792052662
>>> _______________________________________________
>>> ncl-talk mailing list
>>> ncl-talk at mailman.ucar.edu
>>> List instructions, subscriber options, unsubscribe:
>>> https://mailman.ucar.edu/mailman/listinfo/ncl-talk
>>
>>
>>
>> --
>> Adam Phillips
>> Associate Scientist,  Climate and Global Dynamics Laboratory, NCAR
>> www.cgd.ucar.edu/staff/asphilli/   303-497-1726
>>
>> <http://www.cgd.ucar.edu/staff/asphilli>
>>
>
>
> --
> *Md. Jalal Uddin*
> Founder and Director of Research Society (
> https://researchsociety20.org/founder-and-director/)
> PhD candidate - Atmospheric Physics and Atmospheric Environment (NUIST)
> M.Sc. in Applied Meteorology (NUIST)
> B.Sc. in Disaster Management (PSTU)
> Email: dmjalal90 at nuist.edu.cn
> Personal website: https://dmjalal90.weebly.com/
> Research gate profile: https://www.researchgate.net/profile/Md_Uddin125
> YouTube channel: http://www.youtube.com/c/MdJalalUddindmjalal90
> Cell: +8613260859092, +8801792052662
>


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

<http://www.cgd.ucar.edu/staff/asphilli>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://mailman.ucar.edu/pipermail/ncl-talk/attachments/20210127/0ac66aeb/attachment.html>