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_v_m, ONLY: filtreg_v |
20 |
use rotatf_m, only: rotatf |
21 |
|
22 |
IMPLICIT NONE |
23 |
|
24 |
|
25 |
|
26 |
! .......... variables en arguments ............. |
27 |
|
28 |
INTEGER, INTENT (IN) :: klevel |
29 |
REAL rotin(iim + 1, jjm, klevel), rotout(ip1jm, klevel) |
30 |
|
31 |
! .......... variables locales ................ |
32 |
|
33 |
REAL ghy(ip1jm, klevel), ghx(ip1jmp1, klevel) |
34 |
! ........................................................ |
35 |
|
36 |
|
37 |
CALL filtreg_v(rotin, intensive = .true.) |
38 |
|
39 |
CALL nxgrad(klevel, rotin, ghx, ghy) |
40 |
CALL rotatf(klevel, ghx, ghy, rotout) |
41 |
|
42 |
RETURN |
43 |
END SUBROUTINE laplacien_rot |