#!/bin/csh -f
#
# you can ask for grid data at onest30@gmail.com
#subject of the message must be "Help retrieving grid data from your web"

# PAper size
gmtset PAPER_MEDIA letter
#  ---- Projection. You may change it if you know how to
set proj = -JM0/0/18.00c
#  ---- Frame annotations. You may change it if you know how to
set frm = -B2/1WSen
#  ---- Map limits. You may change it if you know how to
set lim = -R-93.7/-78/6.9/17.0666667
# 
#  ---- Longitude annotation style. The +ddd:mm:ss form => [0;360] range 
set deg_form=--PLOT_DEGREE_FORMAT=ddd:mm:ss

set grd = nicaraguamirone2.grd
set cpt = nicrelief3.cpt
set ps = nicrelief3.ps
 
# -------- Start by creating the basemap frame
psbasemap $lim $proj $frm -X2.0c -Y5.5c -P $deg_form -K > $ps

# -------- Compute the illumination grid
grdgradient $grd -E0/30/0.55/0.6/0.4/10 -Gnicrelief3_intens.grd=cf
 
# -------- Plot the the base image using grdimage & illumination
grdimage $grd -R -J -C$cpt -Inicrelief3_intens.grd=cf -O >> $ps
convert $ps gmt32.jpg

    Source: geocities.com/mx/onest302

               ( geocities.com/mx)