GCC Code Coverage Report

Directory: ./
File: phys/ocean_forced_mod.f90
Date: 2022-01-11 19:19:34
Exec Total Coverage
Lines: 72 78 92.3%
Branches: 70 94 74.5%
Line Branch Exec Source
1 !
2 ! $Id: ocean_forced_mod.F90 3974 2021-07-30 14:27:30Z asima $
3 !
4 MODULE ocean_forced_mod
5 !
6 ! This module is used for both the sub-surfaces ocean and sea-ice for the case of a
7 ! forced ocean, "ocean=force".
8 !
9 IMPLICIT NONE
10
11 CONTAINS
12 !
13 !****************************************************************************************
14 !
15 1440 SUBROUTINE ocean_forced_noice( &
16 480 itime, dtime, jour, knon, knindex, &
17 p1lay, cdragh, cdragq, cdragm, precip_rain, precip_snow, &
18 temp_air, spechum, &
19 AcoefH, AcoefQ, BcoefH, BcoefQ, &
20 AcoefU, AcoefV, BcoefU, BcoefV, &
21 ps, u1, v1, gustiness, tsurf_in, &
22 radsol, snow, agesno, &
23 480 qsurf, evap, fluxsens, fluxlat, flux_u1, flux_v1, &
24
2/4
✗ Branch 0 not taken.
✓ Branch 1 taken 480 times.
✗ Branch 2 not taken.
✓ Branch 3 taken 480 times.
 480 tsurf_new, dflux_s, dflux_l, sens_prec_liq, rhoa)
25 !
26 ! This subroutine treats the "open ocean", all grid points that are not entierly covered
27 ! by ice.
28 ! The routine receives data from climatologie file limit.nc and does some calculations at the
29 ! surface.
30 !
31 USE dimphy
32 USE calcul_fluxs_mod
33 USE limit_read_mod
34 USE mod_grid_phy_lmdz
35 USE indice_sol_mod
36 USE phys_output_var_mod, ONLY : sens_prec_liq_o, sens_prec_sol_o, lat_prec_liq_o, lat_prec_sol_o
37 use config_ocean_skin_m, only: activate_ocean_skin
38
39 INCLUDE "YOMCST.h"
40 INCLUDE "clesphys.h"
41 INCLUDE "flux_arp.h"
42
43 ! Input arguments
44 !****************************************************************************************
45 INTEGER, INTENT(IN) :: itime, jour, knon
46 INTEGER, DIMENSION(klon), INTENT(IN) :: knindex
47 REAL, INTENT(IN) :: dtime
48 REAL, DIMENSION(klon), INTENT(IN) :: p1lay
49 REAL, DIMENSION(klon), INTENT(IN) :: cdragh, cdragq, cdragm
50 REAL, DIMENSION(klon), INTENT(IN) :: precip_rain, precip_snow
51 REAL, DIMENSION(klon), INTENT(IN) :: temp_air, spechum
52 REAL, DIMENSION(klon), INTENT(IN) :: AcoefH, AcoefQ, BcoefH, BcoefQ
53 REAL, DIMENSION(klon), INTENT(IN) :: AcoefU, AcoefV, BcoefU, BcoefV
54 REAL, DIMENSION(klon), INTENT(IN) :: ps
55 REAL, DIMENSION(klon), INTENT(IN) :: u1, v1, gustiness
56 REAL, DIMENSION(klon), INTENT(IN) :: tsurf_in
57 real, intent(in):: rhoa(:) ! (knon) density of moist air (kg / m3)
58
59 ! In/Output arguments
60 !****************************************************************************************
61 REAL, DIMENSION(klon), INTENT(INOUT) :: radsol
62 REAL, DIMENSION(klon), INTENT(INOUT) :: snow
63 REAL, DIMENSION(klon), INTENT(INOUT) :: agesno !? put to 0 in ocean
64
65 ! Output arguments
66 !****************************************************************************************
67 REAL, DIMENSION(klon), INTENT(OUT) :: qsurf
68 REAL, DIMENSION(klon), INTENT(OUT) :: evap, fluxsens, fluxlat
69 REAL, DIMENSION(klon), INTENT(OUT) :: flux_u1, flux_v1
70 REAL, DIMENSION(klon), INTENT(OUT) :: tsurf_new
71 REAL, DIMENSION(klon), INTENT(OUT) :: dflux_s, dflux_l
72 REAL, intent(out):: sens_prec_liq(:) ! (knon)
73
74 ! Local variables
75 !****************************************************************************************
76 INTEGER :: i, j
77 960 REAL, DIMENSION(klon) :: cal, beta, dif_grnd
78 960 REAL, DIMENSION(klon) :: alb_neig, tsurf_lim, zx_sl
79 960 REAL, DIMENSION(klon) :: u0, v0
80 960 REAL, DIMENSION(klon) :: u1_lay, v1_lay
81 LOGICAL :: check=.FALSE.
82 960 REAL sens_prec_sol(knon)
83 960 REAL, DIMENSION(klon) :: lat_prec_liq, lat_prec_sol
84
85 !****************************************************************************************
86 ! Start calculation
87 !****************************************************************************************
88
1/2
✗ Branch 0 not taken.
✓ Branch 1 taken 480 times.
 480 IF (check) WRITE(*,*)' Entering ocean_forced_noice'
89
90 !****************************************************************************************
91 ! 1)
92 ! Read sea-surface temperature from file limit.nc
93 !
94 !****************************************************************************************
95 !--sb:
96 !!jyg if (knon.eq.1) then ! single-column model
97
1/2
✗ Branch 0 not taken.
✓ Branch 1 taken 480 times.
 480 if (klon_glo.eq.1) then ! single-column model
98 ! EV: now surface Tin flux_arp.h
99 !CALL read_tsurf1d(knon,tsurf_lim) ! new
100 ✗ DO i = 1, knon
101 ✗ tsurf_lim(i) = tg
102 ENDDO
103
104 else ! GCM
105 480 CALL limit_read_sst(knon,knindex,tsurf_lim)
106 endif ! knon
107 !sb--
108
109 !****************************************************************************************
110 ! 2)
111 ! Flux calculation
112 !
113 !****************************************************************************************
114 ! Set some variables for calcul_fluxs
115 !cal = 0.
116 !beta = 1.
117 !dif_grnd = 0.
118
119
120 ! EV: use calbeta to calculate beta
121 ! Need to initialize qsurf for calbeta but it is not modified by this routine
122
2/2
✓ Branch 0 taken 477120 times.
✓ Branch 1 taken 480 times.
 477600 qsurf(:)=0.
123 480 CALL calbeta(dtime, is_oce, knon, snow, qsurf, beta, cal, dif_grnd)
124
125
126
2/2
✓ Branch 0 taken 480 times.
✓ Branch 1 taken 477120 times.
 477600 alb_neig(:) = 0.
127
2/2
✓ Branch 0 taken 480 times.
✓ Branch 1 taken 477120 times.
 477600 agesno(:) = 0.
128
4/4
✓ Branch 0 taken 480 times.
✓ Branch 1 taken 477120 times.
✓ Branch 2 taken 480 times.
✓ Branch 3 taken 477120 times.
 955200 lat_prec_liq = 0.; lat_prec_sol = 0.
129
130 ! Suppose zero surface speed
131
2/2
✓ Branch 0 taken 480 times.
✓ Branch 1 taken 477120 times.
 477600 u0(:)=0.0
132
2/2
✓ Branch 0 taken 480 times.
✓ Branch 1 taken 477120 times.
 477600 v0(:)=0.0
133
2/2
✓ Branch 0 taken 480 times.
✓ Branch 1 taken 477120 times.
 477600 u1_lay(:) = u1(:) - u0(:)
134
2/2
✓ Branch 0 taken 477120 times.
✓ Branch 1 taken 480 times.
 477600 v1_lay(:) = v1(:) - v0(:)
135
136 ! Calcul de tsurf_new, evap, fluxlat, fluxsens, dflux_s, dflux_l and qsurf
137 CALL calcul_fluxs(knon, is_oce, dtime, &
138 merge(tsurf_in, tsurf_lim, activate_ocean_skin == 2), p1lay, cal, &
139 beta, cdragh, cdragq, ps, &
140 precip_rain, precip_snow, snow, qsurf, &
141 radsol, dif_grnd, temp_air, spechum, u1_lay, v1_lay, gustiness, &
142 f_qsat_oce,AcoefH, AcoefQ, BcoefH, BcoefQ, &
143 tsurf_new, evap, fluxlat, fluxsens, dflux_s, dflux_l, &
144
5/8
✓ Branch 0 taken 480 times.
✗ Branch 1 not taken.
✓ Branch 2 taken 477120 times.
✓ Branch 3 taken 480 times.
✗ Branch 4 not taken.
✓ Branch 5 taken 477120 times.
✗ Branch 8 not taken.
✓ Branch 9 taken 480 times.
 477600 sens_prec_liq, sens_prec_sol, lat_prec_liq, lat_prec_sol, rhoa)
145
1/4
✗ Branch 0 not taken.
✓ Branch 1 taken 480 times.
✗ Branch 2 not taken.
✗ Branch 3 not taken.
 480 if (activate_ocean_skin == 2) tsurf_new = tsurf_lim
146
147
2/2
✓ Branch 0 taken 362977 times.
✓ Branch 1 taken 480 times.
 363457 do j = 1, knon
148 362977 i = knindex(j)
149 362977 sens_prec_liq_o(i,1) = sens_prec_liq(j)
150 362977 sens_prec_sol_o(i,1) = sens_prec_sol(j)
151 362977 lat_prec_liq_o(i,1) = lat_prec_liq(j)
152 363457 lat_prec_sol_o(i,1) = lat_prec_sol(j)
153 enddo
154
155
156 ! - Flux calculation at first modele level for U and V
157 CALL calcul_flux_wind(knon, dtime, &
158 u0, v0, u1, v1, gustiness, cdragm, &
159 AcoefU, AcoefV, BcoefU, BcoefV, &
160 p1lay, temp_air, &
161 480 flux_u1, flux_v1)
162
163 480 END SUBROUTINE ocean_forced_noice
164 !
165 !***************************************************************************************
166 !
167 480 SUBROUTINE ocean_forced_ice( &
168 itime, dtime, jour, knon, knindex, &
169 480 tsurf_in, p1lay, cdragh, cdragm, precip_rain, precip_snow, temp_air, spechum, &
170 AcoefH, AcoefQ, BcoefH, BcoefQ, &
171 AcoefU, AcoefV, BcoefU, BcoefV, &
172 ps, u1, v1, gustiness, &
173 radsol, snow, qsol, agesno, tsoil, &
174 qsurf, alb1_new, alb2_new, evap, fluxsens, fluxlat, flux_u1, flux_v1, &
175
1/2
✗ Branch 0 not taken.
✓ Branch 1 taken 480 times.
 480 tsurf_new, dflux_s, dflux_l, rhoa)
176 !
177 ! This subroutine treats the ocean where there is ice.
178 ! The routine reads data from climatologie file and does flux calculations at the
179 ! surface.
180 !
181 ✗ USE dimphy
182 USE geometry_mod, ONLY: longitude,latitude
183 USE calcul_fluxs_mod
184 USE surface_data, ONLY : calice, calsno
185 USE limit_read_mod
186 USE fonte_neige_mod, ONLY : fonte_neige
187 USE indice_sol_mod
188 USE phys_output_var_mod, ONLY : sens_prec_liq_o, sens_prec_sol_o, lat_prec_liq_o, lat_prec_sol_o
189
190 ! INCLUDE "indicesol.h"
191 INCLUDE "dimsoil.h"
192 INCLUDE "YOMCST.h"
193 INCLUDE "clesphys.h"
194 INCLUDE "flux_arp.h"
195
196 ! Input arguments
197 !****************************************************************************************
198 INTEGER, INTENT(IN) :: itime, jour, knon
199 INTEGER, DIMENSION(klon), INTENT(IN) :: knindex
200 REAL, INTENT(IN) :: dtime
201 REAL, DIMENSION(klon), INTENT(IN) :: tsurf_in
202 REAL, DIMENSION(klon), INTENT(IN) :: p1lay
203 REAL, DIMENSION(klon), INTENT(IN) :: cdragh, cdragm
204 REAL, DIMENSION(klon), INTENT(IN) :: precip_rain, precip_snow
205 REAL, DIMENSION(klon), INTENT(IN) :: temp_air, spechum
206 REAL, DIMENSION(klon), INTENT(IN) :: AcoefH, AcoefQ, BcoefH, BcoefQ
207 REAL, DIMENSION(klon), INTENT(IN) :: AcoefU, AcoefV, BcoefU, BcoefV
208 REAL, DIMENSION(klon), INTENT(IN) :: ps
209 REAL, DIMENSION(klon), INTENT(IN) :: u1, v1, gustiness
210 real, intent(in):: rhoa(:) ! (knon) density of moist air (kg / m3)
211
212 ! In/Output arguments
213 !****************************************************************************************
214 REAL, DIMENSION(klon), INTENT(INOUT) :: radsol
215 REAL, DIMENSION(klon), INTENT(INOUT) :: snow, qsol
216 REAL, DIMENSION(klon), INTENT(INOUT) :: agesno
217 REAL, DIMENSION(klon, nsoilmx), INTENT(INOUT) :: tsoil
218
219 ! Output arguments
220 !****************************************************************************************
221 REAL, DIMENSION(klon), INTENT(OUT) :: qsurf
222 REAL, DIMENSION(klon), INTENT(OUT) :: alb1_new ! new albedo in visible SW interval
223 REAL, DIMENSION(klon), INTENT(OUT) :: alb2_new ! new albedo in near IR interval
224 REAL, DIMENSION(klon), INTENT(OUT) :: evap, fluxsens, fluxlat
225 REAL, DIMENSION(klon), INTENT(OUT) :: flux_u1, flux_v1
226 REAL, DIMENSION(klon), INTENT(OUT) :: tsurf_new
227 REAL, DIMENSION(klon), INTENT(OUT) :: dflux_s, dflux_l
228
229 ! Local variables
230 !****************************************************************************************
231 LOGICAL :: check=.FALSE.
232 INTEGER :: i, j
233 REAL :: zfra
234 REAL, PARAMETER :: t_grnd=271.35
235 960 REAL, DIMENSION(klon) :: cal, beta, dif_grnd, capsol
236 960 REAL, DIMENSION(klon) :: alb_neig, tsurf_tmp
237 960 REAL, DIMENSION(klon) :: soilcap, soilflux
238 960 REAL, DIMENSION(klon) :: u0, v0
239 960 REAL, DIMENSION(klon) :: u1_lay, v1_lay
240 960 REAL sens_prec_liq(knon), sens_prec_sol (knon)
241 960 REAL, DIMENSION(klon) :: lat_prec_liq, lat_prec_sol
242
243
244 !****************************************************************************************
245 ! Start calculation
246 !****************************************************************************************
247
1/2
✗ Branch 0 not taken.
✓ Branch 1 taken 480 times.
 480 IF (check) WRITE(*,*)'Entering surface_seaice, knon=',knon
248
249 !****************************************************************************************
250 ! 1)
251 ! Flux calculation : tsurf_new, evap, fluxlat, fluxsens, flux_u1, flux_v1
252 ! dflux_s, dflux_l and qsurf
253 !****************************************************************************************
254
255
2/2
✓ Branch 0 taken 477120 times.
✓ Branch 1 taken 480 times.
 477600 tsurf_tmp(:) = tsurf_in(:)
256
257 ! calculate the parameters cal, beta, capsol and dif_grnd and then recalculate cal
258 480 CALL calbeta(dtime, is_sic, knon, snow, qsol, beta, capsol, dif_grnd)
259
260
261
1/2
✓ Branch 0 taken 480 times.
✗ Branch 1 not taken.
 480 IF (soil_model) THEN
262 ! update tsoil and calculate soilcap and soilflux
263 CALL soil(dtime, is_sic, knon, snow, tsurf_tmp, qsol, &
264
5/6
✓ Branch 0 taken 480 times.
✗ Branch 1 not taken.
✓ Branch 2 taken 104835 times.
✓ Branch 3 taken 480 times.
✓ Branch 4 taken 104835 times.
✓ Branch 5 taken 480 times.
 210150 & longitude(knindex(1:knon)), latitude(knindex(1:knon)), tsoil,soilcap, soilflux)
265
2/2
✓ Branch 0 taken 480 times.
✓ Branch 1 taken 104835 times.
 105315 cal(1:knon) = RCPD / soilcap(1:knon)
266
2/2
✓ Branch 0 taken 104835 times.
✓ Branch 1 taken 480 times.
 105315 radsol(1:knon) = radsol(1:knon) + soilflux(1:knon)
267
2/2
✓ Branch 0 taken 477120 times.
✓ Branch 1 taken 480 times.
 477600 dif_grnd = 1.0 / tau_gl
268 ELSE
269 ✗ dif_grnd = 1.0 / tau_gl
270 ✗ cal = RCPD * calice
271 ✗ WHERE (snow > 0.0) cal = RCPD * calsno
272 ENDIF
273
274 ! beta = 1.0
275
4/4
✓ Branch 0 taken 480 times.
✓ Branch 1 taken 477120 times.
✓ Branch 2 taken 480 times.
✓ Branch 3 taken 477120 times.
 955200 lat_prec_liq = 0.; lat_prec_sol = 0.
276
277 ! Suppose zero surface speed
278
2/2
✓ Branch 0 taken 480 times.
✓ Branch 1 taken 477120 times.
 477600 u0(:)=0.0
279
2/2
✓ Branch 0 taken 480 times.
✓ Branch 1 taken 477120 times.
 477600 v0(:)=0.0
280
2/2
✓ Branch 0 taken 480 times.
✓ Branch 1 taken 477120 times.
 477600 u1_lay(:) = u1(:) - u0(:)
281
2/2
✓ Branch 0 taken 477120 times.
✓ Branch 1 taken 480 times.
 477600 v1_lay(:) = v1(:) - v0(:)
282 CALL calcul_fluxs(knon, is_sic, dtime, &
283 tsurf_tmp, p1lay, cal, beta, cdragh, cdragh, ps, &
284 precip_rain, precip_snow, snow, qsurf, &
285 radsol, dif_grnd, temp_air, spechum, u1_lay, v1_lay, gustiness, &
286 f_qsat_oce,AcoefH, AcoefQ, BcoefH, BcoefQ, &
287 tsurf_new, evap, fluxlat, fluxsens, dflux_s, dflux_l, &
288 480 sens_prec_liq, sens_prec_sol, lat_prec_liq, lat_prec_sol, rhoa)
289
2/2
✓ Branch 0 taken 104835 times.
✓ Branch 1 taken 480 times.
 105315 do j = 1, knon
290 104835 i = knindex(j)
291 104835 sens_prec_liq_o(i,2) = sens_prec_liq(j)
292 104835 sens_prec_sol_o(i,2) = sens_prec_sol(j)
293 104835 lat_prec_liq_o(i,2) = lat_prec_liq(j)
294 105315 lat_prec_sol_o(i,2) = lat_prec_sol(j)
295 enddo
296
297 ! - Flux calculation at first modele level for U and V
298 CALL calcul_flux_wind(knon, dtime, &
299 u0, v0, u1, v1, gustiness, cdragm, &
300 AcoefU, AcoefV, BcoefU, BcoefV, &
301 p1lay, temp_air, &
302 480 flux_u1, flux_v1)
303
304 !****************************************************************************************
305 ! 2)
306 ! Calculations due to snow and runoff
307 !
308 !****************************************************************************************
309 CALL fonte_neige( knon, is_sic, knindex, dtime, &
310 tsurf_tmp, precip_rain, precip_snow, &
311 480 snow, qsol, tsurf_new, evap)
312
313 ! Calculation of albedo at snow (alb_neig) and update the age of snow (agesno)
314 !
315 480 CALL albsno(klon, knon, dtime, agesno(:), alb_neig(:), precip_snow(:))
316
317
4/4
✓ Branch 0 taken 480 times.
✓ Branch 1 taken 104835 times.
✓ Branch 2 taken 29152 times.
✓ Branch 3 taken 75683 times.
 105315 WHERE (snow(1:knon) .LT. 0.0001) agesno(1:knon) = 0.
318
319
2/2
✓ Branch 0 taken 480 times.
✓ Branch 1 taken 477120 times.
 477600 alb1_new(:) = 0.0
320
2/2
✓ Branch 0 taken 480 times.
✓ Branch 1 taken 104835 times.
 105315 DO i=1, knon
321 104835 zfra = MAX(0.0,MIN(1.0,snow(i)/(snow(i)+10.0)))
322 105315 alb1_new(i) = alb_neig(i) * zfra + 0.6 * (1.0-zfra)
323 ENDDO
324
325
2/2
✓ Branch 0 taken 477120 times.
✓ Branch 1 taken 480 times.
 477600 alb2_new(:) = alb1_new(:)
326
327 480 END SUBROUTINE ocean_forced_ice
328
329 !************************************************************************
330 ! 1D case
331 !************************************************************************
332 ! SUBROUTINE read_tsurf1d(knon,sst_out)
333 !
334 ! This subroutine specifies the surface temperature to be used in 1D simulations
335 !
336 ! USE dimphy, ONLY : klon
337 !
338 ! INTEGER, INTENT(IN) :: knon ! nomber of points on compressed grid
339 ! REAL, DIMENSION(klon), INTENT(OUT) :: sst_out ! tsurf used to force the single-column model
340 !
341 ! INTEGER :: i
342 ! COMMON defined in lmdz1d.F:
343 ! real ts_cur
344 ! common /sst_forcing/ts_cur
345 !
346 ! DO i = 1, knon
347 ! sst_out(i) = ts_cur
348 ! ENDDO
349 !
350 ! END SUBROUTINE read_tsurf1d
351 !
352 !
353 !************************************************************************
354 END MODULE ocean_forced_mod
355
356
357
358
359
360
361