function sphcolor_polar_hatch(field,unc, lat, lon,num)

% field=squeeze(station(1,:,:))+.001;
% unc= uncertainity in the field
% lat=crut_arctic_lat;
% lon=crut_lon;
% num=5;


%plots spherical data on a Polar projection with continents and a contour
%to show regions of significance

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);

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;
%tightmap2;
MS=12;
hold on;

%DO=find((field-unc) > 0); %find gridpoints where the field is significant
%- considers sign and magnitude
DO=find((abs(field)-unc) > 0); %find gridpoints where the field is significant - only cares about magnitude
plotm(LAT(DO), LON(DO), 'm+', 'MarkerSize', MS);
%ht = clabelm(c,h); 
%set(ht,'color','k')
%hidem(gca)

load coast
hold on; plotm(lat,long, 'k');
%tightmap2


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