gwave_llc_eg_1deg/input/gendata.m

snx=90;
sny=90;

t0=zeros(snx,sny);
tb=t0;

x = [1:snx]-.5-snx/2;
y = [1:sny]-.5-sny/2;

[xx,yy]=meshgrid(x,y);

a = 1;
sigma=10;

tb = a*exp(-.5*(xx.^2+yy.^2)/sigma);

% mit_writefield('ssh.bump',tb,'real*8');
% mit_writefield('ssh.flat',t0,'real*8');

fid = fopen('ssh.bump', 'w', 'ieee-be');
fwrite(fid,tb,'real*8',0,'ieee-be');
fclose(fid);

fid = fopen('ssh.flat', 'w', 'ieee-be');
fwrite(fid,t0,'real*8',0,'ieee-be');
fclose(fid);


1	edhill	1.1	snx=90;
2			sny=90;
3
4			t0=zeros(snx,sny);
5			tb=t0;
6
7			x = [1:snx]-.5-snx/2;
8			y = [1:sny]-.5-sny/2;
9
10			[xx,yy]=meshgrid(x,y);
11
12			a = 1;
13			sigma=10;
14
15			tb = aexp(-.5(xx.^2+yy.^2)/sigma);
16
17			% mit_writefield('ssh.bump',tb,'real*8');
18			% mit_writefield('ssh.flat',t0,'real*8');
19
20			fid = fopen('ssh.bump', 'w', 'ieee-be');
21			fwrite(fid,tb,'real*8',0,'ieee-be');
22			fclose(fid);
23
24			fid = fopen('ssh.flat', 'w', 'ieee-be');
25			fwrite(fid,t0,'real*8',0,'ieee-be');
26			fclose(fid);
27
28