Re: Creating 3D grid from scattered data

[<SPILLANE@xxxxxxxxxxxxx>]

Hello,
	I think the only problem with your .jnl file is that you don't
limit the region of the z-axis in the append steps - this is missing
from example 4 too.
	The following demonstrates how three slabs can be written
def axis/x=1:5:1 xax 
def axis/y=10:40:10 yax 
def axis/z=100:300:100 zax 
 
def grid/x=xax/y=yax/z=zax whole 
 
def axis/z=100:100:100 zax1 
def grid/like=whole/z=zax1 slab 
 
file/form=free/var=v1/col=5/g=slab z1.d 
let v=v1[g=whole] 
list v 
save/file=xyz.nc/klimits=1:3 v 
sp ncdump xyz.nc 
can data 1 
 
def axis/z=200:200:100 zax1 
file/form=free/var=v1/col=5/g=slab z2.d 
list v 
save/file=xyz.nc/z=200/app v 
sp ncdump xyz.nc 
can data 1 
 
def axis/z=300:300:100 zax1 
file/form=free/var=v1/col=5/g=slab z3.d 
list v 
save/file=xyz.nc/z=300/app v 
sp ncdump xyz.nc 

The "list v" steps show that, as expected, v holds only one slab at a time.
Without the "/z=200" or equivalently "region/z=200" the missing values in
the other layers overwrite the previously written data.

Good luck,
Mick