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