#!/bin/ksh
IAm=${0##*/}
lf="$IAm.$$"
xmin=-18
ymin=-2
zmin=-6
xmax=-8
ymax=8
zmax=4

totden=$(awk -v xmin="$xmin" -v xmax="$xmax" -v ymin="$ymin" -v ymax="$ymax" \
             -v zmin="$zmin" -v zmax="$zmax" -v lf="$lf" '
BEGIN { printf("scale=20\nx=0\n") # Initialize bc
}
$5!="" && $2<=xmax && $2>=xmin && $3<=ymax && $3>=ymin && $4<=zmax && $4>=zmin {
        if(match($5, /E([+])?/))
                printf("x+=(%s*10^%s)^2\n", substr($5, 1, RSTART - 1),
                        substr($5, RSTART + RLENGTH))
        else    printf("NR=%d: $5 bad format: %s\n", NR, $5) > lf
}
END {   printf("x\nquit\n")
}' "${1:-density.mesh_index}" | bc)
es=$?
if [ -s "$lf" ]
then    printf "%s: Log file (%s) is not empty\n" "$IAm" "$lf" >&2
        cat "$lf" >&2
        if [ "$es" -eq 0 ]
        then    es=1
        fi
fi
printf "totden: %s\n" "$totden"
exit $es

totden: 29418651008124.442882561190769208458247

#!/bin/bash
xmin=-18
ymin=-2
zmin=-6

xmax=-8
ymax=8
zmax=4

totden=0

o=$2
s=$3

for (( e = $o; e <= $s; e++ ))
do
        x=$( sed -n ${e}p $1 | awk '{print $2}' | sed 's/E/\*10\^/' | sed 's/\+//' )
        if [ "$x" == "" ]
        then    continue
        fi
        xbool=$( echo "$x <= $xmax && $x >= $xmin" | bc )
        if [ "$xbool" == 1 ]
        then
                y=$( sed -n ${e}p $1 | awk '{print $3}' | sed 's/E/\*10\^/' | sed 's/\+//' )
                ybool=$( echo "$y <= $ymax && $y >= $ymin" | bc )
                if [ "$ybool" == 1 ]
                then
                        z=$( sed -n ${e}p $1 | awk '{print $4}' | sed 's/E/\*10\^/' | sed 's/\+//' )
                        zbool=$( echo "$z <= $zmax && $z >= $zmin" | bc )
                        if [ "$zbool" == 1 ]
                        then
                                psi=$( sed -n ${e}p $1 | awk '{print $5}' | sed 's/E/\*10\^/' | sed 's/\+//' )
                                dens=$( echo "scale=20; ($psi)^2" | bc )
                                totden=$( echo "scale=20;$totden + $dens" | bc )
                        fi
                fi
        fi
done
echo $totden

29418651008124.442882561190769208458247

scale=20
x=0
x+=(7.357049247944539628087100*10^-005)^2
x+=(8.933437123800284495461800*10^-005)^2
x+=(1.032677570213522605163300*10^-004)^2
 ... ... ...
x+=(1.623299231409858035877900*10^-005)^2
x+=(9.208050322421728367170400*10^-006)^2
x+=(4.405340835827515605220000*10^-006)^2
x
quit