fid=fopen('ecco_oasim_edp_below.bin','r','b');
edp=fread(fid,'float32');
edp=reshape(edp,360,160,13,12);

fid=fopen('ecco_oasim_esp_below.bin','r','b');
esp=fread(fid,'float32');
esp=reshape(esp,360,160,13,12);

etot=edp+esp;

edpn=edp/2;
espn=esp+edp/2;

fid=fopen('ecco_oasim_edpesp_below.bin','w','b');
for it=1:12,
  for il=1:13,
   clear tmp, tmp=etot(:,:,il, it);
   fwrite(fid,tmp,'float32');
  end % for il
end

fid=fopen('ecco_oasim_edp_half_below.bin','w','b');
for it=1:12,
  for il=1:13,
   clear tmp, tmp=edpn(:,:,il, it);
   fwrite(fid,tmp,'float32');
  end % for il
end
fid=fopen('ecco_oasim_esp_halfedp_below.bin','w','b');
for it=1:12,
  for il=1:13,
   clear tmp, tmp=espn(:,:,il, it);
   fwrite(fid,tmp,'float32');
  end % for il
end



figure
for it=1:12
  subplot(4,3,it), pcolor(squeeze(sum(edp(:,:,:,it),3))'); caxis([0 100]);  shading flat, colorbar
end 
figure
for it=1:12
  subplot(4,3,it), pcolor(squeeze(sum(esp(:,:,:,it),3))');  caxis([0 100]); shading flat, colorbar
end
figure
for it=1:12
  subplot(4,3,it), pcolor(squeeze(sum(etot(:,:,:,it),3))');  caxis([0 200]); shading flat, colorbar
end

