save MultiScalarField3D in Array?

Palabos General
#1

Hello everybody,

I need to compare two MultiScalarFields at different timesteps.
For that I testet the possibility to save the field and load it at a later timestep, but that doesn’t work. At timestep “501” the field is full of “0”.
Can anybody help me?

[code=“cpp”]

MultiScalarField3D V1(lattice, velSection);

if (iT%1000 == 0){
V1=(*computeVelocityNorm(lattice, velSection));
pcout << “V1_0=”<< setprecision(4) << V1 << endl;
}

if (iT%1000 == 501){
pcout << “V1_1=”<< setprecision(4) << V1 << endl;
}



I am thankful for any help