REM Ehsan Momeni

​

set filename=Ehsan_Momeni.ps
set input_grid=Total_Dem_Iceland.nc
set reformatted_grid=reformateed_Total_Dem_Iceland.nc
set renge=335/347/63/67
set scale=M6.5i
set cpt=GMT_wysiwyg.cpt

gmt grdreformat %input_grid%=ns/1/0/NaN %reformatted_grid%=ns

gmt grd2cpt %reformatted_grid% -Z -C%cpt% > NewCpt.cpt
gmt psbasemap -R%renge% -J%scale% -Bxa2f1 -Bya2f1 -BWESN+t"Iceland Topography" -X2.5i -Y1.8i -V -K > %filename%

gmt grdgradient %reformatted_grid% -A45/270 -Greformatted_grid.int -Ne0.6 -fg -V
gmt grdimage %reformatted_grid% -R%renge% -J%scale% -CNewCpt.cpt -Ireformatted_grid.int -O -K  -V >> %filename%

gmt pscoast -R%renge% -J%scale% -Df -B -W0.01p,255 -P -O -Sblue -K -V -Lf344/63.5/65/150+u -N1/1.2p,0 -W1p >> %filename%
gmt psscale -CNewCpt.cpt -D8/-1.5/17/0.15h  -Bxf100a1000+l"Elevation" -By+lm -I -V -O >> %filename%
 
 
del *.int *.gmt *.xml gmt.conf gmt.history NewCpt.cpt 
%filename%