OutCurrents.m
View Options

			function [I, pos]=OutCurrents(obj,timestep, subdiv)
			% Computes the Outgoing currens at the simulation axial boundaries at timestep timestep
			% This is simply the surface integral of the axial flux
			if nargin<3
			subdiv=1;
			end
			flux=obj.Axialflux(timestep,[1 obj.nz+1]);
			Iz=squeeze(trapz(obj.rgrid,flux.obj.rgrid)2piobj.qsim/obj.weight);
			Iz(1,:)=-Iz(1,:);
			gamm=obj.Metallicflux(timestep, subdiv);
			Im=zeros(length(gamm.p),length(timestep));
			pos=cell(size(gamm.p));
			for i=1:length(gamm.p)
			p=gamm.p{i};
			pos{i}=p;
			flux=gamm.gamma{i};
			for j=1:length(timestep)
			Im(i,j)=pi/2sum((p(2,1:end-1)+p(2,2:end)).(flux(2:end,j)+flux(1:end-1,j))'...
			.*sqrt((p(1,2:end)-p(1,1:end-1)).^2+(p(2,2:end)-p(2,1:end-1)).^2));
			end
			end
			I=-cat(1,Iz,Im*obj.qsim/obj.weight);
			end

fennecs/matlab/@fennecshdf5/OutCurrents.m81171805e064breakdown