Hi,
I have a number of netcdf files (one per processor of a parallel job)
that I need to combine
to visualize. This has proven to be a painfully slow process and I
wonder if any of you have any suggestions on speeding it up.
Currently I am using the following ncl:
fnames = systemfunc("ls movies/*.nc")
f = addfiles(fnames,"r") ; Open netCDF files.
fcnt = dimsizes(fnames)
vcnt=new((/fcnt,3/),"integer")
do i=0,fcnt-1
vcnt(i,:) = dimsizes(f[i]->Q)
end do
ntimes = vcnt(0,0)
nlev = vcnt(0,1)
nxy = dim_sum(vcnt(:,2))
var = new((/ntimes,nlev,nxy/),"double")
do n=1,ntimes
do k=1,nlev
var(n,k,:)=f[:]->Q(n,k,:)
end do
end do
-- Jim Edwards jedwards@ucar.edu IBM Applications Analyst NCAR SCD BOULDER CO 303-497-1842_______________________________________________ ncl-talk mailing list ncl-talk@ucar.edu http://mailman.ucar.edu/mailman/listinfo/ncl-talk
This archive was generated by hypermail 2b29 : Tue Oct 12 2004 - 17:06:33 MDT