The following works pretty much the way you would expect it to:
;- Lambert Azimuthal Equal Area (Great Lakes)
map = map_proj_init(111, center_lat=45, $
center_lon=-85, limit=[40, -95, 50, -75])
;- Default position
window, /free
plot, map.uv_box[[0, 2]], map.uv_box[[1, 3]], $
xstyle=5, ystyle=5, /nodata
gshhs_plot, 'gshhs_i.b', level=2, map=map, noclip=0
map_grid, /label, /box, map=map
;- Default position (isotropic)
window, /free
plot, map.uv_box[[0, 2]], map.uv_box[[1, 3]], $
xstyle=5, ystyle=5, /nodata, /isotropic
gshhs_plot, 'gshhs_i.b', level=2, map=map, noclip=0
map_grid, /label, /box, map=map
;- Position specified
window, /free
plot, map.uv_box[[0, 2]], map.uv_box[[1, 3]], $
xstyle=5, ystyle=5, /nodata, $
position=[0.1, 0.1, 0.6, 0.6]
gshhs_plot, 'gshhs_i.b', level=2, map=map, noclip=0
map_grid, /label, /box, map=map
;- Position specified (isotropic)
window, /free
plot, map.uv_box[[0, 2]], map.uv_box[[1, 3]], $
xstyle=5, ystyle=5, /nodata, $
position=[0.1, 0.1, 0.6, 0.6], /isotropic
gshhs_plot, 'gshhs_i.b', level=2, map=map, noclip=0
map_grid, /label, /box, map=map
|