| 1 | ! | 
|---|
| 2 | ! $Header$ | 
|---|
| 3 | ! | 
|---|
| 4 | SUBROUTINE laplacien ( klevel, teta, divgra ) | 
|---|
| 5 |   ! | 
|---|
| 6 |   ! P. Le Van | 
|---|
| 7 |   ! | 
|---|
| 8 |   !   ************************************************************ | 
|---|
| 9 |   !    ....     calcul de  (div( grad ))   de   teta  ..... | 
|---|
| 10 |   !   ************************************************************ | 
|---|
| 11 |   ! klevel et teta  sont des arguments  d'entree pour le s-prog | 
|---|
| 12 |   !  divgra     est  un argument  de sortie pour le s-prog | 
|---|
| 13 |   ! | 
|---|
| 14 |   USE comgeom_mod_h | 
|---|
| 15 |   USE dimensions_mod, ONLY: iim, jjm, llm, ndm | 
|---|
| 16 | USE paramet_mod_h | 
|---|
| 17 | IMPLICIT NONE | 
|---|
| 18 |   ! | 
|---|
| 19 |  | 
|---|
| 20 |  | 
|---|
| 21 |  | 
|---|
| 22 |   ! | 
|---|
| 23 |   !    .........      variables  en arguments   .............. | 
|---|
| 24 |   ! | 
|---|
| 25 |   INTEGER :: klevel | 
|---|
| 26 |   REAL :: teta( ip1jmp1,klevel ), divgra( ip1jmp1,klevel ) | 
|---|
| 27 |   ! | 
|---|
| 28 |   !    ............     variables  locales      .............. | 
|---|
| 29 |   ! | 
|---|
| 30 |   REAL :: ghy(ip1jm,llm), ghx(ip1jmp1,llm) | 
|---|
| 31 |   !    ....................................................... | 
|---|
| 32 |  | 
|---|
| 33 |  | 
|---|
| 34 |   ! | 
|---|
| 35 |   CALL SCOPY ( ip1jmp1 * klevel, teta, 1, divgra, 1 ) | 
|---|
| 36 |  | 
|---|
| 37 |   CALL filtreg( divgra,  jjp1, klevel,  2, 1, .TRUE., 1 ) | 
|---|
| 38 |   CALL   grad ( klevel,divgra,   ghx , ghy              ) | 
|---|
| 39 |   CALL  divergf ( klevel, ghx , ghy  , divgra           ) | 
|---|
| 40 |  | 
|---|
| 41 |   RETURN | 
|---|
| 42 | END SUBROUTINE laplacien | 
|---|