Go to the documentation of this file.00001
00002
00003
00004
00005
00006
00007
00008
00009
00010
00011
00012
00013
00014
00015
00016
00017
00018
00019
00020
00021
00022
00023
00024
00025
00026 template<class Type>
00027 void Foam::calcTypes::interpolate::writeInterpolateField
00028 (
00029 const IOobject& header,
00030 const fvMesh& mesh,
00031 bool& processed
00032 )
00033 {
00034 typedef GeometricField<Type, fvPatchField, volMesh> fieldType;
00035 typedef GeometricField<Type, fvsPatchField, surfaceMesh> surfaceFieldType;
00036
00037 if (header.headerClassName() == fieldType::typeName)
00038 {
00039 Info<< " Reading " << header.name() << endl;
00040 fieldType field(header, mesh);
00041
00042 Info<< " Calculating interpolate" << header.name() << endl;
00043 surfaceFieldType interpolateField
00044 (
00045 IOobject
00046 (
00047 "interpolate" + header.name(),
00048 mesh.time().timeName(),
00049 mesh,
00050 IOobject::NO_READ
00051 ),
00052 fvc::interpolate(field)
00053 );
00054 interpolateField.write();
00055
00056 processed = true;
00057 }
00058 }
00059
00060
00061