tracer.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.
 *
 * For a copy of the GNU General Public License see LICENSE or go to
 *     http://www.gnu.org/licenses/gpl.html
 *
 */

/* ---  tracer.h:  Header file for tracer data structures and routinues   */

#ifndef TRACER_H
#define TRACER_H

#include "depth.h"
#include "funwaveC_timestep.h"

typedef enum {
  TRACER_POINT_SOURCE,  TRACER_INITIAL_CONDITION
} tracer_source_t;

typedef struct {
  switch_t  tracer_flag;
  tracer_source_t  source_flag;    /* source type */
  char filename[80];
  int  i_source;               /* i location of the source */
  int  j_source;               /* j location of the source */
  double  point_source_flux;   /* Units of  stuff/(m^2 s) where stuff can = Liters,gallons, mol., etc. */
  double kappa;                /* tracer diffusivity (m^2/s) */ 
  double half_life;            /* tracer half life  (sec) */
  double source_start_time;    /* start time of source in seconds */
  double source_end_time;      /* end time of source in seconds */
} tracerinfo;


extern switch_t tracer_flag;

void init_passive_tracer_variables(int nx,int ny, tracerinfo *T);
void deallocate_passive_tracer_variables();
void passive_tracer_report(tracerinfo *T);
void calc_passive_tracer_rhs(field2D *QT, const field2D *UU, const field2D *VV, const field2D *QQ, const field2D_depth_var *Depth_Var);
void update_passive_tracerAB3(field2D *QQ, field2D_time_derivative_var *DTVARS, const field2D_depth_var *Depth_Var, double dt);
void update_passive_tracerAB2(field2D *QQ, const field2D *QT1, const field2D *QT2, const field2D_depth_var *Depth_Var, double dt);
void update_passive_tracer_euler(field2D *QQ, const field2D *QQT, const field2D_depth_var *Depth_Var, double dt);


#endif    /* TRACER_H */