
  clear;
  close all;
  more  off;

  idir = '/ches/scr10/pierre/out_xsede/run101/';

  for ifil = 1 : 60
    fnam = [ idir 'ocean_avg_' num2str( ifil, '%04u' ) '.nc' ];
    tic  = ncread( fnam, 'TIC' );
    clear fnam;
    tic  = permute( tic, [4, 3, 2, 1] );
    tic  = squeeze( min( squeeze( min( tic ) ) ) );
    tic( find( tic >= 0.    ) ) =   1;
    tic( find( tic <  0.    ) ) = - 1;
    tic( find( isnan( tic ) ) ) =   0;
    tic  = tic * (- 1.);
    figure(1);
      colormap( jet );
      hand = imagesc( tic );
      axis xy; axis equal;
      caxis([- 1, 1]);
      title(['Month ' num2str( ifil, '%02u' ) ', Negative TIC in red']);
    print( '-dpng', [idir 'min_value_of_tic_month_' num2str( ifil, '%02u' ) '.png'] );
  end; clear ifil;
