67 |
% list_param={'XC','YC','RAC','RC'}; |
% list_param={'XC','YC','RAC','RC'}; |
68 |
% grid_load(dirGrid,5,list_param); |
% grid_load(dirGrid,5,list_param); |
69 |
% gcmfaces_bindata; |
% gcmfaces_bindata; |
70 |
MYBASININDEX=convert2array(v4_read_bin(fileBasin,1,0)); |
MYBASININDEX=convert2array(read_bin(fileBasin,1,0)); |
71 |
grid.name=model; |
grid.name=model; |
72 |
grid.mygrid=mygrid; |
grid.mygrid=mygrid; |
73 |
grid.mytri=mytri; |
grid.mytri=mytri; |
76 |
% read T/S Atlas |
% read T/S Atlas |
77 |
fldT=mygrid.mskC; fldT(:)=0; fldS=fldT; |
fldT=mygrid.mskC; fldT(:)=0; fldS=fldT; |
78 |
for tt=1:12; |
for tt=1:12; |
79 |
fldT(:,:,:,tt)=v4_read_bin(fileT,tt).*mygrid.mskC; |
fldT(:,:,:,tt)=read_bin(fileT,tt).*mygrid.mskC; |
80 |
fldS(:,:,:,tt)=v4_read_bin(fileS,tt).*mygrid.mskC; |
fldS(:,:,:,tt)=read_bin(fileS,tt).*mygrid.mskC; |
81 |
end; |
end; |
82 |
atlas.name=model; |
atlas.name=model; |
83 |
atlas.T={convert2array(fldT)}; atlas.S={convert2array(fldS)}; |
atlas.T={convert2array(fldT)}; atlas.S={convert2array(fldS)}; |