[Thread Prev][Thread Next][Index]

Re: ERROR Y coordinate field



Hi Ansley and thanks for your great help,

At 20:41 26/01/2005, Ansley Manke wrote:
Hi Julien,
I don't see anything wrong with your scripts.  You might try closing
the coastline file inside script gg.
I've just searched in the documentation, but don't find anything on closing files, sorry ; what's the command please ?

Somehow it seems that the comands
in gg.jnl  may be trying to read data that does not exist into the variables
v1 and v2 (though this does not happen with a quick test that I tried.)

After the script runs and you see the error, try commands like
SHOW DATA
LIST/I=1:20 v1, v2

to see if you can diagnose what is happening.
!-> REPEAT: L=3
 **ERROR: all data have same value: Y coordinate field
CONTOUR/FILL/i=3:28/j=5:35/axes=1,1,1,1/line/key rain + 2, longicrs, latitcrs
Command file, command group, or REPEAT execution aborted
yes? show data
     currently SET data sets:
    1> /project/aspa_ferret/data/out.nc  (default)
 name     title                             I         J         K         L
 LATITCRS LATITUDE (SOUTH NEGATIVE)        1:30      1:39      ...       ...
 LONGICRS LONGITUDE (WEST NEGATIVE)        1:30      1:39      ...       ...
 RAIN     ACCUMULATED NONCONVECTIVE PRECI  1:30      1:39      ...       1:97
 VNORM    10-meter Speed Norm              1:30      1:39      ...       1:97
 TAU      Optical Thickness                1:30      1:39      ...       1:97

    2> /project/aspa_ferret/data/coastline_france.txt
 name     title                             I         J         K         L
 V1       V1                               1:140829  ...       ...       ...
 V2       V2                               1:140829  ...       ...       ...

yes? LIST/I=1:20 v1, v2
 **ERROR: variable unknown or not in data set: V1

  Also, you can begin
your script with the command
SET MODE DIAG
which will output text lines as the script runs, and you will be able to see
where Ferret is doing read operations, and so on.  That might be helpful.
the result :

yes? g tt
set mode diag
use out.nc
repeat/l=2:3 ( \
fill/i=3:28/j=5:35/axes=1,1,1,1/line/key rain + 2, longicrs, latitcrs ;
\
go gg ; \
frame/file=toto`l`.gif )
!-> REPEAT: L=2
getgrid EX#1 C: 5 I: 1 1 J: 1 1 K: 1 1 L: 1 1
allocate dynamic grid GJJ2 J_CROSS I_CROSS NORMAL TIME
strip limits reconciliation : EX#1
eval EX#1 C: 5 I: 3 28 J: 5 35 K: -111 -111 L: 2 2
found RAIN M: 1 I: 3 28 J: 5 35 K: -111 -111 L: 2 2
constan cnst M: 8 I: -111 -111 J: -111 -111 K: -111 -111 L: -111 -111
-DELETE cnst M: 8 I: -111 -111 J: -111 -111 K: -111 -111 L: -111 -111
getgrid EX#2 C: 6 I: 1 1 J: 1 1 K: 1 1 L: 1 1
allocate dynamic grid GJJ1 J_CROSS I_CROSS NORMAL NORMAL
strip limits reconciliation : EX#2
eval EX#2 C: 6 I: 3 28 J: 5 35 K: -111 -111 L: 2 2
found LONGICRS M: 2 I: 3 28 J: 5 35 K: -111 -111 L: -111 -111
getgrid EX#3 C: 7 I: 1 1 J: 1 1 K: 1 1 L: 1 1
allocate dynamic grid GJJ1 J_CROSS I_CROSS NORMAL NORMAL
strip limits reconciliation : EX#3
eval EX#3 C: 7 I: 3 28 J: 5 35 K: -111 -111 L: 2 2
found LATITCRS M: 4 I: 3 28 J: 5 35 K: -111 -111 L: -111 -111
setting up 2D plot
**ERROR: all data have same value: Y coordinate field
CONTOUR/FILL/i=3:28/j=5:35/axes=1,1,1,1/line/key rain + 2, longicrs, latitcrs
Command file, command group, or REPEAT execution aborted
-DELETE EX#1 M: 8 I: 3 28 J: 5 35 K: -111 -111 L: 2 2


In fact it seems there is a problem in the first loop L=2...


Many thanks in advance,

Julien


Ansley

Julien Demaria wrote:

Hi,

I need some help on a little code, I really don't understand why I've this error "**ERROR: all data have same value: Y coordinate field" :


the script:
-----------
use out.nc
repeat/l=2:3 ( \
fill/i=3:28/j=5:35/axes=1,1,1,1/line/key rain + 2, longicrs, latitcrs ; \
go gg ; \
frame/file=toto`l`.gif )
-----------


the file gg.jnl:
-----------
set data/save
set grid/save
define axis/x/x=1:150000:1 longmax
define grid/x=longmax biggrid
columns/skip=1/grid=biggrid coastline_france.txt
set var/bad=-9999 v1
set var/bad=-9999 v2
plot/over/nolabels/noaxis/vs/line=14 v1, v2
set data/restore
set grid/restore
-----------


ncdump -h out.nc :
-----------
netcdf out {
dimensions:
time = 97 ;
i_cross = 39 ;
j_cross = 30 ;
pressure = 22 ;
variables:
float time(time) ;
time:units = "minutes since 2005-01-11 00:00:00" ;
time:long_name = "time" ;
float latitcrs(i_cross, j_cross) ;
latitcrs:units = "degrees" ;
latitcrs:long_name = "LATITUDE (SOUTH NEGATIVE)" ;
float longicrs(i_cross, j_cross) ;
longicrs:units = "degrees" ;
longicrs:long_name = "LONGITUDE (WEST NEGATIVE)" ;
float rain(time, i_cross, j_cross) ;
rain:units = "cm" ;
rain:long_name = "ACCUMULATED NONCONVECTIVE PRECIPITATION"
;
float vnorm(time, i_cross, j_cross) ;
vnorm:units = "m s{-1}" ;
vnorm:long_name = "10-meter Speed Norm" ;
float tau(time, i_cross, j_cross) ;
tau:units = "-" ;
tau:long_name = "Optical Thickness" ;
}
-----------


The error happens when l=3, on the second loop of the repeat.

If I remove the biggrid use in the gg.jnl file, it works ! Why ?

I've also seen the same error with X on another similar code...


Thanks in advance for help,

Julien






[Thread Prev][Thread Next][Index]

Dept of Commerce / NOAA / OAR / PMEL / TMAP

Contact Us | Privacy Policy | Disclaimer | Accessibility Statement