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