Re: multi-netcdf files for ensemble simulations

[Ansley Manke <ansley.b.manke@xxxxxxxx>]

Hi Sonia,
I am going to quote your script here -- your attachment doesn't show
up in the email archives.
--------------------------
  def axis/x=1:360:1 xaxis
  def axis/y=-90:90:1 yaxis
  def axis/z=1:2:1/DEPTH zaxis
  def grid/x=xaxis/y=yaxis/z=zaxis glace_grid
 
  use "DL_HAD3_ok/HadAM3_ENSR_SIM01_SoilWet.nc"
  def axis/z=1:1:1 zaxis1
  def grid/like=glace_grid/z=zaxis1 glace_grid1
  let/title="soilwetHAD3" sw=soilwet[G=glace_grid1]
  save/file=testglace.nc/klimits=1:2/k=1 sw
  can data/all
 
  use "DL_HAD3_ok/HadAM3_ENSR_SIM02_SoilWet.nc"
  def axis/z=2:2:1 zaxis1
  def grid/like=glace_grid/z=zaxis1 glace_grid1
  let/title="soilwetHAD3" sw=soilwet[G=glace_grid1]
  save/append/file=testglace.nc/k=2 sw
  can data/all
--------------------------

In Example 4  of  Chap. 10 Sec. 2, you'll notice that the first SAVE
command, which writes the NetCDF file header, writes variables that are
defined using the "zaxis" Z axis, which has all the desired output levels.
This way the file is set with the right number of levels. So you want to 
put
your variable sw on the first level of  axis "zaxis" when you first 
write it out.

The example in the Users Guide is also a bit different from what 
you'redoing,
in that it's reading ascii files with the FILE command. That is why the axes
zaxis1, zaxis2, are used.  You do not need those. Do your input files
have Z axes?  If not, then you will need to force them onto a Z axis, by 
adding
to it a variable whose values are zero, at the correct depth, i.e.
    let sw=soilwet + 0*z[gz=zgrid1,k=1].

Your script would then be something like this:

  def axis/z=1:2:1/DEPTH zaxis
 
  use "DL_HAD3_ok/HadAM3_ENSR_SIM01_SoilWet.nc"
  let/title="soilwetHAD3" sw=soilwet + 0*z[gz=zaxis,z=1]
  save/file=testglace.nc/klimits=1:2/k=1 sw
  can data/all
 
  use "DL_HAD3_ok/HadAM3_ENSR_SIM02_SoilWet.nc"
  let/title="soilwetHAD3" sw=soilwet + 0*z[gz=zaxis,z=2]
  save/append/file=testglace.nc/k=2 sw
  can data/all


Sonia Seneviratne wrote:

> Dear Ferreters,
>
> Hopefully someone already found a solution to the following problem:
>
> I have a series of NetCDF files entailing output from GCM ensemble simulations. 
> In order to analyse the data, it would help me to be able to concatenate the 
> NetCDF files on one axis, with each axis "level" corresponding to one member of 
> the ensemble (I thought of using the Z axis for this as I only have surface 
> fields). 
> I did a few tests following the instructions from Chap. 10  Sec. 2, however, the 
> output file does never entail more than one z-level (see attachment). Did anyone 
> have this problem before and manage to find a solution?
>
> Thank you very much!
> Sonia
>
> ---------------------------------------------------------------------
> Sonia Seneviratne (Dr. Sc.)
> NASA Goddard Space Flight Center, Code 900.3, Bldg. 33, Rm A105
> Greenbelt Road, Greenbelt, MD 20771, USA
> Phone: +1 301 614 5930; Fax: +1 301 614 6297
> senevira@janus.gsfc.nasa.gov
> ---------------------------------------------------------------------
>
>