# 2d memecho fit
# 2002 Jan Keith Horne @ St-Andrews
# 2003 Aug KDH @St.A
# 2008 Sep KDH @St.A 

set run = '2d'
	# parameters
set dt = 0.5		# time step (-1 splits median)
set taumin = 0		# min delay
set taumax = 50		# max delay
set cinit = -1		# continuum (-1,0,1 linear, mean, nearest)
set bg = 50		# background (0-1) or (2-99%)

set reload = 'L'	# reload 'L' or not ' '
set chi2 = 1.1		# Chi^2/N
set equal = N		# equalize Chi^2 ? (y/n)
set stifftau = .01	# (increase to stiffen maps vs driver)
set aspect = 1		# (0-inf increase to stiifen vel vs tau)

set its = 100		# iterations per go
set go1 = 'G'
set go2 = 'G'
set go3 = 'G'
set go4 = 'G'
set go5 = 'G'

set noauto = ' '	# autopilot disable 'A' enable ' '
set safe = '0.3'	# safety (shorten steps by factor)

	# echoes and labels
set echo0 = bflx
set lab0 =  'B'
set echo1 = vflx
set lab1 =  'V'

	# lightcurve data
set dir = '.'
set pre = ''
set ext = '.dat'
set dat0 = $dir/$pre$echo0$ext
set dat1 = $dir/$pre$echo1$ext

	# lightcurve fits
set dir = '.'
set pre = '2d_arp151_'
set ext = '.t'
set lcv0 = $dir/$pre$echo0$ext
set lcv1 = $dir/$pre$echo1$ext

	# delay maps
set dir = '.'
set pre = '2d_arp151_'
set ext = '.d'
set map0 = $dir/$pre$echo0$ext
set map1 = $dir/$pre$echo1$ext

	# trailed spectrogram data
set pre = 'arp151'
set trs = ''
set ext = '.trs'
set trsdat = $dir/$pre$trs$ext

	# trailed spectrogram fit
set pre = '2d_arp151'
set ext = '.vt'
set trsfit = $dir/$pre$trs$ext

	# velocity-delay map
set ext = '.vd'
set trsmap = $dir/$pre$trs$ext

	# plots
set plotecho = '1 999'	# which echoes to plot?
set printer = 'k'	# which printer?
set spool = 'n'		# spool plots to printer?
set dir = '.'
set pre = '2d_arp151_'
set plotinit = $dir/$pre'init.ps'
set plotiter = $dir/$pre'iter.ps'
set plotfit = $dir/$pre'fit.ps'

set plotmap = $dir/$pre'map.ps'
set ploteps = $dir/$pre'eps.ps'
set ploteta = $dir/$pre'eta.ps'

alias memecho	$EXE/memecho.exe

memecho << $
$dat0
$lab0
$dat1
$lab1
end
$trsdat
end
$dt
$taumin $taumax
$cinit
$bg

	! load previous fit
$reload
$lcv0
$lcv1
$map1
$trsfit
$trsmap


	! plot initial fit
p
1 3
54
59
64
69
74
85
109
114
119
124
129
140
154
159
164
169
174
200
300
413
418
423
428
433
450
500
550
588
593
598
603
608
685
0 0
p
$plotinit
$printer
p
$spool


	! Chi^2 equalization
c
$equal
	! increase w to stiffen maps
w
$stifftau
	! (0-inf) increase to stiffen vel vs tau
a
$aspect

	! MEM iterations
m


l
$chi2	! chi-2

$noauto
$safe	! autopilot disable
s
$safe	! safety
i
$its	! iterations
$go1
$go2
$go3
$go4
$go5
	! convergence plot
c
p
$plotiter
$printer
l
$spool
	! MEM iterations complete
q

	! plot final fit
p
0 0
p
$plotfit
$printer
p
$spool

	! save results
s
$lcv0
$lcv1
$map1
$trsfit
$trsmap

	! map
i
m
q
q
q
t
Arp 151
q
x
C=$chi2 E=$equal
y
W=$stifftau A=$aspect
q
p
$plotmap
$printer
l
$spool
n
n

	! residuals
i
r
q
q
q
t
Arp 151
q
x
C=$chi2 E=$equal
y
W=$stifftau A=$aspect
q
p
$ploteps
$printer
l
$spool
n
n

	! normalised residuals
i
n
q
q
q
t
Arp 151
q
x
C=$chi2 E=$equal
y
W=$stifftau A=$aspect
q
p
$ploteta
$printer
l
$spool
n
n

	! fini
q
y
$
exit