N_gwflow.h

Go to the documentation of this file.

/*****************************************************************************
*
* MODULE:       Grass PDE Numerical Library
* AUTHOR(S):    Soeren Gebbert, Berlin (GER) Dec 2006
*               soerengebbert <at> gmx <dot> de
*               
* PURPOSE:      groundwater flow in porous media 
*               part of the gpde library
*
* COPYRIGHT:    (C) 2000 by the GRASS Development Team
*
*               This program is free software under the GNU General Public
*               License (>=v2). Read the file COPYING that comes with GRASS
*               for details.
*
*****************************************************************************/

#ifndef _N_GWFLOW_H_
#define _N_GWFLOW_H_
#include "N_pde.h"

#define N_GW_CONFINED 0         /*confined groundwater */
#define N_GW_UNCONFINED 1       /*unconfined groundwater */

#define N_GW_DRY_CELL 0         /*a dry cell */
#define N_GW_SURFACE_CELL 1     /*a surface cell */
#define N_GW_NORMAL_CELL 2      /*a normal wet cell */

typedef struct
{
    N_array_3d *phead;          
    N_array_3d *phead_start;    
    N_array_3d *hc_x;           
    N_array_3d *hc_y;           
    N_array_3d *hc_z;           
    N_array_3d *q;              
    N_array_2d *r;              
    N_array_3d *s;              
    N_array_3d *nf;             
    /*river */
    N_array_3d *river_leak;     
    N_array_3d *river_head;     
    N_array_3d *river_bed;      
    /*drainage */
    N_array_3d *drain_leak;     
    N_array_3d *drain_bed;      
    N_array_3d *status;         
    N_array_3d *drycells;       
    double dt;                  
} N_gwflow_data3d;

typedef struct
{
    N_array_2d *phead;          
    N_array_2d *phead_start;    
    N_array_2d *hc_x;           
    N_array_2d *hc_y;           
    N_array_2d *q;              
    N_array_2d *r;              
    N_array_2d *s;              
    N_array_2d *nf;             
    /*river */
    N_array_2d *river_leak;     
    N_array_2d *river_head;     
    N_array_2d *river_bed;      
    /*drainage */
    N_array_2d *drain_leak;     
    N_array_2d *drain_bed;      
    N_array_2d *top;            
    N_array_2d *bottom;         
    N_array_2d *status;         
    N_array_2d *drycells;       
    double dt;                  
    int gwtype;                 
} N_gwflow_data2d;

extern N_data_star *N_callback_gwflow_3d(void *gwdata, N_geom_data * geom,
                                         int col, int row, int depth);
extern N_data_star *N_callback_gwflow_2d(void *gwdata, N_geom_data * geom,
                                         int col, int row);
extern N_gwflow_data3d *N_alloc_gwflow_data3d(int cols, int rows, int depths,
                                              int river, int drain);
extern N_gwflow_data2d *N_alloc_gwflow_data2d(int cols, int rows, int river,
                                              int drain);
extern void N_free_gwflow_data3d(N_gwflow_data3d * data);

extern void N_free_gwflow_data2d(N_gwflow_data2d * data);
#endif