/[MITgcm]/MITgcm/verification/exp1/input/gendata.m
ViewVC logotype

Diff of /MITgcm/verification/exp1/input/gendata.m

Parent Directory Parent Directory | Revision Log Revision Log | View Revision Graph Revision Graph | View Patch Patch

revision 1.1 by adcroft, Mon Dec 14 23:53:23 1998 UTC revision 1.2 by adcroft, Wed May 5 18:32:36 1999 UTC
# Line 1  Line 1 
1    ieee='b';
2    accuracy='real*8';
3    
4  Ho=2000;  Ho=2000;
5  nx=60;  nx=60;
6  ny=60;  ny=60;
# Line 7  h=-Ho*ones(nx,ny); Line 10  h=-Ho*ones(nx,ny);
10  % Walls  % Walls
11  h(end,:)=0;  h(end,:)=0;
12  h(:,end)=0;  h(:,end)=0;
13  f77write('topog.box',h,'real*8','b');  fid=fopen('topog.box','w',ieee); fwrite(fid,h,accuracy); fclose(fid);
14    
15  % Wind-stress  % Wind-stress
16  tauMax=0.1;  tauMax=0.1;
# Line 15  x=((1:nx)-0.5)/(nx-1); % nx-1 accounts f Line 18  x=((1:nx)-0.5)/(nx-1); % nx-1 accounts f
18  y=((1:ny)-0.5)/(ny-1); % ny-1 accounts for a solid wall  y=((1:ny)-0.5)/(ny-1); % ny-1 accounts for a solid wall
19  [X,Y]=ndgrid(x,y);  [X,Y]=ndgrid(x,y);
20  tau=tauMax*sin(pi*Y);  tau=tauMax*sin(pi*Y);
21  f77write('windx.sin_y',tau,'real*8','b');  fid=fopen('windx.sin_y','w',ieee); fwrite(fid,tau,accuracy); fclose(fid);
Legend:
Removed from v.1.1  
changed lines
  Added in v.1.2
  ViewVC Help
Powered by ViewVC 1.1.22