# bscaling Ts Tp Tis Tip P
bscaling
	bbox 
	bplot 10 10 0 1 128
	bplot 10 100 0 1 128
	bplot 10 1000 0 1 128
	bplot 10 10000 0 1 128
	bplot 10 100000 0 1 128
	relocate 0 0
	ltype 2
	draw 128 128
	ltype 0
bplot 5
	define Ts $1
	define Tp $2
	define Tis $3
	define Tip $4
	define Pmax $5
        define Pmp1 ($Pmax + 1)
	set dimen(R) = $Pmp1
	set dimen(P) = $Pmp1
	do i=1,$Pmax {
		set P[($i)] = ($i)
		set R[($i)] = ($Ts + $Tp)/($Ts + P[($i)]*$Tis + ($Tp/P[($i)]) + $Tip)
		# set R[($i)] = ($Ts + $Tp)/($Ts + P[($i)]*P[($i)]*$Tis + ($Tp/P[($i)]) + $Tip)
	}
	ltype 0
	connect P R
bbox 
	erase
	location 5000 30000 6000 28000
	define TeX_strings 1
	expand 1.0
	limits 0 128 0 128
	box
	expand 1.5
	angle 0
	ylabel S(N) (speedup)
	xlabel N (processors)
	expand 1.0
	ltype 0
ALogplot4 
	ALogbox 0 11 -15 2.5
	data Brown.Ciftan
	identification L = 4 A_x(t) fits for x = R,(\Delta,||),(\Delta,\perp).
	ANdata A.4.1.4427.average.pub 1 4		# A_R fit
	ALogpoints 9 3 0
	ARLogpoints 9 0 1
	relocate 2 1.5
	label L = 4
	relocate 3 1.5
	label \point 9 0 1/\tau_R = $b  (\chi^2 = $chisq)
	ANdata A.4.1.4427.average.pub 2 4		# A_perp fit
	AperpLogpoints 3 0 3
	relocate 3 0.0
	label \point 3 0 1/\tau_{\Delta,\perp} = $b  (\chi^2 = $chisq)
	ANdata A.4.1.4427.average.pub 2 4		# A_par fit
	AparLogpoints 4 0 2
	relocate 3 0.75
	label \point 4 0 1/\tau_{\Delta,||} = $b  (\chi^2 = $chisq)
ALogplot8
	ALogbox 0 21 -10 5
	# ALogbox 0 21 -2.5 -1.5
	data Brown.Ciftan
	expand 1.1
	identification L = 8 A_x(t) fits for x = R,(\Delta,||),(\Delta,\perp).
	relocate 2 4
	label L = 8
	ANdata A.8.1.4427.average.pub 1 4		# A_R fit
	ALogpoints 8 3 0
	ARLogpoints 8 0 1
	relocate 2.0 3.
	expand 0.99
	label \point 8 0 1/\tau_R = $b  (\chi^2 = $chisq)
	expand 1.1
	ANdata A.8.1.4427.average.pub 10 15		# A_perp fit
	AperpLogpoints 3 0 3
	AperpLogline 0
	relocate 2.0 2.0
	expand 0.99
	label \point 3 0 1/\tau_{\Delta,\perp} = $b  (\chi^2 = $chisq)
	expand 1.1
	ANdata A.8.1.4427.average.pub 10 15		# A_par fit
	AparLogpoints 4 0 2
	AparLogline 0
	relocate 2.0 1.0
	expand 0.99
	label \point 4 0 1/\tau_{\Delta,||} = $b  (\chi^2 = $chisq)
	expand 1.1
	expand 1.1
ALogplot16
	ALogbox 0 80 -15 2.5
	# ALogbox 0 50 -3 -2.5
	data Brown.Ciftan
	expand 1.1
	identification L = 16 A_x(t) fits for x = R,(\Delta,||),(\Delta,\perp).
	relocate 10 3
	label L = 16
	ANdata A.16.1.4427.average.pub 1 4		# A_R fit
	ALogpoints 8 3 0
	ARLogpoints 8 0 1
	relocate 5.0 1.5
	expand 0.99
	label \point 8 0 1/\tau_R = $b  (\chi^2 = $chisq)
	expand 1.1
	expand 1.1
	ANdata A.16.1.4427.average.pub 30 60	# A_perp fit
	AperpLogpoints 3 0 3
	AperpLogline 0
	relocate 5.0 0.75
	expand 0.99
	label \point 3 0 1/\tau_{\Delta,\perp} = $b  (\chi^2 = $chisq)
	expand 1.1
	expand 1.1
	ANdata A.16.1.4427.average.pub 30 60	# A_par fit
	AparLogpoints 4 0 2
	AparLogline 0
	relocate 5.0 0.0
	expand 0.99
	label \point 4 0 1/\tau_{\Delta,||} = $b  (\chi^2 = $chisq)
	expand 1.1
	expand 1.1
ALogplot24
	ALogbox 0 105 -10 2.5
	# ALogbox 0 50 -3 -2.5
	data Brown.Ciftan
	expand 1.1
	identification L = 24 A_x(t) fits for x = R,(\Delta,||),(\Delta,\perp).
	relocate 10 1.5
	label L = 24
	ANdata A2.24.1.4427.average 2 10		# A_R fit
	ALogpoints 8 3 0
	ARLogpoints 8 0 1
	relocate 10 1.0
	label \point 8 0 1/\tau_R = $b  (\chi^2 = $chisq)
	ANdata A2.24.1.4427.average 60 100	# A_perp fit
	AperpLogpoints 3 0 3
	AperpLogline 0
	relocate 10 0.5
	label \point 3 0 1/\tau_{\Delta,\perp} = $b  (\chi^2 = $chisq)
	ANdata A2.24.1.4427.average 80 100	# A_par fit
	AparLogpoints 4 0 2
	AparLogline 0
	relocate 10 0.0
	label \point 4 0 1/\tau_{\Delta,||} = $b  (\chi^2 = $chisq)
	expand 1.1
ALogplot32
	ALogbox 0 160 -10 2.5
	data Brown.Ciftan
	identification L = 32 A_x(t) fits for x = R,(\Delta,||),(\Delta,\perp).
	ANdata A.32.1.4427.average.pub 1 4 		# A_R fit
	ALogpoints 9 3 0
	ARLogpoints 9 0 1
	relocate 10 1.5
	label L = 32
	relocate 20 1.5
	label \point 9 0 1/\tau_R = $b  (\chi^2 = $chisq)
	ANdata A.32.1.4427.average.pub 90 130		# A_perp fit
	AperpLogpoints 3 0 3
	relocate 20 0.75
	label \point 3 0 1/\tau_{\Delta,\perp} = $b  (\chi^2 = $chisq)
	ANdata A.32.1.4427.average.pub 90 130		# A_par fit
	AparLogpoints 4 0 2
	relocate 20 0.0
	label \point 4 0 1/\tau_{\Delta,||} = $b  (\chi^2 = $chisq)
	relocate 38 -2.0
	expand 1.1
	label Total Autocorrelation
	draw_arrow 38 -2.0 30 -2.75
	expand 1.01
APerpLogplot
	APerpLogbox 0 160 -10 2.5
	identification A_\perp(L,t) on semilog scale
	ANdata A.4.1.4427.a13559 2 6
	AperpLogpoints 3 0 0
	relocate 50 1.5
	label -1/\tau_{\Delta,\perp}(4) = $b, \chi^2 = $chisq
	ANdata A.8.1.4427.average 15 20 # 08297 15 25
	AperpLogpoints 4 0 0
	relocate 50 0.75
	label -1/\tau(8) = $b, \chi^2 = $chisq
	ANdata A.16.1.4427.average 30 60
	AperpLogpoints 5 0 0
	relocate 50 0.00
	label -1/\tau(16) = $b, \chi^2 = $chisq
	ANdata A.32.1.4427.average 85 150
	AperpLogpoints 6 0 0
	relocate 50 -0.75
	label -1/\tau(32) = $b, \chi^2 = $chisq
AParLogplot
	AParLogbox 0 160 -13 2.5
	identification A_{||}(L,t) on semilog scale
	ANdata A.4.1.4427.a13559 2 4
	AparLogpoints 3 0 0
	relocate 50 1.5
	label -1/\tau_{\Delta,||}(4) = $b, \chi^2 = $chisq
	ANdata A.8.1.4427.average 8 18 # 08297 5 17
	AparLogpoints 4 0 0
	relocate 50 0.75
	label -1/\tau_{||}(8) = $b, \chi^2 = $chisq
	ANdata A.16.1.4427.average 30 60
	AparLogpoints 5 0 0
	relocate 50 0.0
	label -1/\tau_{\Delta,||}(16) = $b, \chi^2 = $chisq
	ANdata A.32.1.4427.average 50 150
	AparLogpoints 6 0 0
	relocate 50 -0.75
	label -1/\tau_{\Delta,||}(32) = $b, \chi^2 = $chisq
ARLogplot
	ARLogbox 0 21 -5 1.0
	identification A_R(L,t) on semilog scale
	ANdata A2.4.1.4427.average 2 10
	ARLogpoints 6 0 0
	relocate 3 0.5
	label L = 4, b_R = $b, \chi^2 = $chisq
	ANdata A2.8.1.4427.average 2 10
	ARLogpoints 6 0 0
	relocate 3 0.0
	label L = 8, b_R = $b, \chi^2 = $chisq
	ANdata A2.16.1.4427.average 2 10
	ARLogpoints 6 0 0
	relocate 3 -0.5
	label L = 16, b_R = $b, \chi^2 = $chisq
	ANdata A.32.1.4427.average 2 10
	ARLogpoints 6 0 0
	relocate 3 -1.0
	label L = 32, b_R = $b, \chi^2 = $chisq
NALogpicture
	ALogbox 0 11 -10 0
	# ANdata A.4.1.4427.a11706
	ANdata A.8.1.4427.a12484  #11986
	ALogpoints 5 0 0
	ARLogpoints 5 0 1
	AparLogpoints 5 1 2
	AperpLogpoints 5 2 3
	AxLogpoints 5 3 4
	ANdata A.8.1.4427.a11986
	ALogpoints 4 0 0
	ARLogpoints 4 0 1
	AparLogpoints 4 1 2
	AperpLogpoints 4 2 3
	AxLogpoints 4 3 4
#	ANdata A.4.1.4427.a11837
	AparLogpoints 4 1 2
#	ALogpoints 4 0 0
	ARLogpoints 4 0 1
	AperpLogpoints 4 2 3
	AxLogpoints 4 3 4
#
	# ASumLogpoints 3 1 0
	# AMRLogpoints 3 3 4
ARLogpicture
	ALogbox 0 250 -5 2
	ARdata A.4.1.4427.a07199
	ARLogpoints 3 1 2
	relocate 50 1.5
	label L = 4, b = $b, \chi^2 \ $chisq
	ARdata A.8.1.4427.a25595
	ARLogpoints 4 1 2
	relocate 50 1
	label L = 8, b = $b, \chi^2 \ $chisq
	ARdata A.16.1.4427.a07780
	ARLogpoints 5 1 2
	relocate 50 0.5
	label L = 16, b = $b, \chi^2 \ $chisq
	ARdata A.32.1.4427.a10215  # 08662
	ARLogpoints 6 1 2
	relocate 50 0
	label L = 32, b = $b, \chi^2 \ $chisq
ARLogtest
	ALogbox 0 20 -2 2
	ARdata A.4.1.4427.a07199.save
	ARLogpoints 3 1 2
#	relocate 5 1.5
	label L = 4, b = $b, \chi^2 = $chisq
	ARdata A.4.1.4427.a00995
	ARLogpoints 3 1 2
	relocate 5 1.25
	label L = 4, b = $b, \chi^2 = $chisq
	ARdata A.4.1.4427.a01034
	ARLogpoints 3 0 2
	relocate 5 1.0
	label L = 4, b = $b, \chi^2 = $chisq
	ARdata A.4.1.4427.a00865
	ARLogpoints 3 0 2
	relocate 5 0.75
	label L = 4, b = $b, \chi^2 = $chisq
	ARdata A.4.1.4427.a00902
	ARLogpoints 3 2 2
	relocate 5 0.5
	label L = 4, b = $b, \chi^2 = $chisq
	ARdata A.4.1.4427.a00925
	ARLogpoints 3 2 2
	relocate 5 0.0
	label L = 4, b = $b, \chi^2 = $chisq
	ARdata A.4.1.4427.a00933
	ARLogpoints 3 2 2
	relocate 5 -0.5
	label L = 4, b = $b, \chi^2 = $chisq
	ARdata A.4.1.4427.a00940
	ARLogpoints 3 2 2
	relocate 5 -1.0
	label L = 4, b = $b, \chi^2 = $chisq
	ARdata A.4.1.4427.a00963
	ARLogpoints 3 2 2
	relocate 5 -1.5
	label L = 4, b = $b, \chi^2 = $chisq
Adata 3
	data $1
	lines 0 0
	read {t 3 A 4 Aerr 5}
	define nskip $2
	define nmax $3
	define ni ($nmax - $nskip - 1)
	set dimen(tF) = $($nmax - $nskip)
	set dimen(AF) = $($nmax - $nskip)
	set dimen(AFerr) = $($nmax - $nskip)
	do i=0,$ni {
		set tF[($i)] = t[($i+$nskip)]
		set AF[($i)] = A[($i+$nskip)]
		set AFerr[($i)] = Aerr[($i+$nskip)]
	}
ARdata 1
	data $1
	lines 0 0
	read {t 3 A 4 Aerr 5 AR 6 ARerr 7}
	set RA=A+AR
	set RAerr=Aerr+ARerr
	define nskip $2
	define nmax $3
	define ni ($nmax - $nskip - 1)
	set dimen(tF) = $($nmax - $nskip)
	set dimen(AF) = $($nmax - $nskip)
	set dimen(AFerr) = $($nmax - $nskip)
	do i=0,$ni {
		set tF[($i)] = t[($i+$nskip)]
		set AF[($i)] = A[($i+$nskip)]
		set AFerr[($i)] = Aerr[($i+$nskip)]
	}
	define nskip 25
	define nmax 50
	define ni ($nmax - $nskip - 1)
	set dimen(tRF) = $($nmax - $nskip)
	set dimen(ARF) = $($nmax - $nskip)
	set dimen(ARFerr) = $($nmax - $nskip)
	do i=0,$ni {
		set tRF[($i)] = t[($i+$nskip)]
		set ARF[($i)] = AR[($i+$nskip)]
		set ARFerr[($i)] = ARerr[($i+$nskip)]
	}
Apicture
	Abox 0 200
	Adata A.4.1.4427.a23076
	Apoints 3 0
	Aline 0
	Adata A.8.1.4427.a06207  # 17489
	Apoints 4 0
	Aline 0
	Adata A.16.1.4427.a06582 # 19580
	Apoints 5 0
	Aline 2
	Adata A.32.1.4427.a26923
	Apoints 6 0
	Aline 0
ANdata 3
	data $1
	lines 0 0
	read {t 1 A 2 Aerr 3 Apar 4 Aparerr 5 Aperp 6 Aperperr 7 Ax 8 Axerr 9 AR 10 ARerr 11}
	# Renormalize Ax, as it is twice the published value
	set Ax2 = 2*Ax
	set Ax2err = sqrt(2)*Ax
	# Subtract the asymptotic behavior of Apar, which turns out
	# to be AR!
	set Apar = Apar - AR
	set Aparerr = Aparerr + ARerr
	define nskip $2
	define nmax ($3+1)
	define ni ($nmax - $nskip - 1)
	set dimen(tF) = $($nmax - $nskip)
	set dimen(ARF) = $($nmax - $nskip)
	set dimen(ARFerr) = $($nmax - $nskip)
	set dimen(AperpF) = $($nmax - $nskip)
	set dimen(AperpFerr) = $($nmax - $nskip)
	set dimen(AparF) = $($nmax - $nskip)
	set dimen(AparFerr) = $($nmax - $nskip)
	do i=0,$ni {
		set tF[($i)] = t[($i+$nskip)]
		set ARF[($i)] = AR[($i+$nskip)]
		set ARFerr[($i)] = ARerr[($i+$nskip)]
		set AperpF[($i)] = Aperp[($i+$nskip)]
		set AperpFerr[($i)] = Aperperr[($i+$nskip)]
		set AparF[($i)] = Apar[($i+$nskip)]
		set AparFerr[($i)] = Aparerr[($i+$nskip)]
	}
Abox 2
	erase
	location 5000 30000 6000 28000
	define TeX_strings 1
	expand 1.0
	limits $1 $2  -0.1 1.1
	box
	relocate $1 0.0
	draw $2 0.0
	expand 1.5
	angle 0
	ylabel <S_i(0)S_i(t)>
	xlabel t
	expand 1.0
	ltype 0
Apoints 2
	ptype $1 $2
	points t A
	error_y t A Aerr
Aline 1
	ltype $1
	connect t A
	ltype 0
ARpoints 2
	ptype $1 $2
	points t AR
	error_y t AR ARerr
ARline 1
	ltype $1
	connect t AR
	ltype 0
ALogpoints 3
	ptype $1 $2
	set LA = ln(abs(A))
	set LAerr = ln(abs(A + Aerr))
	set LAerr = LAerr - LA
	points t LA
	error_y t LA LAerr
	ltype $3
	connect t LA
	ltype 0
ANLogpoints 3
	ptype $1 $2
	set LA = ln(abs(A))
	set LAerr = ln(abs(A + Aerr))
	set LAerr = LAerr - LA
	points t LA
	error_y t LA LAerr
	set tfit={0 1e4}
	mylsq t LA LAerr tfit LAfit
	ltype $3
	connect tfit LAfit
	ltype 0
ASumLogpoints 3
	ptype $1 $2
	set LASum = ln(abs(Apar+Aperp+2*Ax))
	points t LASum
	error_y t LA LAerr
	ltype $3
	connect t LASum
	ltype 0
AMRLogpoints 3
	ptype $1 $2
	set LAMR = ln(abs(A - AR))
	set LAMRerr = ln(abs((A-AR) + (Aerr + ARerr)))
	set LAMRerr = LAMRerr - LAMR
	points t LAMR
	error_y t LAMR LAMRerr
	ltype $3
	connect t LAMR
	ltype 0
AparLogpoints 3
	ptype $1 $2
	set LApar = ln(abs(Apar))
	points t LApar
	lnerr t LApar Aparerr
	set LAparF = ln(abs(AparF))
	set LAparFerr = ln(abs(AparF) - abs(AparFerr))
	set LAparFerr = LAparF - LAparFerr
	set tfit={0 1e4}
	mylsq tF LAparF LAparFerr tfit LAparfit
	ltype $3
	connect tfit LAparfit
	ltype 0
AparLogline 1
	ltype $1
	connect t LApar
	ltype 0
AperpLogpoints 3
	ptype $1 $2
	set LAperp = ln(abs(Aperp))
	points t LAperp
	lnerr t LAperp Aperperr
	set LAperpF = ln(abs(AperpF))
	set LAperpFerr = ln(abs(AperpF) - abs(AperpFerr))
	set LAperpFerr = LAperpF - LAperpFerr
	set tfit={0 1e4}
	mylsq tF LAperpF LAperpFerr tfit LAperpfit
	ltype $3
	connect tfit LAperpfit
	ltype 0
AperpLogline 1
	ltype $1
	connect t LAperp
	ltype 0
AxLogpoints 3
	ptype $1 $2
	set LAx2 = ln(abs(Ax2))
	set LAx2err = ln(abs(Ax2 + Ax2err))
	set LAx2err = LAx2err - LAx2
	points t LAx2
	error_y t LAx2 LAx2err
	ltype $3
	connect t LAx2
	ltype 0
ALogline 1
	ltype $1
	connect t (ln(abs(A)))
	ltype 0
ALogLogpoints 2
	ptype $1 $2
	set Lt = ln(t)
	set LA = ln(abs(A))
	set LAerr = ln(abs(A + Aerr))
	set LAerr = LAerr - LA
	points Lt LA
	# error_y t LA LAerr
ALogLogline 1
	ltype $1
	connect (ln(t)) (ln(abs(A)))
	ltype 0
ARLogpoints 3
	ptype $1 $2
	set LAR = ln(abs(AR))
	points t LAR
	lnerr t LAR ARerr
	set LARF = ln(abs(ARF))
	set LARFerr = ln(abs(ARF) - abs(ARFerr))
	set LARFerr = LARF - LARFerr
	set tfit={0 1e4}
	mylsq tF LARF LARFerr tfit LARfit
	ltype $3
	connect tfit LARfit
	ltype 0
ANRLogpoints 3
	ptype $1 $2
	set LAR = ln(abs(AR))
	set LARerr = ln(abs(AR + ARerr))
	set LARerr = LARerr - LAR
	set tfit={0 1e4}
	mylsq t LAR LARerr tfit LARfit
	points t LAR
	# points t (ln(abs(AR)))
	ltype $3
	connect tfit LARfit
	ltype 0
ARLogline 1
	ltype $1
	connect t (ln(abs(AR)))
	ltype 0
RALogpoints 2
	ptype $1 $2
	set LRA = ln(abs(RA))
	set LRAerr = ln(abs(RA + RAerr))
	set LRAerr = LRAerr - LRA
	points t LRA
	# error_y t LRA LRAerr
RALogline 1
	ltype $1
	connect t (ln(abs(RA)))
	ltype 0
# Afit (slow amplitude) kappa_1 kappa_2 l
Afit 4
	set dimen(tfit)=501
	set dimen(Afit)=501
	do i=0,500 {
		set tfit[($i)] = $i/10.0
		set Afit[($i)] = (1 - $1)*exp(-$2*tfit[($i)]) + $1*exp(-$3*tfit[($i)])
	}
	ltype $4
	connect tfit Afit
	ltype 0
# Alogplotb
Alogplotb
	Alogbox 0 5 -10 0
	Adata A1.dyn
	# Alogpoint 3 0
	Alogline 2
	Adata A2.dyn
	# Alogpoint 4 0
	Alogline 0
	expand 1.0
ALogbox 4
	erase
	location 5000 30000 6000 28000
	define TeX_strings 1
	expand 1.0
	limits $1 $2 $3 $4
	box
	expand 1.5
	angle 0
	ylabel ln(A_x(t))
	xlabel t
	expand 1.0
	ltype 0
ARLogbox 4
	erase
	location 5000 30000 6000 28000
	define TeX_strings 1
	expand 1.0
	limits $1 $2 $3 $4
	box
	expand 1.5
	angle 0
	ylabel ln<\over"1""N" (\sum_i {\bf S}_i(0))\cdot \over"1""N" (\sum_j {\bf S}_j(t))>
	xlabel t
	expand 1.0
	ltype 0
APerpLogbox 4
	erase
	location 5000 30000 6000 28000
	define TeX_strings 1
	expand 1.0
	limits $1 $2 $3 $4
	box
	expand 1.5
	angle 0
	ylabel ln<\over"1""N"(\sum_i {\bf S}_{\perp,i}(0)\cdot{\bf S}_{\perp,i}(t))>
	xlabel t
	expand 1.0
	ltype 0
AParLogbox 4
	erase
	location 5000 30000 6000 28000
	define TeX_strings 1
	expand 1.0
	limits $1 $2 $3 $4
	box
	expand 1.5
	angle 0
	expand 1.5
	ylabel ln<\over"1""N" (\sum_i {\bf S}_{||,i}(0)\cdot{\bf S}_{||,i}(t)) - \over"1""N"(\sum_i {\bf S}_i(0))\cdot\over"1""N"(\sum_j {\bf S}_j(t))>
	expand 1.5
	xlabel t
	expand 1.01
	ltype 0
ALogLogbox 4
	erase
	location 5000 30000 6000 28000

	define TeX_strings 1
	expand 1.0
	limits $1 $2 $3 $4
	box
	relocate $1 0.0
	draw $2 0.0
	expand 1.5
	angle 0
	ylabel ln(<S_i(0)S_i(t)>)
	xlabel ln(t)
	expand 1.0
	ltype 0
ALLTauPerpbox 4
	erase
	location 5000 30000 6000 28000
	define TeX_strings 1
	expand 1.01
	limits $1 $2 $3 $4
	box
	expand 1.5
	angle 0
	ylabel log(\tau_\perp)
	xlabel log(L)
	expand 1.01
	ltype 0
ALLTauParbox 4
	erase
	location 5000 30000 6000 28000
	define TeX_strings 1
	expand 1.01
	limits $1 $2 $3 $4
	box
	expand 1.5
	angle 0
	ylabel log(\tau_{||})
	xlabel log(L)
	expand 1.01
	ltype 0
ARLLTaubox 4
	erase
	location 5000 30000 6000 28000
	define TeX_strings 1
	expand 1.01
	limits $1 $2 $3 $4
	box
	expand 1.5
	angle 0
	ylabel log(\tau_R(L))
	xlabel log(L)
	expand 1.01
	ltype 0
ALLplot
	ALLTaubox 0 2 -0.5 2.5
	data Brown.Ciftan
	identification Preliminary Autocorrelation study.
	data zdata
	read {L 1 Tinv 2 TiErr 3}
	ptype 3 0
	set T=(1/Tinv)
	set LL=lg(L)
	set LT=lg(T)
	set LTE=lg(1/(Tinv-TiErr))
	set LTE=LTE - LT
	set LLF = {0.5 1.6}
	mylsq LL LT LTE LLF LTF
	relocate 0.2 2.1
	expand 1.01
	label  T = 1.4427: \chi^2 = $chisq Z = $b
	expand 1.5
	ptype 3 0
	points LL LT
	error_y LL LT LTE
	ltype 2
	connect LLF LTF
	ltype 0
APerpLLplot 1
	ALLTauPerpbox 0 2 -0.5 2.5
	data Brown.Ciftan
	identification Z_\perp for \tau_\perp(L) obtained from A_\perp(L,t)
	data zdata.$1.perp
	read {L 1 Tinv 2 TiErr 3}
	ptype 3 0
	set T=(1/Tinv)
	set LL=lg(L)
	set LT=lg(T)
	set LTE=lg(1/(Tinv-TiErr))
	set LTE=LTE - LT
	set LLF = {0.5 1.6}
	mylsq LL LT LTE LLF LTF
	relocate 0.2 2.1
	expand 1.01
	label  T = 1.4427: \chi^2 = $chisq Z_\perp = $b
	expand 1.5
	ptype 3 0
	points LL LT
	error_y LL LT LTE
	ltype 2
	connect LLF LTF
	ltype 0
AParLLplot 1
	ALLTauParbox 0 2 -0.5 2.5
	data Brown.Ciftan
	identification Z_{||} for \tau_{||}(L) obtained from \bar{A}_{||}(L,t)
	data zdata.$1.par
	read {L 1 Tinv 2 TiErr 3}
	ptype 3 0
	set T=(1/Tinv)
	set LL=lg(L)
	set LT=lg(T)
	set LTE=lg(1/(Tinv-TiErr))
	set LTE=LTE - LT
	set LLF = {0.5 1.6}
	mylsq LL LT LTE LLF LTF
	relocate 0.2 2.1
	expand 1.01
	label  T = 1.4427: \chi^2 = $chisq Z_{||} = $b
	expand 1.5
	ptype 3 0
	points LL LT
	error_y LL LT LTE
	ltype 2
	connect LLF LTF
	ltype 0
ARLLplot 1
	ARLLTaubox 0 2 0.0 3.5
	data Brown.Ciftan
	identification Z_R for \tau_R(L) obtained from A_R(L,t)
	data zdata.$1.rot
	read {L 1 Tinv 2 TiErr 3}
	ptype 3 0
	set T=(1/Tinv)
	set LL=lg(L)
	set LT=lg(T)
	set LTE=lg(1/(Tinv-TiErr))
	set LTE=LTE - LT
	set LLF = {0.5 1.6}
	mylsq LL LT LTE LLF LTF
	relocate 0.2 3.1
	expand 1.01
	label  T = 1.4427: \chi^2 = $chisq Z_R = $b
	expand 1.5
	ptype 3 0
	points LL LT
	error_y LL LT LTE
	ltype 2
	connect LLF LTF
	ltype 0