Go to the documentation of this file.00001 Info<< "Reading field h\n" << endl;
00002 volScalarField h
00003 (
00004 IOobject
00005 (
00006 "h",
00007 runTime.timeName(),
00008 mesh,
00009 IOobject::MUST_READ,
00010 IOobject::AUTO_WRITE
00011 ),
00012 mesh
00013 );
00014
00015 Info<< "Reading field h0 if present\n" << endl;
00016 volScalarField h0
00017 (
00018 IOobject
00019 (
00020 "h0",
00021 runTime.findInstance("polyMesh", "points"),
00022 mesh,
00023 IOobject::READ_IF_PRESENT
00024 ),
00025 mesh,
00026 dimensionedScalar("h0", dimLength, 0.0)
00027 );
00028
00029 Info<< "Reading field U\n" << endl;
00030 volVectorField U
00031 (
00032 IOobject
00033 (
00034 "U",
00035 runTime.timeName(),
00036 mesh,
00037 IOobject::MUST_READ,
00038 IOobject::AUTO_WRITE
00039 ),
00040 mesh
00041 );
00042
00043 Info<< "Creating field hU\n" << endl;
00044 volVectorField hU
00045 (
00046 IOobject
00047 (
00048 "hU",
00049 runTime.timeName(),
00050 mesh
00051 ),
00052 h*U,
00053 U.boundaryField().types()
00054 );
00055
00056 Info<< "Creating field hTotal for post processing\n" << endl;
00057 volScalarField hTotal
00058 (
00059 IOobject
00060 (
00061 "hTotal",
00062 runTime.timeName(),
00063 mesh,
00064 IOobject::READ_IF_PRESENT,
00065 IOobject::AUTO_WRITE
00066 ),
00067 h+h0
00068 );
00069 hTotal.write();
00070
00071 # include "createPhi.H"
00072
00073 Info<< "Creating Coriolis Force" << endl;
00074 const dimensionedVector F("F", ((2.0*Omega) & gHat)*gHat);