!!ATTENTION HYPOTHESE de FLUX COMPENSATOIRE DESCENDANT ET DONC comme schema amont on va chercher trac au dessus!!!!! !! tentative de prise en compte d'un flux compensatoire montant !!!! si on voulait le prendre en compte on
! Reecriture du schéma explicite avec les notations du schéma implicite ! Calcul de s1
| Type | Intent | Optional | Attributes | Name | ||
|---|---|---|---|---|---|---|
| integer, | intent(in) | :: | ngrid | |||
| integer, | intent(in) | :: | nlay | |||
| real, | intent(in) | :: | ptimestep | |||
| integer, | intent(in), | dimension(ngrid) | :: | lmax | ||
| real, | intent(in), | dimension(ngrid,nlay) | :: | eup | ||
| real, | intent(in), | dimension(ngrid,nlay) | :: | dup | ||
| real, | intent(in), | dimension(ngrid,nlay) | :: | edn | ||
| real, | intent(in), | dimension(ngrid,nlay) | :: | ddn | ||
| real, | intent(in), | dimension(ngrid,nlay) | :: | masse | ||
| real, | intent(in), | dimension(ngrid,nlay) | :: | trac | ||
| real, | intent(out), | dimension(ngrid,nlay) | :: | dtrac |