[ncl-talk] How to get the result for each grid?

[Atul Saini]

Ok, I will check and let you know.


On Mon, Sep 9, 2019, 17:04 Rashed Mahmood <rashidcomsis at gmail.com> wrote:

>
> Are you sure about the data file? Why got no values for x min and max?
> Also why was the min and max same value for X, 5050?
>
> There is no problem with function, I am sure. It's the data file that need
> to be checked thoroughly? If you cannot find issue with the file you would
> need provide the data file to the list.
>
>
>
>
> On Mon, Sep 9, 2019 at 1:26 PM Atul Saini <atulsainimail at gmail.com> wrote:
>
>> I got the following answer-
>>
>>
>>
>>
>>
>>
>>
>>
>>
>>
>>
>>
>>
>>
>>
>>
>>
>>
>>
>>
>>
>>
>>
>>
>>
>>
>>
>>
>>
>>
>>
>> *ncl 6> printVarSummary(x)Variable: xType: floatTotal Size: 6966000
>> bytes            1741500 valuesNumber of Dimensions: 3Dimensions and sizes:
>> [time | 100] x [lat | 129] x [lon | 135]Coordinates:             time: [
>> 1.. 100]            lat: [6.5..38.5]            lon: [66.5..100]Number Of
>> Attributes: 3  long_name : GRIDDED RAINFALL  _FillValue : -999
>> missing_value : -999ncl 7> printVarSummary(X)Variable: XType: floatTotal
>> Size: 69660 bytes            17415 valuesNumber of Dimensions: 2Dimensions
>> and sizes: [129] x [135]Coordinates: Number Of Attributes: 1  _FillValue :
>> -999ncl 8> ncl 9> print(x at _FillValue)(0) -999ncl 10> print(X at _FillValue)(0)
>> -999*
>> Regards,
>>
>>
>>
>>
>>
>> On Mon, Sep 9, 2019 at 4:46 PM Rashed Mahmood <rashidcomsis at gmail.com>
>> wrote:
>>
>>> Why there are no values for min and max of x variable?
>>> I suspect that there is something wrong with your data file, however,
>>> you should do some more checking... e.g.
>>>
>>> printVarSummary(x)
>>> printVarSummary(X)
>>>
>>> print(x at _FillValue)
>>> print(X at _FillValue)
>>>
>>>
>>>
>>>
>>> On Mon, Sep 9, 2019 at 1:11 PM Atul Saini <atulsainimail at gmail.com>
>>> wrote:
>>>
>>>> Rashed sir,
>>>>                   It gives the following answer-
>>>>
>>>>
>>>>
>>>>
>>>>
>>>>
>>>> *ncl 4> printMinMax(x,True)(0) (0) GRIDDED RAINFALLncl 5>
>>>> printMinMax(X,True)(0) (0) min=5050   max=5050*
>>>>
>>>> Regards,
>>>>
>>>>
>>>>
>>>>
>>>>
>>>> On Mon, Sep 9, 2019 at 4:38 PM Rashed Mahmood <rashidcomsis at gmail.com>
>>>> wrote:
>>>>
>>>>> I do not see any problem with the lines you just posted.
>>>>> Can you add following in the above code you wrote and let me know:
>>>>>
>>>>> printMinMax(x,True)
>>>>> printMinMax(X,True)
>>>>>
>>>>>
>>>>> On Mon, Sep 9, 2019 at 12:18 PM Atul Saini <atulsainimail at gmail.com>
>>>>> wrote:
>>>>>
>>>>>> Hello Rashed Sir,
>>>>>>                            I used dim_sum_n but it could not yield
>>>>>> the required result at each grid point. I want to get the sum of the
>>>>>> variable through time for all the grid points and therefore I used -
>>>>>>
>>>>>>
>>>>>>
>>>>>>
>>>>>> *a=addfile("RFIM_1901_1902_days_rf_tym.nc","r")x=a->tym(0:99,:,:)X=dim_sum_n(x,0)*
>>>>>> *print(X)*
>>>>>>
>>>>>> Printing X doesn't provide the correct result. Please guide. Your
>>>>>> reply is highly appreciated.
>>>>>>
>>>>>>
>>>>>>
>>>>>> Regards,
>>>>>>
>>>>>>
>>>>>>
>>>>>>
>>>>>>
>>>>>> On Mon, Sep 9, 2019 at 2:58 PM Rashed Mahmood <rashidcomsis at gmail.com>
>>>>>> wrote:
>>>>>>
>>>>>>> I think you want to use dim_sum_n_wrap() or dim_sum_n() functions
>>>>>>> instead of using "sum".  "sum" would result in total of everything for a
>>>>>>> variable. See the docs at the following link:
>>>>>>>
>>>>>>> https://www.ncl.ucar.edu/Document/Functions/Built-in/dim_sum_n.shtml
>>>>>>>
>>>>>>>
>>>>>>> On Mon, Sep 9, 2019 at 8:55 AM Atul Saini via ncl-talk <
>>>>>>> ncl-talk at ucar.edu> wrote:
>>>>>>>
>>>>>>>> Hi NCL Talk,
>>>>>>>>                     I am trying to calculate the value for each
>>>>>>>> grid using the function as explained in the script, but at the end the
>>>>>>>> result is one value only. I want the value at each grid point.
>>>>>>>> *Script is as follows and datafile is attached herewith.*
>>>>>>>> ;begin
>>>>>>>> a=addfile("Data_File.nc","r")
>>>>>>>> x=a->tym(0:4,:,:)
>>>>>>>> y=a->rf(0:4,:,:)
>>>>>>>> v=a->tym(5:9,:,:)
>>>>>>>> w=a->rf(5:9,:,:)
>>>>>>>> A1=
>>>>>>>> ((sum(y))*(sum(sqrt(x)))-(sum(x))*(sum(x*y)))/(5*(sum(sqrt(x)))-(sqrt(sum(x))))
>>>>>>>> B1= (5*(sum(x*y))-(sum(x)*sum(y))/(5*(sum(sqrt(x)))-(sqrt(sum(x)))))
>>>>>>>> A2=
>>>>>>>> ((sum(w))*(sum(sqrt(v)))-(sum(v))*(sum(v*w)))/(5*(sum(sqrt(v)))-(sqrt(sum(v))))
>>>>>>>> B2= (5*(sum(v*w))-(sum(v)*sum(w))/(5*(sum(sqrt(v)))-(sqrt(sum(v)))))
>>>>>>>> chpOnset = ((A1-A2)/(B2-B1))
>>>>>>>> print(chpOnset)
>>>>>>>> *Please, help me to find out the mistake I am making. Your reply is
>>>>>>>> very important to me. Thanks in advance.*
>>>>>>>> Regards,
>>>>>>>>
>>>>>>>>
>>>>>>>>
>>>>>>>> _______________________________________________
>>>>>>>> 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/20190909/bd8279e4/attachment.html>