Re: writing several .nc files to 1 file

[Emilio Mayorga <emiliom@xxxxxxxxxxxxxxxx>]

Hi Brent,

In your SAVE statement, I imagine ao[d=1] is the variable from one of the
files, as is. Also, I'm guessing you're not working with Ferret's
multi-file netcdf scheme.  I will make a suggestion that I haven't tested,
but it's something that's worked in the time dimension.

It seems you need to first define the common Z axis with DEFINE AXIS. Say
the name is "depth". You have irregularly spaced values, and I can't
remember right now how to create such an axis. After the axis is created,
create a "z stamp" variable:

LET zstamp = Z[GZ=depth] * 0
Now to write the first ao z-slice, define a new ao variable to stamp ao
with its z value:
LET aoz = ao[d=1] + zstamp[K=1]

Now you can write this first value into the file
SAVE/FILE=verttest.nc/KLIMITS=1:17/CLOBBER/K=1 aoz

For subsequent z slices, eg, K=5:
LET aoz = ao[d=1] + zstamp[K=5]
SAVE/FILE=verttest.nc/APPEND/K=5 aoz

Of course, you should add the title and units info to aoz if you want to
preserve it, and you should probably put this together in a script file
and use REPEAT/K=1:17

I hope this works, and good luck.

Emilio Mayorga
________________________________________________________________
University of Washington
School of Oceanography			emiliom@u.washington.edu
Box 357940				ph. (206) 221-6756
Seattle, WA  98195-7940	  USA		FAX (206) 685-3351
http://boto.ocean.washington.edu/staff_pages/emilio

On Fri, 4 May 2001, Brent A McDaniel wrote:

> 
> Hi Ferret Users,
> 
> Hopefully this should be an easy problem for the more experienced users.
> I've got 17 netcdf that I want to combine into 1 netcdf file.  They all
> are of dimension 144x37x1 (no time) of the same variable (ao).  The 144x37
> are the same but they each have a different z-axis. (1000, 925, 850,
> 700....)  I want to combine them into a 144x37x17 file.  I tried the
> simple brute force technique of
> save/file=verttest.nc/klimits=1:17/clobber ao[d=1]
> and got the message
>  LISTing to file verttest.nc
>  **ERROR: invalid command: Exceeds grid limits: klimits=1:17
>           Variable: AO
> 
> For completeness, here's the output of show data/full
> yes? show data/full 1
>      currently SET data sets:
>     1> ./airanom.1000.1000.nc
> 
>  name     title                             I         J         K
> L
>  AO       (TEST97*-1)/3600                 1:144     1:38      1:1
> ...
>                on grid GGP1 with -1.E+34 for missing data
>              X=1.2W(-1.2):1.2W  Y=3.8S:91.2N  Z=999.5:1001
> 
> Can someone help me out please?
> 
> Thanks in advance,
> 
> Brent
> 
>