TMAP ERR: limit on storage for coordinates has been reached

[Rob Scott <rscott@xxxxxxxxxxxxxxxxxx>]

Dear ferret:

I am trying to combine 488 netCDF files, each representing 1 week of data 
on a mercator lat lon grid, into one big file.
However I get the above error while reading the 136th netCDF file.

I have tried following the advice of Steve Hankin:
http://ferret.wrc.noaa.gov/Ferret/Mail_Archives/fu_2000/msg00183.html

yet the problem persists.

In particular, I'm doing the following:

I set up the time axis, and a 3D grid based on the 2D mercator grid that 
"grid_0001" in on, and an empty variable on this grid called "ssha"

DEFINE AXIS/T=15627:19036:7/UNITS=DAYS taxis
DEFINE GRID/LIKE=grid_0001/T=taxis grid_aviso
LET/TITLE="AVISO Sea Surface Height Anomaly"/UNITS=m \
     ssha = (1/0) * x[G=grid_aviso] * y[G=grid_aviso] * t[G=grid_aviso]

SAVE/CLOBBER/FILE="ssha_aviso_AR.cdf"/RIGID/ILIMITS=1:91/JLIMITS=260:364/LLIMITS=1:488 
ssha

LET dummie = x[G=grid_aviso] * y[G=grid_aviso] * t[G=grid_aviso]

REPEAT/k=1:488 (LET a = names[i=`k`]; \
   SET DATA/regulart `a`; \
   LET/TITLE="AVISO Sea Surface Height Anomaly"/UNITS=m ssha = 
RESHAPE(ssha_tmp,dummie[L=`k`]); \
   SAVE/APPEND/FILE="ssha_aviso_AR.cdf"/RIGID/ILIMITS=1:91/JLIMITS=260:364/LLIMITS=1:488 
ssha[L=`k`]; \
   CANCEL DATA `a`)


I'm guessing that the problem arises because the mercator grid is irregular 
in lat.
But I'm cancelling each file after I use it, so I don't see why it's 
accumulating this info.

Thanks for your help,
Rob

PS
I'm running ferret v5.41 on Linux 2.4.3-12smp








Rob Scott
Research Associate
Institute for Geophysics
University of Texas at Austin
4412 Spicewood Springs Rd.
Building 600
Austin, TX
78759

Office: +1-512-471-0375