close all
clear all

%Baseline 
cd('/data/scihub-users/jamrda26/mitgcm/cases/GROUNDING_FULL2/input');



% use streamice generated thickness
 fid = fopen('vfacemask.box','r','b'); vfm=fread(fid,[3 200],'real*8'); fclose(fid);
 
 % use streamice generated thickness
 fid = fopen('ufacemask.box','r','b'); ufm=fread(fid,[3 200],'real*8'); fclose(fid);
 
 % use streamice generated thickness
 fid = fopen('hinit_dynamic.box','r','b'); hin=fread(fid,[3 200],'real*8'); fclose(fid);
 
 % use streamice generated thickness
 fid = fopen('vdirich.box','r','b'); vdi=fread(fid,[3 200],'real*8'); fclose(fid);
 
 % use streamice generated thickness
 fid = fopen('HBCy.box','r','b'); hby=fread(fid,[3 200],'real*8'); fclose(fid);
 
 % use streamice generated thickness
 fid = fopen('hmask.box','r','b'); hma=fread(fid,[3 200],'real*8'); fclose(fid);
 
 figure
 imagesc(vfm)
 colorbar
  figure
 imagesc(ufm)
 colorbar
 figure
 imagesc(hin)
 colorbar
 figure
 imagesc(vdi)
 colorbar
 figure
 imagesc(hby)
 colorbar
 figure
 imagesc(hma)
 colorbar
 
%  HM(1:3,1:200)=0;
%  HM(1:3,1:100)=1;
% 
%  fid = fopen('hmask.box','w','b'); fwrite(fid,HM,'real*8'); fclose(fid);
%