shallowwater1d.h

Go to the documentation of this file.

#include <basic.h>
#include <math.h>
#include <matops.h>

#define _SHALLOW_WATER_1D_ "shallow-water-1d"

/* define ndims and nvars for this model */
#undef _MODEL_NDIMS_
#undef _MODEL_NVARS_

#define _MODEL_NDIMS_ 1

#define _MODEL_NVARS_ 2

/* choices for upwinding schemes */
#define _ROE_     "roe"

#define _LLF_     "llf-char"

/* grid direction */
#undef _XDIR_

#define _XDIR_ 0

#define _ShallowWater1DGetFlowVar_(u,h,v) \
  { \
    h = u[0]; \
    v = u[1] / h; \
  }

#define _ShallowWater1DSetFlux_(f,h,v,g) \
  { \
    f[0] = (h) * (v); \
    f[1] = (h) * (v) * (v) + 0.5 * (g) * (h) * (h); \
  }

#define _ShallowWater1DRoeAverage_(uavg,uL,uR,p) \
  { \
    double h ,v; \
    double hL,vL; \
    double hR,vR; \
    _ShallowWater1DGetFlowVar_(uL,hL,vL); \
    _ShallowWater1DGetFlowVar_(uR,hR,vR); \
    h  = 0.5 * (hL + hR); \
    v  = (sqrt(hL)*vL + sqrt(hR)*vR) / (sqrt(hL) + sqrt(hR)); \
    uavg[0] = h; \
    uavg[1] = h*v; \
  }

#define _ShallowWater1DEigenvalues_(u,D,p,dir) \
  { \
    double h,v,c; \
    _ShallowWater1DGetFlowVar_(u,h,v); \
    c = sqrt(p->g*h); \
    D[0*_MODEL_NVARS_+0] = (v+c);  D[0*_MODEL_NVARS_+1] = 0;     \
    D[1*_MODEL_NVARS_+0] = 0;      D[1*_MODEL_NVARS_+1] = (v-c); \
  }

#define _ShallowWater1DLeftEigenvectors_(u,L,p,dir) \
  { \
    double R[_MODEL_NVARS_*_MODEL_NVARS_]; \
    _ShallowWater1DRightEigenvectors_(u,R,p,dir); \
    _MatrixInvert_(R,L,_MODEL_NVARS_); \
  }

#define _ShallowWater1DRightEigenvectors_(u,R,p,dir) \
  { \
    double h,v,c; \
    _ShallowWater1DGetFlowVar_(u,h,v); \
    c    = sqrt(p->g*h); \
    R[0*_MODEL_NVARS_+0] = 1.0; \
    R[1*_MODEL_NVARS_+0] = v+c; \
    R[0*_MODEL_NVARS_+1] = 1.0; \
    R[1*_MODEL_NVARS_+1] = v-c; \
  }

typedef struct shallowwater1d_parameters {
  double  g;         
  int     bt_type,   
          topo_flag; 
  double  *b;        
  char    upw_choice[_MAX_STRING_SIZE_]; 
  int (*SourceUpwind)(double*,double*,double*,double*,int,void*,double);
} ShallowWater1D;

int    ShallowWater1DInitialize (void*,void*);
int    ShallowWater1DCleanup    (void*);