Hi everyone
I am trying to get some contours from a xyz file. Then I convert to grd file with xyz2grd and since I obtained a very rough contour distribution I re-sample it to a finer resolution with grdsample. However, I am obtaining some gaps with NaN values in the middle of my plot and I cannot figure what is the problem.
I attach the code and the plot.
I will appreciate your help!
#! /bin/bash
ps=contours.ps
E1=xyz-original.dat
lon0=-78.0
lon1=-74.0
lat0=8.0
lat1=12.0
pscoast -JM6i -R$lon0/$lon1/$lat0/$lat1 -P -Ba2f1/a2f1WSen -N1 -N3 -Df -W -K -V -Xc -Yc -Sghostwhite > $ps
psxy $E1 -R -JM -Sc0.1 -Gblack -O -K >> $ps
mesh=xyz_0.25.grd
I=0.25
xyz2grd $E1 -G$mesh -I$I -R$lon0/$lon1/$lat0/$lat1
grdsample $mesh -Gxyz_0.05.grd -I0.05
grdcontour xyz_0.05.grd -JM -R -L-3/0 -C0.025 -W0.5,blue -K -O >> $ps
grdcontour xyz_0.05.grd -JM -R -L0/3 -C0.025 -A0.05+f11p -W0.5,red -K -O >> $ps
grd2xyz xyz_0.05.grd > xyz_0.05.xyz
psxy xyz_0.05.xyz -R -JM -Sc0.05 -Ggreen -O -K >> $ps
psxy -R -JM -SJ -W1.0,black -O << END >> $ps
-75.826 10.202 20.0 135 110
END
psconvert $ps -A -Tt
contours.ps (222.3 KB)
xyz-original.dat (6.0 KB)