I have 24 netcdf files each represents different depths
(0,2,5,10,.......150,200), i tried to append these files by
defining a z-axis and appending along that axis. I tried the
following script in the ferret which was worked previously but
not working now and showing an error message which i can not
understand. I am looking forward for the help. Thanks in
advance. Error message and script as follows:
LISTing to file
HYCOM_temp_apr_oct2014.cdf
***** ERROR IN
CD_GET_GENERIC_BOUNDS *****
**TMAP ERR: error in line
definition
file coords dont match
variable coords on axis DEPTH
LIST/FORMAT=CDF/append/file=HYCOM_temp_apr_oct2014.cdf
TEMP[K=1:24]
Command file, command group, or
REPEAT execution aborted
-------------------------------------
ferret script ----------------------------
let a = {SPAWN: "ls *.nc4"}
define
axis/z/from_data/units=meters/DEPTH/npoints=24
zaxis={0,2,4,6,8,10,12,15,20,25,30,35,40,45,50,60,70,80,90,100,125,150,200,250}
define grid/Z=zaxis gdepth
LET depth_stamp = Z[G=gdepth] * 0
let TEMP = WATER_TEMP + depth_stamp
let zz = Z[gz=gdepth]
let nz = `zz,return=ksize`
list nz
pause
sp rm model_temp_apr_oct2014.cdf
REPEAT/range=1:`nz`/NAME=r ( define
symbol c = `r`; define symbol b = `a[i=($c)]`; \
use
"($b)";\
save/append/file=model_temp_apr_oct2014.cdf TEMP[K=1:24];\
can
da "($b) "; )
--------------------------------------------------------------------------------------------