; b558_analysis.bat
; Computes both the simple Monte Carlo average and the Wiener
; transfer functions.  Then computes a phase 4 predict.  Also
; saves the transfer functions and generates
; some plots.
;
; Author:  Robert Kessel
;          Electro-Optics Technology Section
;          Code 8123
;          Naval Research Laboratory
;          Washington, D.C.  20375-5354
; email:   kessel@ncst.nrl.navy.mil
;
; This code was written for Bill Palya's and Don Walter's
; use at JSU.  However, anyone else who finds it useful
; should feel free to use it, rewrite it, or whatever else,
; save charging any money, goods or services from any other
; person for the code or use thereof.
; 10/23/03

file_path = ':bird_558:'

r_in_file = 'b558_p3_input.txt'
r_in_vari_file = 'b558_p3_input_var.txt'
b_out_file = 'b558_p3_output.txt'
b_out_vari_file = 'b558_p3_output_var.txt'

n_t_funs = 4096
seed = 2866595l

ndata = 256
trial = 2000.

; generate and save Wiener transfer function
monte_wien_t_fun = monte_carlo_wien_t_fun(file_path, r_in_file, r_in_vari_file, $
    b_out_file, b_out_vari_file, n_t_funs, seed)
;
SAVE, monte_wien_t_fun, FILENAME = file_path+'b558_monte_wien_t_fun.dat'
plot_t_fun, ndata, trial, monte_wien_t_fun, $
    'Bird 558, Phase 3 Monte Carlo Wiener tranfer function', file_path, $
    'b558_monte_wien_t_fun.ps'
plot_t_fun_ps, ndata, trial, monte_wien_t_fun, $
    'Bird 558, Phase 3 Monte Carlo Wiener tranfer function', file_path, $
    'b558_monte_wien_t_fun.ps'

; phase 4
r_in_file = 'b558_p4_input.txt'
r_in_vari_file = 'b558_p4_input_var.txt'
b_out_file = 'b558_p4_output.txt'
b_out_vari_file = 'b558_p4_output_var.txt'

; restore the original E42 transfer function for a comparison
e42_conj_t_fun = CONJ(read_256_t_fun('b558_tfun_v1_p3.dat'))
;
predict_from_t_fun, file_path, r_in_file, r_in_vari_file, b_out_file, $
    b_out_vari_file, e42_conj_t_fun, times, b_out, b_out_vari, b_out_pred
b_out_pred_simp = FLOAT(b_out_pred)
;
predict_from_t_fun, file_path, r_in_file, r_in_vari_file, b_out_file, $
    b_out_vari_file, monte_wien_t_fun, times, b_out, b_out_vari, b_out_pred
b_out_pred_wein = FLOAT(b_out_pred)
;
plot_predict_v3, ndata, times, b_out, b_out_vari, b_out_pred_simp, $
    b_out_pred_wein, 558, 4, file_path
plot_predict_v3_ps, ndata, times, b_out, b_out_vari, b_out_pred_simp, $
    b_out_pred_wein, 558, 4, file_path