bin/foamGraphResUVWP

   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 #     foamGraphResUVWP
  27 #
  28 # Description
  29 #
  30 #------------------------------------------------------------------------------
  31 if (! -f "$1" ) then
  32     echo "No file: '$1'"
  33     exit 1
  34 endif
  35
  36 set tmpU = newU$$.dat
  37 cat $1 | grep "Solving for Ux" | grep -v "solution singularity" | sed s/,//g | awk 'BEGIN { NENTRIES = 0 }{NENTRIES++; printf("%d %e\n", NENTRIES, $8)}' - > $tmpU
  38
  39 set tmpV = newV$$.dat
  40 cat $1 | grep "Solving for Uy" | grep -v "solution singularity" | sed s/,//g | awk 'BEGIN { NENTRIES = 0 }{NENTRIES++; printf("%d %e\n", NENTRIES, $8)}' - > $tmpV
  41
  42 set tmpW = newW$$.dat
  43 cat $1 | grep "Solving for Uz" | grep -v "solution singularity" | sed s/,//g | awk 'BEGIN { NENTRIES = 0 }{NENTRIES++; printf("%d %e\n", NENTRIES, $8)}' - > $tmpW
  44
  45 set tmpP = newP$$.dat
  46 cat $1 | grep "Solving for p" | grep -v "solution singularity" | sed s/,//g | awk 'BEGIN { NENTRIES = 0 }{NENTRIES++; printf("%d %e\n", NENTRIES, $8)}' - > $tmpP
  47
  48 set linesU = `cat $tmpU | wc -l`
  49 set linesV = `cat $tmpV | wc -l`
  50 set linesW = `cat $tmpW | wc -l`
  51 set linesP = `cat $tmpP | wc -l`
  52
  53 set ngraphs = 0
  54
  55 if ($linesU > 1) then
  56 @ ngraphs = $ngraphs + 1
  57 endif
  58
  59 if ($linesV > 1) then
  60 @ ngraphs = $ngraphs + 1
  61 endif
  62
  63 if ($linesW > 1) then
  64 @ ngraphs = $ngraphs + 1
  65 endif
  66
  67 if ($linesP > 1) then
  68 @ ngraphs = $ngraphs + 1
  69 endif
  70
  71 cat << EOF > residualUVWP.dat
  72 EOF
  73 #Solver Performance
  74 #number of iterations
  75 #residual
  76 #0 0
  77 #0 0
  78 #$ngraphs
  79 #EOF
  80
  81 if ($linesU > 1) then
  82 cat << EOF >> residualUVWP.dat
  83
  84 EOF
  85 #U
  86 #0
  87 #$linesU
  88 #EOF
  89 cat $tmpU >> residualUVWP.dat
  90 endif
  91
  92 if ($linesV > 1) then
  93 cat << EOF >> residualUVWP.dat
  94
  95 EOF
  96 #V
  97 #0
  98 #$linesV
  99 #EOF
 100 cat $tmpV >> residualUVWP.dat
 101 endif
 102
 103 if ($linesW > 1) then
 104 cat << EOF >> residualUVWP.dat
 105
 106 EOF
 107 #W
 108 #0
 109 #$linesW
 110 #EOF
 111 cat $tmpW >> residualUVWP.dat
 112 endif
 113
 114 if ($linesP > 1) then
 115 cat << EOF >> residualUVWP.dat
 116
 117 EOF
 118 #P
 119 #0
 120 #$linesP
 121 #EOF
 122 cat $tmpP >> residualUVWP.dat
 123 endif
 124
 125 rm $tmpU
 126 rm $tmpV
 127 rm $tmpW
 128 rm $tmpP
 129
 130 #------------------------------------------------------------------------------