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 |
IMPLICIT NONE |
USE comgeom |
19 |
c |
USE filtreg_m, ONLY: filtreg |
20 |
|
IMPLICIT NONE |
21 |
c |
|
22 |
c .......... variables en arguments ............. |
|
23 |
c |
|
24 |
INTEGER klevel |
! .......... variables en arguments ............. |
25 |
REAL rotin( ip1jm,klevel ), rotout( ip1jm,klevel ) |
|
26 |
c |
INTEGER, INTENT (IN) :: klevel |
27 |
c .......... variables locales ................ |
REAL rotin(iim + 1, jjm, klevel), rotout(ip1jm, klevel) |
28 |
c |
|
29 |
REAL ghy(ip1jm,klevel), ghx(ip1jmp1,klevel) |
! .......... variables locales ................ |
30 |
c ........................................................ |
|
31 |
c |
REAL ghy(ip1jm, klevel), ghx(ip1jmp1, klevel) |
32 |
c |
! ........................................................ |
33 |
CALL filtreg ( rotin , jjm, klevel, 2, 1, .FALSE., 1 ) |
|
34 |
|
|
35 |
CALL nxgrad ( klevel, rotin, ghx , ghy ) |
CALL filtreg(rotin, jjm, klevel, 2, 1, .FALSE.) |
36 |
CALL rotatf ( klevel, ghx , ghy , rotout ) |
|
37 |
c |
CALL nxgrad(klevel, rotin, ghx, ghy) |
38 |
RETURN |
CALL rotatf(klevel, ghx, ghy, rotout) |
39 |
END |
|
40 |
|
RETURN |
41 |
|
END SUBROUTINE laplacien_rot |