%cmin=21;
cmin=18;
%cmax=27;
%cmax=30;
cmax=32;
V=[cmin cmax];
caxis('manual')

% make land mask
XCenter=Nx/2;
YCenter=Ny/2;
R = 58 ;
% create single circular mask
mask1=ones(Nx,Ny); 
mask2=zeros(Nx,Ny); 
for i=2:Nx-1
  for j=2:Ny-1
    DispX=XCenter-i;
    DispY=YCenter-j;
      if ( DispX^2 + DispY^2 >= R^2 ) 
        mask1(i,j)=0;
        mask2(i,j)=1;
      end
   end
end
%%        txy=reshape(t(:,:,2),[Nx,Ny]);
%%        txy=reshape(t(:,:,15),[Nx,Ny]);
%        txy=reshape(t(:,:,1),[Nx,Ny]);
%        txy(find(txy==0))=txy(find(txy==0))*NaN;
%%        txy(find(mask==1))=txy(find(txy==1))*0;
%        txy=txy.*mask1+(mask2*19);
%        imagesc(lat,long,txy');shading('flat');
%%        axis image;caxis(V);colorbar;colormap(flipud(gray))
%        axis image;colorbar;

%contourf(t(40:84,40:84,1)',5);caxis(V);axis square;colorbar
%imagesc(t(40:Nx-40,40:Ny-40,1)');caxis(V);axis square;colorbar
imagesc(t(:,:,2)');caxis(V);axis square;colorbar
%imagesc(squeeze(t(:,Nx/2,:))');caxis(V);axis square;colorbar
%contourf(flipud(squeeze(t(:,Nx/2,:))'),10);caxis(V);axis square;colorbar