int iter_num = {argv {4}}

// initialize file bat.out
if (iter_num == 1) // start of batch run

  echo % recall_exc_wt_mult EC2pyr_maxweight expand_recall_time bg1_perc bg2_perc bg3_perc bg4_perc \
	mean_perc_bg SD_perc_bg mean_diff_perc_bg SD_diff_perc_bg recall_accuracy > data/bat.out
end

// Load Hippo simulation
bat_Hippo // only constants.g loaded in here

// reset some global constants defined in constants.g using command line arguments passed from rbat
float	recall_exc_wt_mult = {argv {1}} * 1
float	EC2pyr_maxweight = {argv {2}} * 1e-1
float	expand_recall_time = {argv {3}} * 1e-3

// load rest of simulation files (do not reload constants.g)
bat_hippo_sim

// run sequence storage/recall function
lrs

// save results of run in bat.out
floatformat %.4f

/*
if (expand_recall_time == 1) // new value of parameter
 echo % >> data/bat.out // newline
 echo "% ******************************************************" >> data/bat.out
 echo % >> data/bat.out // newline
end
*/


echo {recall_exc_wt_mult} {EC2pyr_maxweight} " " -n >> data/bat.out
floatformat %.4f
echo {expand_recall_time} " " -n >> data/bat.out
floatformat %.4f
echo {getfield /bg_perc table->table[0]} {getfield /bg_perc table->table[1]} \
  {getfield /bg_perc table->table[2]} {getfield /bg_perc table->table[3]} \
  {mean_perc_bg} {SD_perc_bg} {mean_diff_perc_bg} {SD_diff_perc_bg} {RCL_ACC} >> data/bat.out


// terminate genesis 
quit