%axfix.m
   
if istart==0
	fac=cos(41.5*pi/180);

	latmin=ming(latfix);
	latmax=maxg(latfix);
	lonmin=ming(lonfix);
	lonmax=maxg(lonfix);
	if (lonmax-lonmin)*fac>(latmax-latmin)
	      lonmax=lonmax+.02;
	      lonmin=lonmin-.02;
	      latave=(latmin+latmax)/2;
	      latmin=latave-(lonmax-lonmin)/2*fac;
	      latmax=latave+(lonmax-lonmin)/2*fac;
	else
	      latmax=latmax+.02*fac;
	      latmin=latmin-.02*fac;
	      lonave=(lonmin+lonmax)/2;
	      lonmin=lonave-(latmax-latmin)/2/fac;
	      lonmax=lonave+(latmax-latmin)/2/fac;
	end
     istart=1;
     axis('square')
     load gaybath.sma
else

hold off
axis([lonmin lonmax latmin latmax])
gaysh2
ax=axis;
lonmin=ax(1); lonmax=ax(2); latmin=ax(3); latmax=ax(4);


hold on
xlabel('longitude')
ylabel('latitude')

end
