Re: Plotting more than 495000 points (again)?

[verschel@xxxxxxxxxxxxxxxxxxxxx (Mark A. Verschell)]

Sophie,

	This is a problem I run across fairly often. There are two ways to deal
with this. What you tried to do with the command:
yes? SHADE count[I=1:4001:10,J=1:3501:10]
was to try to regrid the data. This is done through describing a new grid and
regridding before plotting. Something like (ignoring the fact that I don't know
what your underlying grid looks like, so I assume there is a 1 to 1
correspondence between index (i,j) and coordinate (x,y)):

yes? DEFINE AXIS/x=1:4001:10 newx
yes? DEFINE AXIS/y=1:3501:10 newy
yes? DEFINE GRID/x=newx/y=newy newg
yes? SHADE count[g=newg]

Alternatively, if you REALLY want to plot all the points:

yes? stat count ! To get the min and max for setting levs
yes? SHADE/xl=1:4001/lev=(0,0,0) count[i=1:101,j=1:3501]
yes? REPEAT/i=102:3902:100 SHADE/over/nolab/lev=(0,0,0) count[i=`i`:`i+99`]

where the 0,0,0 is replaced by info from the stat. For example if the min and
max you get from a stat of count is 10 and 100, then the /lev could be
/lev=(10,100,10)

I would suggest doing the regridding however!

Mark


-- 
 Mark Verschell                     NASA/Goddard Space Flight Center
 verschell@neptune.gsfc.nasa.gov    Code 970/Lab. for Hydrospheric Proc.
 301-614-5669  Fax: 301-614-5666    Greenbelt, MD 20771

                  USQC Certified: 100% Microsoft Free