GCC Code Coverage Report

Directory: ./
File: phys/surf_landice_mod.f90
Date: 2022-01-11 19:19:34
Exec Total Coverage
Lines: 63 80 78.8%
Branches: 78 118 66.1%
Line Branch Exec Source
1 !
2 MODULE surf_landice_mod
3
4 IMPLICIT NONE
5
6 CONTAINS
7 !
8 !****************************************************************************************
9 !
10 1440 SUBROUTINE surf_landice(itime, dtime, knon, knindex, &
11 rlon, rlat, debut, lafin, &
12 rmu0, lwdownm, albedo, pphi1, &
13 swnet, lwnet, tsurf, p1lay, &
14 cdragh, cdragm, precip_rain, precip_snow, temp_air, spechum, &
15 AcoefH, AcoefQ, BcoefH, BcoefQ, &
16 AcoefU, AcoefV, BcoefU, BcoefV, &
17 ps, u1, v1, gustiness, rugoro, pctsrf, &
18 snow, qsurf, qsol, agesno, &
19 480 tsoil, z0m, z0h, SFRWL, alb_dir, alb_dif, evap, fluxsens, fluxlat, &
20 tsurf_new, dflux_s, dflux_l, &
21 alt, slope, cloudf, &
22 snowhgt, qsnow, to_ice, sissnow, &
23 480 alb3, runoff, &
24 flux_u1, flux_v1)
25
26 USE dimphy
27 USE geometry_mod, ONLY : longitude,latitude
28 USE surface_data, ONLY : type_ocean, calice, calsno, landice_opt, iflag_albcalc
29 USE fonte_neige_mod, ONLY : fonte_neige,run_off_lic,fqcalving_global,ffonte_global,fqfonte_global,runofflic_global
30 USE cpl_mod, ONLY : cpl_send_landice_fields
31 USE calcul_fluxs_mod
32 USE phys_output_var_mod
33 !FC
34 USE ioipsl_getin_p_mod, ONLY : getin_p
35
36
37
38 USE indice_sol_mod
39
40 ! INCLUDE "indicesol.h"
41 INCLUDE "dimsoil.h"
42 INCLUDE "YOMCST.h"
43 INCLUDE "clesphys.h"
44
45 ! Input variables
46 !****************************************************************************************
47 INTEGER, INTENT(IN) :: itime, knon
48 INTEGER, DIMENSION(klon), INTENT(in) :: knindex
49 REAL, INTENT(in) :: dtime
50 REAL, DIMENSION(klon), INTENT(IN) :: swnet ! net shortwave radiance
51 REAL, DIMENSION(klon), INTENT(IN) :: lwnet ! net longwave radiance
52 REAL, DIMENSION(klon), INTENT(IN) :: tsurf
53 REAL, DIMENSION(klon), INTENT(IN) :: p1lay
54 REAL, DIMENSION(klon), INTENT(IN) :: cdragh, cdragm
55 REAL, DIMENSION(klon), INTENT(IN) :: precip_rain, precip_snow
56 REAL, DIMENSION(klon), INTENT(IN) :: temp_air, spechum
57 REAL, DIMENSION(klon), INTENT(IN) :: AcoefH, AcoefQ
58 REAL, DIMENSION(klon), INTENT(IN) :: BcoefH, BcoefQ
59 REAL, DIMENSION(klon), INTENT(IN) :: AcoefU, AcoefV, BcoefU, BcoefV
60 REAL, DIMENSION(klon), INTENT(IN) :: ps
61 REAL, DIMENSION(klon), INTENT(IN) :: u1, v1, gustiness
62 REAL, DIMENSION(klon), INTENT(IN) :: rugoro
63 REAL, DIMENSION(klon,nbsrf), INTENT(IN) :: pctsrf
64
65 LOGICAL, INTENT(IN) :: debut !true if first step
66 LOGICAL, INTENT(IN) :: lafin !true if last step
67 REAL, DIMENSION(klon), INTENT(IN) :: rlon, rlat
68 REAL, DIMENSION(klon), INTENT(IN) :: rmu0
69 REAL, DIMENSION(klon), INTENT(IN) :: lwdownm !ylwdown
70 REAL, DIMENSION(klon), INTENT(IN) :: albedo !mean albedo
71 REAL, DIMENSION(klon), INTENT(IN) :: pphi1
72 REAL, DIMENSION(klon), INTENT(IN) :: alt !mean altitude of the grid box
73 REAL, DIMENSION(klon), INTENT(IN) :: slope !mean slope in grid box
74 REAL, DIMENSION(klon), INTENT(IN) :: cloudf !total cloud fraction
75
76 ! In/Output variables
77 !****************************************************************************************
78 REAL, DIMENSION(klon), INTENT(INOUT) :: snow, qsol
79 REAL, DIMENSION(klon), INTENT(INOUT) :: agesno
80 REAL, DIMENSION(klon, nsoilmx), INTENT(INOUT) :: tsoil
81
82 ! Output variables
83 !****************************************************************************************
84 REAL, DIMENSION(klon), INTENT(OUT) :: qsurf
85 REAL, DIMENSION(klon), INTENT(OUT) :: z0m, z0h
86 !albedo SB >>>
87 ! REAL, DIMENSION(klon), INTENT(OUT) :: alb1 ! new albedo in visible SW interval
88 ! REAL, DIMENSION(klon), INTENT(OUT) :: alb2 ! new albedo in near IR interval
89 REAL, DIMENSION(6), INTENT(IN) :: SFRWL
90 REAL, DIMENSION(klon,nsw), INTENT(OUT) :: alb_dir,alb_dif
91 !albedo SB <<<
92 REAL, DIMENSION(klon), INTENT(OUT) :: evap, fluxsens, fluxlat
93 REAL, DIMENSION(klon), INTENT(OUT) :: tsurf_new
94 REAL, DIMENSION(klon), INTENT(OUT) :: dflux_s, dflux_l
95 REAL, DIMENSION(klon), INTENT(OUT) :: flux_u1, flux_v1
96
97 REAL, DIMENSION(klon), INTENT(OUT) :: alb3
98 REAL, DIMENSION(klon), INTENT(OUT) :: qsnow !column water in snow [kg/m2]
99 REAL, DIMENSION(klon), INTENT(OUT) :: snowhgt !Snow height (m)
100 REAL, DIMENSION(klon), INTENT(OUT) :: to_ice
101 REAL, DIMENSION(klon), INTENT(OUT) :: sissnow
102 REAL, DIMENSION(klon), INTENT(OUT) :: runoff !Land ice runoff
103
104
105 ! Local variables
106 !****************************************************************************************
107 960 REAL, DIMENSION(klon) :: soilcap, soilflux
108 960 REAL, DIMENSION(klon) :: cal, beta, dif_grnd
109 960 REAL, DIMENSION(klon) :: zfra, alb_neig
110 960 REAL, DIMENSION(klon) :: radsol
111 960 REAL, DIMENSION(klon) :: u0, v0, u1_lay, v1_lay, ustar
112 INTEGER :: i,j,nt
113 REAL, DIMENSION(klon) :: fqfonte,ffonte
114 REAL, DIMENSION(klon) :: emis_new !Emissivity
115 REAL, DIMENSION(klon) :: swdown,lwdown
116 REAL, DIMENSION(klon) :: precip_snow_adv, snow_adv !Snow Drift precip./advection (not used in inlandsis)
117 REAL, DIMENSION(klon) :: erod !erosion of surface snow (flux, kg/m2/s like evap)
118 REAL, DIMENSION(klon) :: zsl_height, wind_velo !surface layer height, wind spd
119 REAL, DIMENSION(klon) :: dens_air, snow_cont_air !air density; snow content air
120 REAL, DIMENSION(klon) :: alb_soil !albedo of underlying ice
121 REAL, DIMENSION(klon) :: pexner !Exner potential
122 REAL :: pref
123 REAL, DIMENSION(klon,nsoilmx) :: tsoil0 !modif
124 REAL :: dtis ! subtimestep
125 LOGICAL :: debut_is, lafin_is ! debut and lafin for inlandsis
126
127 CHARACTER (len = 20) :: modname = 'surf_landice'
128 CHARACTER (len = 80) :: abort_message
129
130
131 960 REAL,DIMENSION(klon) :: alb1,alb2
132 960 REAL, DIMENSION (klon,6) :: alb6
133 ! End definition
134 !****************************************************************************************
135 !FC
136 !FC
137 REAL,SAVE :: alb_vis_sno_lic
138 !$OMP THREADPRIVATE(alb_vis_sno_lic)
139 REAL,SAVE :: alb_nir_sno_lic
140 !$OMP THREADPRIVATE(alb_nir_sno_lic)
141 LOGICAL, SAVE :: firstcall = .TRUE.
142 !$OMP THREADPRIVATE(firstcall)
143
144
145 !FC firtscall initializations
146 !******************************************************************************************
147
2/2
✓ Branch 0 taken 1 times.
✓ Branch 1 taken 479 times.
 480 IF (firstcall) THEN
148 1 alb_vis_sno_lic=0.77
149 1 CALL getin_p('alb_vis_sno_lic',alb_vis_sno_lic)
150 1 PRINT*, 'alb_vis_sno_lic',alb_vis_sno_lic
151 1 alb_nir_sno_lic=0.77
152 1 CALL getin_p('alb_nir_sno_lic',alb_nir_sno_lic)
153 1 PRINT*, 'alb_nir_sno_lic',alb_nir_sno_lic
154
155 ! z0m=1.e-3
156 ! z0h = z0m
157 1 firstcall=.false.
158 ENDIF
159 !******************************************************************************************
160
161 ! Initialize output variables
162
2/2
✓ Branch 0 taken 477120 times.
✓ Branch 1 taken 480 times.
 477600 alb3(:) = 999999.
163
2/2
✓ Branch 0 taken 477120 times.
✓ Branch 1 taken 480 times.
 477600 alb2(:) = 999999.
164
2/2
✓ Branch 0 taken 477120 times.
✓ Branch 1 taken 480 times.
 477600 alb1(:) = 999999.
165
166
2/2
✓ Branch 0 taken 477120 times.
✓ Branch 1 taken 480 times.
 477600 runoff(:) = 0.
167 !****************************************************************************************
168 ! Calculate total absorbed radiance at surface
169 !
170 !****************************************************************************************
171
2/2
✓ Branch 0 taken 477120 times.
✓ Branch 1 taken 480 times.
 477600 radsol(:) = 0.0
172
2/2
✓ Branch 0 taken 72960 times.
✓ Branch 1 taken 480 times.
 73440 radsol(1:knon) = swnet(1:knon) + lwnet(1:knon)
173
174 !****************************************************************************************
175 ! landice_opt = 0 : soil_model, calcul_flux, fonte_neige, ...
176 ! landice_opt = 1 : prepare and call INterace Lmdz SISvat (INLANDSIS)
177 !****************************************************************************************
178
179
180
1/2
✗ Branch 0 not taken.
✓ Branch 1 taken 480 times.
 480 IF (landice_opt .EQ. 1) THEN
181
182 !****************************************************************************************
183 ! CALL to INLANDSIS interface
184 !****************************************************************************************
185 ✗ abort_message='Pb de coherence: landice_opt = 1 mais CPP_INLANDSIS = .false.'
186 ✗ CALL abort_physic(modname,abort_message,1)
187
188
189
190 ELSE
191
192 !****************************************************************************************
193 ! Soil calculations
194 !
195 !****************************************************************************************
196
197 ! EV: use calbeta
198 480 CALL calbeta(dtime, is_lic, knon, snow, qsol, beta, cal, dif_grnd)
199
200
201 ! use soil model and recalculate properly cal
202
1/2
✓ Branch 0 taken 480 times.
✗ Branch 1 not taken.
 480 IF (soil_model) THEN
203 CALL soil(dtime, is_lic, knon, snow, tsurf, qsol, &
204
5/6
✓ Branch 0 taken 480 times.
✗ Branch 1 not taken.
✓ Branch 2 taken 72960 times.
✓ Branch 3 taken 480 times.
✓ Branch 4 taken 72960 times.
✓ Branch 5 taken 480 times.
 146400 & longitude(knindex(1:knon)), latitude(knindex(1:knon)), tsoil, soilcap, soilflux)
205
2/2
✓ Branch 0 taken 72960 times.
✓ Branch 1 taken 480 times.
 73440 cal(1:knon) = RCPD / soilcap(1:knon)
206
2/2
✓ Branch 0 taken 72960 times.
✓ Branch 1 taken 480 times.
 73440 radsol(1:knon) = radsol(1:knon) + soilflux(1:knon)
207 ELSE
208 ✗ cal = RCPD * calice
209 ✗ WHERE (snow > 0.0) cal = RCPD * calsno
210 ENDIF
211
212
213 !****************************************************************************************
214 ! Calulate fluxes
215 !
216 !****************************************************************************************
217 ! beta(:) = 1.0
218 ! dif_grnd(:) = 0.0
219
220 ! Suppose zero surface speed
221
2/2
✓ Branch 0 taken 477120 times.
✓ Branch 1 taken 480 times.
 477600 u0(:)=0.0
222
2/2
✓ Branch 0 taken 477120 times.
✓ Branch 1 taken 480 times.
 477600 v0(:)=0.0
223
2/2
✓ Branch 0 taken 477120 times.
✓ Branch 1 taken 480 times.
 477600 u1_lay(:) = u1(:) - u0(:)
224
2/2
✓ Branch 0 taken 477120 times.
✓ Branch 1 taken 480 times.
 477600 v1_lay(:) = v1(:) - v0(:)
225
226 CALL calcul_fluxs(knon, is_lic, dtime, &
227 tsurf, p1lay, cal, beta, cdragh, cdragh, ps, &
228 precip_rain, precip_snow, snow, qsurf, &
229 radsol, dif_grnd, temp_air, spechum, u1_lay, v1_lay, gustiness, &
230 1.,AcoefH, AcoefQ, BcoefH, BcoefQ, &
231 480 tsurf_new, evap, fluxlat, fluxsens, dflux_s, dflux_l)
232
233 CALL calcul_flux_wind(knon, dtime, &
234 u0, v0, u1, v1, gustiness, cdragm, &
235 AcoefU, AcoefV, BcoefU, BcoefV, &
236 p1lay, temp_air, &
237 480 flux_u1, flux_v1)
238
239 !****************************************************************************************
240 ! Calculate snow height, age, run-off,..
241 !
242 !****************************************************************************************
243 CALL fonte_neige(knon, is_lic, knindex, dtime, &
244 tsurf, precip_rain, precip_snow, &
245 480 snow, qsol, tsurf_new, evap)
246
247
248 !****************************************************************************************
249 ! Calculate albedo
250 !
251 !****************************************************************************************
252 480 CALL albsno(klon,knon,dtime,agesno(:),alb_neig(:), precip_snow(:))
253
254
4/4
✓ Branch 0 taken 72960 times.
✓ Branch 1 taken 480 times.
✓ Branch 2 taken 20003 times.
✓ Branch 3 taken 52957 times.
 73440 WHERE (snow(1 : knon) .LT. 0.0001) agesno(1 : knon) = 0.
255
2/2
✓ Branch 0 taken 72960 times.
✓ Branch 1 taken 480 times.
 73440 zfra(1:knon) = MAX(0.0,MIN(1.0,snow(1:knon)/(snow(1:knon)+10.0)))
256 alb1(1:knon) = alb_neig(1:knon)*zfra(1:knon) + &
257
2/2
✓ Branch 0 taken 72960 times.
✓ Branch 1 taken 480 times.
 73440 0.6 * (1.0-zfra(1:knon))
258 !
259 !IM: plusieurs choix/tests sur l'albedo des "glaciers continentaux"
260 ! alb1(1 : knon) = 0.6 !IM cf FH/GK
261 ! alb1(1 : knon) = 0.82
262 ! alb1(1 : knon) = 0.77 !211003 Ksta0.77
263 ! alb1(1 : knon) = 0.8 !KstaTER0.8 & LMD_ARMIP5
264 !IM: KstaTER0.77 & LMD_ARMIP6
265
266 ! Attantion: alb1 and alb2 are not the same!
267
2/2
✓ Branch 0 taken 72960 times.
✓ Branch 1 taken 480 times.
 73440 alb1(1:knon) = alb_vis_sno_lic
268
2/2
✓ Branch 0 taken 72960 times.
✓ Branch 1 taken 480 times.
 73440 alb2(1:knon) = alb_nir_sno_lic
269
270
271 !****************************************************************************************
272 ! Rugosity
273 !
274 !****************************************************************************************
275
2/2
✓ Branch 0 taken 477120 times.
✓ Branch 1 taken 480 times.
 477600 z0m=1.e-3
276
2/2
✓ Branch 0 taken 477120 times.
✓ Branch 1 taken 480 times.
 477600 z0h = z0m
277
2/2
✓ Branch 0 taken 477120 times.
✓ Branch 1 taken 480 times.
 477600 z0m = SQRT(z0m**2+rugoro**2)
278
279
280
281 END IF ! landice_opt
282
283
284 !****************************************************************************************
285 ! Send run-off on land-ice to coupler if coupled ocean.
286 ! run_off_lic has been calculated in fonte_neige or surf_inlandsis
287 !
288 !****************************************************************************************
289
1/2
✗ Branch 0 not taken.
✓ Branch 1 taken 480 times.
 480 IF (type_ocean=='couple') THEN
290 ✗ CALL cpl_send_landice_fields(itime, knon, knindex, run_off_lic)
291 ENDIF
292
293 ! transfer runoff rate [kg/m2/s](!) to physiq for output
294
2/2
✓ Branch 0 taken 72960 times.
✓ Branch 1 taken 480 times.
 73440 runoff(1:knon)=run_off_lic(1:knon)/dtime
295
296
297 !****************************************************************************************
298 ! Etienne: comment these lines because of duplication just below
299 ! snow_o=0.
300 ! zfra_o = 0.
301 ! DO j = 1, knon
302 ! i = knindex(j)
303 ! snow_o(i) = snow(j)
304 ! zfra_o(i) = zfra(j)
305 ! ENDDO
306 !
307 !****************************************************************************************
308
2/2
✓ Branch 0 taken 477120 times.
✓ Branch 1 taken 480 times.
 477600 snow_o=0.
309
2/2
✓ Branch 0 taken 477120 times.
✓ Branch 1 taken 480 times.
 477600 zfra_o = 0.
310
2/2
✓ Branch 0 taken 72960 times.
✓ Branch 1 taken 480 times.
 73440 DO j = 1, knon
311 72960 i = knindex(j)
312 72960 snow_o(i) = snow(j)
313 73440 zfra_o(i) = zfra(j)
314 ENDDO
315
316
317 !albedo SB >>>
318 select case(NSW)
319 case(2)
320 ✗ alb_dir(1:knon,1)=alb1(1:knon)
321 ✗ alb_dir(1:knon,2)=alb2(1:knon)
322 case(4)
323 ✗ alb_dir(1:knon,1)=alb1(1:knon)
324 ✗ alb_dir(1:knon,2)=alb2(1:knon)
325 ✗ alb_dir(1:knon,3)=alb2(1:knon)
326 ✗ alb_dir(1:knon,4)=alb2(1:knon)
327 case(6)
328
2/2
✓ Branch 0 taken 72960 times.
✓ Branch 1 taken 480 times.
 73440 alb_dir(1:knon,1)=alb1(1:knon)
329
2/2
✓ Branch 0 taken 72960 times.
✓ Branch 1 taken 480 times.
 73440 alb_dir(1:knon,2)=alb1(1:knon)
330
2/2
✓ Branch 0 taken 72960 times.
✓ Branch 1 taken 480 times.
 73440 alb_dir(1:knon,3)=alb1(1:knon)
331
2/2
✓ Branch 0 taken 72960 times.
✓ Branch 1 taken 480 times.
 73440 alb_dir(1:knon,4)=alb2(1:knon)
332
2/2
✓ Branch 0 taken 72960 times.
✓ Branch 1 taken 480 times.
 73440 alb_dir(1:knon,5)=alb2(1:knon)
333
2/2
✓ Branch 0 taken 72960 times.
✓ Branch 1 taken 480 times.
 73440 alb_dir(1:knon,6)=alb2(1:knon)
334
335
2/8
✗ Branch 0 not taken.
✗ Branch 1 not taken.
✗ Branch 2 not taken.
✓ Branch 3 taken 480 times.
✓ Branch 4 taken 480 times.
✗ Branch 5 not taken.
✗ Branch 6 not taken.
✗ Branch 7 not taken.
 960 IF ((landice_opt .EQ. 1) .AND. (iflag_albcalc .EQ. 2)) THEN
336 ✗ alb_dir(1:knon,1)=alb6(1:knon,1)
337 ✗ alb_dir(1:knon,2)=alb6(1:knon,2)
338 ✗ alb_dir(1:knon,3)=alb6(1:knon,3)
339 ✗ alb_dir(1:knon,4)=alb6(1:knon,4)
340 ✗ alb_dir(1:knon,5)=alb6(1:knon,5)
341 ✗ alb_dir(1:knon,6)=alb6(1:knon,6)
342 ENDIF
343
344 end select
345
4/4
✓ Branch 0 taken 2880 times.
✓ Branch 1 taken 480 times.
✓ Branch 2 taken 2862720 times.
✓ Branch 3 taken 2880 times.
 2866080 alb_dif=alb_dir
346 !albedo SB <<<
347
348
349
350
351 480 END SUBROUTINE surf_landice
352 !
353 !****************************************************************************************
354 !
355 END MODULE surf_landice_mod
356
357
358
359