# Last edited on 2014-03-24 02:01:15 by stolfilocal
# gawk 'BEGIN{ Q = 0.901; DB = 19; } /^[ ]*[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; }'

 19 3746.0 0.869 1.0   1.00000
 20 3646.0 0.694 1.0   0.90100
 21 3547.0 0.979 1.0   0.81180
 22 3462.0 0.958 1.0   0.73143
 23 3469.0 0.573 1.0   0.65902