%********** INPUTS ***********************************
%SELECT GRID SIZE: N=NUMBER OF GRID NODES
N=128; %For efficiency (FFTs) N should be a power of 2
%SELECT GRID SPACING, CORRELATION LENGTH,
%AND SURFACE HEIGHT VARIANCE, ALL IN METERS
%(NOTE: delx should be smaller than lc for )
%( good surface representation)
lambda=10./3.; %radiation wavelength (m)
lc =1.000*lambda; %transverse correlation length (m)
zdev=0.100*lambda; %surface height variance (m)
delx=lambda/7.; %cell size (grid spacing) (m)
%********** END INPUTS ********************************
%CREATE 2-D GAUSSIAN FILTER
%normalize correlation length to cel1 size
lc=lc/delx;
%x and y are dimensionless for now
x=linspace(-(N-1)/2,(N-1)/2,N);y=x;
%x and y coords.