vermix/input/gendata.m

% This is a matlab script that generates the input data
% variable x resolution

prec='real*8';
ieee='b';

% Dimensions of grid
nx=1;
ny=1;
nz=26;
% Vertical grid (meters)
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];
sprintf('delR = %d * %7.6g,',nz,dz)
zi = [0 cumsum(dz)];
z = 0.5*(zi(1:end-1)+zi(2:end));

% Initial temperature
gravity=9.81;
talpha=2.0e-4;
N2=2.e-5;
Tz=N2/(gravity*talpha);

Tref=-Tz*(z - zi(end))+2;
[sprintf('Tref =') sprintf(' %8.6g,',Tref)];
fid=fopen('T_26.init','w',ieee); fwrite(fid,Tref,prec); fclose(fid);

% Flux
Q=zeros(72,1);
Q(1:18) = 350;
Q(25:66) = -350*3/7;
fid=fopen('Qnet_72.forcing','w',ieee); fwrite(fid,Q,prec); fclose(fid);

%
taux = 0.1*ones([nx,ny],72);
fid=fopen('taux_72.forcing','w',ieee); fwrite(fid,taux,prec); fclose(fid);

%%%%%
tke = 1e-6*ones(1,nz);
tke(1,1) = 1e-3;
fid=fopen('TKE.init','w',ieee); fwrite(fid,tke,prec); fclose(fid);

% double diffusive initial conditions: 
% salt fingers in upper layers (dt>0,ds>0,Rrho=talpha*dt/sbeta*ds > 1)
% diffusive convection in lower layers (dt<0,ds<0,Rrho=talpha*dt/sbeta*ds < 1)
sbeta = 7.e-4;
tdd = max(Tref,8-Tref);
Rrho = 0*tdd+1.1;
Rrho(8:end) = 0.5;
Sz2=talpha*diff(tdd)/sbeta./Rrho(1:end-1);
sdd = cumsum([34,Sz2]);

fid=fopen('T.doublediff','w',ieee); fwrite(fid,tdd,prec); fclose(fid);

fid=fopen('S.doublediff','w',ieee); fwrite(fid,sdd,prec); fclose(fid);
1	% This is a matlab script that generates the input data
2	% variable x resolution
3
4	prec='real*8';
5	ieee='b';
6
7	% Dimensions of grid
8	nx=1;
9	ny=1;
10	nz=26;
11	% 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];
13	sprintf('delR = %d * %7.6g,',nz,dz)
14	zi = [0 cumsum(dz)];
15	z = 0.5*(zi(1:end-1)+zi(2:end));
16
17	% Initial temperature
18	gravity=9.81;
19	talpha=2.0e-4;
20	N2=2.e-5;
21	Tz=N2/(gravity*talpha);
22
23	Tref=-Tz*(z - zi(end))+2;
24	[sprintf('Tref =') sprintf(' %8.6g,',Tref)];
25	fid=fopen('T_26.init','w',ieee); fwrite(fid,Tref,prec); fclose(fid);
26
27	% Flux
28	Q=zeros(72,1);
29	Q(1:18) = 350;
30	Q(25:66) = -350*3/7;
31	fid=fopen('Qnet_72.forcing','w',ieee); fwrite(fid,Q,prec); fclose(fid);
32
33	%
34	taux = 0.1*ones([nx,ny],72);
35	fid=fopen('taux_72.forcing','w',ieee); fwrite(fid,taux,prec); fclose(fid);
36
37	%%%%%
38	tke = 1e-6*ones(1,nz);
39	tke(1,1) = 1e-3;
40	fid=fopen('TKE.init','w',ieee); fwrite(fid,tke,prec); fclose(fid);
41
42	% double diffusive initial conditions:
43	% salt fingers in upper layers (dt>0,ds>0,Rrho=talphadt/sbetads > 1)
44	% diffusive convection in lower layers (dt<0,ds<0,Rrho=talphadt/sbetads < 1)
45	sbeta = 7.e-4;
46	tdd = max(Tref,8-Tref);
47	Rrho = 0*tdd+1.1;
48	Rrho(8:end) = 0.5;
49	Sz2=talpha*diff(tdd)/sbeta./Rrho(1:end-1);
50	sdd = cumsum([34,Sz2]);
51
52	fid=fopen('T.doublediff','w',ieee); fwrite(fid,tdd,prec); fclose(fid);
53
54	fid=fopen('S.doublediff','w',ieee); fwrite(fid,sdd,prec); fclose(fid);