number=$( awk '{print $0}' number.dat )
gmx trjconv -f prod.xtc -s prod.tpr -o prod.mol.xtc -pbc mol<<EOF
0
EOF
for i in {1..5}
do
true
done
while IFS= read -r line
do
echo "$line"
done < file.txt
string="blabla"
sub-string=${string:3:0-1}${string:0-1}
mpirun -np 16 lmp_mpi -in input.lammps -var f 0.01
newline='variable = '$a
oldline=$(cat myfile| grep 'variable = ')
sed -i '/'"$oldline"'/c\'"$newline" myfile