breaking.h

来自「波浪数值模拟」· C头文件 代码 · 共 62 行

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



/* breaking.h   Falk Feddersen */       

#ifndef BREAKING_H
#define BREAKING_H

#include "bdefs.h"
#include "depth.h"

typedef enum {
  ZELT_BR, FUNWAVE_BR, FALK_BR
} breaking_model_t;

typedef struct {
  switch_t breaking_flag;     /* on or off */
  breaking_model_t  brmodel;  /* the type of breaking model: for default lets make it Zelt Coastal Eng. (1991) */
  double delta_b ;            /* the \delta_b term in Kennedy et al & Funwave manual = 1.2 */
  double c_eta_tI;            /* the coefficient for \eta_t^{(I)}: set to 0.65*(gh)^(1/2) as default */
  double c_eta_tF;            /* the coefficient for \eta_t^{(F)}: set to 0.15*(gh)^(1/2) as default */
  double cT;                  /* the coefficient for T^* = cT * sqrt(g/H) - default cT=5.0 */
  double Tstar;               /* T^* : set to 5 (h/g)^(1/2)                      */
  double frequency;
  double wavelength;
  double omega;    // 2*pi*f  f = wave frequency
} breaking_info_t;


void  init_breaking(breaking_info_t *BR, field2D_depth_var *DD);
void  deallocate_breaking();
void  breaking_report();                 /* prints out breaking information */
field2D*  breaking_eddy_viscosity();     /* This is for output purposes - returns a pointer to the eddy viscosity */

void  rhs_uvmom_breaking_add(field2D *UTBR, field2D *VTBR, const field2D *U, const field2D *V, const field2D *ETAT);





#endif /* BREAKING_H */