Звоню его от поглазеть и скрипты bash, заговоры только первый участок

Итак, вот в чем дело. Я начинаю свою магистерскую диссертацию по вычислительному материаловедению и пытаюсь собрать несколько сценариев, чтобы помочь подготовиться к анализу данных.

Я готовил сценарий GAWK, который в основном возьмет некоторые данные (расположенные в 4 Столбцах) и захватит два из них и построит их в GNUPLOT. Чтобы достичь этой цели, я читаю в файле данных, который имеет несколько временных шагов и связанных с ними данных в нем, разделить файл на отдельные .файлы dat для каждого шаг во времени.

Оттуда я просто создаю базовый входной скрипт для GNUPLOT и строю каждый временной шаг, как они происходят в файле данных.

Проблема в том, что по какой-то причине все созданные участки являются точно такими же участками (в этом случае всегда первый временной шаг), но они сохраняются как правильный временной шаг.

Я уже просмотрел и отследил каждую переменную / имя файла на протяжении всего сценария и, наконец, определил, что проблема каким-то образом связана с вызовом GNUPLOT из a скрипт. Я достал свою системную команду и написал короткий сценарий bash, который вызывает gnuplot из цикла for:

#!/bin/bash
for file in ./*gnu
do
   gnuplot $file
done

И это по-прежнему вызывает ту же проблему, когда все сюжеты одинаковы. Затем я просто выполнил команду gnuplot *gnu из командной строки в каталоге, содержащем .файлы gnu, и это сработало.

Я думаю, мне просто интересно, есть ли какой-то буфер, который мне нужно смыть, или я просто что-то упускаю?

Сценарий GAWK приведен ниже. Я новичок в этом деле. это все еще так, поэтому, если вы хотите прокомментировать сценарий с некоторой конструктивной критикой, я был бы очень признателен и за это.

#!/opt/local/bin/gawk -v inputf=$1 -f                                                   

# Write gnuplot files and plot RDF data                                                 
function plot_rdf(timestep, Load_RDF_dat)
{
# Set number of digits in filenames to 6 so data is organized                           
    if (timestep < 10){
        pad_timestep="00000"timestep;
    }
    else if (timestep < 100){
        pad_timestep="0000"timestep;
    }
    else if (timestep < 1000){
        pad_timestep="000"timestep;
    }
    else if (timestep < 10000){
        pad_timestep="00"timestep;
    }
    else if (timestep < 100000){
        pad_timestep="0"timestep;
    }
    else{
        pad_timestep=timestep;
    }

# Give output filenames                                                                 
       gnu_file="plot_RDF_"pad_timestep".gnu";
       png_file="RDF_"pad_timestep".png";

# Create input files for gnuplot                                                        
       print "set output ""png_file""" >> gnu_file;
       print "set terminal png" >> gnu_file;
       print "plot './"Load_RDF_dat"' u 1:2" >> gnu_file;
       close(gnu_file);
       system("gnuplot "gnu_file);
}


# Main part of script                                                                   
{
# Parse the RDF data and save it to GNUPLOT readable files                              
    while(getline < inputf){
       if ($1 == "#"){
           # skips the three commented header lines                                     
           next;
       }
       else if (NF == 2){
           timestep=$1;
           bin_num=$2;
           print "Reading timestep "timestep;
           RDF_dat="RDF_"timestep".dat";
           next;
       }
       else if (NF == 4){
           print $2" "$3 >> RDF_dat;
           if ($1 == bin_num){
               plot_rdf(timestep, RDF_dat);
               close(RDF_dat);
           }
           next;
       }
    }
    close(inputf);
    close(RDF_dat);
 }

Фрагмент файла данных, который я читаю:

# Time-averaged data for fix rdf
# TimeStep Number-of-rows
# Row c_allrdf[1] c_allrdf[2] c_allrdf[3]
500 100
1 0.005 0 0
2 0.015 0 0
3 0.025 0 0
4 0.035 0 0
5 0.045 0 0
6 0.055 1.16597 0.00133333
7 0.065 2.08865 0.00466667
8 0.075 1.56958 0.008
9 0.085 0.733433 0.01
10 0.095 0.587288 0.012
600 100
1 0.005 0 0
2 0.015 0 0
3 0.025 2.79219 0.000666667
4 0.035 2.86766 0.002
5 0.045 0 0.002
6 0.055 0.582985 0.00266667
7 0.065 2.08865 0.006
8 0.075 0.62783 0.00733333
9 0.085 0.488955 0.00866667
10 0.095 1.17458 0.0126667

Обычно в каждом разделе timestep содержится 100 наборов данных, но я решил сократить их здесь, чтобы вы поняли идею.