/Users/huttone/Devel/sedflux-new/sedflux/trunk/ew/hydrotrend/hydrosetparams.c

Go to the documentation of this file.
00001 /*
00002  *  HydroSetParams.c
00003  *
00004  *  Author1:    S.D. Peckman (december 2001)
00005  *  Author2:    A.J. Kettner (september 2002)
00006  *
00007  *
00008  *  In hydroreadinput.c, units are converted as follows:
00009  *  Multiply totalarea by 1e6.  (km^2 -> m^2)
00010  *  Multiply basinlength by 1000.  (km -> m)
00011  *  Divide Ko by 1000.  (mm/day -> m/day)
00012  *  Divide Lapse rate by 1000.  (mm/day -> m/day)
00013  *  Some others are done by this routine.
00014  *
00015  */
00016 
00017 #include "hydroparams.h"
00018 #include "hydroclimate.h"
00019 
00020 /*---------------------------
00021  *  Start of HydroSetParams
00022  *---------------------------*/
00023 int hydrosetparams() {
00024 
00025 int ii, err;
00026 
00027 /*----------------------------------------
00028  *  Hardwire parameters for complete run
00029  *----------------------------------------*/
00030 err = 0;
00031 webflag = 0;            /*0 = weboption is off. In HydroTrend.c, this option can be turned on */
00032 
00033 /*---------------------------------------
00034  *  Hardwired Parameters for all Epochs
00035  *---------------------------------------*/
00036 for (ii=0; ii<maxepoch; ii++) {
00037     rhowater[ii]     = rho_water;    /* 1000 */
00038     rhosed[ii]       = rho_sed;      /* 2670 */
00039     alphac[ii]       = 0.98;    /* saturation excess coeff */
00040     betac[ii]        = 1.00;    /* saturation excess exponent */
00041     alphag[ii]       = -0.0001;    /* groundwater precip offset (m/day) */
00042     betag[ii]        = 0.85;    /* groundwater precip slope */
00043     alphagwe[ii]     = 0.0020;     /* groundwater evap coeff (m/day)*/
00044     betagwe[ii]      = 1.0;     /* groundwater evap exponent */
00045     pmax[ii]         = 0.400;   /* precip need to reach max cond. (m/day) */
00046     Meltrate[ii]     = 0.003;   /* (m/degC) */
00047     percentgw[ii]    = 0.15;    /* percent of nival&ice as groundwater */
00048     pcr[ii]          = 0.010;     /* crit. precip for infilt. excess (m/day) */
00049         bethaexpo                = 1.38;        /* volume-surface area exponent glaciers */
00050         bethaglacier     = 31.11716; /* volume-surface area multiplier glaciers */
00051     /* Meltrate[ii]  = 0.002;   (m/degC) */
00052     /* percentgw[ii] = 0.06;    percent of nival&ice as groundwater */
00053     /* pcr[ii]        = 2.0;    crit. precip for infilt. excess (mm/day) */
00054 
00055     /* widcof[ii]        = ???;   */
00056     /* widpow[ii]   = 0.5;   */
00057     /* velcof[ii]   = ???;   */
00058     /* velpow[ii]   = 0.1;   */
00059 }
00060 
00061 return(err);
00062 }  /* end of HydroSetParams */

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