17 |
USE comgeom, ONLY: apoln, apols, cuvsurcv, cvusurcu, unsaire |
USE comgeom, ONLY: apoln, apols, cuvsurcv, cvusurcu, unsaire |
18 |
USE filtreg_m, ONLY: filtreg |
USE filtreg_m, ONLY: filtreg |
19 |
|
|
|
! div est un argument de sortie pour le s-prog |
|
|
|
|
|
! variables en arguments |
|
|
|
|
20 |
INTEGER, intent(in):: klevel |
INTEGER, intent(in):: klevel |
21 |
REAL, intent(in):: x(ip1jmp1, klevel), y(ip1jm, klevel) |
REAL, intent(in):: x(ip1jmp1, klevel), y(ip1jm, klevel) |
22 |
real div(ip1jmp1, klevel) |
real, intent(out):: div(ip1jmp1, klevel) ! in (unit of x, y) m-2 |
23 |
|
|
24 |
! variables locales |
! Variables locales : |
25 |
|
|
26 |
INTEGER l, ij |
INTEGER l, ij |
27 |
REAL aiy1(iip1) , aiy2(iip1) |
REAL aiy1(iim) , aiy2(iim) |
28 |
REAL sumypn, sumyps |
REAL sumypn, sumyps |
29 |
|
|
|
REAL SSUM |
|
|
|
|
30 |
!------------------------------------------------------------ |
!------------------------------------------------------------ |
31 |
|
|
32 |
DO l = 1, klevel |
DO l = 1, klevel |
46 |
aiy1(ij) = cuvsurcv(ij) * y(ij , l) |
aiy1(ij) = cuvsurcv(ij) * y(ij , l) |
47 |
aiy2(ij) = cuvsurcv(ij+ ip1jmi1) * y(ij+ ip1jmi1, l) |
aiy2(ij) = cuvsurcv(ij+ ip1jmi1) * y(ij+ ip1jmi1, l) |
48 |
ENDDO |
ENDDO |
49 |
sumypn = SSUM (iim, aiy1, 1) / apoln |
sumypn = SUM(aiy1) / apoln |
50 |
sumyps = SSUM (iim, aiy2, 1) / apols |
sumyps = SUM(aiy2) / apols |
51 |
|
|
52 |
DO ij = 1, iip1 |
DO ij = 1, iip1 |
53 |
div(ij , l) = - sumypn |
div(ij , l) = - sumypn |