私はプログラマーではありませんが、ファイルを切り取り、そのファイルから多次元配列を生成するスクリプトを作成する必要があります。

私は化学者なので、配列は次のようになります。

[pe][pH][element][concentration].

これにより、個々の要素の濃度変化をグラフ化できます。pHそして着る。

だから私は入力ファイルとインデックスを切り取るためにこのスクリプトを書いています。着る、pHそして鐘、しかし今は詰まっています。私はそれを達成する方法を知りません。

    #!/bin/bash
    inputFile='test.inp'

    awk '/Distribution of species/ {f=1} /Saturation indices/ {f=0} f' $inputFile > distriPack.inout
   csplit -z -f distri_ distriPack.inout /Distribution/ {*}
   sed -i '1,5d; $d' distri_*
   sed -i 's/^ *//' distri_*
   cut -d ' ' -f1 distri_* | sort | uniq > especes.inp
   grep -E "pH  =   " $inputFile | cut -d '=' -f 2 | cut -d ' ' -f 4 > ph.inp
   grep -E "pe  =   " $inputFile | cut -d '=' -f 2 | cut -d ' ' -f 4 > pe.inp

   declare -a tableVariable=(ipe ipH ies)

  ipe=0 
  tablepe=0
  for pe in cat "$pe.inp"
  do 
      if ([ $ipe -eq 0 ] || pePrev=pe)
          then 
           tablepe=$((tablepe+1))
           ipe=$((ipe+1))
           pePrev=$pe
           fi
             ipH=0
             tablepH=0
             for pH in cat "$pH.inp"
             do 
             if ([ $ipH -eq 0 ] || pHPrev=pH )
             then 
             tablepH=$((tablepH+1)) > tablepH.inp
             ipH=$((ipH+1))
             pHPrev=$pH
             fi
                ies=0
                tablees=0
                for espece in cat "$espece.inp"
                do
                    if ([ $ies -eq 0 ] || especePrev=espece )
                    then 
            #                       indEspece=0
            #                       for espece in distri_$pe
            #                       do 
            #                          if ([ $indEspece -eq 0 ])
                       maVariable= grep -E "espece" < distri_10  | cut -d ' ' -f 1 

                       ${tableVariable[ $((ipe)) $((ipH)) $((ies)) ]}=$((maVariable))
          #                     done            
                    ies=$((ies+1))
                    tablees=$((tablees+1))
                    especePrev=$espece
                    fi  
                    done
                   done

           done 

val=${monTableau[${i_pe} ${i_pH} ${es} ${pa}]}>monTableau.inp

私は次のようなものが欲しい：for [pe][pH][specie]=[0][12][15]

インデックステーブル = 濃度。

しかし、bashとDistri_（pe / pH index value）ファイルの2番目の列で、インデックスkの種類に対応for pe = i, pH = jする行を表現する方法がわかりません。私の入力ファイルの「一部」行は次のとおりです。specie=kdo indexTable[i][j][k]=value

    >Initial solution 1. \newline
    >Description of solution 
    >pH = 0.0 
    >pe = 0.0 
    >Distribution of species 
    >Species Molality log Activity 
    >H+       1.1e+00 1.0e+00 
    >OH-      1.5e-14 9.5e-15 
    >Am+2     0.0e+00 0.0e+00


    >Initial solution 2. pe 0 pH 0.5 
    >Description of solution 
    >pH = 0.5 
    >pe = 0.0 
    >Distribution of species 
    >Species Molality Log Activity 
    >H+ 4.1e-01 3.1e-01 
    >OH- 4.5e-14 3.1e-14 
    >Am+2 0.0e+00 0.0e+00 

私は次のようなものを得たいと思います： outputfile_[pe value].out、ここでデータは次のようにソートされます

    >Column 1    C2 .....      Cn 
    >Specie\pH   0             n 
    >H+          [H+] at pH=0  [H+] at pH=n 

など