Hello,
I wanted to point out that the file has an error in some of the date
strings, 221-JUL-2012, which you will want to fix before
proceeding. Is it 22-Jul or 21-Jul?
If there were exactly one profile for each time, then we could read
the data as a grid. However because some days have several profiles
listed at the same times, and some days have no data, I think the
best solution is to use a scatter-to-grid function. Here are the
options:
yes? show func scat2grid*zt
SCAT2GRIDGAUSS_ZT(ZPTS,TPTS,F,ZAXPTS,TAXPTS,ZSCALE,TSCALE,CUTOFF,0)
Use Gaussian weighting to grid scattered data to a ZT
grid.
ZPTS: Z coordinates of scattered ZT locations
TPTS: T coordinates of scattered ZT locations
F: F(Z,T) data at scattered ZT locations. May also be fcn
of X,Y,E,F
ZAXPTS: Z axis coordinates of a regular output grid
TAXPTS: T axis coordinates of a regular output grid
ZSCALE: Mapping scale in Z direction, in data units (e.g.
m or km)
TSCALE: Mapping scale in T direction, in data units (e.g.
hrs, days)
CUTOFF: Weight cutoff: use pts within cutoff*scale
0: Not used; included for backwards compatibility
SCAT2GRIDLAPLACE_ZT(ZPTS,TPTS,F,ZAXPTS,TAXPTS,KAY,NRNG)
Use Laplace weighting to grid scattered data to a ZT grid.
ZPTS: Z coordinates of scattered ZT locations
TPTS: T coordinates of scattered ZT locations
F: F(Z,T) data at scattered ZT locations. May also be fcn
of X,Y,E,F
ZAXPTS: Z axis coordinates of a regular output grid
TAXPTS: T axis coordinates of a regular output grid
KAY: Laplace/spline mixing parameter KAY; 0.0=pure
Laplace; inf=pure spline
NRNG: Effect parameter NRNG; grid vals undef if > NRNG
grid steps from pts
These will interpolate your data to a grid in depth and time. The
COLUMNS command translates date strings to numbers which are days
since the start of 1990, so for the function you should define a
time axis in those units. So it might look something like this:
COLUMNS/SKIP=1/delim=";"/var="Latitude,Longitude,Depth,TIME,Chl_a"/type="num,num,num,eurodate,num"
TSC.dat
define axis/t=14-Jul-2012:3-Aug-2012:1/units=days/t0=1-jan-1990
timeaxis
define axis/z=0:120:1/depth/units=meters zaxis
let zaxpts = z[gz=zaxis]
let taxpts = t[gt=taxis]
let cvar =
SCAT2GRIDGAUSS_ZT(DEPTH,TIME,ZSEQUENCE(chl_a),zaxpts,taxpts,1,1,5,0)
shade cvar ! or fill cvar
You will want to read in the documentation about the parameters in
the gridding function that you choose, and adjust them according to
your judgement of how the gridding should be done. My choices above
are just some workable values, but not necessarily the best choices.
You could use the /units= or /title= on the "let cvar" command, to
give your variable those attributes.
On 4/2/2013 9:12 PM, Karnan C wrote:
Dear sir,
221-Jul-2012 is replaced by 21-jul-2012.
I can understand ascii data file reading.
Here I have a problem with irregular time (15-07-2012 then
20-07-2012 to 02-08-2012).
I need to show a small white space between 15 and 20-jul-2012
data in the vertical section of fill.
(I need a fill image not plot).
|