最初の列に基づいて行を Grep

Question 1

^CC入力行の先頭ではなく、前にスペースがあるため機能しません。grep '^ C' "$1"あなたが欲しいものをする必要があります。

（ところで、これをgrep | wc -l使用できますgrep -c。ああ、行の引用符がecho少し奇妙です。変数を引用符の中に入れるだけです。）

Answer

^CC入力行の先頭ではなく、前にスペースがあるため機能しません。grep '^ C' "$1"あなたが欲しいものをする必要があります。

（ところで、これをgrep | wc -l使用できますgrep -c。ああ、行の引用符がecho少し奇妙です。変数を引用符の中に入れるだけです。）

Question 2

awkで完全に試してみてください。たとえば、次のスクリプトは2つの配列を使用して、指定されたatoms原子の各入力行を記憶し、count各原子の数を保持します。入力ファイル全体を読み込んだ後、目的の形式でデータを出力します。

/^ [[:alpha:]]/ {
  if (count[$1] == 0) {
    atoms[$1]=$0;
  } else {
    atoms[$1]=atoms[$1] "\n" $0;
  }
    count[$1]++
}

END {
  atypes = length(count);

  print "BASIS\nco2\n\n"

  print "Atomtypes=" atypes " Generators=0 Integrals=1.00D-15 Angstrom"

  print "Charge=6.0 Atoms=" count["C"]
  print atoms["C"]

  print "Charge=1.0 Atoms=" count["H"]
  print atoms["H"]

  if (count["O"] > 0) {
    print "Charge=8.0 Atoms=" count["O"]
    print atoms["O"]
  }
}

各原子に関連するチャージルックアップテーブルがあれば、大幅に改善され、通常の変換スクリプトに置き換えることができますが、必ずしも必要ではありません。変換ユーティリティはすでに存在します。obabelから呼び出されます。Open Babel: オープンソース化学ツールボックスxyzプロジェクトを変換し、次のような複数の化学ファイル形式間で変換できますdalmol。

obabel -i xyz input.xyz -o dalmol -O output.dalmol

olabel -L formats | lessサポートされている形式の完全なリストを取得するには実行してください。

DebianまたはDebian派生製品（Ubuntu、Mintなど）を実行している場合はインストールできますapt-get install openbabel。梱包手順は次のとおりです。

Package: openbabel
Version: 3.1.1+dfsg-6
Installed-Size: 630
Maintainer: Debichem Team <[email protected]>
Architecture: amd64
Depends: libc6 (>= 2.14), libgcc-s1 (>= 3.0), libopenbabel7 (>= 3.1.1+dfsg), libstdc++6 (>= 5.2)
Description-en: Chemical toolbox utilities (cli)
 Open Babel is a chemical toolbox designed to speak the many languages of
 chemical data. It allows one to search, convert, analyze, or store data from
 molecular modeling, chemistry, solid-state materials, biochemistry, or related
 areas.  Features include:
 .
  * Hydrogen addition and deletion
  * Support for Molecular Mechanics
  * Support for SMARTS molecular matching syntax
  * Automatic feature perception (rings, bonds, hybridization, aromaticity)
  * Flexible atom typer and perception of multiple bonds from atomic coordinates
  * Gasteiger-Marsili partial charge calculation
 .
 File formats Open Babel supports include PDB, XYZ, CIF, CML, SMILES, MDL
 Molfile, ChemDraw, Gaussian, GAMESS, MOPAC and MPQC.
 .
 This package includes the following utilities:
  * obabel: Convert between various chemical file formats
  * obenergy: Calculate the energy for a molecule
  * obminimize: Optimize the geometry, minimize the energy for a molecule
  * obgrep: Molecular search program using SMARTS pattern
  * obgen: Generate 3D coordinates for a molecule
  * obprop: Print standard molecular properties
  * obfit: Superimpose two molecules based on a pattern
  * obrotamer: Generate conformer/rotamer coordinates
  * obconformer: Generate low-energy conformers
  * obchiral: Print molecular chirality information
  * obrotate: Rotate dihedral angle of molecules in batch mode
  * obprobe: Create electrostatic probe grid

Answer

awkで完全に試してみてください。たとえば、次のスクリプトは2つの配列を使用して、指定されたatoms原子の各入力行を記憶し、count各原子の数を保持します。入力ファイル全体を読み込んだ後、目的の形式でデータを出力します。

/^ [[:alpha:]]/ {
  if (count[$1] == 0) {
    atoms[$1]=$0;
  } else {
    atoms[$1]=atoms[$1] "\n" $0;
  }
    count[$1]++
}

END {
  atypes = length(count);

  print "BASIS\nco2\n\n"

  print "Atomtypes=" atypes " Generators=0 Integrals=1.00D-15 Angstrom"

  print "Charge=6.0 Atoms=" count["C"]
  print atoms["C"]

  print "Charge=1.0 Atoms=" count["H"]
  print atoms["H"]

  if (count["O"] > 0) {
    print "Charge=8.0 Atoms=" count["O"]
    print atoms["O"]
  }
}

各原子に関連するチャージルックアップテーブルがあれば、大幅に改善され、通常の変換スクリプトに置き換えることができますが、必ずしも必要ではありません。変換ユーティリティはすでに存在します。obabelから呼び出されます。Open Babel: オープンソース化学ツールボックスxyzプロジェクトを変換し、次のような複数の化学ファイル形式間で変換できますdalmol。

obabel -i xyz input.xyz -o dalmol -O output.dalmol

olabel -L formats | lessサポートされている形式の完全なリストを取得するには実行してください。

DebianまたはDebian派生製品（Ubuntu、Mintなど）を実行している場合はインストールできますapt-get install openbabel。梱包手順は次のとおりです。

Package: openbabel
Version: 3.1.1+dfsg-6
Installed-Size: 630
Maintainer: Debichem Team <[email protected]>
Architecture: amd64
Depends: libc6 (>= 2.14), libgcc-s1 (>= 3.0), libopenbabel7 (>= 3.1.1+dfsg), libstdc++6 (>= 5.2)
Description-en: Chemical toolbox utilities (cli)
 Open Babel is a chemical toolbox designed to speak the many languages of
 chemical data. It allows one to search, convert, analyze, or store data from
 molecular modeling, chemistry, solid-state materials, biochemistry, or related
 areas.  Features include:
 .
  * Hydrogen addition and deletion
  * Support for Molecular Mechanics
  * Support for SMARTS molecular matching syntax
  * Automatic feature perception (rings, bonds, hybridization, aromaticity)
  * Flexible atom typer and perception of multiple bonds from atomic coordinates
  * Gasteiger-Marsili partial charge calculation
 .
 File formats Open Babel supports include PDB, XYZ, CIF, CML, SMILES, MDL
 Molfile, ChemDraw, Gaussian, GAMESS, MOPAC and MPQC.
 .
 This package includes the following utilities:
  * obabel: Convert between various chemical file formats
  * obenergy: Calculate the energy for a molecule
  * obminimize: Optimize the geometry, minimize the energy for a molecule
  * obgrep: Molecular search program using SMARTS pattern
  * obgen: Generate 3D coordinates for a molecule
  * obprop: Print standard molecular properties
  * obfit: Superimpose two molecules based on a pattern
  * obrotamer: Generate conformer/rotamer coordinates
  * obconformer: Generate low-energy conformers
  * obchiral: Print molecular chirality information
  * obrotate: Rotate dihedral angle of molecules in batch mode
  * obprobe: Create electrostatic probe grid

最初の列に基づいて行を Grep

答え1

答え2

関連情報