从查找表的特定行替换

Question

我正在尝试获取一个脚本来使用 GAMESS 包自动执行一些任务，我希望以后可以从中推断出更复杂的情况。唉，看来我的 Unix 编程技能达不到标准。

我有一个通用的 GAMESS 输入文件 'ion.inp' 的形式：

 $CONTRL SCFTYP=<tag4> ICHARG=<tag5> MULT=<tag6> ISPHER=1 NPRINT=-5 $END
 $BASIS GBASIS=<tag> $END
 $DATA
 <tag1> energy
 Dnh 2

 <tag2> <tag3>  .0  .0  .0
 $END

我（作为 MWE）有一个“ion.inp”参数的查找表，例如“table.dat”，其中 取自表的每一行。

<tag1>          | <tag2> | <tag3> | <tag4> | <tag5> | <tag6>
Hidrogen        | H      | 1.0    | ROHF   | 0      | 2
Hidrogen cation | H      | 1.0    | RHF    | 1      | 1

为了可移植性，我想获得一个使用 POSIX sh、sed 或 awk 的解决方案，但经过一些试验（使用 sh 或 sed，我对 awk 一点也不熟悉，尽管我知道它是一种潜力在这种情况下的解决方案）我无法让它工作。

文件 'ion.inp' 可以就地编辑，因为它将在 sh 循环中运行。除了这个所谓的简单替换之外，我已经完成了所有其他工作。

任何帮助将不胜感激！

Answer 1

这是一个使用sed 和awk 的示例。该脚本名为script.awk。一切都在标准输出上输出，但您可以使用 AWK 脚本中的> 将其重定向到文件。请参阅下面的第二种解决方案。

我们的想法是使用awk 和包含进行替换的数据的table.dat 文件来驱动进程，然后对于每批替换（文件的每一行），我们使用sed 来执行一旦我们有了每个标签和它的值，实际的替换。

• BEGIN { FS = "\\s*\\|\\s*" } 将字段分隔符更改为使用“可选空格后跟 | 然后再跟可选空格”。这意味着$1、$2、...将为我们提供编号为 1、2、...的标签的值
• NR == 1 { next } 用于跳过无用的第一行，因为标签从 1 开始排序，没有任何间隙。如果不是这样，我们将不得不调整 AWK 脚本。
• { ... } 为每一行我们构建sed 命令并执行它。 sed 的输出成为该特定行的 awk 的输出。

BEGIN { FS = "\\s*\\|\\s*" }
NR == 1 { next }
{
  s = ""
  for (i=1; i <= NF; i++)
    s = s sprintf(";s/%s/%s/g", "<tag" i ">", $i)
  system("sed '" s "' ion.inp")
}

$ cat table.dat
<tag1>          | <tag2> | <tag3> | <tag4> | <tag5> | <tag6>
Hidrogen        | H      | 1.0    | ROHF   | 0      | 2
Hidrogen cation | H      | 1.0    | RHF    | 1      | 1

$ cat ion.inp
$CONTRL SCFTYP=<tag4> ICHARG=<tag5> MULT=<tag6> ISPHER=1 NPRINT=-5 $END
$BASIS GBASIS=<tag> $END
$DATA
<tag1> energy
Dnh 2

<tag2> <tag3>  .0  .0  .0
$END

$ awk -f script.awk table.dat
$CONTRL SCFTYP=ROHF ICHARG=0 MULT=2 ISPHER=1 NPRINT=-5 $END
$BASIS GBASIS=<tag> $END
$DATA
Hidrogen energy
Dnh 2

H 1.0  .0  .0  .0
$END
$CONTRL SCFTYP=RHF ICHARG=1 MULT=1 ISPHER=1 NPRINT=-5 $END
$BASIS GBASIS=<tag> $END
$DATA
Hidrogen cation energy
Dnh 2

H 1.0  .0  .0  .0
$END

使用命令awk -f script2.awk table.dat 重定向到文件。脚本script2.awk是：

BEGIN { FS = "\\s*\\|\\s*" }
NR == 1 { next }
{
  s = ""
  for (i=1; i <= NF; i++)
    s = s sprintf(";s/%s/%s/g", "<tag" i ">", $i)
  system("sed '" s "' ion.inp > " sprintf("output%02d.txt", NR-1))
}

$ cat output01.txt
$CONTRL SCFTYP=ROHF ICHARG=0 MULT=2 ISPHER=1 NPRINT=-5 $END
$BASIS GBASIS=<tag> $END
$DATA
Hidrogen energy
Dnh 2

H 1.0  .0  .0  .0
$END

$ cat output02.txt
$CONTRL SCFTYP=RHF ICHARG=1 MULT=1 ISPHER=1 NPRINT=-5 $END
$BASIS GBASIS=<tag> $END
$DATA
Hidrogen cation energy
Dnh 2

H 1.0  .0  .0  .0
$END