applications/utilities/postProcessing/dataConversion/foamToGMV/gmvOutputLagrangian.H

   1 gmvFile << "tracers " << particles.size() << nl;
   2 forAllConstIter(Cloud<passiveParticle>, particles, iter)
   3 {
   4     gmvFile << iter().position().x() << ' ';
   5 }
   6 gmvFile << nl;
   7
   8 forAllConstIter(Cloud<passiveParticle>, particles, iter)
   9 {
  10     gmvFile << iter().position().y() << ' ';
  11 }
  12 gmvFile << nl;
  13
  14 forAllConstIter(Cloud<passiveParticle>, particles, iter)
  15 {
  16     gmvFile << iter().position().z() << ' ';
  17 }
  18 gmvFile << nl;
  19
  20 forAll(lagrangianScalarNames, i)
  21 {
  22     const word& name = lagrangianScalarNames[i];
  23
  24     IOField<scalar> fld
  25     (
  26         IOobject
  27         (
  28             name,
  29             runTime.timeName(),
  30             cloud::prefix,
  31             mesh,
  32             IOobject::MUST_READ,
  33             IOobject::NO_WRITE
  34         )
  35     );
  36
  37     if (fld.size())
  38     {
  39         gmvFile << name  << nl;
  40
  41         forAll(fld, n)
  42         {
  43             gmvFile << fld[n] << token::SPACE;
  44         }
  45         gmvFile << nl;
  46     }
  47
  48
  49 }
  50
  51 forAll(lagrangianVectorNames, i)
  52 {
  53     const word& name = lagrangianVectorNames[i];
  54
  55     IOField<vector> fld
  56     (
  57         IOobject
  58         (
  59             name,
  60             runTime.timeName(),
  61             cloud::prefix,
  62             mesh,
  63             IOobject::MUST_READ,
  64             IOobject::NO_WRITE
  65         )
  66     );
  67
  68     if (fld.size())
  69     {
  70         gmvFile << name + "x" << nl;
  71
  72         forAll(fld, n)
  73         {
  74             gmvFile << fld[n].x() << token::SPACE;
  75         }
  76         gmvFile << nl;
  77
  78         gmvFile << name + "y" << nl;
  79
  80         forAll(fld, n)
  81         {
  82             gmvFile << fld[n].y() << token::SPACE;
  83         }
  84         gmvFile << nl;
  85
  86         gmvFile << name + "z" << nl;
  87
  88         forAll(fld, n)
  89         {
  90             gmvFile << fld[n].z() << token::SPACE;
  91         }
  92         gmvFile << nl;
  93     }
  94
  95 }
  96
  97
  98 gmvFile << "endtrace"<< nl;