Hi,
I'd like to obtain with gnuplot a graph like this:
using a 3D dataset like that shown by the red points.
Now I obtain a contour graph with the following option:
/set pm3d interpolate 8,8
set dgrid3d 8,45,20
set palette maxcolors 100
set view map
splot "KH.dat" using 1:2:3 with pm3d;
/
I have two question:
* A) How can I add the real points (not the grid!!) on the graph?
because if I use this option (splot "KH.dat" using 1:2:3 with
pm3d,\ "KH.dat" using 1:2:3 w points palette ps 1 pt 4;) I draw
the grid points on the contour graph, but not the real points...
* If possible, how can I clip the grid3d? I'd like to obtain a
contour graph only when I have the real data (see the figure1) and
not "extrapolate" anything.
Thanks for the attention!
Roberto
Here attached, a 3d example dataset:
# x z Density of current[Re]
45 -15 -3.182
60 -15 -19.381
75 -15 43.782
90 -15 22.177
105 -15 -31.759
120 -15 7.625
135 -15 -22.805
150 -15 -42.79
165 -15 38.199
180 -15 52.525
195 -15 60.093
210 -15 21.255
225 -15 -27.971
240 -15 23.938
255 -15 -74.669
270 -15 -98.24
285 -15 -9.867
300 -15 17.152
315 -15 30.979
330 -15 7.646
345 -15 18.718
360 -15 7.106
375 -15 17.715
390 -15 -0.75
405 -15 -2.039
420 -15 -8.993
435 -15 -25.078
90 -30 15.666
105 -30 37.538
120 -30 -40.039
135 -30 -38.69
150 -30 12.286
165 -30 23.028
180 -30 70.035
195 -30 52.557
210 -30 30.41
225 -30 31.778
240 -30 -76.756
255 -30 -63.216
270 -30 -67.121
285 -30 -67.721
300 -30 16.801
315 -30 35.805
330 -30 30.623
345 -30 -5.963
360 -30 22.293
375 -30 20.595
390 -30 25.933
135 -45 -31.461
150 -45 6.237
165 -45 48.476
180 -45 23.709
195 -45 71.678
210 -45 97.939
225 -45 -30.632
240 -45 -66.418
255 -45 -92.299
270 -45 -42.569
285 -45 -34.905
300 -45 -55.369
315 -45 31.355
330 -45 44.209
345 -45 56.547
180 -60 42.678
195 -60 48.503
210 -60 8.012
225 -60 -4.586
240 -60 -32.172
255 -60 -30.978
270 -60 -74.704
285 -60 -38.047
300 -60 -29.84
225 -75 -14.471
240 -75 31.208
255 -75 -0.483
--
Roberto Marzocchi
PhD student in Hydrogeology
SUPSI - Institute of Earth Sciences
Via Trevano C.P. 72, 6952 Canobbio, Switzerland
Phone: +41 58 6666 212 Fax: +41 58 6666 209
E-mail: _roberto.marzocchi@supsi.ch_ <mailto:rob...@su...>
*Please consider the environment before printing this mail note!*
* *
|
