GCC Code Coverage Report

Directory: ./
File: phys/cvltr_noscav.f90
Date: 2022-01-11 19:19:34
Exec Total Coverage
Lines: 61 61 100.0%
Branches: 68 68 100.0%
Line Branch Exec Source
1 !
2 ! $Id $
3 !
4 1920 SUBROUTINE cvltr_noscav(it,pdtime,da, phi, mp,wght_cvfd,paprs,pplay,x,upd,dnd,dx)
5 USE dimphy
6 USE infotrac_phy, ONLY : nbtr
7 IMPLICIT NONE
8 !=====================================================================
9 ! Objet : convection des traceurs / KE
10 ! Auteurs: M-A Filiberti and J-Y Grandpeix
11 !=====================================================================
12 include "YOMCST.h"
13 include "YOECUMF.h"
14
15 ! Entree
16 REAL,INTENT(IN) :: pdtime
17 INTEGER, INTENT(IN) :: it
18 REAL,DIMENSION(klon,klev),INTENT(IN) :: da
19 REAL,DIMENSION(klon,klev,klev),INTENT(IN) :: phi
20 REAL,DIMENSION(klon,klev),INTENT(IN) :: mp
21 REAL,DIMENSION(klon,klev),INTENT(IN) :: wght_cvfd ! weights of the layers feeding convection
22 REAL,DIMENSION(klon,klev+1),INTENT(IN) :: paprs ! pression aux 1/2 couches (bas en haut)
23 REAL,DIMENSION(klon,klev),INTENT(IN) :: pplay ! pression pour le milieu de chaque couche
24 REAL,DIMENSION(klon,klev,nbtr),INTENT(IN) :: x ! q de traceur (bas en haut)
25 REAL,DIMENSION(klon,klev),INTENT(IN) :: upd ! saturated updraft mass flux
26 REAL,DIMENSION(klon,klev),INTENT(IN) :: dnd ! saturated downdraft mass flux
27
28 ! Sortie
29 REAL,DIMENSION(klon,klev,nbtr),INTENT(OUT) :: dx ! tendance de traceur (bas en haut)
30
31 ! Variables locales
32 ! REAL,DIMENSION(klon,klev) :: zed
33 1920 REAL,DIMENSION(klon,klev,klev) :: zmd
34 1920 REAL,DIMENSION(klon,klev,klev) :: za
35 1920 REAL,DIMENSION(klon,klev) :: zmfd,zmfa
36 1920 REAL,DIMENSION(klon,klev) :: zmfp,zmfu
37 1920 REAL,DIMENSION(klon,nbtr) :: qfeed ! tracer concentration feeding convection
38 1920 REAL,DIMENSION(klon,klev) :: deltap
39 INTEGER :: i,k,j
40 REAL :: pdtimeRG
41 REAL :: smallest_mp
42 real conserv
43 real smfd
44 real smfu
45 real smfa
46 real smfp
47 ! =========================================
48 ! calcul des tendances liees au downdraft
49 ! =========================================
50 !
51 smallest_mp = tiny(mp(1,1))
52 !cdir collapse
53
2/2
✓ Branch 0 taken 960 times.
✓ Branch 1 taken 954240 times.
 955200 qfeed(:,it) = 0.
54
2/2
✓ Branch 0 taken 37440 times.
✓ Branch 1 taken 960 times.
 38400 DO j=1,klev
55
2/2
✓ Branch 0 taken 37215360 times.
✓ Branch 1 taken 37440 times.
 37253760 DO i=1,klon
56 ! zed(i,j)=0.
57 37215360 zmfd(i,j)=0.
58 37215360 zmfa(i,j)=0.
59 37215360 zmfu(i,j)=0.
60 37252800 zmfp(i,j)=0.
61 END DO
62 END DO
63 !cdir collapse
64
2/2
✓ Branch 0 taken 37440 times.
✓ Branch 1 taken 960 times.
 38400 DO k=1,klev
65
2/2
✓ Branch 0 taken 1460160 times.
✓ Branch 1 taken 37440 times.
 1498560 DO j=1,klev
66
2/2
✓ Branch 0 taken 1451399040 times.
✓ Branch 1 taken 1460160 times.
 1452896640 DO i=1,klon
67 1451399040 zmd(i,j,k)=0.
68 1452859200 za (i,j,k)=0.
69 END DO
70 END DO
71 END DO
72 ! entrainement
73 ! DO k=1,klev-1
74 ! DO i=1,klon
75 ! zed(i,k)=max(0.,mp(i,k)-mp(i,k+1))
76 ! END DO
77 ! END DO
78
79 ! calcul de la matrice d echange
80 ! matrice de distribution de la masse entrainee en k
81
82
2/2
✓ Branch 0 taken 36480 times.
✓ Branch 1 taken 960 times.
 37440 DO k=1,klev-1
83
2/2
✓ Branch 0 taken 36261120 times.
✓ Branch 1 taken 36480 times.
 36298560 DO i=1,klon
84 36297600 zmd(i,k,k)=max(0.,mp(i,k)-mp(i,k+1))
85 END DO
86 END DO
87
2/2
✓ Branch 0 taken 960 times.
✓ Branch 1 taken 36480 times.
 37440 DO k=2,klev
88
2/2
✓ Branch 0 taken 711360 times.
✓ Branch 1 taken 36480 times.
 748800 DO j=k-1,1,-1
89
2/2
✓ Branch 0 taken 707091840 times.
✓ Branch 1 taken 711360 times.
 707839680 DO i=1,klon
90 !! if(mp(i,j+1).ne.0) then
91 !! zmd(i,j,k)=zmd(i,j+1,k)*min(1.,mp(i,j)/mp(i,j+1))
92 !! ENDif
93 707803200 zmd(i,j,k)=zmd(i,j+1,k)*mp(i,j)/max(mp(i,j),mp(i,j+1),smallest_mp)
94 END DO
95 END DO
96 END DO
97
2/2
✓ Branch 0 taken 37440 times.
✓ Branch 1 taken 960 times.
 38400 DO k=1,klev
98
2/2
✓ Branch 0 taken 1422720 times.
✓ Branch 1 taken 37440 times.
 1461120 DO j=1,klev-1
99
2/2
✓ Branch 0 taken 1414183680 times.
✓ Branch 1 taken 1422720 times.
 1415643840 DO i=1,klon
100 1415606400 za(i,j,k)=max(0.,zmd(i,j+1,k)-zmd(i,j,k))
101 END DO
102 END DO
103 END DO
104 !
105 ! rajout du terme lie a l ascendance induite
106 !
107
2/2
✓ Branch 0 taken 36480 times.
✓ Branch 1 taken 960 times.
 37440 DO j=2,klev
108
2/2
✓ Branch 0 taken 36261120 times.
✓ Branch 1 taken 36480 times.
 36298560 DO i=1,klon
109 36297600 za(i,j,j-1)=za(i,j,j-1)+mp(i,j)
110 END DO
111 END DO
112 !
113 ! tendances
114 !
115
2/2
✓ Branch 0 taken 37440 times.
✓ Branch 1 taken 960 times.
 38400 DO k=1,klev
116
2/2
✓ Branch 0 taken 1460160 times.
✓ Branch 1 taken 37440 times.
 1498560 DO j=1,klev
117
2/2
✓ Branch 0 taken 1451399040 times.
✓ Branch 1 taken 1460160 times.
 1452896640 DO i=1,klon
118 1452859200 zmfd(i,j)=zmfd(i,j)+za(i,j,k)*(x(i,k,it)-x(i,j,it))
119 END DO
120 END DO
121 END DO
122 !
123 ! =========================================
124 ! calcul des tendances liees aux flux satures
125 ! =========================================
126 !RL
127 ! Feeding concentrations
128
2/2
✓ Branch 0 taken 37440 times.
✓ Branch 1 taken 960 times.
 38400 DO j=1,klev
129
2/2
✓ Branch 0 taken 37215360 times.
✓ Branch 1 taken 37440 times.
 37253760 DO i=1,klon
130 37252800 qfeed(i,it)=qfeed(i,it)+wght_cvfd(i,j)*x(i,j,it)
131 END DO
132 END DO
133 !RL
134 !
135
2/2
✓ Branch 0 taken 37440 times.
✓ Branch 1 taken 960 times.
 38400 DO j=1,klev
136
2/2
✓ Branch 0 taken 37215360 times.
✓ Branch 1 taken 37440 times.
 37253760 DO i=1,klon
137 !RL
138 !! zmfa(i,j,it)=da(i,j)*(x(i,1,it)-x(i,j,it)) ! da
139 37252800 zmfa(i,j)=da(i,j)*(qfeed(i,it)-x(i,j,it)) ! da
140 !RL
141 END DO
142 END DO
143 !
144 !! print *,'it, qfeed(1,it), x(1,1,it) ', it, qfeed(1,it), x(1,1,it) !jyg
145 !! print *,'wght_cvfd ', (j, wght_cvfd(1,j), j=1,5) !jyg
146 !
147
2/2
✓ Branch 0 taken 37440 times.
✓ Branch 1 taken 960 times.
 38400 DO k=1,klev
148
2/2
✓ Branch 0 taken 1460160 times.
✓ Branch 1 taken 37440 times.
 1498560 DO j=1,klev
149
2/2
✓ Branch 0 taken 1451399040 times.
✓ Branch 1 taken 1460160 times.
 1452896640 DO i=1,klon
150 1452859200 zmfp(i,j)=zmfp(i,j)+phi(i,j,k)*(x(i,k,it)-x(i,j,it))
151 END DO
152 END DO
153 END DO
154
2/2
✓ Branch 0 taken 36480 times.
✓ Branch 1 taken 960 times.
 37440 DO j=1,klev-1
155
2/2
✓ Branch 0 taken 36261120 times.
✓ Branch 1 taken 36480 times.
 36298560 DO i=1,klon
156 36297600 zmfu(i,j)=max(0.,upd(i,j+1)+dnd(i,j+1))*(x(i,j+1,it)-x(i,j,it))
157 END DO
158 END DO
159
2/2
✓ Branch 0 taken 36480 times.
✓ Branch 1 taken 960 times.
 37440 DO j=2,klev
160
2/2
✓ Branch 0 taken 36261120 times.
✓ Branch 1 taken 36480 times.
 36298560 DO i=1,klon
161 36297600 zmfu(i,j)=zmfu(i,j)+min(0.,upd(i,j)+dnd(i,j))*(x(i,j,it)-x(i,j-1,it))
162 END DO
163 END DO
164
165 ! =========================================
166 ! calcul final des tendances
167 ! =========================================
168
2/2
✓ Branch 0 taken 37440 times.
✓ Branch 1 taken 960 times.
 38400 DO k=1, klev
169
2/2
✓ Branch 0 taken 37215360 times.
✓ Branch 1 taken 37440 times.
 37253760 DO i=1, klon
170 37252800 deltap(i,k)=paprs(i,k)-paprs(i,k+1)
171 ENDDO
172 ENDDO
173 960 pdtimeRG=pdtime*RG
174 !cdir collapse
175
2/2
✓ Branch 0 taken 37440 times.
✓ Branch 1 taken 960 times.
 38400 DO k=1, klev
176
2/2
✓ Branch 0 taken 37215360 times.
✓ Branch 1 taken 37440 times.
 37253760 DO i=1, klon
177 dx(i,k,it)=(zmfd(i,k)+zmfu(i,k) &
178 37252800 +zmfa(i,k)+zmfp(i,k))*pdtimeRG/deltap(i,k)
179 ENDDO
180 ENDDO
181
182 !! test de conservation du traceur
183 conserv=0.
184 smfd = 0.
185 smfu = 0.
186 smfa = 0.
187 smfp = 0.
188 DO k=1, klev
189 DO i=1, klon
190 conserv=conserv+dx(i,k,it)* &
191 deltap(i,k)/RG
192 smfd = smfd + zmfd(i,k)*pdtime
193 smfu = smfu + zmfu(i,k)*pdtime
194 smfa = smfa + zmfa(i,k)*pdtime
195 smfp = smfp + zmfp(i,k)*pdtime
196 ENDDO
197 ENDDO
198 !! print *,'it',it,'cvltr_noscav conserv, smfd, smfu, smfa, smfp ',conserv, &
199 !! smfd, smfu, smfa, smfp
200
201 960 END SUBROUTINE cvltr_noscav
202