void FileIO::flush(Timing timing, double dt)
{
if(timing.check(dataFileFlushTiming, dt)) H5Fflush(file, H5F_SCOPE_GLOBAL);
}
int Analysis::writeData(Timing timing, double dt)
{
if (timing.check(dataOutputMoments, dt) ) {
FA_Mom_Tp->write(getTemperatureParallel().data());
FA_Mom_HeatFlux->write(getHeatFlux().data());
FA_Mom_Density->write(getNumberDensity().data());
FA_Mom_Time->write(&timing);
writeMessage("Data I/O : Moments output");
}
if (timing.check(dataOutputStatistics, dt) ) {
// Ugly and error-prone
getPowerSpectrum();
Array2d pSpecX(Range(1, plasma->nfields), Range(0, Nx/2)); pSpecX(Range(1, plasma->nfields), Range(0, Nx/2)) = pSpec((int) DIR_X, Range(1, plasma->nfields), Range(0, Nx/2));
Array2d pSpecY(Range(1, plasma->nfields), Range(0, Nky)); pSpecY(Range(1, plasma->nfields), Range(0, Nky)) = pSpec((int) DIR_Y, Range(1, plasma->nfields), Range(0, Nky));
Array2d pPhaseX(Range(1, plasma->nfields), Range(0, Nx/2)); pPhaseX(Range(1, plasma->nfields), Range(0, Nx/2)) = pPhase((int) DIR_X, Range(1, plasma->nfields), Range(0, Nx/2));
Array2d pPhaseY(Range(1, plasma->nfields), Range(0, Nky)) ; pPhaseY(Range(1, plasma->nfields), Range(0, Nky)) = pPhase((int) DIR_Y, Range(1, plasma->nfields), Range(0, Nky));
FA_grow_x->write( pSpecX.data()); FA_grow_y->write( pSpecY.data()); FA_grow_t->write(&timing);
FA_freq_x->write(pPhaseX.data()); FA_freq_y->write(pPhaseY.data()); FA_freq_t->write(&timing);
// Heat Flux
Array3d heatKy; heatKy.reference(getHeatFluxKy());
FA_heatKy->write(heatKy.data());
Array3d particleKy; particleKy.reference(getParticleFluxKy());
FA_particleKy->write(particleKy.data());
ScalarValues scalarValues;
// calculate kineic Energy first, need for initial_e ! sum over sumdomains
scalarValues.timestep = timing.step;
scalarValues.time = timing.time;
getFieldEnergy(scalarValues.phiEnergy, scalarValues.ApEnergy, scalarValues.BpEnergy);
// Get scalar Values for every species
for(int s = NsGlD; s <= NsGuD; s++) {
scalarValues.particle_number[s-1] = getParticelNumber(s) ;
scalarValues.entropy [s-1] = getEntropy(s) ;
scalarValues.kinetic_energy [s-1] = getKineticEnergy(s) ;
scalarValues.particle_flux [s-1] = getTotalParticleFlux(s) ;
scalarValues.heat_flux [s-1] = getTotalHeatFlux(s) ;
}
SVTable->append(&scalarValues);
// write out to Terminal/File
std::stringstream messageStream;
messageStream << "Step : " << scalarValues.timestep << " Time " << scalarValues.time
<< " Field : (phi)" << scalarValues.phiEnergy << " (Ap)" << scalarValues.ApEnergy << " (Bp) " << scalarValues.BpEnergy << std::endl;
double charge = 0., kinetic_energy=0.;
for(int s = NsGlD; s <= NsGuD; s++) {
messageStream << plasma->species(s).name << " N :" << scalarValues.particle_number[s-1] << " Kinetic Energy : " << scalarValues.kinetic_energy[s-1] ;
messageStream << " Particle Flux :" << scalarValues.particle_flux[s-1] << " Heat Flux : " << scalarValues.heat_flux[s-1] << std::endl;
charge += plasma->species(s).q * scalarValues.particle_number[s-1];
kinetic_energy += scalarValues.kinetic_energy[s-1];
}
messageStream << std::endl << "------------------------------------------------------------------" <<
std::endl << "Total Energy " << kinetic_energy+scalarValues.phiEnergy + scalarValues.ApEnergy + scalarValues.BpEnergy << " Total Charge = " << ((plasma->species(0).n0 != 0.) ? 0. : charge) << std::endl;
parallel->print(messageStream);
}
return HELIOS_SUCCESS;
}
