adjustment.128x64x1/input/gendata.m

aea29c8517 Alis*0001 nx=128;ny=64;                                   % Grid dimensions
                0002
                0003 radi=30;                                        % Radius of patch
                0004 deltaPs=100;                                    % Amplitude
                0005
                0006 x=(0.5:nx)/nx*360;                              %  0  < x < 360
                0007 y=(0.5:ny)/ny*180-90;                           % -90 < y < +90
                0008 [X,Y]=ndgrid(x,y);
                0009
                0010 h=deltaPs*0.5*(1-cos(max(Y-90+radi,0*Y)/radi*pi));      % Horizontal structure
                0011 h=h*deltaPs/max(max(abs(h)));           % Re-scale to give max val
                0012
                0013 fid=fopen('ps.init','w','b');
                0014 fwrite(fid,h,'real*8');
                0015 fclose(fid);