Re: read from 239 netCDF files and write to only one netCDF file

• To: Qiaozhen Mu <qiaozhen@xxxxxxxxxxxxxxxxxx>
• Subject: Re: read from 239 netCDF files and write to only one netCDF file
• From: Phil Rasch <pjr@xxxxxxxx>
• Date: Wed, 9 Jan 2002 08:26:44 -0700

Try the NCO operator suite for this type of task. It is much more
suited for it. See
http://nco.sourceforge.net 


On Tue, Jan 08, 2002 at 05:10:04PM -0600, Qiaozhen Mu wrote:
> Dear all hands,
> 
> I try to get the field CLOUD out from 239 history files in netCDF format. 
> There 
> are 18 levels, 128 longitudes and 64 latitudes for each CLOUD field in each 
> history file. So when I ran the ncl, it reminded me that 
> 
> """""
> fatal:NclMalloc Failed:[errno=11]:Resource temporarily unavailable
> fatal:_NclBuildConcatArray : Memory allocation failure
> 
> fatal:Execute: Error occurred at or near line 6
> """""
> 
> There isn't enough memory to deal with it. What can I do? The following is 
> the 
> CLOUD.ncl I used to get the CLOUD out.
> 
> 
>           begin
>             diri = "/disk/cj2/ccm/ccm-mqz/run.amip2/"  ; directory where fil
> es reside
>             fils = systemfunc ("ls "+diri+"19??-??.nc")
>             f    = addfiles (fils, "r")   ; note the "s" of addfile
> 
>             ListSetType (f, "cat")        ; concatenate or "merge"
>             T    = f[:]->CLOUD
>             print (dimsizes(T))           ; T will have 60 time steps
>             T!0  = "time"
>             T!1  = "lev"
>             T!2  = "lat"
>             T!3  = "lon"
>             T&time = f[:]->time           ; time coord variable
>             T&lev  = f[0]->lev            ; get lev from the 1st file
>             T&lat  = f[0]->lat            ; get lat from the 1st file
>             T&lon  = f[0]->lon            ; get lon from the 1st file
> 
>             filo = "./nc/CLOUD.nc"
>             system ("/usr/bin/rm "+filo)  ; remove any pre-existing file
>             fo   = addfile(filo , "c")    ; open output file
>             fo@title = "CLOUD from model files"
>             
>             fo->CLOUD = T                    ; write CLOUD to a file
>           end
> 
> Anyone knows how to deal with this using any netCDF softwares? I can't 
> insatall  
> nco on my SUN computer now. Thanks in advance.
> 
> Best regards,
> Qiaozhen
> 
> ***************************************
> *The University of Texas at Austin    *
> *Institute for Geophysics             *
> *4412 Spicewood Springs Rd., Bldg. 600*
> *Austin, Texas 78759-8500             *
> *phone:(512) 471-0462                 *
> ***************************************

-- 
Phil Rasch, Climate Modeling Section, National Center for Atmospheric Research
Mail     --> P.O. Box 3000, Boulder CO 80307  
Shipping --> 1850 Table Mesa Dr, Boulder, CO 80305
email: pjr@xxxxxxxx, Web: http://www.cgd.ucar.edu/cms/pjr Phone: 303-497-1368, 
FAX: 303-497-1324 

• References:
  • read from 239 netCDF files and write to only one netCDF file
    • From: Qiaozhen Mu