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