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