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

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

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

revision 1.1 by adcroft, Thu Jun 29 19:29:04 2000 UTC revision 1.2 by adcroft, Thu May 22 19:59:28 2003 UTC
# Line 6  ieee='b'; Line 6  ieee='b';
6    
7  % Dimensions of grid  % Dimensions of grid
8  nx=60;  nx=60;
9  ny=3;  ny=1;
10  nz=20;  nz=20;
11  % Nominal depth of model (meters)  % Nominal depth of model (meters)
12  H=200.0;  H=200.0;
# Line 49  end Line 49  end
49  fid=fopen('T.init','w',ieee); fwrite(fid,t,prec); fclose(fid);  fid=fopen('T.init','w',ieee); fwrite(fid,t,prec); fclose(fid);
50    
51  % Sloping channel  % Sloping channel
52  slope=0.0  slope=0.03
53  offset=2.5e3;  offset=2.5e3;
54  dmax=-40.0;  dmax=-40.0;
55  d=0.0*rand([nx,ny]);  d=0.0*rand([nx,ny]);
# Line 65  d(i,j) = dmax; Line 65  d(i,j) = dmax;
65  end  end
66  end  end
67  end  end
68  %d(nx,:)=0.0;  d(nx,:)=0.0;
69  fid=fopen('topog.slope','w',ieee); fwrite(fid,d,prec); fclose(fid);  fid=fopen('topog.slope','w',ieee); fwrite(fid,d,prec); fclose(fid);
70  plot(x,d(:,2))  plot(x,d(:,1))
71    
72  fid=fopen('delXvar','w',ieee); fwrite(fid,dx,prec); fclose(fid);  fid=fopen('delXvar','w',ieee); fwrite(fid,dx,prec); fclose(fid);
73    
# Line 98  end Line 98  end
98  end  end
99  %d1(end-1:end,:)=d1(1:2,:); % debug by aja  %d1(end-1:end,:)=d1(1:2,:); % debug by aja
100  fid=fopen('topog.convex','w',ieee); fwrite(fid,d1,prec); fclose(fid);  fid=fopen('topog.convex','w',ieee); fwrite(fid,d1,prec); fclose(fid);
101  hold on  plot(x,d1(:,1),'g')
 plot(x,d1(:,2),'g')  
102    
103  %convex slope  %convex slope
104  d2=zeros(nx,ny);  d2=zeros(nx,ny);
# Line 126  end Line 125  end
125  %d2(end-1:end,:)=d2(1:2,:); % debug by aja  %d2(end-1:end,:)=d2(1:2,:); % debug by aja
126  fid=fopen('topog.concave','w',ieee); fwrite(fid,d2,prec); fclose(fid);  fid=fopen('topog.concave','w',ieee); fwrite(fid,d2,prec); fclose(fid);
127  hold on  hold on
128  plot(x,d2(:,2),'r')  plot(x,d2(:,1),'r')
129    hold off
130    
131    
132  fid=fopen('delXvar','w',ieee); fwrite(fid,dx,prec); fclose(fid);  fid=fopen('delXvar','w',ieee); fwrite(fid,dx,prec); fclose(fid);
Legend:
Removed from v.1.1  
changed lines
  Added in v.1.2
  ViewVC Help
Powered by ViewVC 1.1.22