splitrayden.awk

一个地震层析成像的MGT工具包

# awk -script to split ray-density tensor output file from
# velsort (raydenxx.xx.xyz) into class-files for plotting
#
# classes are a combination of ray distr.(=qual)
# and total weighted ray length (=Nrayf)
#

BEGIN{lnfac = 0.4343	# define 1/ln(10) for calculating log from ln (log in awk is ln)
      vlen  = 0.3}	# define plot-vector length
{
 if ($1 == 1) {		# plot circles, no further classification 
    Nray = log($3) * lnfac           # (circle diam. in col. 3 ofoutput)
    Nrayf = int(Nray)
    if ((Nrayf==1) or (Nrayf==2)) print $4,$5,Nrayf >"qual1_1.out"
    if (Nrayf==3) print $4, $5, Nrayf >"class2_1"
    }
 if ($1 == 2) {         # qual==2
    Nray = log($3) * lnfac
    Nrayf = int(Nray)
    if (Nrayf == 3) {
     print ">" >"qual1.out"
     print $4, $5 >"qual1.out"
     print $6, $7 >"qual1.out"
     }
    if (Nrayf == 2) {
     print ">" >"qual2.out"
     print $4, $5 >"qual2.out"
     print $6, $7 >"qual2.out"
     }
    if (Nrayf <= 1) {
     print ">" >"qual3.out"
     print $4, $5 >"qual3.out"
     print $6, $7 >"qual3.out"
     }
 }
 if ($1 == 3) {         # qual==3
    Nray = log($3) * lnfac
    Nrayf = int(Nray)
    if (Nrayf == 3) {
     print ">" >"qual2.out"
     print $4, $5 >"qual2.out"
     print $6, $7 >"qual2.out"
     }
    if (Nrayf == 2) {
     print ">" >"qual3.out"
     print $4, $5 >"qual3.out"
     print $6, $7 >"qual3.out"
     }
    if (Nrayf <= 1) {
     print ">" >"qual4.out"
     print $4, $5 >"qual4.out"
     print $6, $7 >"qual4.out"
     }
 }
 if ($1 == 4) {         # plot vector sized by Nrayf, colored by qual for $1 = [2,3,4]
    Nray = log($3) * lnfac
    Nrayf = int(Nray)
    if (Nrayf == 3) {
     print ">" >"qual3.out"
     print $4, $5 >"qual3.out"
     print $6, $7 >"qual3.out"
     }
    if (Nrayf == 2) {
     print ">" >"qual4.out"
     print $4, $5 >"qual4.out"
     print $6, $7 >"qual4.out"
     }
    if (Nrayf <= 1) {
     print ">" >"qual4.out"
     print $4, $5 >"qual4.out"
     print $6, $7 >"qual4.out"
     }
 }
}
END{ }