fft_espic2d.m
6945cf16c0c4geom_test

fft_espic2d.m
View Options

			file='testfinefem2.h5';
			filedata=dir(file);
			timestamp=filedata.date;
			clear M
			if (~exist('M','var') \| ~strcmp(M.file,file) \| timestamp~= M.timestamp)
			M=load_espic2d(file);
			end
			close all hidden;

			figure;
			ax1=gca;
			pcolor(ax1,M.zgrid,M.rgrid,M.N(:,:,end)/double(M.nplasma));
			xlim(ax1,[M.zgrid(1) M.zgrid(end)])
			ylim(ax1,[M.rgrid(1) M.rgrid(end)])
			xlabel(ax1,'Z [m]')
			ylabel(ax1,'R [m]')
			c = colorbar(ax1);
			c.Label.String= '\Phi [V]';
			view(ax1,2)
			%set(ax1,'colorscale','log')

			stepmin=800;

			[x,y]=ginput(1);
			nz=find(x>M.zgrid,1,'last');
			nr=find(y>M.rgrid,1,'last');

			x=(M.zgrid(nz-1)+M.zgrid(nz))/2;
			y=(M.rgrid(nr-1)+M.rgrid(nr))/2;


			dt=M.dt*double(M.it1);

			FS=1/dt;
			L=size(M.N,3)-stepmin;
			f=FS*double(0:L/2)/L/1e9;

			S=M.N(nz,nr,stepmin:end)/double(M.nplasma);
			S=M.pot(nz,nr,stepmin:end)*M.phinorm;

			P2=abs(fft(S(:))/L);
			P1=P2(1:L/2+1);
			P1(2:end-1) = 2*P1(2:end-1);

			tmin=stepmin*dt;
			tmax=dt*size(M.N,3);


			figure
			plot(f(2:end),P1(2:end))
			title(sprintf('Spectrum at z=%3.2e,r=%3.2e, between [%3.2e,%3.2e]s',x,y,tmin,tmax))
			xlabel('f (GHz)')
			ylabel('\|\Phi(f)\|')
			%set(gca,'yscale','log')