function sphcolor_polar_south(field, lat, lon,num)

% field=squeeze(station(1,:,:))+.001;
% %field=rand(NN,N);
% lat=crut_arctic_lat;
% lon=crut_lon;
% num=5;


%plots spherical data on a robinson projection with continents

if max(lon)<=4
    error('input longitude should be in degrees')
end

if max(abs(lat))>90
    error('input order is wrong - input latitude vector before longitude vector')
end

field=squeeze(field);
a=size(field);
NN=a(1,1); N=a(1,2);

if length(lat)<NN | length(lat)>NN
    error('field size and latitude input must be same size')
end

if length(lon)<N | length(lon)>N
    error(' field size and longitude input must be same size')
end

minlat=min(lat);
minlat=30;

LAT=lat;
load coast
figure(num); clf;
%axesm('stereo','MapLatLimit',[min(LAT) max(LAT)],'Origin',[90 0],'FLatLimit', [60 90])
axesm('stereo','Origin',[-90 0], 'FLonLimit', [-180 180])
axesm('stereo','Origin',[-90 -30], 'MapLatLimit',[minlat 90],'FLonLimit', [-180 180])
latdo=find(lat >= minlat);
%framem; 
plotm(lat,long, 'k');
lat=LAT; clear LAT



[LON,LAT]=meshgrid(lon,lat);
hold on;
color_trans; colormap(cmap)
caxis([-max(max(abs(field))) max(max(abs(field)))])
pcolorm(LAT, LON, field)

colorbar;
%tightmap;
shading interp;
tightmap2;

%fill_landmap('over', 0.5*[1 1 1]);  %makes land grey