Functions
This file is part of the Snoopy code Snoopy code is free either version of the	or (at your option) any later version.%Snoopy code is distributed in the hope that it will be useful
	corr (i)
compute transport over time	steps (allow statistical computations for%transport vs Romega) figure(1) subplot(2
compute transport over time plot(t, log10(vel.e)) title('Velocity field	subplot (2, 1, 2)
compute transport over time plot(t, log10(vel.e)) title('Velocity field	plot (t, log10(abs(vel.vxvy))) figure(2) subplot(3
compute transport over time plot(t, log10(vel.e)) title('Velocity field	plot (t, vel.vxmax, t, vel.vxmin)
	title ('Velocity field;extrema vx, vy, vz')
	subplot (3, 1, 2)
	subplot (3, 1, 3)
	figure (3) subplot(2
plot(t, log10(mag.e)) title('Magnetic field	plot (t, log10(abs(mag.bxby))) figure(4) subplot(3
plot(t, log10(mag.e)) title('Magnetic field	plot (t, mag.bxmax, t, mag.bxmin)
	title ('Magnetic field:extrema bx, by, bz')
figure(5) plot(t(mean_start	legend ('Reynolds Stress','Maxwell Stress','alpha coefficient')
	title ('Transport Cumulative average')
	figure (6) subplot(3
semilogy(t, vel.e) title('Energy	semilogy (t, mag.e) subplot(3
semilogy(t, vel.e) title('Energy	semilogy (t, abs(budget.alpha)) figure(7) plot(t
	legend ('alpha coefficient','Maxwell Stress','Reynolds Stress')
	title ('Turbulent Transport') figure(8) plot(t
budget mag	xlabel ('t','fontsize', 16) ylabel('\alpha'
budget mag	set (gca,'fontsize', 16) figure(9) semilogy(t
budget mag vel mag	legend ('Enstrophy','J^2') xlabel('t'
budget mag vel mag mag	title ('Magnetic helicity') xlabel('t'
Variables
This file is part of the Snoopy code Snoopy code is free	software
This file is part of the Snoopy code Snoopy code is free either version of the	License
This file is part of the Snoopy code Snoopy code is free either version of the but WITHOUT ANY	WARRANTY
without even the implied warranty of MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE See the GNU General Public License for more details You should have received a copy of the GNU General Public License along with Snoopy code If	not
without even the implied warranty of MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE See the GNU General Public License for more details You should have received a copy of the GNU General Public License along with Snoopy code If see< http:%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%Example File:how to read the timevar file in Matlab%This example reads all the variable in the timevar file%and display the time history of several quantities...rep='../';mean_start=100;%Where do we start our time-averages?ncorr=20;timevar=(transpose(importdata([rep,'timevar'])));t=timevar(1,:);vel.e=timevar(2,:);mag.e=timevar(3,:);vel.vxmax=timevar(4,:);vel.vxmin=timevar(5,:);vel.vymax=timevar(6,:);vel.vymin=timevar(7,:);vel.vzmax=timevar(8,:);vel.vzmin=timevar(9,:);vel.vxvy=timevar(10,:);mag.bxmax=timevar(11,:);mag.bxmin=timevar(12,:);mag.bymax=timevar(13,:);mag.bymin=timevar(14,:);mag.bzmax=timevar(15,:);mag.bzmin=timevar(16,:);mag.bxby=timevar(17,:);vel.thmax=timevar(18,:);vel.thmin=timevar(19,:);vel.vort=timevar(20,:);mag.curr=timevar(21,:);mag.helicity=timevar(22,:);budget.alpha=vel.vxvy-mag.bxby;clear mean_transport;if length(vel.vxvy) > mean_start	mean_vel_transport = cumsum(vel.vxvy(mean_start:end))./(1:length(vel.vxvy(mean_start:end)))
	mean_mag_transport = cumsum(mag.bxby(mean_start:end))./(1:length(mag.bxby(mean_start:end)))
	mean_tot_transport = cumsum(budget.alpha(mean_start:end))./(1:length(budget.alpha(mean_start:end)))
	e
	alpha_b = -mean(mag.bxby)
	alpha
for	i
end	corr_time
semilogy(t, vel.e) title('Energy budget	t
semilogy(t, vel.e) title('Energy budget mag	bxby
semilogy(t, vel.e) title('Energy budget mag vel	vxvy
budget mag	linewidth
budget mag	fontsize
budget mag vel	vort
budget mag vel mag	curr
budget mag vel mag mag	helicity

Function Documentation

corr ( i )

figure ( 3 )

figure ( 6 )

legend	(	'alpha coefficient'	,
		'Maxwell Stress'	,
		'Reynolds Stress'
	)

budget mag vel mag legend	(	'Enstrophy'	,
		'J^2'
	)

figure (5) plot(t(mean_start legend	(	'Reynolds Stress'	,
		'Maxwell Stress'	,
		'alpha coefficient'
	)

This file is part of the Snoopy code Snoopy code is free either version of the or ( at your option )

plot (t,log10(mag.e)) title('Magnetic field plot	(	t	,
		log10(abs(mag.bxby))
	)

compute transport over time plot (t,log10(vel.e)) title('Velocity field plot	(	t	,
		log10(abs(vel.vxvy))
	)

plot	(	t	,
		mag.	bxmax,
		t	,
		mag.	bxmin
	)

plot	(	t	,
		vel.	vxmax,
		t	,
		vel.	vxmin
	)

semilogy (t,vel.e) title('Energy semilogy	(	t	,
		mag.	e
	)

semilogy (t,vel.e) title('Energy semilogy	(	t	,
		abs(budget.alpha)
	)

end if set	(	gca	,
		'fontsize'	,
		16
	)

compute transport over time steps ( allow statistical computations for%transport vs Romega )

subplot	(	3	,
		1	,
		2
	)

compute transport over time plot (t,log10(vel.e)) title('Velocity field subplot	(	2	,
		1	,
		2
	)

subplot	(	3	,
		1	,
		3
	)

budget mag vel mag mag title ( 'Magnetic helicity' )

title ( 'Transport Cumulative average' )

title	(	'Velocity field;extrema	vx,
		vy	,
		vz'
	)

title ( 'Turbulent Transport' )

title	(	'Magnetic field:extrema	bx,
		by	,
		bz'
	)

budget mag xlabel	(	't'	,
		'fontsize'	,
		16
	)

Variable Documentation

budget alpha

Initial value:

alpha_v-alpha_b
end

% Compute the correlation time of alpha
clear corr

alphanorm=budget.alpha-mean(budget.alpha)

alpha_b = -mean(mag.bxby)

budget mag bxby

end corr_time

Initial value:

sum(corr)*(t(2)-t(1))

plot(corr)

budget mag vel mag curr

e

Initial value:

mean(vel.e(mean_start:end))
    alpha_v=mean_vel_transport(end)
    alpha_b=mean_mag_transport(end)
    alpha=alpha_v-alpha_b
else
    disp(['Caution, we are not removing initial conditions'])

end if fontsize

budget mag vel mag mag helicity

for i

Initial value:

1:ncorr
    j=i-1

This file is part of the Snoopy code Snoopy code is free either version of the License

end if linewidth

mean_mag_transport = cumsum(mag.bxby(mean_start:end))./(1:length(mag.bxby(mean_start:end)))

mean_tot_transport = cumsum(budget.alpha(mean_start:end))./(1:length(budget.alpha(mean_start:end)))

without even the implied warranty of MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE See the GNU General Public License for more details You should have received a copy of the GNU General Public License along with Snoopy code If see<http:%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% Example File: how to read the timevar file in Matlab% This example reads all the variable in the timevar file% and display the time history of several quantities...rep='../';mean_start=100; % Where do we start our time-averages?ncorr=20;timevar=(transpose(importdata([rep,'timevar'])));t=timevar(1,:);vel.e=timevar(2,:);mag.e=timevar(3,:);vel.vxmax=timevar(4,:);vel.vxmin=timevar(5,:);vel.vymax=timevar(6,:);vel.vymin=timevar(7,:);vel.vzmax=timevar(8,:);vel.vzmin=timevar(9,:);vel.vxvy=timevar(10,:);mag.bxmax=timevar(11,:);mag.bxmin=timevar(12,:);mag.bymax=timevar(13,:);mag.bymin=timevar(14,:);mag.bzmax=timevar(15,:);mag.bzmin=timevar(16,:);mag.bxby=timevar(17,:);vel.thmax=timevar(18,:);vel.thmin=timevar(19,:);vel.vort=timevar(20,:);mag.curr=timevar(21,:);mag.helicity=timevar(22,:);budget.alpha=vel.vxvy-mag.bxby;clear mean_transport;if length(vel.vxvy) > mean_start mean_vel_transport = cumsum(vel.vxvy(mean_start:end))./(1:length(vel.vxvy(mean_start:end)))

without even the implied warranty of MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE See the GNU General Public License for more details You should have received a copy of the GNU General Public License along with Snoopy code If not

This file is part of the Snoopy code Snoopy code is free software

end if t

budget mag vel vort

semilogy (t,vel.e) title('Energy budget mag vel vxvy

This file is part of the Snoopy code Snoopy code is free either version of the but WITHOUT ANY WARRANTY

Functions

Variables

Function Documentation

Variable Documentation