GCC Code Coverage Report


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