ShallowWater1DComputeCFL.c File Reference

Contains the function to compute maximum CFL over the domain for the 1D shallow water equations physical model. More...

#include <stdlib.h>
#include <basic.h>
#include <math.h>
#include <arrayfunctions.h>
#include <mathfunctions.h>
#include <physicalmodels/shallowwater1d.h>
#include <hypar.h>

Go to the source code of this file.

Functions
double	ShallowWater1DComputeCFL (void *s, void *m, double dt, double t)

Detailed Description

Contains the function to compute maximum CFL over the domain for the 1D shallow water equations physical model.

Author

Debojyoti Ghosh

Definition in file ShallowWater1DComputeCFL.c.

Function Documentation

ShallowWater1DComputeCFL()

double ShallowWater1DComputeCFL	(	void *	s,
		void *	m,
		double	dt,
		double	t
	)

Computes the maximum CFL number over the domain. Note that the CFL is computed over the local domain on this processor only.

Parameters

s	Solver object of type HyPar
m	MPI object of type MPIVariables
dt	Time step size for which to compute the CFL
t	Time

Definition at line 17 of file ShallowWater1DComputeCFL.c.

{
  HyPar           *solver = (HyPar*)   s;
  ShallowWater1D  *param  = (ShallowWater1D*) solver->physics;

  int *dim    = solver->dim_local;
  int ghosts  = solver->ghosts;
  int ndims   = solver->ndims;
  int index[ndims];
  double *u   = solver->u;

  double max_cfl = 0;
  int done = 0; _ArraySetValue_(index,ndims,0);
  while (!done) {
    int p; _ArrayIndex1D_(ndims,dim,index,ghosts,p);
    double h, v, c, dxinv, local_cfl;
    _ShallowWater1DGetFlowVar_((u+_MODEL_NVARS_*p),h,v);
    _GetCoordinate_(0,index[0],dim,ghosts,solver->dxinv,dxinv); /* 1/dx */
    c = sqrt(param->g*h); /* speed of gravity waves */
    local_cfl = (absolute(v)+c)*dt*dxinv; /* local cfl for this grid point */
    if (local_cfl > max_cfl) max_cfl = local_cfl;
    _ArrayIncrementIndex_(ndims,dim,index,done);
  }

  return(max_cfl);
}