FILE2 = ...
['MLD2001_f_01.mat';...
'MLD2001_f_02.mat';...
'MLD2001_f_03.mat';...
'MLD2001_f_04.mat';...
'MLD2001_f_05.mat';...
'MLD2001_f_06.mat';...
'MLD2001_f_07.mat';...
'MLD2001_f_08.mat';...
'MLD2001_f_09.mat';...
'MLD2001_f_10.mat'];

for i = 1:10
  eval(['load ' FILE2(i,:) ])
  ic = find(x(1,:)>=39.8672 & x(1,:)<=100.1602); ir = find(y(:,1)>=-30 & y(:,1)<=30); 
  x_plot=[x(ir,ic)];
  y_plot=[y(ir,ic)];
  new_mld = replace(mixeddepth(ir,ic),1e10,NaN);
  mld_mean(i) = mean(new_mld(isfinite(new_mld)));
  subplot(5,2,i)
  contourf(x_plot,y_plot,new_mld)
end

