source: LMDZ6/trunk/libf/phylmd/ocean_forced_mod.F90 @ 3947

Last change on this file since 3947 was 3815, checked in by lguez, 4 years ago

Merge Ocean_skin branch back into trunk
  • Property copyright set to
    Name of program: LMDZ
    Creation date: 1984
    Version: LMDZ5
    License: CeCILL version 2
    Holder: Laboratoire de m\'et\'eorologie dynamique, CNRS, UMR 8539
    See the license file in the root directory
  • Property svn:eol-style set to native
  • Property svn:keywords set to Id
File size: 13.4 KB
Line 
1!
2! $Id: ocean_forced_mod.F90 3815 2021-02-01 14:30:57Z emillour $
3!
4MODULE 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
11CONTAINS
12!
13!****************************************************************************************
14!
15  SUBROUTINE ocean_forced_noice( &
16       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       qsurf, evap, fluxsens, fluxlat, flux_u1, flux_v1, &
24       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    REAL, DIMENSION(klon)       :: cal, beta, dif_grnd
78    REAL, DIMENSION(klon)       :: alb_neig, tsurf_lim, zx_sl
79    REAL, DIMENSION(klon)       :: u0, v0
80    REAL, DIMENSION(klon)       :: u1_lay, v1_lay
81    LOGICAL                     :: check=.FALSE.
82    REAL sens_prec_sol(knon)
83    REAL, DIMENSION(klon) :: lat_prec_liq, lat_prec_sol   
84
85!****************************************************************************************
86! Start calculation
87!****************************************************************************************
88    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    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      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    qsurf(:)=0.
123    CALL calbeta(dtime, is_oce, knon, snow, qsurf, beta, cal, dif_grnd)
124
125
126    alb_neig(:) = 0.
127    agesno(:) = 0.
128    lat_prec_liq = 0.; lat_prec_sol = 0.
129
130! Suppose zero surface speed
131    u0(:)=0.0
132    v0(:)=0.0
133    u1_lay(:) = u1(:) - u0(:)
134    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         sens_prec_liq, sens_prec_sol, lat_prec_liq, lat_prec_sol, rhoa)
145    if (activate_ocean_skin == 2) tsurf_new = tsurf_lim
146
147    do j = 1, knon
148      i = knindex(j)
149      sens_prec_liq_o(i,1) = sens_prec_liq(j)
150      sens_prec_sol_o(i,1) = sens_prec_sol(j)
151      lat_prec_liq_o(i,1) = lat_prec_liq(j)
152      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         flux_u1, flux_v1) 
162
163  END SUBROUTINE ocean_forced_noice
164!
165!***************************************************************************************
166!
167  SUBROUTINE ocean_forced_ice( &
168       itime, dtime, jour, knon, knindex, &
169       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       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 calcul_fluxs_mod
183    USE surface_data,     ONLY : calice, calsno
184    USE limit_read_mod
185    USE fonte_neige_mod,  ONLY : fonte_neige
186    USE indice_sol_mod
187    USE phys_output_var_mod, ONLY : sens_prec_liq_o, sens_prec_sol_o, lat_prec_liq_o, lat_prec_sol_o
188
189!   INCLUDE "indicesol.h"
190    INCLUDE "dimsoil.h"
191    INCLUDE "YOMCST.h"
192    INCLUDE "clesphys.h"
193    INCLUDE "flux_arp.h"
194
195! Input arguments
196!****************************************************************************************
197    INTEGER, INTENT(IN)                  :: itime, jour, knon
198    INTEGER, DIMENSION(klon), INTENT(IN) :: knindex
199    REAL, INTENT(IN)                     :: dtime
200    REAL, DIMENSION(klon), INTENT(IN)    :: tsurf_in
201    REAL, DIMENSION(klon), INTENT(IN)    :: p1lay
202    REAL, DIMENSION(klon), INTENT(IN)    :: cdragh, cdragm
203    REAL, DIMENSION(klon), INTENT(IN)    :: precip_rain, precip_snow
204    REAL, DIMENSION(klon), INTENT(IN)    :: temp_air, spechum
205    REAL, DIMENSION(klon), INTENT(IN)    :: AcoefH, AcoefQ, BcoefH, BcoefQ
206    REAL, DIMENSION(klon), INTENT(IN)    :: AcoefU, AcoefV, BcoefU, BcoefV
207    REAL, DIMENSION(klon), INTENT(IN)    :: ps
208    REAL, DIMENSION(klon), INTENT(IN)    :: u1, v1, gustiness
209    real, intent(in):: rhoa(:) ! (knon) density of moist air  (kg / m3)
210
211! In/Output arguments
212!****************************************************************************************
213    REAL, DIMENSION(klon), INTENT(INOUT)          :: radsol
214    REAL, DIMENSION(klon), INTENT(INOUT)          :: snow, qsol
215    REAL, DIMENSION(klon), INTENT(INOUT)          :: agesno
216    REAL, DIMENSION(klon, nsoilmx), INTENT(INOUT) :: tsoil
217
218! Output arguments
219!****************************************************************************************
220    REAL, DIMENSION(klon), INTENT(OUT)            :: qsurf
221    REAL, DIMENSION(klon), INTENT(OUT)            :: alb1_new  ! new albedo in visible SW interval
222    REAL, DIMENSION(klon), INTENT(OUT)            :: alb2_new  ! new albedo in near IR interval
223    REAL, DIMENSION(klon), INTENT(OUT)            :: evap, fluxsens, fluxlat
224    REAL, DIMENSION(klon), INTENT(OUT)            :: flux_u1, flux_v1
225    REAL, DIMENSION(klon), INTENT(OUT)            :: tsurf_new
226    REAL, DIMENSION(klon), INTENT(OUT)            :: dflux_s, dflux_l     
227
228! Local variables
229!****************************************************************************************
230    LOGICAL                     :: check=.FALSE.
231    INTEGER                     :: i, j
232    REAL                        :: zfra
233    REAL, PARAMETER             :: t_grnd=271.35
234    REAL, DIMENSION(klon)       :: cal, beta, dif_grnd, capsol
235    REAL, DIMENSION(klon)       :: alb_neig, tsurf_tmp
236    REAL, DIMENSION(klon)       :: soilcap, soilflux
237    REAL, DIMENSION(klon)       :: u0, v0
238    REAL, DIMENSION(klon)       :: u1_lay, v1_lay
239    REAL sens_prec_liq(knon), sens_prec_sol (knon)   
240    REAL, DIMENSION(klon)       :: lat_prec_liq, lat_prec_sol   
241
242
243!****************************************************************************************
244! Start calculation
245!****************************************************************************************
246    IF (check) WRITE(*,*)'Entering surface_seaice, knon=',knon
247
248!****************************************************************************************
249! 1)
250! Flux calculation : tsurf_new, evap, fluxlat, fluxsens, flux_u1, flux_v1
251!                    dflux_s, dflux_l and qsurf
252!****************************************************************************************
253
254    tsurf_tmp(:) = tsurf_in(:)
255
256! calculate the parameters cal, beta, capsol and dif_grnd and then recalculate cal
257    CALL calbeta(dtime, is_sic, knon, snow, qsol, beta, capsol, dif_grnd)
258
259   
260    IF (soil_model) THEN
261! update tsoil and calculate soilcap and soilflux
262       CALL soil(dtime, is_sic, knon, snow, tsurf_tmp, tsoil,soilcap, soilflux)
263       cal(1:knon) = RCPD / soilcap(1:knon)
264       radsol(1:knon) = radsol(1:knon)  + soilflux(1:knon)
265       dif_grnd = 1.0 / tau_gl
266    ELSE
267       dif_grnd = 1.0 / tau_gl
268       cal = RCPD * calice
269       WHERE (snow > 0.0) cal = RCPD * calsno
270    ENDIF
271
272!    beta = 1.0
273    lat_prec_liq = 0.; lat_prec_sol = 0.
274
275! Suppose zero surface speed
276    u0(:)=0.0
277    v0(:)=0.0
278    u1_lay(:) = u1(:) - u0(:)
279    v1_lay(:) = v1(:) - v0(:)
280    CALL calcul_fluxs(knon, is_sic, dtime, &
281         tsurf_tmp, p1lay, cal, beta, cdragh, cdragh, ps, &
282         precip_rain, precip_snow, snow, qsurf,  &
283         radsol, dif_grnd, temp_air, spechum, u1_lay, v1_lay, gustiness, &
284         f_qsat_oce,AcoefH, AcoefQ, BcoefH, BcoefQ, &
285         tsurf_new, evap, fluxlat, fluxsens, dflux_s, dflux_l, &
286         sens_prec_liq, sens_prec_sol, lat_prec_liq, lat_prec_sol, rhoa)
287    do j = 1, knon
288      i = knindex(j)
289      sens_prec_liq_o(i,2) = sens_prec_liq(j)
290      sens_prec_sol_o(i,2) = sens_prec_sol(j)
291      lat_prec_liq_o(i,2) = lat_prec_liq(j)
292      lat_prec_sol_o(i,2) = lat_prec_sol(j)
293    enddo
294
295! - Flux calculation at first modele level for U and V
296    CALL calcul_flux_wind(knon, dtime, &
297         u0, v0, u1, v1, gustiness, cdragm, &
298         AcoefU, AcoefV, BcoefU, BcoefV, &
299         p1lay, temp_air, &
300         flux_u1, flux_v1) 
301
302!****************************************************************************************
303! 2)
304! Calculations due to snow and runoff
305!
306!****************************************************************************************
307    CALL fonte_neige( knon, is_sic, knindex, dtime, &
308         tsurf_tmp, precip_rain, precip_snow, &
309         snow, qsol, tsurf_new, evap)
310   
311! Calculation of albedo at snow (alb_neig) and update the age of snow (agesno)
312!
313    CALL albsno(klon, knon, dtime, agesno(:), alb_neig(:), precip_snow(:)) 
314
315    WHERE (snow(1:knon) .LT. 0.0001) agesno(1:knon) = 0.
316
317    alb1_new(:) = 0.0
318    DO i=1, knon
319       zfra = MAX(0.0,MIN(1.0,snow(i)/(snow(i)+10.0)))
320       alb1_new(i) = alb_neig(i) * zfra +  0.6 * (1.0-zfra)
321    ENDDO
322
323    alb2_new(:) = alb1_new(:)
324
325  END SUBROUTINE ocean_forced_ice
326
327!************************************************************************
328! 1D case
329!************************************************************************
330!  SUBROUTINE read_tsurf1d(knon,sst_out)
331!
332! This subroutine specifies the surface temperature to be used in 1D simulations
333!
334!      USE dimphy, ONLY : klon
335!
336!      INTEGER, INTENT(IN)                  :: knon     ! nomber of points on compressed grid
337!      REAL, DIMENSION(klon), INTENT(OUT)   :: sst_out  ! tsurf used to force the single-column model
338!
339!       INTEGER :: i
340! COMMON defined in lmdz1d.F:
341!       real ts_cur
342!       common /sst_forcing/ts_cur
343!
344!       DO i = 1, knon
345!        sst_out(i) = ts_cur
346!       ENDDO
347!
348!      END SUBROUTINE read_tsurf1d
349!
350!
351!************************************************************************
352END MODULE ocean_forced_mod
353
354
355
356
357
358
Note: See TracBrowser for help on using the repository browser.