[ncl-talk] geolocation_circle and gc_inout
Md. Jalal Uddin
dmjalal90 at gmail.com
Tue Apr 9 02:31:34 MDT 2019
Thank you, sir. I have solved the problem from your previous code (April 6).
[image: polyg_mar.png]
On Tue, Apr 9, 2019 at 12:13 PM Dennis Shea <shea at ucar.edu> wrote:
>
> [1] The variable looks like:
>
> double pre(time, longitude, latitude) ;
> pre:long_name = "precipitation" ;
> pre:units = "mm" ;
> pre:*_FillValue = NaN *; *<=== NaN not recognized by NC*L
>
>
> [2] The precipitation *does* have values of 0.0. In fact: **MANY** of
> them
>
> *%>* *ncl_filedump* all_mar.nc -v pre | less
>
> This is **your** responsibility to look-at (explore) the variable.
>
> [3] **You** set values 0.0 to _FillValue
>
> nonzerototmar = *where*(totmar.eq.0.0, prc2_reorder at _FillValue,
> totmar)
>
> By default: _FillValue are set to background (no color)
>
> [4] The polyg geolocation scripts [29-31] examples had the plotting in the
> following loop.
>
> do i =0,nLoc-1
> do j = 0,nRad-1
>
> [SNIP]
>
> --
> You had no plotting whatsoever.
>
> Attached is a subset of what I *think* you intended.
> I'm done!
>
>
>
>
> On Fri, Apr 5, 2019 at 6:49 PM Md. Jalal Uddin <dmjalal90 at gmail.com>
> wrote:
>
>> Sir,
>>
>> As I don't have 0.0 value within 500 km. What should I do to get
>> the precipitation value within circles?
>>
>> It would be a great help to me.
>>
>> Best Regards,
>> Jalal
>>
>>
>> On Tue, Apr 2, 2019 at 6:59 AM Md. Jalal Uddin <dmjalal90 at gmail.com>
>> wrote:
>>
>>> Sorry, sir. I don't have the 0.0 value within 500 km. Please look at the
>>> attached figures.
>>>
>>> On Mon, Apr 1, 2019 at 10:14 PM Dennis Shea <shea at ucar.edu> wrote:
>>>
>>>> [1] Save time: change *N=180* to (say) *N=90*
>>>>
>>>> *gc_inout* <http://www.ncl.ucar.edu/Document/Functions/Built-in/gc_inout.shtml> while written in "C" , it is not particularly 'fast'. It is best to reduce the number of points to check.
>>>>
>>>> [2] * *YOU*** should put timing checks into your code. This will tell you where the time is being spent. EG:
>>>>
>>>> beg_script = *get_cpu_time*()
>>>> :
>>>> :
>>>> beg_geo_circle = *get_cpu_time*()
>>>> circle = *geolocation_circle*(slat, slon, srad, sradu, N, optgeo) ; circle is type list
>>>> * print*("geolocation_circle: " + (*get_cpu_time*() - beg_geo_circle)+" seconds")
>>>> :
>>>> :
>>>> beg_gc_inout = *get_cpu_time*()
>>>> latlon_circle := *gc_inout*(lat2d,lon2d, poly_lat,poly_lon)
>>>> print("gc_inout: " + (*get_cpu_time*() - beg_gc_inout)+" seconds")
>>>> :
>>>> :
>>>>
>>>> : at end of script
>>>>
>>>> print("Script: Total time: " + (*get_cpu_time*() - beg_script)+" seconds")
>>>>
>>>> [3]
>>>> It seems you have 0.0 in your data as both a valid value and as a missing value [_FillValue]?
>>>> In your case, you [possibly] have 0.0 values *within* and* outside* 500km. Is that correct?
>>>> Currently, values outside the circle are set to False by gc_inout.
>>>> You have also, reset value 0.0 to _FillValue.
>>>>
>>>> If True, you will have to add code to draw a circle.
>>>>
>>>>
>>>>
>>>>
>>>>
>>>> On Fri, Mar 29, 2019 at 5:21 AM Md. Jalal Uddin <dmjalal90 at gmail.com>
>>>> wrote:
>>>>
>>>>> In addition, the do loop is too slow. How can I make it faster?
>>>>>
>>>>> On Fri, Mar 29, 2019 at 5:04 PM Md. Jalal Uddin <dmjalal90 at gmail.com>
>>>>> wrote:
>>>>>
>>>>>> Thank you, sir, for your kind help. I was able to mask data for
>>>>>> multiple cyclone's centres, but the figure did not show the whole circle
>>>>>> data (circle_mar.png). I attached the output figure and code in the
>>>>>> attachment. Could I have your kind suggestion?
>>>>>>
>>>>>>
>>>>>> On Fri, Mar 29, 2019 at 11:35 AM Dennis Shea <shea at ucar.edu> wrote:
>>>>>>
>>>>>>> Please see examples 29 and 30 at:
>>>>>>>
>>>>>>> *https://www.ncl.ucar.edu/Applications/polyg.shtml*
>>>>>>> <https://www.ncl.ucar.edu/Applications/polyg.shtml>
>>>>>>>
>>>>>>> On Thu, Mar 28, 2019 at 2:54 AM Md. Jalal Uddin <dmjalal90 at gmail.com>
>>>>>>> wrote:
>>>>>>>
>>>>>>>> Hi all,
>>>>>>>>
>>>>>>>> Could anyone give me the idea to mask rainfall data for multiple
>>>>>>>> circles? geolocation_circle produce ''circle_lat and lon'' that is 3
>>>>>>>> dimensions, but the data lat and lon 1 dimension.
>>>>>>>>
>>>>>>>> How can be used gc_inout to mask the data inside multiple circles?
>>>>>>>> Please look at the multiple circle's figures in the attached file.
>>>>>>>>
>>>>>>>> Any help would be appreciated.
>>>>>>>>
>>>>>>>> Variable: circle_lat
>>>>>>>> Type: float
>>>>>>>> Total Size: 10800 bytes
>>>>>>>> 2700 values
>>>>>>>> Number of Dimensions: 3
>>>>>>>> Dimensions and sizes: [location | 15] x [radii | 1] x [circle |
>>>>>>>> 180]
>>>>>>>> Coordinates:
>>>>>>>> Number Of Attributes: 7
>>>>>>>> _FillValue : 9.96921e+36
>>>>>>>> radii_input_units : km
>>>>>>>> radii_input_values : 500
>>>>>>>> radii_converted_units : great_circle_degrees
>>>>>>>> radii_converted_values : 4.522022
>>>>>>>> long_name : circle latitudes
>>>>>>>> units : degrees_north
>>>>>>>>
>>>>>>>> Variable: circle_lon
>>>>>>>> Type: float
>>>>>>>> Total Size: 10800 bytes
>>>>>>>> 2700 values
>>>>>>>> Number of Dimensions: 3
>>>>>>>> Dimensions and sizes: [location | 15] x [radii | 1] x [circle |
>>>>>>>> 180]
>>>>>>>> Coordinates:
>>>>>>>> Number Of Attributes: 7
>>>>>>>> _FillValue : 9.96921e+36
>>>>>>>> radii_converted_values : 4.522022
>>>>>>>> radii_converted_units : great_circle_degrees
>>>>>>>> radii_input_values : 500
>>>>>>>> radii_input_units : km
>>>>>>>> long_name : circle longitudes
>>>>>>>> units : degrees_east
>>>>>>>>
>>>>>>>> --
>>>>>>>> *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
>>>>>>
>>>>>
>>>>>
>>>>> --
>>>>> *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
>>>>>
>>>>
>>>
>>> --
>>> *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
>>>
>>
>>
>> --
>> *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
>>
>
--
*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/20190409/f3f25b90/attachment.html>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: polyg_mar.png
Type: image/png
Size: 111859 bytes
Desc: not available
URL: <http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20190409/f3f25b90/attachment.png>
More information about the ncl-talk
mailing list