VlasovWriteEFieldAndPotential.c

Go to the documentation of this file.

#include <stdlib.h>
#include <string.h>
#include <basic.h>
#include <common.h>
#include <arrayfunctions.h>
#include <mathfunctions.h>
#include <io.h>
#include <mpivars.h>
#include <hypar.h>
#include <physicalmodels/vlasov.h>

int VlasovWriteSpatialField( void*, void*, double*, char* );
int VlasovPlotSpatialField( void*, void*, double*, double, char* );

int VlasovWriteEFieldAndPotential(  void*   s,    
                                    void*   m,    
                                    double  a_t    )
{
  HyPar         *solver = (HyPar*)        s;
  MPIVariables  *mpi    = (MPIVariables*) m;
  Vlasov        *param  = (Vlasov*)       solver->physics;

  {
    char fname_root[_MAX_STRING_SIZE_] = "efield";
    VlasovWriteSpatialField( solver, mpi, param->e_field, fname_root );
    if (!strcmp(solver->plot_solution,"yes")) {
      VlasovPlotSpatialField( solver, mpi, param->e_field, a_t, fname_root );
    }
  }

  if (param->self_consistent_electric_field) {
    char fname_root[_MAX_STRING_SIZE_] = "potential";
    VlasovWriteSpatialField( solver, mpi, param->potential, fname_root );
    if (!strcmp(solver->plot_solution,"yes")) {
      VlasovPlotSpatialField( solver, mpi, param->potential, a_t, fname_root );
    }
  }

  return 0;
}