[rucky96]

Hi all,

I am trying to save some data from some variables that I am working out every time step with:

        energyFile << runTime.timeName() << "\t\t"
                   << Em.weightedAverage(mesh.V()).value() << "\t\t" //Em is volScalarField
                   << Ek.weightedAverage(mesh.V()).value() << "\t\t"//Ek is volScalarField
                   << Et.weightedAverage(mesh.V()).value() << "\t\t"//Et is volScalarField
                   << max(mag(J)).value() << "\t\t"//J is volVectorField
                   << a.weightedAverage(mesh.V()).value()+b.weightedAverage(mesh.V()).value() << "\t\t" //a & b are volScalarField
                   << pR << "\t\t" ////pR is scalar
                   << endl;

When I run the simulation I obtain:

0.01            0.271417                0.000693869             0.272111                3.13241         0.0263245               5.1935
00.02           0.271321                0.000464667             0.271786                3.44374         0.0238673               4.8994
00.03           0.271221                0.00033977              0.271561                3.48565         0.0224566               5.15044 000.04          0.271117                0.000268453             0.271386                3.47999         0.021592                4.63075 000.05          0.271011                0.000229149             0.27124         3.41899         0.0210318               4.6224100.06            0.270903                0.000210638             0.271114                3.29874         0.0206501               4.8617
000.07          0.270795                0.000206592             0.271001                3.14365         0.0203772               4.84191 00.08           0.270686                0.000213266             0.270899                2.99081         0.0201733               4.81369

which is awfully saved and I cannot plot it.

However, if i remove the scalar "pR"

        energyFile << runTime.timeName() << "\t\t"
                   << Em.weightedAverage(mesh.V()).value() << "\t\t" //Em is volScalarField
                   << Ek.weightedAverage(mesh.V()).value() << "\t\t"//Ek is volScalarField
                   << Et.weightedAverage(mesh.V()).value() << "\t\t"//Et is volScalarField
                   << max(mag(J)).value() << "\t\t"//J is volVectorField
                   << a.weightedAverage(mesh.V()).value()+b.weightedAverage(mesh.V()).value() << "\t\t" //a & b are volScalarField
                   //<< pR << "\t\t" ////pR is scalar
                   << endl;

the data is properly saved:

0.01            0.271417                0.000693869             0.272111                3.13241         0.0263245
0.02            0.271321                0.000464667             0.271786                3.44374         0.0238673
0.03            0.271221                0.00033977              0.271561                3.48565         0.0224566
0.04            0.271117                0.000268453             0.271386                3.47999         0.021592
0.05            0.271011                0.000229149             0.27124                 3.41899         0.0210318
0.06            0.270903                0.000210638             0.271114                3.29874         0.0206501
0.07            0.270795                0.000206592             0.271001                3.14365         0.0203772
0.08            0.270686                0.000213266             0.270899                2.99081         0.0201733

I am compiling in a Linux server and my laptop has a Windows 10 Pro OS in case it has something to do.

Does anyone known what is happening and how can I save a scalar?

Regards.

[Danubi]

Are you having the same problem running the case in parallel and in serial?


In case you run in parallel you should use:

if(Pstream::master())

and the printing commands within { }.


Hope this helps

[rucky96]

Hi Danubi, thanks for answering.

I tried the following:

        if(Pstream::master())
        {
                energyFile << runTime.timeName() << "\t\t" 
                        << Em.weightedAverage(mesh.V()).value() << "\t\t"
                        << Ek.weightedAverage(mesh.V()).value() << "\t\t"
                        << Et.weightedAverage(mesh.V()).value() << "\t\t"
                        << max(mag(J)).value() << "\t\t"
                        << a.weightedAverage(mesh.V()).value()+b.weightedAverage(mesh.V()).value() << "\t\t"
                        << pR << "\t\t"
                        << endl;
        }

But when the simulation ends the first time step, it is "frozen". Neither print the values ​​nor jump to the second time step nor returns FATAL ERROR. Do I have to add something else?

PS: The problem also disappears when I run it in serial

[rucky96]

In the end I have solved it, I took the idea of ​​Danubi and continued testing combinations until it has come out:

	energyFile << runTime.timeName() << "\t\t" 
		<< Em.weightedAverage(mesh.V()).value() << "\t\t"
                << Ek.weightedAverage(mesh.V()).value() << "\t\t"
                << Et.weightedAverage(mesh.V()).value() << "\t\t"
                << max(mag(J)).value() << "\t\t"
                << a.weightedAverage(mesh.V()).value()+b.weightedAverage(mesh.V()).value() << "\t\t";
	if(Pstream::master()){energyFile << pR << "\t\t" << endl;}

Don't ask me why, I don't know what I did but works.