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

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

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

revision 1.2 by adcroft, Mon Mar 22 16:23:28 1999 UTC revision 1.3 by adcroft, Wed May 5 18:32:39 1999 UTC
# Line 36  z=-dz/2:-dz:-H; Line 36  z=-dz/2:-dz:-H;
36  % Temperature profile  % Temperature profile
37  [sprintf('Tref =') sprintf(' %8.6g,',Tz*z-mean(Tz*z))]  [sprintf('Tref =') sprintf(' %8.6g,',Tz*z-mean(Tz*z))]
38    
39    ieee='b';
40    accuracy='real*8';
41    
42  % Gaussian bump  % Gaussian bump
43  h=-H+dh*exp( -(X.^2+Y.^2)/(2*(L^2)) );  h=-H+dh*exp( -(X.^2+Y.^2)/(2*(L^2)) );
44  f77write('topog.bump',h,'real*8','b');  fid=fopen('topog.bump','w',ieee); fwrite(fid,h,accuracy); fclose(fid);
45    
46  % Side walls + bump  % Side walls + bump
47  h(:,1)=0;  h(:,1)=0;
48  h(:,ny)=0;  h(:,ny)=0;
49  f77write('topog.bumpchannel',h,'real*8','b');  fid=fopen('topog.bumpchannel','w',ieee); fwrite(fid,h,accuracy); fclose(fid);
50    
51  % Simple channel  % Simple channel
52  h(:,1)=0;  h(:,1)=0;
53  h(:,2:ny-1)=-H;  h(:,2:ny-1)=-H;
54  h(:,ny)=0;  h(:,ny)=0;
55  f77write('topog.channel',h,'real*8','b');  fid=fopen('topog.channel','w',ieee); fwrite(fid,h,accuracy); fclose(fid);
Legend:
Removed from v.1.2  
changed lines
  Added in v.1.3
  ViewVC Help
Powered by ViewVC 1.1.22