function sphcolor_polar_cont_south(field,FIELD, mlat, mlon,num)

%field is the color;
%FILD is the contour
% 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(mlon)<=4
    error('input longitude should be in degrees')
end

if max(abs(mlat))>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(mlat)<NN | length(mlat)>NN
    error('field size and latitude input must be same size')
end

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

minlat=min(mlat);
minlat=30;

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]);
[LON,LAT]=meshgrid(mlon,mlat); hold on;


latdo=find(mlat >= minlat);
%framem; 
plotm(lat,long, 'k');


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

colorbar;
%tightmap;
shading interp;
tightmap2;

LWc=1.8;
FSl=18; %labeling fontsize
FN='Arial';

hold on;
vals1=[-10 -8 -6 -4 -2]/4;
vals2=[2 4 6 8 10]/4;
[c,h]=contourm(LAT, LON, FIELD, vals1, 'm--', 'LineWidth', LWc);
ht = clabelm(c,h); set(ht,'color','m', 'FontSize', FSl, 'FontName', FN );
hold on; 
[c,h]=contourm(LAT, LON, FIELD, vals2, 'g--', 'LineWidth', LWc);
ht = clabelm(c,h); set(ht,'color','g','FontSize', FSl, 'FontName', FN );

hidem(gca);
plotm(lat,long, 'k', 'LineWidth', 1.9);
Box On;