# Last edited on 2014-03-19 19:12:41 by stolfilocal
# gawk 'BEGIN{ Q = 0.400; DB = 16; } /^[ ]*[0-9]/{ d = $1; p = $2; w = $3; z = exp(log(Q)*(d-DB)); printf "%3d %6.1f %5.3f %3.1f %9.5f\n", d, p, w, 1, z; next; } //{ print; }'

 16 3842.0 0.580 1.0   1.00000
 17 3758.0 0.957 1.0   0.40000
 18 3718.0 0.432 1.0   0.16000
 19 3746.0 0.865 1.0   0.06400