View of /trunk/gadgetMerluza2015/model/hke.predict.lt.sh

#!/bin/sh

# shell script to run long term projections (50 years)
# use params.pred, hke.rec, hke.fleet, hke.amount.pred

for i in 0.001 0.01 0.02 0.03 0.04 0.05 0.06 0.07 0.08 0.09 0.10 0.11 0.12 0.13 0.14 0.15 0.16 0.17 0.18 0.19 0.20 0.21 0.22 0.23 0.24 0.25 0.26 0.27 0.28 0.29 0.30 0.31 0.32 0.33 0.34 0.35 0.36 0.37 0.38 0.39 0.40 0.41 0.42 0.43 0.44 0.45 0.46 0.47 0.48 0.49 0.50 0.51 0.52 0.53 0.54 0.55 0.56 0.57 0.58 0.59 0.60 0.61 0.62 0.63 0.64 0.65 0.66 0.67 0.68 0.69 0.70 0.71 0.72 0.73 0.74 0.75 0.76 0.77 0.78 0.79 0.80 0.81 0.82 0.83 0.84 0.85 0.86 0.87 0.88 0.89 0.90 0.91 0.92 0.93 0.94 0.95 0.96 0.97 0.98 0.99 1 1.01 1.1 1.2 1.3 1.4 1.5 1.6 1.7 1.8 1.9 2

#for i in 0.10 0.20 0.30 0.40 0.50 0.60 0.70 0.80 0.90 1

do
    echo "New run: F = $i"
    sed -e "s/XX/$i/g" < params.prj> ../out/out.lt.prj/params.lt.$i    
    gadget -s -main hke.main.lt -i ../out/out.lt.prj/params.lt.$i 
    rm ../out/out.lt.prj/params.lt.$i 

    mv ../out/out.lt.prj/catch.lt ../out/out.lt.prj/catch.lt.$i
    mv ../out/out.lt.prj/F1_3.lt ../out/out.lt.prj/F1_3.lt.$i
    mv ../out/out.lt.prj/Flen.lt ../out/out.lt.prj/Flen.lt.$i
    mv ../out/out.lt.prj/land.lt ../out/out.lt.prj/land.lt.$i
    mv ../out/out.lt.prj/lenN.lt ../out/out.lt.prj/lenN.lt.$i
done

    rm ../out/out.lt.prj/error.lt 
# borramos el archivo creado en printfiles.lt para solventar problema (último printfile necesita acabar con allyears y no con 2050)