#!/bin/bash

N=10
beta=40
nom=200

saverage.pl ../sig.inp.1[789].1 > Sig.average

grep -v '#' Sig.average > Sig.work

for x in `seq 1 $N`; do strans Sig.work '#1' "#$((2*x))" "#$((2*x+1))" "#$((2*(x+$N)))" "#$((2*(x+$N)+1))" > Sig.$x; done

for x in `seq 1 $N`; do mkdir n.$x; done

for x in `seq 1 $N`; do cp Sig.$x pmesh.py ancont.py chi2f.so n.$x; done

for x in `seq 1 $N`; do cd n.$x; ./ancont.py -sig Sig.$x -nom $nom -beta $beta -wexp 1.15 -Ng 60 -FermiLiquid True -L0 20 -lcut 0.0 -b 0.9 -alpha4 1. -p0 0.1 -vunity 0.01 ; cd ../; done


strans n.1/Sig.out n.2/Sig.out n.3/Sig.out n.4/Sig.out n.5/Sig.out n.6/Sig.out n.7/Sig.out n.8/Sig.out n.9/Sig.out n.10/Sig.out '#1' '#1:2' '#1:3' '#2:2' '#2:3' '#3:2' '#3:3' '#4:2' '#4:3' '#5:2' '#5:3' '#6:2' '#6:3' '#7:2' '#7:3' '#8:2' '#8:3' '#9:2' '#9:3' '#10:2' '#10:3'  > sig.inp
 
#cp sig.inp onreal

#cd onreal; ./iterate.py ham.dat; cd ../