leachdly.c File Reference

#include "swconst.h"

Include dependency graph for leachdly.c:

Go to the source code of this file.

Functions
void	leachdly (double wfluxout[], int numlyrs, double nitrate[], double critflow, double frlechd[], double stream[], double basef, double stormf, double *inorglch)

Function Documentation

void leachdly	(	double	wfluxout[],
		int	numlyrs,
		double	nitrate[],
		double	critflow,
		double	frlechd[],
		double	stream[],
		double	basef,
		double	stormf,
		double *	inorglch
	)

    {

      int lyr, nxt;
      int iel = 0;  /* nitrogen */
      double amtlea[MAXLYR];
      double strm_N;
      double base_N;
      double linten;

      strm_N = 0.0;
      base_N = 0.0;
      *inorglch = 0.0;

      for(lyr=0; lyr < numlyrs; lyr++) {
        amtlea[lyr] = 0.0;
        nxt = lyr + 1;

        /* wfluxout > 0. indicates a saturated water flow out of layer lyr */

        if ((wfluxout[lyr] >  0.0) && (nitrate[lyr] > 1.0E-10)) {
          linten = min(1.0 - (critflow - wfluxout[lyr]) / critflow, 1.0);
          linten = max(linten, 0.0);
          amtlea[lyr] = (frlechd[iel]) * nitrate[lyr] * linten;

          /* If you are at the bottom layer, compute storm flow. */
          if (lyr == numlyrs-1) {
            strm_N = amtlea[lyr] * stormf;
          }
          if (amtlea[lyr] < nitrate[lyr]) {
            nitrate[lyr] = nitrate[lyr] - amtlea[lyr];
            nitrate[nxt] = nitrate[nxt] + (amtlea[lyr] - strm_N);
          }
        }
      }

      /* Compute base flow and mineral stream flows. */
      /* Check for underflow. -mdh */

      if (nitrate[nxt] >= 1.0E-10) {
        base_N = nitrate[nxt] * basef;
        nitrate[nxt] = nitrate[nxt] - base_N;
      }

      /* Note: stream flow indices differ from mineral element */ 
      /*       indices by 1 (eg  stream(2) is stream flow for nitrogen). */

      stream[iel+1] = stream[iel+1] + (strm_N + base_N);
      *inorglch = strm_N + base_N;

      return;
    }