Re: Read large volume data from modis files

Do you mean that for each of several housand modis level 2 data files
you want to interpolate via (say) linear interpolation to a fixed
0.1 degree grid? Are these swath data or grided data?

There should be one loop over all the files

  dirModis = "./"
  filModis = systemfunc("cd "+dirModis " ; ls ????*")
  nfil = dimsizes(filModis)

  do nf=0,nfil-1
     f = addfile (dirModis_filModis(nf), "r")
     x = f->
     y = f->
     z = f->

     xnew = ...interpolate...
     ynew = ...interpolate...
     znew = ...interpolate...

; write out a netCDF with time unlimited

  end do

After it is over, use the netCDF operator 'ncrcat' to
merge all the individual files

On 8/21/12 11:06 AM, Xin Xi wrote:
> Hi there,
> I am working on several thousand modis level 2 data files over a few
> years and a large domain (30 degX40 deg). I want to read in several
> variables from these modis files into a fixed grid domain at 0.1 degree
> resolution for each day. I expect quite a few looping procedures, and it
> is very likely that there will be multiple data points in the same grid,
> but the count will vary, depending on the overlapping of modis passes
> and retrieval quality. So I am wondering how I can deal with the data
> array (for each cell) with unknown sizes/lengths? do I set it to
> "unlimited"? Also any helpful advice are appreciated. Thanks!
> Xin
>
>
> _______________________________________________
> ncl-talk mailing list
> List instructions, subscriber options, unsubscribe:
> http://mailman.ucar.edu/mailman/listinfo/ncl-talk
>
_______________________________________________
ncl-talk mailing list
List instructions, subscriber options, unsubscribe:
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
Received on Wed Aug 22 14:58:19 2012