--- MITgcm/verification/exp1/input/gendata.m	1998/12/14 23:53:23	1.1
+++ MITgcm/verification/exp1/input/gendata.m	1999/05/05 18:32:36	1.2
@@ -1,3 +1,6 @@
+ieee='b';
+accuracy='real*8';
+
 Ho=2000;
 nx=60;
 ny=60;
@@ -7,7 +10,7 @@
 % Walls
 h(end,:)=0;
 h(:,end)=0;
-f77write('topog.box',h,'real*8','b');
+fid=fopen('topog.box','w',ieee); fwrite(fid,h,accuracy); fclose(fid);
 
 % Wind-stress
 tauMax=0.1;
@@ -15,4 +18,4 @@
 y=((1:ny)-0.5)/(ny-1); % ny-1 accounts for a solid wall
 [X,Y]=ndgrid(x,y);
 tau=tauMax*sin(pi*Y);
-f77write('windx.sin_y',tau,'real*8','b');
+fid=fopen('windx.sin_y','w',ieee); fwrite(fid,tau,accuracy); fclose(fid);