Directory: | ./ |
---|---|
File: | dyn3d_common/flumass.f90 |
Date: | 2022-01-11 19:19:34 |
Exec | Total | Coverage | |
---|---|---|---|
Lines: | 31 | 31 | 100.0% |
Branches: | 38 | 38 | 100.0% |
Line | Branch | Exec | Source |
---|---|---|---|
1 | 2881 | SUBROUTINE flumass (massebx,masseby, vcont, ucont, pbaru, pbarv ) | |
2 | ! | ||
3 | !------------------------------------------------------------------------------- | ||
4 | ! Authors: P. Le Van , Fr. Hourdin. | ||
5 | !------------------------------------------------------------------------------- | ||
6 | ! Purpose: Compute mass flux at s levels. | ||
7 | IMPLICIT NONE | ||
8 | include "dimensions.h" | ||
9 | include "paramet.h" | ||
10 | include "comgeom.h" | ||
11 | !=============================================================================== | ||
12 | ! Arguments: | ||
13 | REAL, INTENT(IN) :: massebx(ip1jmp1,llm) | ||
14 | REAL, INTENT(IN) :: masseby(ip1jm ,llm) | ||
15 | REAL, INTENT(IN) :: vcont (ip1jm ,llm) | ||
16 | REAL, INTENT(IN) :: ucont (ip1jmp1,llm) | ||
17 | REAL, INTENT(OUT) :: pbaru (ip1jmp1,llm) | ||
18 | REAL, INTENT(OUT) :: pbarv (ip1jm ,llm) | ||
19 | !=============================================================================== | ||
20 | ! Method used: A 2 equations system is solved. | ||
21 | ! * 1st one describes divergence computation at pole point nr. i (i=1 to im): | ||
22 | ! (0.5*(pbaru(i)-pbaru(i-1))-pbarv(i))/aire(i) = - SUM(pbarv(n))/aire pole | ||
23 | ! * 2nd one specifies that mean mass flux at pole is equal to 0: | ||
24 | ! SUM(pbaru(n)*local_area(n))=0 | ||
25 | ! This way, we determine additive constant common to pbary elements representing | ||
26 | ! pbaru(0,j,l) in divergence computation equation for point i=1. (i=1 to im) | ||
27 | !=============================================================================== | ||
28 | ! Local variables: | ||
29 | REAL :: sairen, saireun, ctn, ctn0, apbarun(iip1) | ||
30 | REAL :: saires, saireus, cts, cts0, apbarus(iip1) | ||
31 | INTEGER :: l, i | ||
32 | !=============================================================================== | ||
33 |
2/2✓ Branch 0 taken 2881 times.
✓ Branch 1 taken 112359 times.
|
115240 | DO l=1,llm |
34 |
2/2✓ Branch 0 taken 114943257 times.
✓ Branch 1 taken 112359 times.
|
115055616 | pbaru(iip2:ip1jm,l)=massebx(iip2:ip1jm,l)*ucont(iip2:ip1jm,l) |
35 |
2/2✓ Branch 0 taken 118651104 times.
✓ Branch 1 taken 112359 times.
|
118766344 | pbarv( 1:ip1jm,l)=masseby( 1:ip1jm,l)*vcont( 1:ip1jm,l) |
36 | END DO | ||
37 | |||
38 | !--- NORTH POLE | ||
39 |
2/2✓ Branch 0 taken 2881 times.
✓ Branch 1 taken 92192 times.
|
95073 | sairen =SUM(aire (1:iim)) |
40 |
2/2✓ Branch 0 taken 2881 times.
✓ Branch 1 taken 92192 times.
|
95073 | saireun=SUM(aireu(1:iim)) |
41 |
2/2✓ Branch 0 taken 2881 times.
✓ Branch 1 taken 112359 times.
|
115240 | DO l = 1,llm |
42 |
2/2✓ Branch 0 taken 3595488 times.
✓ Branch 1 taken 112359 times.
|
3707847 | ctn=SUM(pbarv(1:iim,l))/sairen |
43 | 112359 | pbaru(1,l)= pbarv(1,l)-ctn*aire(1) | |
44 |
2/2✓ Branch 0 taken 112359 times.
✓ Branch 1 taken 3483129 times.
|
3595488 | DO i=2,iim |
45 | 3595488 | pbaru(i,l)=pbaru(i-1,l)+pbarv(i,l)-ctn*aire(i) | |
46 | END DO | ||
47 |
2/2✓ Branch 0 taken 112359 times.
✓ Branch 1 taken 3595488 times.
|
3707847 | DO i=1,iim |
48 | 3707847 | apbarun(i)=aireu(i)*pbaru(i,l) | |
49 | END DO | ||
50 |
2/2✓ Branch 0 taken 3595488 times.
✓ Branch 1 taken 112359 times.
|
3707847 | ctn0 = -SUM(apbarun(1:iim))/saireun |
51 |
2/2✓ Branch 0 taken 3595488 times.
✓ Branch 1 taken 112359 times.
|
3707847 | DO i = 1,iim |
52 | 3707847 | pbaru(i,l)=2.*(pbaru(i,l)+ctn0) | |
53 | END DO | ||
54 | 115240 | pbaru(iip1,l)=pbaru(1,l) | |
55 | END DO | ||
56 | |||
57 | !--- SOUTH POLE | ||
58 |
2/2✓ Branch 0 taken 2881 times.
✓ Branch 1 taken 92192 times.
|
95073 | saires =SUM(aire (ip1jm+1:ip1jmp1-1)) |
59 |
2/2✓ Branch 0 taken 2881 times.
✓ Branch 1 taken 92192 times.
|
95073 | saireus=SUM(aireu(ip1jm+1:ip1jmp1-1)) |
60 |
2/2✓ Branch 0 taken 112359 times.
✓ Branch 1 taken 2881 times.
|
115240 | DO l = 1,llm |
61 |
2/2✓ Branch 0 taken 3595488 times.
✓ Branch 1 taken 112359 times.
|
3707847 | cts=SUM(pbarv(ip1jmi1+1:ip1jm-1,l))/saires |
62 | 112359 | pbaru(1+ip1jm,l)=-pbarv(1+ip1jmi1,l)+cts*aire(1+ip1jm) | |
63 |
2/2✓ Branch 0 taken 112359 times.
✓ Branch 1 taken 3483129 times.
|
3595488 | DO i=2,iim |
64 | 3595488 | pbaru(i+ip1jm,l)=pbaru(i-1+ip1jm,l)-pbarv(i+ip1jmi1,l)+cts*aire(i+ip1jm) | |
65 | END DO | ||
66 |
2/2✓ Branch 0 taken 112359 times.
✓ Branch 1 taken 3595488 times.
|
3707847 | DO i=1,iim |
67 | 3707847 | apbarus(i)=aireu(i+ip1jm)*pbaru(i+ip1jm,l) | |
68 | END DO | ||
69 |
2/2✓ Branch 0 taken 3595488 times.
✓ Branch 1 taken 112359 times.
|
3707847 | cts0 = -SUM(apbarus(1:iim))/saireus |
70 |
2/2✓ Branch 0 taken 3595488 times.
✓ Branch 1 taken 112359 times.
|
3707847 | DO i = 1,iim |
71 | 3707847 | pbaru(i+ip1jm,l)=2.*(pbaru(i+ip1jm,l)+cts0) | |
72 | END DO | ||
73 | 115240 | pbaru(ip1jmp1,l)=pbaru(1+ip1jm,l) | |
74 | END DO | ||
75 | |||
76 | 2881 | END SUBROUTINE flumass | |
77 |