Re: trouble with saving 3D data

[Mick Spillane <spillane@xxxxxxxxxxxxx>]

>

Hi Muyin,
    I think what you need to do is to define the "combined" temperature variable on
a
new grid.  Then you can either do the integration right away on that variable, or
perhaps write the new variable to a netcdf file if it would simplify other tasks.
    Here is a demo of what I have in mind
 ! NOAA/PMEL TMAP
 ! Program FERRET
 ! Version 5.22 - 07/27/00
 !  7-Nov-00 08:29

def axis/x=1000:5000:1000 xax ; def axis/y=100:400:100 yax
def axis/z=10:30:10 zax3 ; def axis/t=1:2:1 tax
def grid/x=xax/y=yax/z=zax3/t=tax gxyzt

! in the above zax3 and the grid gxyzt is for the non-surface data
! which we define in this example as

let temp=x[g=gxyzt]+y[g=gxyzt]+z[g=gxyzt]+t[g=gxyzt]

!-----------------------------------------------------------------

def grid/x=xax/y=yax/t=tax gxyt

! the grid gxyt above is for the surface data defined in this example as

let surf=x[g=gxyt]+y[g=gxyt]+t[g=gxyt]

!-----------------------------------------------------------------

! now introduce a z-axis and grid that combines surface and non-surface

def axis/z=0:30:10 zax4
def grid/x=xax/y=yax/z=zax4/t=tax grd

! and define a combined variable

let combo=if(z eq 0)then surf[g=grd] else temp[g=grd]

list combo

! looks good, so save it to a netCDF file

save/file=combo.nc combo



--
|--****--****-*---*---***--***--|____spillane@pmel.noaa.gov____|
|-*__---*-----*--*-*--*--*-*--*-|_SCIENCE APPLICATIONS SUPPORT_|
|--***--*-----*-*---*-***--***--|____EPIC/Ferret/PlotPlus______|
|-----*-*-----*-*****-*----*----|__Room 2070 Bldg#3 NOAA/PMEL__|
|-****---****-*-*---*-*----*----|____Phone_:_(206)526-6780_____|