TMAP ERR: limit on storage for coordinates has been reached

To: ferret <ferret_users@xxxxxxxxxxxxxxxxxxx>
Subject: TMAP ERR: limit on storage for coordinates has been reached
From: "Paulo B. Oliveira" <pbo@xxxxxxxxx>
Date: Wed, 17 Sep 2003 15:25:23 +0100
Organization: IPIMAR
Sender: owner-ferret_users@xxxxxxxxxxxxxxxxxxxx

Hi all,

After looking the mail archive concerning this subject,
I believe that I'm facing a different type of ferret 
behaviour related with the use of data irregularly spaced
in time. 

The idea is to plot a series of lon/lat vectors, unevenly
spaced on time, over 2D fields. In my case, SST fronts extracted
from synoptic satellite images, over monthly SST averages. 

The attached scripts reproduce the ferret behaviour using
abstract data because the data files I'm using are too big
to attach:

yes? go vecs_gen  ! generate abstract fronts

and,

yes? go vecs_overlay ! now, plot the fronts over the monthly field

the script runs OK for the first time steps and then continues
with the error message:
 **TMAP ERR: limit on storage for coordinates has been reached
             MAX= 250000
but the vectors are not plotted.

Its interesting to note that if the script is interrupted and 
restarted afterwards, the error message comes right after the 
first time step.

It seems that some coordinate info is acumulating into 
memory at each time step and may be there is a command
to clear it. Any hints will be appreciated.

Regards, 

Paulo.

 ! NOAA/PMEL TMAP
 ! FERRET v5.51  
 ! Linux 2.4.3-12smp - 03/05/03
 ! 17-Sep-03 10:56     

!----------------- create abstract fronts
file/var=myt irreg_t.dat

def ax/z=1:500:1 zv
def ax/t/from_data/name=tv/unit=day/t0=1-jan-1900 myt
def gr/z=zv/t=tv gv

let/tit=longitude/unit=longitude lon = -20+k[g=gv]/50+0*t[g=gv]
let/tit=latitude/unit=latitude   lat = 32+k[g=gv]/50+0*t[g=gv]

set mem/size=30
save/file=vecs.cdf/rigid/head=enhanced/llim=1:7200/clobb lon,lat

quit

 ! NOAA/PMEL TMAP
 ! FERRET v5.51  
 ! Linux 2.4.3-12smp - 03/05/03
 ! 17-Sep-03 10:56     

!---------------------------------------- generate an abstract field
!                                         and overlay the vectors
def ax/x=-30:-6/np=274/unit=longitude x0
def ax/y=32:45/np=149/unit=latitude  y0
def ax/t="15-jan-1990":"15-dec-1999":`365.25/12`/unit=day/t0="1-jan-1900" t0

def gr/x=x0/y=y0/t=t0 g0

let/tit="my monthly field" f0 = x[g=g0]+y[g=g0]+t[g=g0]

canc mode verif
repeat/l=1:120 ( say `l`; shade f0; go vecs_overlay_sub )

quit

def sym date_str `f0,return=tstart,precision=3`
ppl SET mon      $ELEMENT(2,-,date_str)
ppl SET yyyy0    $ELEMENT(3,-,date_str)
ppl SET yyyy     $EDIT(yyyy0,collapse)
set data/save
use vecs.cdf
 def sym ti  1-($mon)-($yyyy)
 def sym tf 28-($mon)-($yyyy)
 set reg/t=($ti):($tf)
 plot/ov/vs/lin/nolab/col=bla lon,lat 
 canc reg/t
 canc data vecs.cdf
set data/restore
quit

Previous by thread: How to make ferret handle a variable number of arguments?
Next by thread: User-Specified Key range

[Thread Prev][Thread Next][Index]

Dept of Commerce / NOAA / OAR / PMEL / TMAP