40 int main(
int argc,
char *argv[])
49 scalar fullCycle = 360.0;
54 while (runTime.theta() > ca0)
60 while (
mag(runTime.theta() - ca0) > eps)
62 scalar t0 = runTime.userTimeToTime(ca0 - runTime.theta());
63 runTime.setDeltaT(t0);
65 Info<<
"CA = " << runTime.theta() <<
endl;
69 scalar Vmax =
sum(
mesh.V().field());
71 while (
mag(runTime.theta()-ca1) > eps)
73 scalar t1 = runTime.userTimeToTime(ca1-runTime.theta());
74 runTime.setDeltaT(t1);
76 Info<<
"CA = " << runTime.theta() <<
endl;
80 scalar Vmin =
sum(
mesh.V().field());
82 Info<<
"\nVmax = " << Vmax;
84 Info<<
"Vmax/Vmin = " << Vmax/Vmin <<
endl;
Ostream & endl(Ostream &os)
Add newline and flush stream.
dimensioned< Type > sum(const DimensionedField< Type, GeoMesh > &df)
dimensioned< scalar > mag(const dimensioned< Type > &)