[ncl-talk] r,theta (radar) plotting
Dennis Shea
shea at ucar.edu
Wed Jul 10 08:11:47 MDT 2019
re: How can I plot for multiple centers?
Use a do loop:
do nc_plot=0,nctr-1 ; loop over centers
:
end do
==============
re: How can I set x tick and y tick (-300,300) like radar_2 example?
radar_2 uses cartesian coordinates not lat/lon
A minor minor addition to place a grid over the middle figure.
*https://www.ncl.ucar.edu/Applications/regrid.shtml*
<https://www.ncl.ucar.edu/Applications/regrid.shtml>
Example 15
If you want to convert to cartesian, you will have to do that yourself.
On Tue, Jul 9, 2019 at 3:00 AM Md. Jalal Uddin <dmjalal90 at gmail.com> wrote:
> Thank you, sir. regrid_15 example works very well for one center (attached
> figure).
>
> How can I plot for multiple centers?
> How can I set x tick and y tick (-300,300) like radar_2 example? (
> https://www.ncl.ucar.edu/Applications/Scripts/radar_2.ncl)
>
> On Tue, Jul 9, 2019 at 5:42 AM Dennis Shea <shea at ucar.edu> wrote:
>
>> Not sure what you are asking: What does "my TRMM rainfall data" look like.
>> ===
>>
>> Generally: TRMM is on a 0.25 global rectilinear grid [400x1440] spanning
>> 49.875S to 49.875N
>>
>> ===
>>
>> If you have the latitudes and longitudes corresponding to some r,theta
>> [say: plat(:,:), plon(:,:) ] then use
>>
>> res at sfXArray = plon
>> res at sfYArray = plat
>>
>> and use appropriate contour function.
>> ====
>>
>> If you want to specify some central location [clat, clon] then maybe
>>
>> *http://www.ncl.ucar.edu/Applications/regrid.shtml*
>> <http://www.ncl.ucar.edu/Applications/regrid.shtml>
>> Example 15
>>
>> If you do not want the plot over a map: eliminate 'mp' related resources
>> and use (say) gsn_csm_contour instead of gsn_csm_contour_map. See the
>> leftmost and center figures.
>>
>>
>>
>>
>>
>>
>>
>>
>>
>>
>>
>> On Sun, Jul 7, 2019 at 11:53 PM Md. Jalal Uddin via ncl-talk <
>> ncl-talk at ucar.edu> wrote:
>>
>>> Hi all,
>>>
>>> How can I plot my TRMM rainfall data as a similar to radar_2 example (
>>> https://www.ncl.ucar.edu/Applications/Scripts/radar_2.ncl)?
>>>
>>> How can I modify here?
>>>
>>> angles = f->azimuth
>>> angles(0:63) = angles(0:63)-360 ; fix to make angles monotonic
>>>
>>> DEGTORAD = 0.017453292519943
>>> xcenter = 0.0
>>> ycenter = 0.0
>>> radius = 960 * 0.25 ; this is radius in kilometers
>>>
>>> ; Create 2D coordinate arrays.
>>> inc = radius / (dsizes(1) - 1)
>>> x = inc * ispan(0,dsizes(1)-1,1)
>>> angles2d = conform(dbz,angles,0)
>>> x2d = conform(dbz,x,1)
>>> xarr = xcenter + x2d * cos(DEGTORAD * angles2d)
>>> yarr = ycenter + x2d * sin(DEGTORAD * angles2d)
>>>
>>>
>>> Best wishes,
>>> 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
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20190710/7c7502cb/attachment.html>
More information about the ncl-talk
mailing list