bdefs.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
 *
 */


/* 
   bdefs.h - definitions header file for bousinessq wave model

   Falk Feddersen 

*/

#ifndef BDEFS_H  /* BDEFS_H */
#define BDEFS_H


#include <stdio.h>
#include <glib.h>

#include "field.h"
#include "funwaveC_config.h"


#define AB3  /* now AB3 time stepping is the default */


/* pick one of the 5 dynamics choices */

typedef enum {
  LINEAR_DYNAMICS, NSWE_DYNAMICS, PEREGRINE_DYNAMICS, NWOGU_DYNAMICS, WEI_KIRBY_DYNAMICS 
} dynamics_t;

//#define NSWE
//#define WEI_KIRBY
//#define NWOGU 
//#define PEREGRINE
//#define LINEAR



/* Check to see which boundary conditions are defined.  The default is SLIP */

#ifndef NOSLIP
#define SLIP
#endif


#define SQR(a)  ((a)*(a))

#define JPM(j)  (((j)==ny-1) ? 0 : (j)+1)
#define JMM(j)  (((j)==0) ? ny-1 : (j)-1)

//#define MAX(a,b) ( (a)>(b) ) ? (a) : (b) )
//#define MIN(a,b) ( (a)<(b) ) ? (a) : (b) )

typedef enum {
  ON, OFF
} switch_t;

/*
typedef struct {
  int nx,ny,nq;
  int pnx,p2nx,psize;
  double dx,dy,idx,idy;
  double g;
  double mu, nu_bi, nu_newt;
  double c_d;
} Param;
*/

void funwaveC_error_message(char *);
void funwaveC_error_message_file_line(char *, char *, int );


/* External Global Variable Declarations */

extern int nx;
extern int ny;

extern double iny;

extern double idx;
extern double idy;
extern double idx2;
extern double idy2;


extern double dx;
extern double dy;

extern const double f;
extern const double gravity;

extern double global_time;   // this is the global time index


extern double c_d;         // Quadratic drag coefficient (non-d)
//extern double nu_bi;      // this is bi-harmonic  m^4/s  
//extern double nu_newt;    // Newtonian friction term m^2/s

extern const double pi;

extern field2D *UU;
extern field2D *VV;
extern field2D *ETA;

extern field2D *FXX, *FYY;

extern field2D *TAUBX, *TAUBY;

extern field2D *H;    /* this is the depth field2D */

//extern double  *HBAR;

/* ---  PASSIVE TRACER VARIABLES --- */

//extern switch_t tracer_flag;

extern field2D *QQ;
extern field2D *QQT1, *QQT2, *QQT3;
extern field2D *QSRC;

/* ---- Floats variables ----- */



//extern field2D *DISS_U,  *DISS_V;


#ifdef NOSLIP
extern double HEXT0;
extern double HEXTN;
#endif

#endif  /* BDEFS_H */