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

Diff of /MITgcm/verification/front_relax/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 Apr 10 16:14:07 2000 UTC revision 1.3 by jmc, Thu Jun 21 01:07:17 2007 UTC
# Line 52  for k=1:nz, Line 52  for k=1:nz,
52               *exp(-(3*zc(k)/Ho)^2);               *exp(-(3*zc(k)/Ho)^2);
53  end  end
54  T(:,end,:)=0;  T(:,end,:)=0;
55  fid=fopen('theta.bin','w',ieee); fwrite(fid,T,prec); fclose(fid);  %- add 10 more levels:
56    t25=zeros(nx,ny,nz+10); t25(:,:,[1:nz]=T;
57    fid=fopen('theta.bin','w',ieee); fwrite(fid,t25,prec); fclose(fid);
58    
59  % Salt (passive tracer)  % Salt (passive tracer)
60  for k=1:nz,  for k=1:nz,
# Line 62  for k=1:nz, Line 64  for k=1:nz,
64   S(:,:,k)=exp(-(2*Y/Ly).^2);           % Exponential with y   S(:,:,k)=exp(-(2*Y/Ly).^2);           % Exponential with y
65  end  end
66  S(:,end,:)=0;  S(:,end,:)=0;
67  fid=fopen('salt.bin','w',ieee); fwrite(fid,S,prec); fclose(fid);  %- add 10 more levels:
68    s25=zeros(nx,ny,nz+10); s25(:,:,[1:nz]=S;
69    fid=fopen('salt.bin','w',ieee); fwrite(fid,s25,prec); fclose(fid);
Legend:
Removed from v.1.2  
changed lines
  Added in v.1.3
  ViewVC Help
Powered by ViewVC 1.1.22