Diffusion schemes -- grid faces (Source File: diff_face.F90)

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