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

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

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

revision 1.3 by dfer, Tue Nov 10 02:46:39 2009 UTC revision 1.4 by dfer, Thu Nov 12 16:08:42 2009 UTC
# Line 9  nx=1; Line 9  nx=1;
9  ny=1;  ny=1;
10  nz=26;  nz=26;
11  % Vertical grid (meters)  % Vertical grid (meters)
12  dz =[10 10 10 10 10 12 12 15 15 18 18 22 22 26 32 38 46 46 46 46 46 46 46 46 46 46];  dz =[10 10 10 10 10 11 12 14 16 18 21 24 27 31 35 40 40 40 40 40 40 40 40 40 40 40];
13  sprintf('delR = %d * %7.6g,',nz,dz)  sprintf('delR = %d * %7.6g,',nz,dz)
14  zi = [0 cumsum(dz)];  zi = [0 cumsum(dz)];
15  z = 0.5*(zi(1:end-1)+zi(2:end));  z = 0.5*(zi(1:end-1)+zi(2:end));
# Line 21  N2=2.e-5; Line 21  N2=2.e-5;
21  Tz=N2/(gravity*talpha);  Tz=N2/(gravity*talpha);
22    
23  Tref=-Tz*(z - zi(end))+2;  Tref=-Tz*(z - zi(end))+2;
24  [sprintf('Tref =') sprintf(' %8.6g,',Tref)];  sprintf('Tref ='), sprintf(' %8.6g,',Tref)
25  fid=fopen('T_26.init','w',ieee); fwrite(fid,Tref,prec); fclose(fid);  fid=fopen('T_26.init','w',ieee); fwrite(fid,Tref,prec); fclose(fid);
26    
27  % Flux  % Flux
Legend:
Removed from v.1.3  
changed lines
  Added in v.1.4
  ViewVC Help
Powered by ViewVC 1.1.22