LinearADRWriteAdvField.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 <plotfuncs.h>
#include <mpivars.h>
#include <hypar.h>
#include <physicalmodels/linearadr.h>

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

  if (param->constant_advection == 0) {

    char fname_root[_MAX_STRING_SIZE_] = "advection_field";
    if (solver->nsims > 1) {
      char index[_MAX_STRING_SIZE_];
      GetStringFromInteger(solver->my_idx, index, (int)log10(solver->nsims)+1);
      strcat(fname_root, "_");
      strcat(fname_root, index);
    }

    int adv_nvar = solver->ndims * solver->nvars;
    WriteArray(  solver->ndims,
                 adv_nvar,
                 solver->dim_global,
                 solver->dim_local,
                 solver->ghosts,
                 solver->x,
                 param->a,
                 solver,mpi,
                 fname_root );

    if (!strcmp(solver->plot_solution, "yes")) {
      PlotArray( solver->ndims,
                 adv_nvar,
                 solver->dim_global,
                 solver->dim_local,
                 solver->ghosts,
                 solver->x,
                 param->a,
                 a_t,
                 solver,mpi,
                 fname_root );
    }

  }

  return(0);
}