[jzlam]

I calculated field U and I need Ux,Uy,Uz in each point of mesh (together with mesh coordinates) like
1 x y z Ux Uy Uz
2 x y z Ux Uy Uz

to be able to import field to my program.

Can anyone give me code to do this? I am not succesfull to do it myself and I am little bit confused from C++ because I am programming in FORTRAN.

Thank for any advice

Jakub Zlamal

[nlc]

ave you tried something like this !

int main(int argc, char *argv[]) 
{ 

# include "setRootCase.H" 
# include "createTime.H" 
# include "createMesh.H" 
# include "createFields.H" 

  // Read  U field
  Info<< "Reading U field \n" << endl;
    volVectorField U
    (
        IOobject
        (
            "U",
            runTime.timeName(),
            mesh,
            IOobject::MUST_READ,
            IOobject::AUTO_WRITE
        ),
        mesh
    );

  pointField C = mesh.C()
  
  forAll(cell,U)
  {
  Info <<  C[cell] << U[Cell] << endl;
  }

  return(0); 
}

[nlc]

or

  forAll(cell,U)
   {
   Info <<  cell 
        <<  C[cell].x() 
        <<  C[cell].y() 
        <<  C[cell].z() 
        << U[cell].x() 
        << U[cell].y() 
        << U[cell].z() 
        << endl;
   }

Regards