#!/bin/sh
#------------------------------------------------------------------------------
# =========                 |
# \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
#  \\    /   O peration     |
#   \\  /    A nd           | Copyright (C) 2011-2017 OpenFOAM Foundation
#    \\/     M anipulation  |
#-------------------------------------------------------------------------------
# License
#     This file is part of OpenFOAM.
#
#     OpenFOAM is free software: you can redistribute it and/or modify it
#     under the terms of the GNU General Public License as published by
#     the Free Software Foundation, either version 3 of the License, or
#     (at your option) any later version.
#
#     OpenFOAM is distributed in the hope that it will be useful, but WITHOUT
#     ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or
#     FITNESS FOR A PARTICULAR PURPOSE.  See the GNU General Public License
#     for more details.
#
#     You should have received a copy of the GNU General Public License
#     along with OpenFOAM.  If not, see <http://www.gnu.org/licenses/>.
#
# Script
#     foamGraphResUVWP
#
# Description
#     Deprecated script extract initial velocity and pressure residuals for each
#     time-step from a log file for graphing.
#
#     Superseded by the more general foamLog script.
#------------------------------------------------------------------------------
Script=${0##*/}

if [ ! -f "$1" ]; then
    echo "$Script: $1 is not a valid filename"
    exit 1
fi

scan()
{
    tmpFile=new$1$$.dat
    cat $2 | \
        grep "Solving for $1" | \
        grep -v "solution singularity" | \
        sed s/,//g | \
        awk 'BEGIN { NENTRIES = 0 }{NENTRIES++; printf("%d %e\n", NENTRIES, $8)}' - > $tmpFile

    if [ "$(cat $tmpFile | wc -l)" -gt 1 ]; then
        echo >> residualUVWP.dat
        cat $tmpFile >> residualUVWP.dat
    fi

    rm $tmpFile
}

rm -f residualUVWP.dat
scan "Ux" $1
scan "Uy" $1
scan "Uz" $1
scan "p" $1

#------------------------------------------------------------------------------
