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
   $Log: diff_face.F90,v $
   Revision 1.3  2005-11-03 20:56:55  hb
   Source term linearisation now fully explicit again, reversion to old method
   Revision 1.2  2005/09/16 13:54:02  lars
   added missing IMPLICIT NONE
   Revision 1.1  2005/06/27 10:54:33  kbk
   new files needed