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) 1991-2009 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 the
  14 #     Free Software Foundation; either version 2 of the License, or (at your
  15 #     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, write to the Free Software Foundation,
  24 #     Inc., 51 Franklin St, Fifth Floor, Boston, MA 02110-1301 USA
  25 #
  26 # Script
  27 #     foamGraphResKE
  28 #
  29 # Description
  30 #
  31 #------------------------------------------------------------------------------
  32
  33 set tmpK = newK$$.dat
  34 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
  35
  36 set tmpEpsilon = newEpsilon$$.dat
  37 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
  38
  39 set linesK = `cat $tmpK | wc -l`
  40 set linesEpsilon = `cat $tmpEpsilon | wc -l`
  41
  42 set ngraphs = 0
  43
  44 if ($linesK > 1) then
  45 @ ngraphs = $ngraphs + 1
  46 endif
  47
  48 if ($linesEpsilon > 1) then
  49 @ ngraphs = $ngraphs + 1
  50 endif
  51
  52 cat << EOF > residualKE.dat
  53 Solver Performance
  54 number of iterations
  55 residual
  56 0 0
  57 0 0
  58 $ngraphs
  59 EOF
  60
  61 if ($linesK > 1) then
  62 cat << EOF >> residualKE.dat
  63 k
  64 0
  65 $linesK
  66 EOF
  67 cat $tmpK >> residualKE.dat
  68 endif
  69
  70 if ($linesEpsilon > 1) then
  71 cat << EOF >> residualKE.dat
  72 epsilon
  73 0
  74 $linesEpsilon
  75 EOF
  76 cat $tmpEpsilon >> residualKE.dat
  77 endif
  78
  79 rm $tmpK
  80 rm $tmpEpsilon
  81
  82 #------------------------------------------------------------------------------