#pgf90 calculate_dd.f90
#./a.out
NGUONG="0.1 1 5 10 15 20 25 30 35 40 45 50 55 60 65 70 75 80 85 90 95 100 120 150 200 300"
#NGUONG="10"
# Neu chay 4GRID can chon tung vung va thay doi gia tri NT (cal_HSS_4GRID.f90)
#KKL:  time=1148
#TC:   time=1078
vung="TC"
ntime="1078"
for ng in $NGUONG ; do
#for reg in $Region ; do
##==== chon quy mo GRID hay region
cp HSS.f90 tmp.f90
sed 's/NNNN/'$ng'/g' -i tmp.f90
sed 's/AAAA/'$ntime'/g' -i tmp.f90
sed 's/VVVV/'$vung'/g' -i tmp.f90
#
pgf90 tmp.f90
./a.out
rm a.out
rm tmp.f90
##==== chon quy mo GRID hay region
cat KQ_${vung}.txt >> ALL_KQ_${vung}.txt
#rm tmp.f90
rm  KQ_${vung}.txt
done
#done
