Homec4science
Diffusion FENNECS: Finite Element Non-Neutral Electron Cloud Simulator (geom_testmumps) Browse

fennecs/matlab/load_espic2d.m
bde8e520bdf2geom_testmumps

Last Change
Raw File

load_espic2d.m
View Options

function Data = load_espic2d( file, M, type)
% Load Simulation data from espic2d simulation
if(nargin==1)
type='all';
else
filedata=dir(file);
timestamp=filedata.date;
try
if (strcmp(M.file,file) && ~max(timestamp > M.timestamp) && strcmp(M.type,'all'))
Data=M;
return
end
catch
end
if(nargin~=3)
type='all';
end
end
Data.type=type;
Data.file=file;
% Physical constants
Data.vlight=299792458;
Data.qe=1.60217662E-19;
Data.me=9.109383E-31;
Data.eps_0=8.85418781762E-12;
Data.job_time = hdf5read(file,'/data/input.00/','job_time');
Data.extra_time = hdf5read(file,'/data/input.00/','extra_time');
Data.dt = hdf5read(file,'/data/input.00/','dt');
Data.tmax= hdf5read(file,'/data/input.00/','tmax');
Data.nrun = hdf5read(file,'/data/input.00/','nrun');
Data.nlres = hdf5read(file,'/data/input.00/','nlres');
Data.nlsave = hdf5read(file,'/data/input.00/','nlsave');
Data.nlres = hdf5read(file,'/data/input.00/','nlres');
Data.nz = hdf5read(file,'/data/input.00/','nz');
Data.lz = hdf5read(file,'/data/input.00/lz');
Data.nplasma = hdf5read(file,'/data/input.00/','nplasma');
Data.potinn = hdf5read(file,'/data/input.00/','potinn');
Data.potout = hdf5read(file,'/data/input.00/','potout');
Data.B0 = hdf5read(file,'/data/input.00/','B0');
Data.Rcurv = hdf5read(file,'/data/input.00/','Rcurv');
Data.width = hdf5read(file,'/data/input.00/','width');
Data.n0 = hdf5read(file,'/data/input.00/','n0');
Data.temp = hdf5read(file,'/data/input.00/','temp');
try
Data.it0 = h5readatt(file,'/data/input.00/','it0d');
Data.it1 = h5readatt(file,'/data/input.00/','it2d');
Data.it2 = h5readatt(file,'/data/input.00/','itparts');
catch
Data.it0 = h5readatt(file,'/data/input.00/','it0');
Data.it1 = h5readatt(file,'/data/input.00/','it1');
Data.it1 = h5readatt(file,'/data/input.00/','it2');
end
Data.omepe=sqrt(abs(Data.n0)*Data.qe^2/(Data.me*Data.eps_0));
Data.omece=Data.qe*Data.B0/Data.me;
% Normalizations
Data.tnorm=1/Data.omepe;
Data.rnorm=Data.vlight*Data.tnorm;
Data.bnorm=Data.B0;
Data.enorm=Data.vlight*Data.bnorm;
Data.phinorm=Data.enorm*Data.rnorm;
Data.vnorm=Data.vlight;
% Grid data
Data.rgrid= hdf5read(file, '/data/var1d/rgrid')*Data.rnorm;
Data.zgrid= hdf5read(file, '/data/var1d/zgrid')*Data.rnorm;
Data.dz=(Data.zgrid(end)-Data.zgrid(1))/double(Data.nz);
Br = hdf5read(file,'/data/fields/Br')*Data.bnorm;
Data.Br= reshape(Br,length(Data.zgrid),length(Data.rgrid));
Bz = hdf5read(file,'/data/fields/Bz')*Data.bnorm;
Data.Bz= reshape(Bz,length(Data.zgrid),length(Data.rgrid));
Data.B=sqrt(Data.Bz.^2+Data.Br.^2);
clear Br Bz
Data.femorder = hdf5read(file,'/data/input.00/femorder');
Data.ngauss = hdf5read(file,'/data/input.00/ngauss');
Data.nnr = hdf5read(file,'/data/input.00/femorder');
Data.radii = hdf5read(file,'/data/input.00/radii');
Data.epot = hdf5read(file,'/data/var0d/epot');
Data.ekin = hdf5read(file,'/data/var0d/ekin');
Data.etot = hdf5read(file,'/data/var0d/etot');
Data.eerr = Data.etot-Data.etot(2);
if(strcmp(type,'all') || strcmp(type,'fields'))
pot = hdf5read(file,'/data/fields/pot')*Data.phinorm;
Data.pot= reshape(pot,length(Data.zgrid),length(Data.rgrid),size(pot,2));
clear pot
Er = hdf5read(file,'/data/fields/Er')*Data.enorm;
Data.Er= reshape(Er,length(Data.zgrid),length(Data.rgrid),size(Er,2));
clear Er
Ez = hdf5read(file,'/data/fields/Ez')*Data.enorm;
Data.Ez= reshape(Ez,length(Data.zgrid),length(Data.rgrid),size(Ez,2));
clear Ez
N = hdf5read(file,'/data/fields/partdensity');
Data.N=reshape(N,length(Data.zgrid)-1,length(Data.rgrid)-1,[]);
Data.N=permute(Data.N,[2,1,3]);
Data.N=padarray(Data.N,[1,1],0,'post');
clear N
Data.t2d=Data.dt*(0:size(Data.N,3)-1)*double(Data.it1);
end
if(strcmp(type,'all') || strcmp(type,'parts'))
Data.R = hdf5read(file,'/data/part/R')*Data.rnorm;
Data.Z = hdf5read(file,'/data/part/Z')*Data.rnorm;
Data.UR = hdf5read(file,'/data/part/UR');
Data.UZ = hdf5read(file,'/data/part/UZ');
Data.UTHET= hdf5read(file,'/data/part/UTHET');
Data.GAMM = sqrt(1+Data.UR.^2+Data.UZ.^2+Data.UTHET.^2);
Data.VR = Data.UR./Data.GAMM*Data.vnorm;
clear Data.UR
Data.VZ = Data.UZ./Data.GAMM*Data.vnorm;
clear Data.UZ
Data.VTHET = Data.UTHET./Data.GAMM*Data.vnorm;
clear Data.UTHET
Data.VPERP = sqrt(Data.VR.*Data.VR+Data.VTHET.*Data.VTHET);
Data.tpart=Data.dt*(0:size(Data.R,2)-1)*double(Data.it2);
end
Data.t=Data.dt.*double(0:length(Data.epot)-1);
Data.t0d=Data.dt.*double(0:length(Data.epot)-1);
filedata=dir(file);
Data.timestamp=filedata.date;
disp('M reloaded')
end
c4science · Help