Homec4science
Diffusion GYACOMO (OBSOLETE: LATEST VERSION ON gitlab.epfl.ch) (master) Browse

gyacomo/matlab/load_params.m
0d685e3b68a2master

Last Change
Raw File

load_params.m
View Options

DATA.CO = h5readatt(filename,'/data/input','CO');
% DATA.K_N = h5readatt(filename,'/data/input','eta_n');
% DATA.K_T = h5readatt(filename,'/data/input','eta_T');
DATA.K_N = h5readatt(filename,'/data/input','K_n');
DATA.K_T = h5readatt(filename,'/data/input','K_T');
DATA.K_E = h5readatt(filename,'/data/input','K_E');
DATA.Q0 = h5readatt(filename,'/data/input','q0');
DATA.SHEAR = h5readatt(filename,'/data/input','shear');
DATA.EPS = h5readatt(filename,'/data/input','eps');
DATA.PMAXI = h5readatt(filename,'/data/input','pmaxi');
DATA.JMAXI = h5readatt(filename,'/data/input','jmaxi');
DATA.PMAXE = h5readatt(filename,'/data/input','pmaxe');
DATA.JMAXE = h5readatt(filename,'/data/input','jmaxe');
DATA.NON_LIN = h5readatt(filename,'/data/input','NON_LIN');
DATA.NU = h5readatt(filename,'/data/input','nu');
DATA.Nx = h5readatt(filename,'/data/input','Nx');
DATA.Ny = h5readatt(filename,'/data/input','Ny');
DATA.L = h5readatt(filename,'/data/input','Lx');
DATA.CLOS = h5readatt(filename,'/data/input','CLOS');
DATA.DT_SIM = h5readatt(filename,'/data/input','dt');
DATA.MU = h5readatt(filename,'/data/input','mu');
% MU = str2num(filename(end-18:end-14)); %bad...
DATA.W_GAMMA = h5readatt(filename,'/data/input','write_gamma') == 'y';
DATA.W_PHI = h5readatt(filename,'/data/input','write_phi') == 'y';
DATA.W_NA00 = h5readatt(filename,'/data/input','write_Na00') == 'y';
DATA.W_NAPJ = h5readatt(filename,'/data/input','write_Napj') == 'y';
DATA.W_SAPJ = h5readatt(filename,'/data/input','write_Sapj') == 'y';
if DATA.NON_LIN == 'y'
DATA.NON_LIN = 1;
else
DATA.NON_LIN = 0;
end
switch abs(DATA.CO)
case 0; DATA.CONAME = 'LB';
case 1; DATA.CONAME = 'DG';
case 2; DATA.CONAME = 'SG';
case 3; DATA.CONAME = 'PA';
case 4; DATA.CONAME = 'FC';
otherwise; DATA.CONAME ='UK';
end
if (DATA.CO <= 0); DATA.CONAME = [DATA.CONAME,'DK'];
else; DATA.CONAME = [DATA.CONAME,'GK'];
end
if (DATA.CLOS == 0); DATA.CLOSNAME = 'Trunc.';
elseif(DATA.CLOS == 1); DATA.CLOSNAME = 'Clos. 1';
elseif(DATA.CLOS == 2); DATA.CLOSNAME = 'Clos. 2';
end
if (DATA.PMAXE == DATA.PMAXI) && (DATA.JMAXE == DATA.JMAXI)
degngrad = ['P_',num2str(DATA.PMAXE),'_J_',num2str(DATA.JMAXE)];
else
degngrad = ['Pe_',num2str(DATA.PMAXE),'_Je_',num2str(DATA.JMAXE),...
'_Pi_',num2str(DATA.PMAXI),'_Ji_',num2str(DATA.JMAXI)];
end
degngrad = [degngrad,'_Kn_%1.1f_nu_%0.0e_',...
DATA.CONAME,'_CLOS_',num2str(DATA.CLOS),'_mu_%0.0e'];
degngrad = sprintf(degngrad,[DATA.K_N,DATA.NU,DATA.MU]);
if ~DATA.NON_LIN; degngrad = ['lin_',degngrad]; end
resolution = [num2str(DATA.Nx),'x',num2str(DATA.Ny),'_'];
gridname = ['L_',num2str(DATA.L),'_'];
DATA.PARAMS = [resolution,gridname,degngrad];
c4science · Help