/[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.1 by mlosch, Mon Oct 18 15:15:05 2004 UTC revision 1.2 by mlosch, Fri Sep 18 12:33:36 2009 UTC
# Line 80  tke = 1e-6*ones([nx,ny,nz]); Line 80  tke = 1e-6*ones([nx,ny,nz]);
80  tke(2,1,1) = 1e-3;  tke(2,1,1) = 1e-3;
81  fid=fopen('TKE.init','w',ieee); fwrite(fid,tke,prec); fclose(fid);  fid=fopen('TKE.init','w',ieee); fwrite(fid,tke,prec); fclose(fid);
82    
83    % double diffusive initial conditions:
84    % salt fingers in upper layers (dt>0,ds>0,Rrho=talpha*dt/sbeta*ds > 1)
85    % diffusive convection in lower layers (dt<0,ds<0,Rrho=talpha*dt/sbeta*ds < 1)
86    sbeta = 7.e-4;
87    tdd = -min(Tref,2-Tref);
88    Rrho = 0*tdd+1.1;
89    Rrho(8:end) = 0.5;
90    Sz2=talpha*diff(tdd)/sbeta./Rrho(1:end-1);
91    sdd = cumsum([34,Sz2]);
92    
93    t=zeros([nx,ny,nz]);
94    for k=1:nz
95     t(:,:,k) = tdd(k);
96    end
97    fid=fopen('T.doublediff','w',ieee); fwrite(fid,t,prec); fclose(fid);
98    
99    s=zeros([nx,ny,nz]);
100    for k=1:nz
101     s(:,:,k) = sdd(k);
102    end
103    fid=fopen('S.doublediff','w',ieee); fwrite(fid,s,prec); fclose(fid);
104    
105    
Legend:
Removed from v.1.1  
changed lines
  Added in v.1.2
  ViewVC Help
Powered by ViewVC 1.1.22