[ncl-talk] loop_query

Mary Haley haley at ucar.edu
Fri Mar 31 13:58:38 MDT 2017 

Kunal,

The error in your code is this line:\

          data = a[:]->bcfire(0,:,:)

What this is doing is reading "bcfire" across all files, but then you are
subsetting it with (0,:,:), so you are effectively only using the first day
of data.

In the example I sent you before, I instructed you to use this:
  data = a[:]->bcfire

This will read "bcfire" across all files.  You can verify that it looks
correct by adding a printVarSummary immediately after the above statement:

  printVarSummary(data)


--Mary

On Fri, Mar 31, 2017 at 12:14 PM, Kunal Bali <kunal.bali9 at gmail.com> wrote:

> I apologize from the deep of my heart. I didn't mean to bother you.
>
> actually
> The script is reading all the 30 days files.
>
> Variable: data
> Type: float
> Total Size: 11880000 bytes
>             2970000 values
> Number of Dimensions: 3
> Dimensions and sizes:    [time | 30] x [latitude | 330] x [longitude | 300]
> Coordinates:
>             time: [913824..914520]
>             latitude: [37.95..5.05]
>             longitude: [68.05..97.95]
> Number Of Attributes: 3
>   long_name :    Wildfire flux of Black Carbon
>
>
> It showing an error
>
> fatal:conform: the dimensions sizes of the second argument do not match
> those indicated by the third argument
> fatal:["Execute.c":8578]:Execute: Error occurred at or near line 146 in
> file kunal.ncl
>
> The error line are
>   lat1d = ndtooned(conform(data,data&latitude,0))
>   lon1d = ndtooned(conform(data,data&longitude,1))
>
> I tried to find to find this answer but I could not find it. I am sorry if
> i am bothering you again.
>
>
>
>
>
>
>
> Kunal Bali
> Research Scholar
> Radio & Atmospheric Science Division
> CSIR - National Physical Laboratory
> New Delhi - 110012
> India
>
>
>
>
>
> On Fri, Mar 31, 2017 at 11:20 PM, Dennis Shea <shea at ucar.edu> wrote:
>
>> I still find your question and the code confusing.
>>
>> It takes time to look at code and determine what is happening. It is your
>> responsibility to pose clear questions and, if appropriate, a clean,
>> minimal script that illustrates the issue(s). We answer questions on a
>> volunteer basis. Our time is valuable too! You must help us help you!
>>
>> Please read the addfiles and adfiles documentation.
>>
>> http://www.ncl.ucar.edu/Document/Functions/Built-in/addfiles.shtml
>> http://www.ncl.ucar.edu/Document/Functions/Built-in/addfile.shtml
>>
>> With addfiles ... please look at examples 1 and 2.
>>
>> Note the difference between 'cat' (the default) and 'join' (alternate
>> approach)
>>
>>
>> ---
>> Create a script that contains *only* this part of your code.
>>
>>
>> ;---Read var1
>>       dir    = "/home/kunal/"
>>       fnames = systemfunc("ls " + dir + "DAY_*.nc")
>>
>>       print(fnames)
>>       nfil = dimsizes(fnames)
>>       print("nfil="+nfil)
>>       print("---")
>>
>>        a      = addfiles(fnames,"r")
>>        data = a[:]->bcfire(0,:,:)        ; only one timestep
>>        printVarSummary(data)        ; *****look at the structure ****
>> and sizes
>>
>>
>>    DATA = a[:]->bcfire
>>    printVarSummary(data)  ; *****look at the structure ****
>>
>> DATA should have something like (time,lat,lon) where
>>
>> dimsizes(time) >1
>>
>>
>>
>>
>>
>>
>>
>>
>> On Fri, Mar 31, 2017 at 11:01 AM, Kunal Bali <kunal.bali9 at gmail.com>
>> wrote:
>>
>>> Sorry for creating the confusion.
>>>
>>> I changed the script as you suggested me and then run the attached
>>> script. But after running this, I am getting the file only of the single
>>> day instead of all 30 files separately.
>>>
>>>
>>>
>>> Kunal Bali
>>> Research Scholar
>>> Radio & Atmospheric Science Division
>>> CSIR - National Physical Laboratory
>>> New Delhi - 110012
>>> India
>>>
>>>
>>>
>>>
>>>
>>> On Fri, Mar 31, 2017 at 10:15 PM, Mary Haley <haley at ucar.edu> wrote:
>>>
>>>> Kunal,
>>>>
>>>> I'm a little confused by what you are asking for.  You said you wanted
>>>> to extract all 30 days of data instantly, which is what the previous
>>>> "addfiles" example was meant to show.  You don't need a loop to extract all
>>>> the data.
>>>>
>>>> Did you actually try my suggestion of changing your script to use
>>>> "addfiles" instead of "addfile"?  If so, and you are still having a
>>>> problem, then it would help to see the new version of your script.
>>>>
>>>> Maybe my confusion is with what you want to do with the data after you
>>>> extract all 30 days of data. Please be more specific if possible.
>>>>
>>>> Thanks,
>>>>
>>>> --Mary
>>>>
>>>>
>>>> On Thu, Mar 30, 2017 at 10:19 PM, Kunal Bali <kunal.bali9 at gmail.com>
>>>> wrote:
>>>>
>>>>> Dear Mary,
>>>>>
>>>>> Thanks for the suggestions.
>>>>> But still, it's extracting one file at a time, not extracting all the
>>>>> 30-day files separately of UK and NEP (as mentioned in the script).
>>>>>
>>>>>
>>>>>
>>>>> Kunal Bali
>>>>>
>>>>>
>>>>>
>>>>>
>>>>>
>>>>> On Fri, Mar 31, 2017 at 3:09 AM, Mary Haley <haley at ucar.edu> wrote:
>>>>>
>>>>>> Kunal,
>>>>>>
>>>>>> You should be able to use "addfiles" to open all your files, and then
>>>>>> read your data using the special [:] syntax:
>>>>>>
>>>>>>   dir    = "/home/kunal/"
>>>>>>   fnames = systemfunc("ls " + dir + "DAY_*.nc")
>>>>>>   a      = addfiles(fnames,"r")
>>>>>>   data = a[:]->bcfire
>>>>>>
>>>>>> For more information see our addfiles examples page:
>>>>>>
>>>>>> http://www.ncl.ucar.edu/Applications/addfiles.shtml
>>>>>>
>>>>>> --Mary
>>>>>>
>>>>>>
>>>>>>
>>>>>>
>>>>>>
>>>>>>
>>>>>> On Thu, Mar 30, 2017 at 12:42 PM, Kunal Bali <kunal.bali9 at gmail.com>
>>>>>> wrote:
>>>>>>
>>>>>>> Dear NCL users
>>>>>>>
>>>>>>>
>>>>>>> I have the data of 30 days. I want to extract this data with the
>>>>>>> loop.The attached script is based on to extract the data of a single day.
>>>>>>> So it's a tedious job to extract all the 30 days data one by one. So could
>>>>>>> anyone tell me that how to start the loop in the attached file so the I can
>>>>>>> get 30 days extracted data instantly.
>>>>>>>
>>>>>>> Thank You
>>>>>>>
>>>>>>> Regards
>>>>>>> Kunal Bali
>>>>>>>
>>>>>>>
>>>>>>>
>>>>>>>
>>>>>>>
>>>>>>> _______________________________________________
>>>>>>> ncl-talk mailing list
>>>>>>> ncl-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
>>>>>
>>>>>
>>>>
>>>
>>> _______________________________________________
>>> ncl-talk mailing list
>>> ncl-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/20170331/130b967c/attachment.html

More information about the ncl-talk mailing list