#!/bin/bash
#!/bin/bash

#	
#	Purpose: 
#	Calls the following shell scripts:
#		iter_c2r.com
#		make_greenfct_segs.com
#		prep_load_grd.com
#	Calls the following programs:
#		flex2dc (F77)
#		grd_interp (C)
#		FITPROG (F90)

if (( "$#" != 1 )); then
  echo "Usage: run_trench_region.com # of iterations "
  exit	
fi


FITPROG=fitflexl1
SEGLIST=cchile_100k_list.txt
LLCLON=263.0
LLCLAT=-40.0039553203
RIGINIT=rigid0_xy.grd
BMOSCALE=1.0E14
VLDSCALE=1.0E12
TOPODAT=topo_cmfd.dat
GRAVDAT=grav_mfd.dat
TOPOWGT=1.0
GRAVWGT=1.0

COUNTER=0
let NITER=$1
echo $NITER

echo "Solver: $FITPROG"


while [ $COUNTER -lt $NITER ]; do
	let COUNTER=COUNTER+1	 
	echo "Parameter Estimation Iteration # $COUNTER"
	echo "Now Generating Green's functions"
	make_greenfct_segs.com $SEGLIST $LLCLON $LLCLAT $RIGINIT $COUNTER $TOPODAT $GRAVDAT $BMOSCALE $VLDSCALE

done

cp w_inv_xy_prev.grd w_inv${6}_xy.grd
cp curv_inv_xy_prev.grd curv_inv${6}_xy.grd
cp grav_inv_xy_prev.grd grav_inv${6}_xy.grd
cp rigid_inv_prev.grd rigid_inv${6}_xy.grd

echo "Generating Green's functions"

make_greenfct_segs.com $SEGLIST $LLCLON $LLCLAT $RIGINIT 0 $TOPODAT $GRAVDAT $BMOSCALE $VLDSCALE
 
ls gg_seg*_iter1.dat > grav1.fnct
ls gb_seg*_iter1.dat > bath1.fnct

echo "Estimating values for bending moment and vertical load"

$FITPROG topo_cmfd.dat bath1.fnct 1.0 grav_mfd.dat grav1.fnct 1.0

prep_load_grd.com 1 1.0 lq.dat $SEGLIST

iter_c2r.com 3 fz_inv1_xy.grd $RIGINIT rigid.grd_lookup age_xy.grd 1

echo "Generating Green's functions"

make_greenfct_segs.com $SEGLIST $LLCLON $LLCLAT rigid_inv1_xy.grd 1 $TOPODAT $GRAVDAT $BMOSCALE $VLDSCALE

ls gg_seg*_iter2.dat > grav2.fnct
ls gb_seg*_iter2.dat > bath2.fnct

echo "Estimating values for bending moment and vertical load"

$FITPROG topo_cmfd.dat bath2.fnct 1.0 grav_mfd.dat grav2.fnct 1.0

prep_load_grd.com 2 1.0 lq.dat $SEGLIST

iter_c2r.com 3 fz_inv2_xy.grd rigid_inv1_xy.grd rigid.grd_lookup age_xy.grd 2

echo "Generating Green's functions"

make_greenfct_segs.com $SEGLIST $LLCLON $LLCLAT rigid_inv2_xy.grd 2 $TOPODAT $GRAVDAT $BMOSCALE $VLDSCALE

ls gg_seg*_iter3.dat > grav3.fnct
ls gb_seg*_iter3.dat > bath3.fnct

echo "Estimating values for bending moment and vertical load"

$FITPROG topo_cmfd.dat bath3.fnct 1.0 grav_mfd.dat grav3.fnct 1.0

prep_load_grd.com 3 1.0 lq.dat $SEGLIST

iter_c2r.com 3 fz_inv2_xy.grd rigid_inv1_xy.grd rigid.grd_lookup age_xy.grd 3

make_greenfct_segs.com $SEGLIST $LLCLON $LLCLAT rigid_inv3_xy.grd 3 $TOPODAT $GRAVDAT $BMOSCALE $VLDSCALE

ls gg_seg*_iter4.dat > grav4.fnct
ls gb_seg*_iter4.dat > bath4.fnct

$FITPROG topo_cmfd.dat bath4.fnct 1.0 grav_mfd.dat grav4.fnct 1.0

prep_load_grd.com 4 1.0 lq.dat $SEGLIST

flex2dc fz_inv4_xy.grd rigid_inv3_xy.grd 1 w_inv4_xy.grd curv_inv4_xy.grd grav_inv4_xy.grd

grd_interp rigid.grd_lookup age_xy.grd curv_inv4_xy.grd rigid_inv4_xy.grd
grd_interp moment.grd_lookup age_xy.grd curv_inv4_xy.grd moment_inv4_xy.grd
grd_interp yld.grd_lookup age_xy.grd curv_inv4_xy.grd yld_inv4_xy.grd