Hello everyone
I am using surface command to interpolate my data and get a smoother contour distribution, but for some reason it doesn’t work well for part of my data. Also, some of the contours are outside of a specific area (red rectangle) and it shouldn’t be like that.
I would appreciate your help 
I attach the input file and the plot. This is the script to create the grid file from xyz and to plot the contours.
ps=slip_dist-interpolate.ps
gmt blockmean File4Interpola.txt -R-76.25/-74.25/8.55/10.8 -I0.01m/0.01m > File4Interpola-bm.xyz
gmt surface File4Interpola-bm.xyz -GFile4Interpola.grd -T0 -I0.01m/0.01m -R -Ll0 -Lu1.0 -Am
gmt pscoast -R-77.5/-73.5/7.5/12 -JM6i -Ba1 -BWSne -N1 -N3 -Df -W -Xc -Yc -Sghostwhite -K > $ps
gmt grdcontour File4Interpola.grd -R -J -C0.1 -S0.01 -A0.2+f10p -W1p,black -O -K >> $ps
gmt psxy -R -JM -SJ -W1.0,red -O << END >> $ps
-75.20112505742531993968 9.57695945571200752958 20.0 200 175
END
gmt psconvert $ps -A -Tt
File4Interpola.txt (43.1 KB)
slip_dist-interpolate.tif (208.7 KB)