I am trying to make a figure similar to this one:
which I made using the “old” style GMT (not using the new subplot functionality). I am trying to reproduce this type of plot with the subplot framework and have 3 columns. Here is my attempt:
#!/usr/bin/env bash
EQMERIDIAN="180" # center equatorial view on 180 deg (12 MLT)
SUBPLOT_WIDTH="2i"
SUBPLOT_HEIGHT="1i"
POLAR_HEIGHT="0.5i"
MIN_LAT="60"
gmt set MAP_FRAME_TYPE plain
eq_proj="-JW${EQMERIDIAN}/${SUBPLOT_WIDTH} -Rg"
np_proj="-JG0/90/${POLAR_HEIGHT} -R-180/180/${MIN_LAT}/90"
sp_proj="-JG180/-90/${POLAR_HEIGHT} -R/-180/180/-90/-${MIN_LAT}"
gmt begin figure png
gmt subplot begin 1x3 -Fs${SUBPLOT_WIDTH}/${SUBPLOT_HEIGHT}
gmt subplot set 0,0
gmt basemap ${eq_proj} -B0
gmt basemap ${np_proj} -B0 -Ya1i
#gmt basemap ${sp_proj} -B0 -Xa2i
gmt subplot set 0,1
gmt basemap ${eq_proj} -B0
gmt subplot set 0,2
gmt basemap ${eq_proj} -B0
gmt subplot end
rm -f $XYZFILE $GRDFILE
gmt end show
but unfortunately I get this error:
basemap [ERROR]: Option -Y given more than once
basemap [ERROR]: Option -Y parsing failure. Correct syntax:
-X[a|c|f|r]<xshift> -Y[a|c|f|r]<yshift>
Shift origin of plot to (<xshift>, <yshift>). Prepend r for shift relative to current point [Default],
prepend a for temporary adjustment of origin, prepend f to position relative to lower left corner of
page, prepend c for offset of center of plot to center of page. For overlays (-O), the default setting
is [r0], otherwise [f2.54c].
basemap [ERROR]: Offending option -Ya0i
When I remove the -Y option from the second basemap command, the result is this:
Can anyone suggest how to do this correctly?