adjustment.cs-32x32x1/input.nlfs/gendata.m

% This is a matlab script that generates the input data

% $Header:  $
% $Name:  $

nx=32;ny=nx;ntx=6;

h=zeros(nx,ntx,ny);

x=0.25*((0.5:nx)/nx-0.5);
[X,Y]=ndgrid(x,x);
R=sqrt(X.^2+Y.^2);


tileno=1;

o=1;
i=round(nx*1/2);
j=round(ny*1/2);
h(i:i+o,tileno,j:j+o)=1;

h(:,1,:)=0.5+0.5*cos(   pi*min(R,0*R+0.04)/0.04 );

%- PS anomaly = 100.mb (10^4 Pa):
h=h*1.e+4;

fid=fopen('ps100mb.bin','w','b');
fwrite(fid,h,'real*8');
fclose(fid);
1	jmc	1.1	% This is a matlab script that generates the input data
2
3			% $Header: $
4			% $Name: $
5
6			nx=32;ny=nx;ntx=6;
7
8			h=zeros(nx,ntx,ny);
9
10			x=0.25*((0.5:nx)/nx-0.5);
11			[X,Y]=ndgrid(x,x);
12			R=sqrt(X.^2+Y.^2);
13
14
15			tileno=1;
16
17			o=1;
18			i=round(nx*1/2);
19			j=round(ny*1/2);
20			h(i:i+o,tileno,j:j+o)=1;
21
22			h(:,1,:)=0.5+0.5cos( pimin(R,0*R+0.04)/0.04 );
23
24			%- PS anomaly = 100.mb (10^4 Pa):
25			h=h*1.e+4;
26
27			fid=fopen('ps100mb.bin','w','b');
28			fwrite(fid,h,'real*8');
29			fclose(fid);