# Last edited on 2014-03-06 19:47:27 by stolfilocal # gawk 'BEGIN{ Q = 0.464; DB = 32; } /^[ ]*[0-9]/{ d = $1; p = $2; z = exp(log(Q)*(d-DB)); print d, p, 1, z; }' 32 4040 1 1 33 3956 1 0.464 34 3917 1 0.215296