densprofiletrex.m
46050b96d3admaster

densprofiletrex.m
View Options

			%% Input for the nonlinear simulation

			% Get the radial density profile from FENNECS
			h=findobj(gcf,'Type','Line');
			rsim=get(h,'XData')/1e3;
			nsim=get(h,'YData');

			% Eliminate data outside the cloud region
			a=1.0e-2;
			b=2.212e-2;
			nsim(rsim<(a+b)/2)=0;
			nsim(rsim>=b)=0;

			% Eliminate lower density values
			rsim=rsim(nsim>max(nsim)/100);
			nsim=nsim(nsim>max(nsim)/100);

			% Interpolation on a new grid
			rm=min(rsim);
			rp=max(rsim);
			nbr=256;
			rinit=linspace(rm,rp,nbr);
			ninit=interp1(rsim,nsim,rinit,'pchip');


			%% Result from the nonlinear simulation
			resfile='/misc/nonlinear/results_trex_shift.h5';

			M=nonlinhdf5(resfile);

			r=M.rgrid;
			dens=-M.rho(:,10,1)'/1.602e-19;

			figure
			plot(rinit,ninit,r,dens)