function mkfig_gridbox

addpath /Users/patrick/ARTS/atmlab/arts/various

DEG2RAD = constants( 'DEG2RAD' );

clf
lat_grid = [85 95];
z_field  = [10 10.1;11 10.9];
%
arts_plot_atmgrids( 2, lat_grid, [], z_field, [0 0]', [] );
hold on

r   = z_field(1,1);
lat = lat_grid(1);
atmplot_pol2cart( r, lat, 'o' );
[x,y] = pol2cart( DEG2RAD * lat, r );
text( x+0.08, y, '(r_{[l+1]}+z_{[k,l+1]},\alpha_{[l+1]})' );

r   = z_field(2,1);
lat = lat_grid(1);
atmplot_pol2cart( r, lat, 'o' );
[x,y] = pol2cart( DEG2RAD * lat, r );
text( x+0.07, y, '(r_{[l+1]}+z_{[k+1,l+1]},\alpha_{[l+1]})' );

r   = z_field(1,2);
lat = lat_grid(2);
atmplot_pol2cart( r, lat, 'o' );
[x,y] = pol2cart( DEG2RAD * lat, r );
h = text( x-0.07, y, '(r_{[l]}+z_{[k,l]},\alpha_{[l]})' );
set( h, 'HorizontalAlignment', 'right')

r   = z_field(2,2);
lat = lat_grid(2);
atmplot_pol2cart( r, lat, 'o' );
[x,y] = pol2cart( DEG2RAD * lat, r );
h = text( x-0.07, y, '(r_{[l]}+z_{[k+1,l]},\alpha_{[l]})' );
set( h, 'HorizontalAlignment', 'right')


axis equal
axis([ -1.9 2.2 9.8  11.15 ])
axes_frame( gca, 'off' )

hold off

return

if yes_or_no('Print figures')
  print ppath_3Dsearch.eps -depsc
  ! epstopdf ppath_3Dsearch.eps
  ! rm ppath_3Dsearch.eps
  figure(3)
  print ppath_ex3.eps -depsc
  ! epstopdf ppath_ex3.eps
  ! rm ppath_ex3.eps
  figure(2)
  print ppath_ex2.eps -depsc
  ! epstopdf ppath_ex2.eps
  ! rm ppath_ex2.eps
  figure(1)
  print ppath_ex1.eps -depsc
  ! epstopdf ppath_ex1.eps
  ! rm ppath_ex1.eps
end