bousinessq_constants.h

波浪数值模拟

/*
 * Copyright (c) 2001-2005 Falk Feddersen
 *
 * This program is free software; you can redistribute it and/or modify
 * it under the terms of the GNU General Public License as published by
 * the Free Software Foundation; either version 2 of the License, or
 * (at your option) any later version.
 *
 * This program is distributed in the hope that it will be useful,
 * but WITHOUT ANY WARRANTY; without even the implied warranty of
 * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
 * GNU General Public License for more details.
 *
 * You should have received a copy of the GNU General Public License
 * along with this program; if not, write to the Free Software
 * Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA  02111-1307  USA
 *
 */

/* ---  bousinessq_constants.h  

This is where the constants for the particular bousinessq wave model are defined.
It assumes that somewhere else  the file #include "bdefs.h" has been included.
So far only the files:
                         bousinessq_dynamics.c   and tridiagonal.c  need this file

*/

#ifdef BOUSINESSQ_CONSTANTS_H
#define BOUSINESSQ_CONSTANTS_H


#ifdef LINEAR
const double beta = -0.531;         /*  betea =  za/h */
const double a1 = 0; //0.5*beta*beta - (1.0/6.0);   /* a1 and a2 are bousinessq parameters in the d(eta)/dt = ... equation */
const double a2 = 0; //beta + 0.5;
const double b1 = 0.0; //0.5* beta*beta;        /* b1 and b2 are bousinessq parameters in the d(u,v)/dt = ... equation */
const double b2 = 0.0; //beta;
const double gamma = 0.0;
#endif

#ifdef NSWE
const double beta = -0.531;         /*  betea =  za/h */
const double a1 = 0; //0.5*beta*beta - (1.0/6.0);   /* a1 and a2 are bousinessq parameters in the d(eta)/dt = ... equation */
const double a2 = 0; //beta + 0.5;
const double b1 = 0.0; //0.5* beta*beta;        /* b1 and b2 are bousinessq parameters in the d(u,v)/dt = ... equation */
const double b2 = 0.0; //beta;
const double gamma = 0.0;
#endif

#ifdef PEREGRINE   /* peregrine depth-averaged dynamics */
const double beta = -0.5;
const double a1 = 0.0; 
const double a2 = 0.0; 
const double b1 = 1.0/6.0; /* not sure why b1 is what it is */
const double b2 = -0.5;    /* or b2 but these are the peregrine JFM (1967) values */
const double gamma = 0.0;
#endif  /* PEREGRINE */

#ifdef WEI_KIRBY
const double beta = -0.531;         /*  betea =  za/h */
const double a1 =  -2.568616666666665e-02;  /* a1 = 0.5*beta_h *beta_h - (1.0/6.0);  a1 and a2 are bous parameters in d(eta)/dt =  */
const double a2 = -0.031;                     /* a2 = beta + 0.5;  */
const double b1 =   1.409805000000000e-01;    /* b1 = 0.5* beta_h*beta_h;   b1 and b2 are bousinessq parameters in the d(u,v)/dt = */
const double b2 =  -0.531;
const double gamma = 1.0;
#endif

#ifdef NWOGU
const double beta = -0.531;         /*  betea =  za/h */
const double a1 =  -2.568616666666665e-02;    /* a1 = 0.5*beta_h *beta_h - (1.0/6.0);  a1 and a2 are bous parameters in d(eta)/dt = ... equation */
const double a2 = -0.031;                     /* a2 = beta + 0.5;  */
const double b1 =   1.409805000000000e-01;    /* b1 = 0.5* beta_h*beta_h;   b1 and b2 are bousinessq parameters in the d(u,v)/dt = ... equation */
const double b2 =  -0.531;
const double gamma = 0.0;
#endif


#endif /* BOUSINESSQ_CONSTANTS_H */