forked from fmihpc/vlasiator
-
Notifications
You must be signed in to change notification settings - Fork 0
/
iowrite.h
99 lines (84 loc) · 4.44 KB
/
iowrite.h
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
/*
* This file is part of Vlasiator.
* Copyright 2010-2016 Finnish Meteorological Institute
*
* For details of usage, see the COPYING file and read the "Rules of the Road"
* at http://www.physics.helsinki.fi/vlasiator/
*
* This program is free software; you can redistribute it and/or modify
* it under the terms of the GNU General Public License as published by
* the Free Software Foundation; either version 2 of the License, or
* (at your option) any later version.
* This program is distributed in the hope that it will be useful,
* but WITHOUT ANY WARRANTY; without even the implied warranty of
* MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
* GNU General Public License for more details.
* You should have received a copy of the GNU General Public License along
* with this program; if not, write to the Free Software Foundation, Inc.,
* 51 Franklin Street, Fifth Floor, Boston, MA 02110-1301 USA.
*/
#ifndef IOWRITE_H
#define IOWRITE_H
#include "mpi.h"
#include <dccrg.hpp>
#include <dccrg_cartesian_geometry.hpp>
#include <string>
#include <vector>
#include <vlsv_writer.h>
#include "spatial_cell.hpp"
#include "datareduction/datareducer.h"
/*!
\brief Write out system into a vlsv file
\param mpiGrid The DCCRG grid with spatial cells
\param dataReducer Contains datareductionoperators that are used to compute data that is added into file
\param index Index to call the correct member of the various parameter vectors
\param writeGhosts Write ghost zones
*/
bool writeGrid(dccrg::Dccrg<SpatialCell,dccrg::Cartesian_Geometry>& mpiGrid,
FsGrid< std::array<Real, fsgrids::bfield::N_BFIELD>, 2>& perBGrid,
FsGrid< std::array<Real, fsgrids::efield::N_EFIELD>, 2>& EGrid,
FsGrid< std::array<Real, fsgrids::ehall::N_EHALL>, 2>& EHallGrid,
FsGrid< std::array<Real, fsgrids::egradpe::N_EGRADPE>, 2>& EGradPeGrid,
FsGrid< std::array<Real, fsgrids::moments::N_MOMENTS>, 2>& momentsGrid,
FsGrid< std::array<Real, fsgrids::dperb::N_DPERB>, 2>& dPerBGrid,
FsGrid< std::array<Real, fsgrids::dmoments::N_DMOMENTS>, 2>& dMomentsGrid,
FsGrid< std::array<Real, fsgrids::bgbfield::N_BGB>, 2>& BgBGrid,
FsGrid< std::array<Real, fsgrids::volfields::N_VOL>, 2>& volGrid,
FsGrid< fsgrids::technical, 2>& technicalGrid,
DataReducer* dataReducer,
const uint& index,
const bool writeGhosts = true
);
/*!
\brief Write out a restart of the simulation into a vlsv file. All block data in remote cells will be reset.
\param mpiGrid The DCCRG grid with spatial cells
\param dataReducer Contains datareductionoperators that are used to compute data that is added into file
\param name File name prefix, file will be called "name.index.vlsv"
\param fileIndex File index, file will be called "name.index.vlsv"
*/
bool writeRestart(dccrg::Dccrg<SpatialCell,dccrg::Cartesian_Geometry>& mpiGrid,
FsGrid< std::array<Real, fsgrids::bfield::N_BFIELD>, 2>& perBGrid,
FsGrid< std::array<Real, fsgrids::efield::N_EFIELD>, 2>& EGrid,
FsGrid< std::array<Real, fsgrids::ehall::N_EHALL>, 2>& EHallGrid,
FsGrid< std::array<Real, fsgrids::egradpe::N_EGRADPE>, 2>& EGradPeGrid,
FsGrid< std::array<Real, fsgrids::moments::N_MOMENTS>, 2>& momentsGrid,
FsGrid< std::array<Real, fsgrids::dperb::N_DPERB>, 2>& dPerBGrid,
FsGrid< std::array<Real, fsgrids::dmoments::N_DMOMENTS>, 2>& dMomentsGrid,
FsGrid< std::array<Real, fsgrids::bgbfield::N_BGB>, 2>& BgBGrid,
FsGrid< std::array<Real, fsgrids::volfields::N_VOL>, 2>& volGrid,
FsGrid< fsgrids::technical, 2>& technicalGrid,
DataReducer& dataReducer,
const std::string& name,
const uint& fileIndex,
const int& stripe);
/*!
\brief Write out simulation diagnostics into diagnostic.txt
@param mpiGrid The DCCRG grid with spatial cells
@param dataReducer Contains datareductionoperators that are used to compute diagnostic data
*/
bool writeDiagnostic(const dccrg::Dccrg<SpatialCell,dccrg::Cartesian_Geometry>& mpiGrid,DataReducer& dataReducer);
bool writeVelocitySpace(dccrg::Dccrg<SpatialCell,dccrg::Cartesian_Geometry>& mpiGrid,
vlsv::Writer& vlsvWriter,int index,const std::vector<uint64_t>& cells);
bool writeVelocityDistributionData(vlsv::Writer& vlsvWriter,dccrg::Dccrg<SpatialCell,dccrg::Cartesian_Geometry>& mpiGrid,
const std::vector<uint64_t>& cells,MPI_Comm comm);
#endif