[ncl-talk] Binning multilevel MODIS data using bin_sum

[AJILESH PP]

Hi,
Thank You for the prompt reply. Sorry for the mistake, Now loop is working
output has 20 levels, but all the levels have the same data values (all
levels are same).

On Thu, 21 Apr 2016 at 20:56 Rick Brownrigg <brownrig at ucar.edu> wrote:

> Hi,
>
> I'm not totally sure, but I wonder if "level" in this line should be "lev"?
>
>
> bin_sum(GBIN(level,:,:),GKNT(level,:,:),lon,lat \
>
> In any case, a few strategic print() statements on your indexing variables
> will show where/when thing go arwy.
>
> HTH...
> Rick
>
>
> On Thu, Apr 21, 2016 at 7:10 AM, AJILESH PP <ajileshparolla at gmail.com>
> wrote:
>
>> Dear all
>> I want to bin MODIS level 2 products for a period of time. I modified an
>> example script given in ncl HDF examples and it is giving me only output of
>> a single level, where the desired output is for 20 levels. I tried to add a
>> level loop and level dimension also as shown in the script below.
>> --------------------------------------
>>   GBIN  = new ( (/nlev,nlat,mlon/), float )
>>   GKNT  = new ( (/nlev,nlat,mlon/), integer )
>>
>> snip
>>
>>  do nf=0,nfil-1
>>    do lev=0,nlev-1
>>      print(nf+"   "+fili(nf))
>>      f       = addfile(diri+fili(nf), "r")
>>                                             ; read data
>>      x          = short2flt_hdf( f->$vNam$ )
>>      lat2d1     = f->Latitude
>>      lon2d1     = f->Longitude
>>      lat2d      = conform(x,lat2d1,(/1,2/))
>>      lon2d      = conform(x,lon2d1,(/1,2/))
>>      level        =
>> (/1000,950,920,850,700,620,500,400,300,250,200,100,150,100,70,50,30,20,10,5/)
>>      x!0          = "level"
>>      x&level      = level
>>      nx           = product(dimsizes(x))
>>     bin_sum(GBIN(level,:,:),GKNT(level,:,:),lon,lat \
>>            ,ndtooned(lon2d),ndtooned(lat2d),ndtooned(x) )
>> printVarSummary(x)
>> ---------------------------------------
>> but it is giving this error:
>> fatal:Subscript out of range, error in subscript #0
>> fatal:An error occurred reading GBIN
>>
>> Can anybody help me to resolve this issue?
>> Please find my script and data as follows
>> Data:   MOD07_L2.A2012166.0*.hdf
>> Script:  bin.ncl
>>
>> Any help would be appreciated
>> Thank You!
>>
>> _______________________________________________
>> 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/20160422/ae4b79dc/attachment.html