= 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 = Global ICQ resolution to bin sizes [degrees] || ICQ Q Value || Vertices || Bin Size [deg] || || 512 || 1,579,014 || 0.041 || || 256 || 396,294 || 0.164 || || 128 || 99,846 || 0.649 || || 64 || 25,350 || 2.556 || || 32 || 6,534 || 9.917 || == Topography == {{{ /bin/cp -f SHAPE.TXT tmp vi tmp 1,$s/D/E/g ZZ convertLatLon tmp > radius.ll }}} == Sigmas == {{{ convertLatLon tmp 1 > sigma.ll }}} == Number of Images == == Slope == == Albedo == {{{ /bin/cp -f SHAPEA.TXT tmp vi tmp 1,$s/D/E/g ZZ convertLatLon tmp 1 > albedo.ll }}} -------- = Regional = Run backplanesGMT (which is in ORExSPCsupport) (or backplanesISIS) == Topography == == Sigmas == == Number of Images == == Slope == == Albedo == == Photometric Data == This includes emission angle, incidence angle and phase angle. Run phasei (located in ORExSPCsupport) -------- = 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 }}}