Diffusion schemes -- grid faces


INTERFACE:

   subroutine diff_face(N,dt,cnpar,h,Bcup,Bcdw,Yup,Ydw,nuY,Lsour,Qsour,Y)
DESCRIPTION:


USES:

   use util,          only  : Dirichlet, Neumann
   use mtridiagonal
 
   IMPLICIT NONE
INPUT PARAMETERS:
 
   number of vertical layers
   integer,  intent(in)                :: N
 
   time step (s)
   REALTYPE, intent(in)                :: dt
 
   "implicitness" parameter
   REALTYPE, intent(in)                :: cnpar
 
   layer thickness (m)
   REALTYPE, intent(in)                :: h(0:N)
 
   type of upper BC
   integer,  intent(in)                :: Bcup
 
   type of lower BC
   integer,  intent(in)                :: Bcdw
 
   value of upper BC
   REALTYPE, intent(in)                :: Yup
 
   value of lower BC
   REALTYPE, intent(in)                :: Ydw
 
   diffusivity of Y
   REALTYPE, intent(in)                :: nuY(0:N)
 
   linear source term
   (treated implicitly)
   REALTYPE, intent(in)                :: Lsour(0:N)
 
   constant source term
   (treated explicitly)
   REALTYPE, intent(in)                :: Qsour(0:N)
INPUT/OUTPUT PARAMETERS:
   REALTYPE                            :: Y(0:N)
REVISION HISTORY:
   Original author(s): Lars Umlauf



Karsten Bolding 2012-12-28