Directory: | ./ |
---|---|
File: | phys/surf_ocean_mod.f90 |
Date: | 2022-01-11 19:19:34 |
Exec | Total | Coverage | |
---|---|---|---|
Lines: | 38 | 74 | 51.4% |
Branches: | 48 | 140 | 34.3% |
Line | Branch | Exec | Source |
---|---|---|---|
1 | ! | ||
2 | ! $Id: surf_ocean_mod.F90 3906 2021-05-19 10:35:18Z jyg $ | ||
3 | ! | ||
4 | MODULE surf_ocean_mod | ||
5 | |||
6 | IMPLICIT NONE | ||
7 | |||
8 | CONTAINS | ||
9 | ! | ||
10 | !****************************************************************************** | ||
11 | ! | ||
12 | 365377 | SUBROUTINE surf_ocean(rlon, rlat, swnet, lwnet, alb1, & | |
13 | windsp, rmu0, fder, tsurf_in, & | ||
14 | itime, dtime, jour, knon, knindex, & | ||
15 | p1lay, z1lay, cdragh, cdragm, precip_rain, precip_snow, temp_air, spechum, & | ||
16 | AcoefH, AcoefQ, BcoefH, BcoefQ, & | ||
17 | AcoefU, AcoefV, BcoefU, BcoefV, & | ||
18 | 480 | ps, u1, v1, gustiness, rugoro, pctsrf, & | |
19 | snow, qsurf, agesno, & | ||
20 | 480 | z0m, z0h, SFRWL, alb_dir_new, alb_dif_new, evap, fluxsens, fluxlat, & | |
21 | tsurf_new, dflux_s, dflux_l, lmt_bils, & | ||
22 |
6/12✗ Branch 0 not taken.
✓ Branch 1 taken 480 times.
✗ Branch 2 not taken.
✓ Branch 3 taken 480 times.
✗ Branch 4 not taken.
✓ Branch 5 taken 480 times.
✗ Branch 6 not taken.
✓ Branch 7 taken 480 times.
✗ Branch 8 not taken.
✓ Branch 9 taken 480 times.
✗ Branch 10 not taken.
✓ Branch 11 taken 480 times.
|
960 | flux_u1, flux_v1, delta_sst, delta_sal, ds_ns, dt_ns, dter, dser, tkt, tks, & |
23 |
4/8✗ Branch 0 not taken.
✓ Branch 1 taken 480 times.
✗ Branch 2 not taken.
✓ Branch 3 taken 480 times.
✗ Branch 4 not taken.
✓ Branch 5 taken 480 times.
✗ Branch 6 not taken.
✓ Branch 7 taken 480 times.
|
960 | taur, sss) |
24 | |||
25 | use albedo, only: alboc, alboc_cd | ||
26 | use bulk_flux_m, only: bulk_flux | ||
27 | USE dimphy, ONLY: klon, zmasq | ||
28 | USE surface_data, ONLY : type_ocean | ||
29 | USE ocean_forced_mod, ONLY : ocean_forced_noice | ||
30 | USE ocean_slab_mod, ONLY : ocean_slab_noice | ||
31 | USE ocean_cpl_mod, ONLY : ocean_cpl_noice | ||
32 | USE indice_sol_mod, ONLY : nbsrf, is_oce | ||
33 | USE limit_read_mod | ||
34 | use config_ocean_skin_m, only: activate_ocean_skin | ||
35 | ! | ||
36 | ! This subroutine will make a call to ocean_XXX_noice according to the ocean mode (force, | ||
37 | ! slab or couple). The calculations of albedo and rugosity for the ocean surface are | ||
38 | ! done in here because they are identical for the different modes of ocean. | ||
39 | |||
40 | |||
41 | INCLUDE "YOMCST.h" | ||
42 | |||
43 | include "clesphys.h" | ||
44 | ! for cycle_diurne and for iflag_z0_oce==-1 (prescribed z0) | ||
45 | |||
46 | ! Input variables | ||
47 | !****************************************************************************** | ||
48 | INTEGER, INTENT(IN) :: itime, jour, knon | ||
49 | INTEGER, DIMENSION(klon), INTENT(IN) :: knindex | ||
50 | REAL, INTENT(IN) :: dtime | ||
51 | REAL, DIMENSION(klon), INTENT(IN) :: rlon, rlat | ||
52 | REAL, DIMENSION(klon), INTENT(IN) :: swnet ! net shortwave radiation at surface | ||
53 | REAL, DIMENSION(klon), INTENT(IN) :: lwnet ! net longwave radiation at surface | ||
54 | REAL, DIMENSION(klon), INTENT(IN) :: alb1 ! albedo in visible SW interval | ||
55 | REAL, DIMENSION(klon), INTENT(IN) :: windsp ! wind at 10 m, in m s-1 | ||
56 | REAL, DIMENSION(klon), INTENT(IN) :: rmu0 | ||
57 | REAL, DIMENSION(klon), INTENT(IN) :: fder | ||
58 | REAL, DIMENSION(klon), INTENT(IN) :: tsurf_in ! defined only for subscripts 1:knon | ||
59 | REAL, DIMENSION(klon), INTENT(IN) :: p1lay,z1lay ! pression (Pa) et altitude (m) du premier niveau | ||
60 | REAL, DIMENSION(klon), INTENT(IN) :: cdragh | ||
61 | REAL, DIMENSION(klon), INTENT(IN) :: cdragm | ||
62 | REAL, DIMENSION(klon), INTENT(IN) :: precip_rain, precip_snow | ||
63 | REAL, DIMENSION(klon), INTENT(IN) :: temp_air, spechum | ||
64 | REAL, DIMENSION(klon), INTENT(IN) :: AcoefH, AcoefQ, BcoefH, BcoefQ | ||
65 | REAL, DIMENSION(klon), INTENT(IN) :: AcoefU, AcoefV, BcoefU, BcoefV | ||
66 | REAL, DIMENSION(klon), INTENT(IN) :: ps | ||
67 | REAL, DIMENSION(klon), INTENT(IN) :: u1, v1, gustiness | ||
68 | REAL, DIMENSION(klon), INTENT(IN) :: rugoro | ||
69 | REAL, DIMENSION(klon,nbsrf), INTENT(IN) :: pctsrf | ||
70 | |||
71 | ! In/Output variables | ||
72 | !****************************************************************************** | ||
73 | REAL, DIMENSION(klon), INTENT(INOUT) :: snow | ||
74 | REAL, DIMENSION(klon), INTENT(INOUT) :: qsurf | ||
75 | REAL, DIMENSION(klon), INTENT(INOUT) :: agesno | ||
76 | REAL, DIMENSION(klon), INTENT(inOUT) :: z0h | ||
77 | |||
78 | REAL, intent(inout):: delta_sst(:) ! (knon) | ||
79 | ! Ocean-air interface temperature minus bulk SST, in K. Defined | ||
80 | ! only if activate_ocean_skin >= 1. | ||
81 | |||
82 | real, intent(inout):: delta_sal(:) ! (knon) | ||
83 | ! Ocean-air interface salinity minus bulk salinity, in ppt. Defined | ||
84 | ! only if activate_ocean_skin >= 1. | ||
85 | |||
86 | REAL, intent(inout):: ds_ns(:) ! (knon) | ||
87 | ! "delta salinity near surface". Salinity variation in the | ||
88 | ! near-surface turbulent layer. That is subskin salinity minus | ||
89 | ! foundation salinity. In ppt. | ||
90 | |||
91 | REAL, intent(inout):: dt_ns(:) ! (knon) | ||
92 | ! "delta temperature near surface". Temperature variation in the | ||
93 | ! near-surface turbulent layer. That is subskin temperature | ||
94 | ! minus foundation temperature. (Can be negative.) In K. | ||
95 | |||
96 | ! Output variables | ||
97 | !************************************************************************** | ||
98 | REAL, DIMENSION(klon), INTENT(OUT) :: z0m | ||
99 | !albedo SB >>> | ||
100 | ! REAL, DIMENSION(klon), INTENT(OUT) :: alb1_new ! new albedo in visible SW interval | ||
101 | ! REAL, DIMENSION(klon), INTENT(OUT) :: alb2_new ! new albedo in near IR interval | ||
102 | REAL, DIMENSION(6), INTENT(IN) :: SFRWL | ||
103 | REAL, DIMENSION(klon,nsw), INTENT(OUT) :: alb_dir_new,alb_dif_new | ||
104 | !albedo SB <<< | ||
105 | REAL, DIMENSION(klon), INTENT(OUT) :: evap, fluxsens, fluxlat | ||
106 | REAL, DIMENSION(klon), INTENT(OUT) :: tsurf_new ! sea surface temperature, in K | ||
107 | REAL, DIMENSION(klon), INTENT(OUT) :: dflux_s, dflux_l | ||
108 | REAL, DIMENSION(klon), INTENT(OUT) :: lmt_bils | ||
109 | REAL, DIMENSION(klon), INTENT(OUT) :: flux_u1, flux_v1 | ||
110 | |||
111 | REAL, intent(out):: dter(:) ! (knon) | ||
112 | ! Temperature variation in the diffusive microlayer, that is | ||
113 | ! ocean-air interface temperature minus subskin temperature. In | ||
114 | ! K. | ||
115 | |||
116 | REAL, intent(out):: dser(:) ! (knon) | ||
117 | ! Salinity variation in the diffusive microlayer, that is | ||
118 | ! ocean-air interface salinity minus subskin salinity. In ppt. | ||
119 | |||
120 | REAL, intent(out):: tkt(:) ! (knon) | ||
121 | ! �paisseur (m) de la couche de diffusion thermique (microlayer) | ||
122 | ! cool skin thickness | ||
123 | |||
124 | REAL, intent(out):: tks(:) ! (knon) | ||
125 | ! �paisseur (m) de la couche de diffusion de masse (microlayer) | ||
126 | |||
127 | REAL, intent(out):: taur(:) ! (knon) | ||
128 | ! momentum flux due to rain, in Pa | ||
129 | |||
130 | real, intent(out):: sss(:) ! (klon) | ||
131 | ! Bulk salinity of the surface layer of the ocean, in ppt. (Only | ||
132 | ! defined for subscripts 1:knon, but we have to declare it with | ||
133 | ! size klon because of the coupling machinery.) | ||
134 | |||
135 | ! Local variables | ||
136 | !************************************************************************* | ||
137 | INTEGER :: i, k | ||
138 | REAL :: tmp | ||
139 | REAL, PARAMETER :: cepdu2=(0.1)**2 | ||
140 | 960 | REAL, DIMENSION(klon) :: alb_eau, z0_lim | |
141 | 960 | REAL, DIMENSION(klon) :: radsol | |
142 | 960 | REAL, DIMENSION(klon) :: cdragq ! Cdrag pour l'evaporation | |
143 | CHARACTER(len=20),PARAMETER :: modname="surf_ocean" | ||
144 | 960 | real rhoa(knon) ! density of moist air (kg / m3) | |
145 | 960 | REAL sens_prec_liq(knon) | |
146 | |||
147 | 960 | REAL t_int(knon) ! ocean-air interface temperature, in K | |
148 | 480 | real s_int(knon) ! ocean-air interface salinity, in ppt | |
149 | |||
150 | !************************************************************************** | ||
151 | |||
152 | |||
153 | !****************************************************************************** | ||
154 | ! Calculate total net radiance at surface | ||
155 | ! | ||
156 | !****************************************************************************** | ||
157 |
2/2✓ Branch 0 taken 477120 times.
✓ Branch 1 taken 480 times.
|
477600 | radsol(1:klon) = 0.0 ! initialisation a priori inutile |
158 |
2/2✓ Branch 0 taken 362977 times.
✓ Branch 1 taken 480 times.
|
363457 | radsol(1:knon) = swnet(1:knon) + lwnet(1:knon) |
159 | |||
160 | !****************************************************************************** | ||
161 | ! Cdragq computed from cdrag | ||
162 | ! The difference comes only from a factor (f_z0qh_oce) on z0, so that | ||
163 | ! it can be computed inside surf_ocean | ||
164 | ! More complicated appraches may require the propagation through | ||
165 | ! pbl_surface of an independant cdragq variable. | ||
166 | !****************************************************************************** | ||
167 | |||
168 |
1/2✓ Branch 0 taken 480 times.
✗ Branch 1 not taken.
|
480 | IF ( f_z0qh_oce .ne. 1.) THEN |
169 | ! Si on suit les formulations par exemple de Tessel, on | ||
170 | ! a z0h=0.4*nu/u*, z0q=0.62*nu/u*, d'ou f_z0qh_oce=0.62/0.4=1.55 | ||
171 | cdragq(1:knon)=cdragh(1:knon)* & | ||
172 |
2/2✓ Branch 0 taken 362977 times.
✓ Branch 1 taken 480 times.
|
363457 | log(z1lay(1:knon)/z0h(1:knon))/log(z1lay(1:knon)/(f_z0qh_oce*z0h(1:knon))) |
173 | ELSE | ||
174 | ✗ | cdragq(1:knon)=cdragh(1:knon) | |
175 | ENDIF | ||
176 | |||
177 |
2/2✓ Branch 0 taken 362977 times.
✓ Branch 1 taken 480 times.
|
363457 | rhoa = PS(:KNON) / (Rd * temp_air(:knon) * (1. + retv * spechum(:knon))) |
178 | !****************************************************************************** | ||
179 | ! Switch according to type of ocean (couple, slab or forced) | ||
180 | !****************************************************************************** | ||
181 | ✗ | SELECT CASE(type_ocean) | |
182 | CASE('couple') | ||
183 | CALL ocean_cpl_noice( & | ||
184 | swnet, lwnet, alb1, & | ||
185 | windsp, fder, & | ||
186 | itime, dtime, knon, knindex, & | ||
187 | p1lay, cdragh, cdragq, cdragm, precip_rain, precip_snow,temp_air,spechum,& | ||
188 | AcoefH, AcoefQ, BcoefH, BcoefQ, & | ||
189 | AcoefU, AcoefV, BcoefU, BcoefV, & | ||
190 | ps, u1, v1, gustiness, tsurf_in, & | ||
191 | radsol, snow, agesno, & | ||
192 | qsurf, evap, fluxsens, fluxlat, flux_u1, flux_v1, & | ||
193 | tsurf_new, dflux_s, dflux_l, sens_prec_liq, sss, delta_sal, rhoa, & | ||
194 | ✗ | delta_sst) | |
195 | |||
196 | CASE('slab') | ||
197 | CALL ocean_slab_noice( & | ||
198 | itime, dtime, jour, knon, knindex, & | ||
199 | p1lay, cdragh, cdragq, cdragm, precip_rain, precip_snow, temp_air, spechum,& | ||
200 | AcoefH, AcoefQ, BcoefH, BcoefQ, & | ||
201 | AcoefU, AcoefV, BcoefU, BcoefV, & | ||
202 | ps, u1, v1, gustiness, tsurf_in, & | ||
203 | radsol, snow, & | ||
204 | qsurf, evap, fluxsens, fluxlat, flux_u1, flux_v1, & | ||
205 | ✗ | tsurf_new, dflux_s, dflux_l, lmt_bils) | |
206 | |||
207 | CASE('force') | ||
208 | CALL ocean_forced_noice( & | ||
209 | itime, dtime, jour, knon, knindex, & | ||
210 | p1lay, cdragh, cdragq, cdragm, precip_rain, precip_snow, & | ||
211 | temp_air, spechum, & | ||
212 | AcoefH, AcoefQ, BcoefH, BcoefQ, & | ||
213 | AcoefU, AcoefV, BcoefU, BcoefV, & | ||
214 | ps, u1, v1, gustiness, tsurf_in, & | ||
215 | radsol, snow, agesno, & | ||
216 | qsurf, evap, fluxsens, fluxlat, flux_u1, flux_v1, & | ||
217 |
1/4✗ Branch 0 not taken.
✗ Branch 1 not taken.
✓ Branch 2 taken 480 times.
✗ Branch 3 not taken.
|
480 | tsurf_new, dflux_s, dflux_l, sens_prec_liq, rhoa) |
218 | END SELECT | ||
219 | |||
220 | !****************************************************************************** | ||
221 | ! fcodron: compute lmt_bils forced case (same as wfbils_oce / 1.-contfracatm) | ||
222 | !****************************************************************************** | ||
223 |
1/2✓ Branch 0 taken 480 times.
✗ Branch 1 not taken.
|
480 | IF (type_ocean.NE.'slab') THEN |
224 |
2/2✓ Branch 0 taken 477120 times.
✓ Branch 1 taken 480 times.
|
477600 | lmt_bils(1:klon)=0. |
225 |
2/2✓ Branch 0 taken 362977 times.
✓ Branch 1 taken 480 times.
|
363457 | DO i=1,knon |
226 | lmt_bils(knindex(i))=(swnet(i)+lwnet(i)+fluxsens(i)+fluxlat(i)) & | ||
227 | 363457 | *pctsrf(knindex(i),is_oce)/(1.-zmasq(knindex(i))) | |
228 | END DO | ||
229 | END IF | ||
230 | |||
231 | !****************************************************************************** | ||
232 | ! Calculate ocean surface albedo | ||
233 | !****************************************************************************** | ||
234 | !albedo SB >>> | ||
235 |
1/2✗ Branch 0 not taken.
✓ Branch 1 taken 480 times.
|
480 | IF (iflag_albedo==0) THEN |
236 | !--old parametrizations of ocean surface albedo | ||
237 | ! | ||
238 | ✗ | IF (iflag_cycle_diurne.GE.1) THEN | |
239 | ! | ||
240 | ✗ | CALL alboc_cd(rmu0,alb_eau) | |
241 | ! | ||
242 | !--ad-hoc correction for model radiative balance tuning | ||
243 | !--now outside alboc_cd routine | ||
244 | ✗ | alb_eau(1:klon) = fmagic*alb_eau(1:klon) + pmagic | |
245 | ✗ | alb_eau(1:klon)=MIN(MAX(alb_eau(1:klon),0.0),1.0) | |
246 | ! | ||
247 | ELSE | ||
248 | ! | ||
249 | ✗ | CALL alboc(REAL(jour),rlat,alb_eau) | |
250 | !--ad-hoc correction for model radiative balance tuning | ||
251 | !--now outside alboc routine | ||
252 | ✗ | alb_eau(1:klon) = fmagic*alb_eau(1:klon) + pmagic | |
253 | ✗ | alb_eau(1:klon)=MIN(MAX(alb_eau(1:klon),0.04),0.60) | |
254 | ! | ||
255 | ENDIF | ||
256 | ! | ||
257 | ✗ | DO i =1, knon | |
258 | ✗ | DO k=1,nsw | |
259 | ✗ | alb_dir_new(i,k) = alb_eau(knindex(i)) | |
260 | ENDDO | ||
261 | ENDDO | ||
262 | !IM 09122015 next line corresponds to the old way of doing in LMDZ5A/IPSLCM5A versions | ||
263 | !albedo for diffuse radiation is taken the same as for direct radiation | ||
264 | ✗ | alb_dif_new(1:knon,:)=alb_dir_new(1:knon,:) | |
265 | !IM 09122015 end | ||
266 | ! | ||
267 |
1/2✓ Branch 0 taken 480 times.
✗ Branch 1 not taken.
|
480 | ELSE IF (iflag_albedo==1) THEN |
268 | !--new parametrization of ocean surface albedo by Sunghye Baek | ||
269 | !--albedo for direct and diffuse radiation are different | ||
270 | ! | ||
271 | 480 | CALL ocean_albedo(knon,rmu0,knindex,windsp,SFRWL,alb_dir_new,alb_dif_new) | |
272 | ! | ||
273 | !--ad-hoc correction for model radiative balance tuning | ||
274 |
4/4✓ Branch 0 taken 2880 times.
✓ Branch 1 taken 480 times.
✓ Branch 2 taken 2177862 times.
✓ Branch 3 taken 2880 times.
|
2181222 | alb_dir_new(1:knon,:) = fmagic*alb_dir_new(1:knon,:) + pmagic |
275 |
4/4✓ Branch 0 taken 2880 times.
✓ Branch 1 taken 480 times.
✓ Branch 2 taken 2177862 times.
✓ Branch 3 taken 2880 times.
|
2181222 | alb_dif_new(1:knon,:) = fmagic*alb_dif_new(1:knon,:) + pmagic |
276 |
4/4✓ Branch 0 taken 2880 times.
✓ Branch 1 taken 480 times.
✓ Branch 2 taken 2177862 times.
✓ Branch 3 taken 2880 times.
|
2181222 | alb_dir_new(1:knon,:)=MIN(MAX(alb_dir_new(1:knon,:),0.0),1.0) |
277 |
4/4✓ Branch 0 taken 2880 times.
✓ Branch 1 taken 480 times.
✓ Branch 2 taken 2177862 times.
✓ Branch 3 taken 2880 times.
|
2181222 | alb_dif_new(1:knon,:)=MIN(MAX(alb_dif_new(1:knon,:),0.0),1.0) |
278 | ! | ||
279 | ✗ | ELSE IF (iflag_albedo==2) THEN | |
280 | ! F. Codron albedo read from limit.nc | ||
281 | CALL limit_read_rug_alb(itime, dtime, jour,& | ||
282 | ✗ | knon, knindex, z0_lim, alb_eau) | |
283 | ✗ | DO i =1, knon | |
284 | ✗ | DO k=1,nsw | |
285 | ✗ | alb_dir_new(i,k) = alb_eau(i) | |
286 | ENDDO | ||
287 | ENDDO | ||
288 | ✗ | alb_dif_new=alb_dir_new | |
289 | ENDIF | ||
290 | !albedo SB <<< | ||
291 | |||
292 | !****************************************************************************** | ||
293 | ! Calculate the rugosity | ||
294 | !****************************************************************************** | ||
295 |
1/2✗ Branch 0 not taken.
✓ Branch 1 taken 480 times.
|
480 | IF (iflag_z0_oce==0) THEN |
296 | ✗ | DO i = 1, knon | |
297 | ✗ | tmp = MAX(cepdu2,gustiness(i)+u1(i)**2+v1(i)**2) | |
298 | z0m(i) = 0.018*cdragm(i) * (gustiness(i)+u1(i)**2+v1(i)**2)/RG & | ||
299 | ✗ | + 0.11*14e-6 / SQRT(cdragm(i) * tmp) | |
300 | ✗ | z0m(i) = MAX(1.5e-05,z0m(i)) | |
301 | ENDDO | ||
302 | ✗ | z0h(1:knon)=z0m(1:knon) ! En attendant mieux | |
303 | |||
304 |
1/2✓ Branch 0 taken 480 times.
✗ Branch 1 not taken.
|
480 | ELSE IF (iflag_z0_oce==1) THEN |
305 |
2/2✓ Branch 0 taken 480 times.
✓ Branch 1 taken 362977 times.
|
363457 | DO i = 1, knon |
306 | 362977 | tmp = MAX(cepdu2,gustiness(i)+u1(i)**2+v1(i)**2) | |
307 | z0m(i) = 0.018*cdragm(i) * (gustiness(i)+u1(i)**2+v1(i)**2)/RG & | ||
308 | 362977 | + 0.11*14e-6 / SQRT(cdragm(i) * tmp) | |
309 | 362977 | z0m(i) = MAX(1.5e-05,z0m(i)) | |
310 | 363457 | z0h(i)=0.4*14e-6 / SQRT(cdragm(i) * tmp) | |
311 | ENDDO | ||
312 | ✗ | ELSE IF (iflag_z0_oce==-1) THEN | |
313 | ✗ | DO i = 1, knon | |
314 | ✗ | z0m(i) = z0min | |
315 | ✗ | z0h(i) = z0min | |
316 | ENDDO | ||
317 | ELSE | ||
318 | ✗ | CALL abort_physic(modname,'version non prevue',1) | |
319 | ENDIF | ||
320 | |||
321 |
1/2✗ Branch 0 not taken.
✓ Branch 1 taken 480 times.
|
480 | if (activate_ocean_skin >= 1) then |
322 | ✗ | if (type_ocean /= 'couple') sss(:knon) = 35. | |
323 | call bulk_flux(tkt, tks, taur, dter, dser, t_int, s_int, ds_ns, dt_ns, & | ||
324 | u = windsp(:knon), t_ocean_1 = tsurf_new(:knon), s1 = sss(:knon), & | ||
325 | rain = precip_rain(:knon) + precip_snow(:knon), & | ||
326 | hf = - fluxsens(:knon), hlb = - fluxlat(:knon), & | ||
327 | rnl = - lwnet(:knon), & | ||
328 | tau = sqrt(flux_u1(:knon)**2 + flux_v1(:knon)**2), rhoa = rhoa, & | ||
329 | xlv = [(rlvtt, i = 1, knon)], rf = - sens_prec_liq, dtime = dtime, & | ||
330 | ✗ | rns = swnet(:knon)) | |
331 | ✗ | delta_sst = t_int - tsurf_new(:knon) | |
332 | ✗ | delta_sal = s_int - sss(:knon) | |
333 | ✗ | if (activate_ocean_skin >= 2) tsurf_new(:knon) = t_int | |
334 | end if | ||
335 | |||
336 | 480 | END SUBROUTINE surf_ocean | |
337 | !**************************************************************************** | ||
338 | ! | ||
339 | END MODULE surf_ocean_mod | ||
340 |