257 |
end |
end |
258 |
end |
end |
259 |
|
|
260 |
|
ground(1:nx,1:ny)=-990; |
261 |
icetopo2(:,1)=-1000; |
icetopo2(:,1)=-1000; |
262 |
|
fid = fopen('shelftopo.ground.bin','w','b'); fwrite(fid,ground,'real*8'); fclose(fid); |
263 |
fid = fopen('shelftopo.round.bin','w','b'); fwrite(fid,icetopo2,'real*8'); fclose(fid); |
fid = fopen('shelftopo.round.bin','w','b'); fwrite(fid,icetopo2,'real*8'); fclose(fid); |
264 |
fid = fopen('etainit.round.bin','w','b'); fwrite(fid,etainit,'real*8'); fclose(fid); |
fid = fopen('etainit.round.bin','w','b'); fwrite(fid,etainit,'real*8'); fclose(fid); |