Discovery解释结果可以输出Landmark single格式的层位,但由于discovery不支持fault flag(就是single格式文件中最后一列的F),因此损失了这个重要的信息。在一些软件中,同一条测线上层位就会全部连接在一起,很是不好。
层位解释中,一般会有几道处于断面上而无层位解释(正断层)。这点正好可以用来弥补single格式层位输出的不足。
做法:
1. 在层位输出时,把NULL值选项选中,并保留缺省的1.0E37
2. 利用下面简单的awk程序稍作转换即可,awk是UNIX, Linux里面的脚本程序。
3. 用法: awk -f single.awk disc_output_single > single_with_faultflag
够简单吧
single.awk:
#! /bin/awk -f
#
#
#
BEGIN{nullAgain="YES"; flag=" "}
{
if( $5 ~ /1.0E37/ )
nullAgain = "YES"
else {
if (nullAgain ~ /YES/) {
flag = " F"
nullAgain = "NO"
}
else flag = " "
printf "%30.2f%10.2f%11.2f%12.2f%12.4f%2s\n", $1,$2,$3,$4,$5,flag
}
}