ShallowWater2DWriteTopography.c

Go to the documentation of this file.

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

int ShallowWater2DWriteTopography(  void*   s,  
                                    void*   m,  
                                    double  a_t  )
{
  HyPar           *solver = (HyPar*)          s;
  MPIVariables    *mpi    = (MPIVariables*)   m;
  ShallowWater2D  *params = (ShallowWater2D*) solver->physics;
  _DECLARE_IERR_;

  if (params->topo_flag) {

    char fname_root[_MAX_STRING_SIZE_] = "topography";
    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);
      strcat(fname_root, "_");
    }

    WriteArray(  solver->ndims,
                 1,
                 solver->dim_global,
                 solver->dim_local,
                 solver->ghosts,
                 solver->x,
                 params->b,
                 solver,
                 mpi,
                 fname_root );

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

  }

  return(0);
}