#!/bin/awk -f
## plot isoplanatic angle, seeing, time constant
## Use as: all.awk <mass-file>
## Oct 13, 2004: correct the seeing bias, roughly
BEGIN{
    system("rm -f tmp[123]")
}
{
  if ($1=="I") {A = $4; B = $6; GA = log(1.+A)/A; GB = log(1.+B)/B;
#  corr = sqrt(GA*GB)
  corr = 1./(1. + 0.7*A)^0.6
 }
    if ($1=="A") {
    ut=substr($3,1,2)+substr($3,4,2)/60+substr($3,7,2)/3600 
      if (ut>20) {ut=ut-24}
#      print ut, $16 >"tmp1"
#      print ut, $20 >"tmp2"
#    print ut, corr >"tmp1"
#      print ut, $4 >"tmp2"
#      print ut, $4*corr >"tmpMASS"
      print ut, $4*corr 
    }
}
# END{ system ("xmgrace   -param all.par tmp[123]")}