/Users/huttone/Devel/sedflux-new/sedflux/trunk/ew/inflow/inflow.c File Reference

#include <stdlib.h>
#include <stdio.h>
#include <math.h>
#include <limits.h>
#include "inflow.h"
#include <sed/sed_sedflux.h>
#include <utils/utils.h>

Go to the source code of this file.

Functions

gboolean inflow (double day, double x[], double slope[], double width[], int n_nodes, double dx, double x_dep, double river_width, double river_vel, double river_depth, double dc, double *gzF0, double *grain_dia, double *lambda, double *rho_sed, double *rho_grain, int n_grains, double rho_rw, double rho_flow, Inflow_const_st *c, double **deposited, double **eroded, FILE *fpout)


Function Documentation

gboolean inflow ( double  day,
double  x[],
double  slope[],
double  width[],
int  n_nodes,
double  dx,
double  x_dep,
double  river_width,
double  river_vel,
double  river_depth,
double  dc,
double *  gzF0,
double *  grain_dia,
double *  lambda,
double *  rho_sed,
double *  rho_grain,
int  n_grains,
double  rho_rw,
double  rho_flow,
Inflow_const_st c,
double **  deposited,
double **  eroded,
FILE *  fpout 
)

inflow turbidity current model.

A steady state turbidity current model based on the Mulder model.

Mulder, T., Savoye, B., and Syvitski, J.P.M., 1997. Numerical modelling of a mid-sized gravity flow: the 1979 Nice turbidity current (dynamics, processes, sediment budget and seafloor impact). Sedimentology, v. 44, pp. 305-326.

Parameters:
day duration of the flow (s). Typically this is one day.
x x coordinate of each node (m).
slope sea-floor slope (rads). a negative slope means downward dipping.
width width of the channel (m).
n_nodes number of nodes in the domain.
dx distance between nodes (m).
x_dep distance (m) from the river mouth before which no deposition is allowed.
river_width width of the river mouth (m).
river_vel velocity of the river at the river mouth (m/s).
river_depth depth of the river at the river mouth (m).
dc water discharge at the river mouth (m^3/s).
gzF0 fraction of each grain type in the flow.
grain_dia grain diameter of each grain type in the flow (m).
lambda removal rate for each grain type (1/s).
rho_sed bulk density of each grain type (kg/m^3).
rho_grain grain density of each grain type (kg/m^3).
n_grains number of grain types.
rho_rw density of the river water (kg/m^3). This is the river water without any sediment.
rho_flow density of the flow (kg/m^3). This is the river water plus the sediment.
c some constants that are required by the turbidity current model.
deposited 2d array for storing the deposition rates for each node. The fast dimension is over grid node number and the slow dimension is over grain type.
eroded 2d array for storing the erosion rates for each node. The fast dimension is over grid node number and the slow dimension is over grain type.
fpout The output file
Returns:
TRUE on success and FALSE if a problem is found.

Definition at line 75 of file inflow.c.

References Inflow_const_st::c_drag, Inflow_phe_query_st::dx, Inflow_const_st::e_a, Inflow_const_st::e_b, eh_dbl_array_mean_weighted(), eh_free, eh_free_2, eh_new, eh_new_2, Inflow_phe_query_st::erode_depth, Inflow_const_st::get_phe, Inflow_const_st::get_phe_data, Inflow_const_st::mu_water, Inflow_phe_query_st::phe, Inflow_const_st::rho_sea_water, S_DAYS_PER_SECOND, sed_gravity(), Inflow_const_st::sua, Inflow_const_st::sub, Inflow_const_st::tan_phi, and Inflow_phe_query_st::x.

Referenced by inflow_wrapper(), sakura(), and sakura_array_set_bc().


Generated on Fri Jan 4 18:04:17 2008 for sedflux by  doxygen 1.5.2