- Timestamp:
- Sep 11, 2023, 9:56:55 PM (15 months ago)
- Location:
- LMDZ6/trunk/libf
- Files:
-
- 5 edited
- 6 moved
Legend:
- Unmodified
- Added
- Removed
-
LMDZ6/trunk/libf/phylmd/cdrag_mod.F90
r4481 r4687 23 23 USE print_control_mod, ONLY: lunout, prt_level 24 24 USE ioipsl_getin_p_mod, ONLY : getin_p 25 USE atke_turbulence_ini_mod, ONLY : ric, cinf, cepsilon, pr_slope, pr_asym, pr_neut25 USE lmdz_atke_turbulence_ini, ONLY : ric, cinf, cepsilon, pr_slope, pr_asym, pr_neut 26 26 27 27 IMPLICIT NONE -
LMDZ6/trunk/libf/phylmd/lmdz_atke_exchange_coeff.F90
r4686 r4687 1 module atke_exchange_coeff_mod1 module lmdz_atke_exchange_coeff 2 2 3 3 implicit none … … 26 26 27 27 28 USE atke_turbulence_ini_mod, ONLY : iflag_atke, kappa, l0, ric, cinf, rpi, rcpd, atke_ok_virtual29 USE atke_turbulence_ini_mod, ONLY : cepsilon, pr_slope, pr_asym, pr_neut, ctkes,rg, rd, rv, atke_ok_vdiff30 USE atke_turbulence_ini_mod, ONLY : viscom, viscoh, clmix, clmixshear, iflag_atke_lmix, lmin, smmin28 USE lmdz_atke_turbulence_ini, ONLY : iflag_atke, kappa, l0, ric, cinf, rpi, rcpd, atke_ok_virtual 29 USE lmdz_atke_turbulence_ini, ONLY : cepsilon, pr_slope, pr_asym, pr_neut, ctkes,rg, rd, rv, atke_ok_vdiff 30 USE lmdz_atke_turbulence_ini, ONLY : viscom, viscoh, clmix, clmixshear, iflag_atke_lmix, lmin, smmin 31 31 32 32 implicit none … … 182 182 lstrat=clmix*sqrt(tke(igrid,ilay))/sqrt(N2(igrid,ilay)) 183 183 lstrat=max(lstrat,lmin) 184 !Monin-Obukhov consistent interpolation, Van de Wiel et al. 2010 185 l_exchange(igrid,ilay)=(1./(2.*l_exchange(igrid,ilay))+sqrt(1./(4.*l_exchange(igrid,ilay) & 186 *l_exchange(igrid,ilay))+1./(2.*lstrat*lstrat)))**(-1.0) 184 !Inverse quadratic interpolation, Van de Wiel et al. 2010 185 l_exchange(igrid,ilay)=(1./(l_exchange(igrid,ilay))+1./(lstrat))**(-1.0) 187 186 ENDIF 188 187 ENDDO … … 198 197 clmixshear*sqrt(tke(igrid,ilay))/sqrt(shear2(igrid,ilay))) 199 198 lstrat=max(lstrat,lmin) 200 !Monin-Obukhov consistent interpolation, Van de Wiel et al. 2010 201 l_exchange(igrid,ilay)=(1./(2.*l_exchange(igrid,ilay))+sqrt(1./(4.*l_exchange(igrid,ilay) & 202 *l_exchange(igrid,ilay))+1./(2.*lstrat*lstrat)))**(-1.0) 203 199 !Inverse quadratic interpolation, Van de Wiel et al. 2010 200 l_exchange(igrid,ilay)=(1./(l_exchange(igrid,ilay))+1./(lstrat))**(-1.0) 204 201 ENDIF 205 202 ENDDO … … 356 353 ! E Vignon, July 2023 357 354 358 USE atke_turbulence_ini_mod, ONLY : rd, cke, viscom355 USE lmdz_atke_turbulence_ini, ONLY : rd, cke, viscom 359 356 360 357 … … 449 446 450 447 451 end module atke_exchange_coeff_mod448 end module lmdz_atke_exchange_coeff -
LMDZ6/trunk/libf/phylmd/lmdz_atke_turbulence_ini.F90
r4686 r4687 1 MODULE atke_turbulence_ini_mod1 MODULE lmdz_atke_turbulence_ini 2 2 3 3 implicit none … … 139 139 END SUBROUTINE atke_ini 140 140 141 END MODULE atke_turbulence_ini_mod141 END MODULE lmdz_atke_turbulence_ini -
LMDZ6/trunk/libf/phylmd/lmdz_call_atke.F90
r4686 r4687 1 module call_atke_mod1 module lmdz_call_atke 2 2 3 USE atke_exchange_coeff_mod, ONLY : atke_compute_km_kh3 USE lmdz_atke_exchange_coeff, ONLY : atke_compute_km_kh 4 4 5 5 implicit none … … 15 15 16 16 17 USE atke_turbulence_ini_mod, ONLY : iflag_num_atke, rg, rd17 USE lmdz_atke_turbulence_ini, ONLY : iflag_num_atke, rg, rd 18 18 19 19 implicit none … … 167 167 168 168 169 end module call_atke_mod169 end module lmdz_call_atke -
LMDZ6/trunk/libf/phylmd/pbl_surface_mod.F90
r4653 r4687 24 24 USE climb_wind_mod, ONLY : climb_wind_down, climb_wind_up 25 25 USE coef_diff_turb_mod, ONLY : coef_diff_turb 26 USE call_atke_mod, ONLY :call_atke26 USE lmdz_call_atke, ONLY : call_atke 27 27 USE ioipsl_getin_p_mod, ONLY : getin_p 28 28 USE cdrag_mod -
LMDZ6/trunk/libf/phylmd/physiq_mod.F90
r4683 r4687 83 83 USE lmdz_wake_ini, ONLY : wake_ini 84 84 USE yamada_ini_mod, ONLY : yamada_ini 85 USE atke_turbulence_ini_mod, ONLY : atke_ini85 USE lmdz_atke_turbulence_ini, ONLY : atke_ini 86 86 USE lmdz_thermcell_ini, ONLY : thermcell_ini 87 87 USE lmdz_thermcell_dtke, ONLY : thermcell_dtke -
LMDZ6/trunk/libf/phylmdiso/lmdz_atke_exchange_coeff.F90
r4686 r4687 1 link ../phylmd/ atke_exchange_coeff_mod.F901 link ../phylmd/lmdz_atke_exchange_coeff.F90 -
LMDZ6/trunk/libf/phylmdiso/lmdz_atke_turbulence_ini.F90
r4686 r4687 1 link ../phylmd/ atke_turbulence_ini_mod.F901 link ../phylmd/lmdz_atke_turbulence_ini.F90 -
LMDZ6/trunk/libf/phylmdiso/lmdz_call_atke.F90
r4686 r4687 1 link ../phylmd/ call_atke_mod.F901 link ../phylmd/lmdz_call_atke.F90 -
LMDZ6/trunk/libf/phylmdiso/pbl_surface_mod.F90
r4653 r4687 24 24 USE climb_wind_mod, ONLY : climb_wind_down, climb_wind_up 25 25 USE coef_diff_turb_mod, ONLY : coef_diff_turb 26 USE call_atke_mod,ONLY : call_atke26 USE lmdz_call_atke, ONLY : call_atke 27 27 USE ioipsl_getin_p_mod, ONLY : getin_p 28 28 USE cdrag_mod -
LMDZ6/trunk/libf/phylmdiso/physiq_mod.F90
r4685 r4687 82 82 USE lmdz_wake_ini, ONLY : wake_ini 83 83 USE yamada_ini_mod, ONLY : yamada_ini 84 USE atke_turbulence_ini_mod, ONLY : atke_ini84 USE lmdz_atke_turbulence_ini, ONLY : atke_ini 85 85 USE lmdz_thermcell_ini, ONLY : thermcell_ini 86 86 USE lmdz_thermcell_dtke, ONLY : thermcell_dtke
Note: See TracChangeset
for help on using the changeset viewer.