1 |
! |
|
2 |
! $Header: /home/cvsroot/LMDZ4/libf/dyn3d/laplacien_rot.F,v 1.1.1.1 2004/05/19 12:53:07 lmdzadmin Exp $ |
! $Header: /home/cvsroot/LMDZ4/libf/dyn3d/laplacien_rot.F,v 1.1.1.1 |
3 |
! |
! 2004/05/19 12:53:07 lmdzadmin Exp $ |
4 |
SUBROUTINE laplacien_rot ( klevel, rotin, rotout,ghx,ghy ) |
|
5 |
c |
SUBROUTINE laplacien_rot(klevel, rotin, rotout, ghx, ghy) |
6 |
c P. Le Van |
|
7 |
c |
! P. Le Van |
8 |
c ************************************************************ |
|
9 |
c ... calcul de ( rotat x nxgrad ) du rotationnel rotin . |
! ************************************************************ |
10 |
c ************************************************************ |
! ... calcul de ( rotat x nxgrad ) du rotationnel rotin . |
11 |
c |
! ************************************************************ |
12 |
c klevel et rotin sont des arguments d'entree pour le s-prog |
|
13 |
c rotout est un argument de sortie pour le s-prog |
! klevel et rotin sont des arguments d'entree pour le s-prog |
14 |
c |
! rotout est un argument de sortie pour le s-prog |
15 |
use dimens_m |
|
16 |
use paramet_m |
USE dimens_m |
17 |
use comgeom |
USE paramet_m |
18 |
use filtreg_m, only: filtreg |
USE comgeom |
19 |
IMPLICIT NONE |
USE filtreg_m, ONLY: filtreg |
20 |
c |
use rotatf_m, only: rotatf |
21 |
|
|
22 |
c |
IMPLICIT NONE |
23 |
c .......... variables en arguments ............. |
|
24 |
c |
|
25 |
INTEGER klevel |
|
26 |
REAL rotin( ip1jm,klevel ), rotout( ip1jm,klevel ) |
! .......... variables en arguments ............. |
27 |
c |
|
28 |
c .......... variables locales ................ |
INTEGER, INTENT (IN) :: klevel |
29 |
c |
REAL rotin(iim + 1, jjm, klevel), rotout(ip1jm, klevel) |
30 |
REAL ghy(ip1jm,klevel), ghx(ip1jmp1,klevel) |
|
31 |
c ........................................................ |
! .......... variables locales ................ |
32 |
c |
|
33 |
c |
REAL ghy(ip1jm, klevel), ghx(ip1jmp1, klevel) |
34 |
CALL filtreg ( rotin , jjm, klevel, 2, 1, .FALSE., 1 ) |
! ........................................................ |
35 |
|
|
36 |
CALL nxgrad ( klevel, rotin, ghx , ghy ) |
|
37 |
CALL rotatf ( klevel, ghx , ghy , rotout ) |
CALL filtreg(rotin, direct = .true., intensive = .true.) |
38 |
c |
|
39 |
RETURN |
CALL nxgrad(klevel, rotin, ghx, ghy) |
40 |
END |
CALL rotatf(klevel, ghx, ghy, rotout) |
41 |
|
|
42 |
|
RETURN |
43 |
|
END SUBROUTINE laplacien_rot |