Grid to GeoTIF
var=radius.ll
new=big3
argR="-R280/287/-9/-3"
argR="-R0/360/-90/90"
argI="-I0.04"
date
echo "Blockmean"
gmt blockmean $var $argR $argI > out
echo "sphereinterpolate"
cat out | gmt sphinterpolate $argR $argI -Gtest.nc
echo "gdal translate"
gdal_translate -of GTiff -b 1 -a_srs 60300.prj NETCDF:test.nc $new.tif
Global
Topography
Sigmas
Number of Images
Slope
Albedo
Regional
Topography
Sigmas
Number of Images
Slope
Albedo
Make Pretties
gmt begin GMT_cont
gmt set GMT_THEME cookbook
gmt grdcontour test.nc
gmt end show
#gmt grdcontour test.nc -C10 -A50
gmt begin GMT_img
gmt set GMT_THEME cookbook
#gmt makecpt -Crainbow
gmt grdimage test.nc -JM6i -B -BWSnE
gmt colorbar -DJTC -Bxa -By+lm
gmt end show
gmt begin GMT_img
gmt makecpt -Crainbow
gmt set GMT_THEME cookbook
gmt grdimage test.nc $argR -JM6i -B -BWSnE
#gmt colorbar -DJTC -I0.4 -Bxa -By+lm
gmt colorbar -DJTC -Bxa -By+lm
gmt end show