ctlf=rnt_ctl('ias20-sst.nc','sst_time');
sst=rnt_loadvar(ctlf,1:length(ctlf.time),'SST');

[ssta, sstm]=RemoveSeas(sst,ctlf.month);

for imon = 1:12
   in=find(ctlf.month == imon);
   sstSTD(:,:,imon) = stdNaN(ssta(:,:,in),3);
end

djf = [1 2 12]';
mam = [3 4 5]';
amj = [6 7 8];
aut = [9 10 11];

subplot(2,2,1)
  rnt_plcm(mean(sstm(:,:,djf),3),grd)
  title 'win'
subplot(2,2,2)
  rnt_plcm(mean(sstm(:,:,mam),3),grd)
  title 'spr'
subplot(2,2,3)
  rnt_plcm(mean(sstm(:,:,amj),3),grd)
  title 'sum'
subplot(2,2,4)
  rnt_plcm(mean(sstm(:,:,aut),3),grd)
  title 'aut'
rnt_font

subplot(2,2,1)
  rnt_plcm(mean(sstSTD(:,:,djf),3),grd)
  title 'win'
subplot(2,2,2)
  rnt_plcm(mean(sstSTD(:,:,mam),3),grd)
  title 'spr'
subplot(2,2,3)
  rnt_plcm(mean(sstSTD(:,:,amj),3),grd)
  title 'sum'
subplot(2,2,4)
  rnt_plcm(mean(sstSTD(:,:,aut),3),grd)
  title 'aut'
rnt_font