Back to home page

MITgcm

 
 

    


Warning, /verification/adjustment.128x64x1/input/gendata.m is written in an unsupported language. File is not indexed.

view on githubraw file Latest commit aea29c85 on 2001-05-29 14:02:00 UTC
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);