bin/foamGraphResKE

   1 #!/bin/csh -f
   2 #------------------------------------------------------------------------------
   3 # =========                 |
   4 # \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
   5 #  \\    /   O peration     |
   6 #   \\  /    A nd           | Copyright (C) 2004-2010 OpenCFD Ltd.
   7 #    \\/     M anipulation  |
   8 #-------------------------------------------------------------------------------
   9 # License
  10 #     This file is part of OpenFOAM.
  11 #
  12 #     OpenFOAM is free software: you can redistribute it and/or modify it
  13 #     under the terms of the GNU General Public License as published by
  14 #     the Free Software Foundation, either version 3 of the License, or
  15 #     (at your option) any later version.
  16 #
  17 #     OpenFOAM is distributed in the hope that it will be useful, but WITHOUT
  18 #     ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or
  19 #     FITNESS FOR A PARTICULAR PURPOSE.  See the GNU General Public License
  20 #     for more details.
  21 #
  22 #     You should have received a copy of the GNU General Public License
  23 #     along with OpenFOAM.  If not, see <http://www.gnu.org/licenses/>.
  24 #
  25 # Script
  26 #     foamGraphResKE
  27 #
  28 # Description
  29 #
  30 #------------------------------------------------------------------------------
  31 if (! -f "$1" ) then
  32     echo "No file: '$1'"
  33     exit 1
  34 endif
  35
  36 set tmpK = newK$$.dat
  37 cat $1 | grep "Solving for k" | grep -v "solution singularity" | sed s/,//g | awk 'BEGIN { NENTRIES = 0 }{NENTRIES++; printf("%f %e\n", NENTRIES, $8)}' - > $tmpK
  38
  39 set tmpEpsilon = newEpsilon$$.dat
  40 cat $1 | grep "Solving for epsilon" | grep -v "solution singularity" | sed s/,//g | awk 'BEGIN { NENTRIES = 0 }{NENTRIES++; printf("%f %e\n", NENTRIES, $8)}' - > $tmpEpsilon
  41
  42 set linesK = `cat $tmpK | wc -l`
  43 set linesEpsilon = `cat $tmpEpsilon | wc -l`
  44
  45 set ngraphs = 0
  46
  47 if ($linesK > 1) then
  48 @ ngraphs = $ngraphs + 1
  49 endif
  50
  51 if ($linesEpsilon > 1) then
  52 @ ngraphs = $ngraphs + 1
  53 endif
  54
  55 cat << EOF > residualKE.dat
  56 Solver Performance
  57 number of iterations
  58 residual
  59 0 0
  60 0 0
  61 $ngraphs
  62 EOF
  63
  64 if ($linesK > 1) then
  65 cat << EOF >> residualKE.dat
  66 k
  67 0
  68 $linesK
  69 EOF
  70 cat $tmpK >> residualKE.dat
  71 endif
  72
  73 if ($linesEpsilon > 1) then
  74 cat << EOF >> residualKE.dat
  75 epsilon
  76 0
  77 $linesEpsilon
  78 EOF
  79 cat $tmpEpsilon >> residualKE.dat
  80 endif
  81
  82 rm $tmpK
  83 rm $tmpEpsilon
  84
  85 #------------------------------------------------------------------------------