[Thread Prev][Thread Next][Index]

Re: TMAP ERR: limit on storage for coordinates has been reached

Hi Rob,
I think you're right that the trouble is due to the irregular grid for
the input files, and that storage space for irregular axes is being
used up as you read in the different files. The message from the Users
List that you are following deals with the time axes, and you do need
that part of the solution, but you also have an irregular x axis being defined
for each input data set.  Canceling the data sets does not recover this
storage space (perhaps it should!).

Starting with Ferret version 5.41, there is a new option "cancel axis/all"
which will cancel the irregular axes, and make the coordinate storage
available for re-use.  If you are running v5.41, add a "cancel axis/all"
with each "cancel data"  command and that should do the trick.

Ansley Manke

Rob Scott wrote:

> 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

Ansley Manke  Pacific Marine Environmental Laboratory  Seattle WA  (206)526-6246

[Thread Prev][Thread Next][Index]

Dept of Commerce / NOAA / OAR / PMEL / TMAP

Contact Us | Privacy Policy | Disclaimer | Accessibility Statement