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) "; )
--------------------------------------------------------------------------------------------