Changeset 4449


Ignore:
Timestamp:
Mar 2, 2023, 4:17:32 PM (2 years ago)
Author:
evignon
Message:

commission du nouveau schema de turbulence developpe
dans le cadre de l'atelier tke

Location:
LMDZ6/trunk/libf
Files:
4 added
4 edited

Legend:

Unmodified
Added
Removed

  • LMDZ6/trunk/libf/phylmd/pbl_surface_mod.F90

    r4370 r4449  
    2323  USE climb_wind_mod,      ONLY : climb_wind_down, climb_wind_up
    2424  USE coef_diff_turb_mod,  ONLY : coef_diff_turb
     25  USE atke_exchange_coeff_mod, ONLY :  atke_compute_km_kh
    2526  USE ioipsl_getin_p_mod,  ONLY : getin_p
    2627  USE cdrag_mod
    … …  
    16561657!!! jyg le 07/02/2012
    16571658       IF (iflag_split .eq.0) THEN
    1658 !!!
    16591659!!! nrlmd & jyg les 02/05/2011, 13/06/2011, 05/02/2012
    16601660      IF (prt_level >=10) THEN
    … …  
    16691669      print *,' args coef_diff_turb: ycdragh ', ycdragh
    16701670      print *,' args coef_diff_turb: ytke ', ytke
    1671 
    16721671       ENDIF
     1672
     1673        IF (iflag_pbl>=50) THEN
     1674
     1675        CALL atke_compute_km_kh(knon,klev,yu,yv,yt, &
     1676             ypplay,ypaprs,ytke,ycoefm, ycoefh)
     1677
     1678        ELSE
     1679
    16731680        CALL coef_diff_turb(dtime, nsrf, knon, ni,  &
    16741681            ypaprs, ypplay, yu, yv, yq, yt, yts, yqsurf, ycdragm, &
    … …  
    16861693       ENDDO
    16871694       ENDIF
     1695
     1696       ENDIF ! iflag_pbl >= 50
     1697
    16881698        IF (prt_level >=10) print *,'coef_diff_turb -> ycoefh ',ycoefh
    1689 !
     1699
     1700
    16901701       ELSE  !(iflag_split .eq.0)
     1702
     1703     
    16911704      IF (prt_level >=10) THEN
    16921705      print *,' args coef_diff_turb: yu_x ',  yu_x 
    … …  
    16991712      print *,' args coef_diff_turb: ycdragh_x ', ycdragh_x
    17001713      print *,' args coef_diff_turb: ytke_x ', ytke_x
    1701 
    1702        ENDIF
     1714      ENDIF
     1715
     1716
     1717        IF (iflag_pbl>=50) THEN
     1718     
     1719        CALL atke_compute_km_kh(knon,klev,yu_x,yv_x,yt_x, &
     1720             ypplay,ypaprs,ytke_x,ycoefm_x, ycoefh_x)
     1721
     1722        ELSE
     1723
    17031724        CALL coef_diff_turb(dtime, nsrf, knon, ni,  &
    17041725            ypaprs, ypplay, yu_x, yv_x, yq_x, yt_x, yts_x, yqsurf_x, ycdragm_x, &
    … …  
    17161737       ENDDO
    17171738       ENDIF
     1739
     1740        ENDIF ! iflag_pbl >= 50
     1741
    17181742        IF (prt_level >=10) print *,'coef_diff_turb -> ycoefh_x ',ycoefh_x
    17191743!
    … …  
    17281752      print *,' args coef_diff_turb: ycdragh_w ', ycdragh_w
    17291753      print *,' args coef_diff_turb: ytke_w ', ytke_w
    1730        ENDIF
     1754      ENDIF
     1755     
     1756        IF (iflag_pbl>=50) THEN
     1757       
     1758        CALL atke_compute_km_kh(knon,klev,yu_w,yv_w,yt_w, &
     1759             ypplay,ypaprs,ytke_w,ycoefm_w, ycoefh_w)
     1760
     1761        ELSE
     1762
    17311763        CALL coef_diff_turb(dtime, nsrf, knon, ni,  &
    17321764            ypaprs, ypplay, yu_w, yv_w, yq_w, yt_w, yts_w, yqsurf_w, ycdragm_w, &
    … …  
    17431775       ENDDO
    17441776       ENDIF
     1777
     1778       ENDIF ! iflag_pbl >= 50
     1779
     1780
    17451781        IF (prt_level >=10) print *,'coef_diff_turb -> ycoefh_w ',ycoefh_w
    1746 !
     1782
    17471783!!!jyg le 10/04/2013
    17481784!!   En attendant de traiter le transport des traceurs dans les poches froides, formule
    … …  
    17541790        ENDDO
    17551791      ENDDO
    1756 !!!
     1792
     1793
    17571794       ENDIF  ! (iflag_split .eq.0)
    1758 !!!
     1795
    17591796       
    17601797!****************************************************************************************

  • LMDZ6/trunk/libf/phylmd/physiq_mod.F90

    r4448 r4449  
    8282    USE wake_ini_mod, ONLY : wake_ini
    8383    USE yamada_ini_mod, ONLY : yamada_ini
    84 
     84    USE atke_turbulence_ini_mod, ONLY : atke_ini
    8585    USE thermcell_ini_mod, ONLY : thermcell_ini
    8686    USE lscp_ini_mod, ONLY : lscp_ini
    … …  
    17481748       CALL wake_ini(rg,rd,rv,prt_level)
    17491749       CALL yamada_ini(klon,lunout,prt_level)
     1750       CALL atke_ini(prt_level, lunout, RG, RD, RPI)
    17501751       CALL thermcell_ini(iflag_thermals,prt_level,tau_thermals,lunout, &
    17511752   &    RG,RD,RCPD,RKAPPA,RLVTT,RETV)

  • LMDZ6/trunk/libf/phylmdiso/pbl_surface_mod.F90

    r4374 r4449  
    2323  USE climb_wind_mod,      ONLY : climb_wind_down, climb_wind_up
    2424  USE coef_diff_turb_mod,  ONLY : coef_diff_turb
     25  USE atke_exchange_coeff_mod, ONLY :  atke_compute_km_kh
    2526  USE ioipsl_getin_p_mod,  ONLY : getin_p
    2627  USE cdrag_mod
    … …  
    19511952
    19521953       ENDIF
     1954
     1955        IF (iflag_pbl>=50) THEN
     1956
     1957        CALL atke_compute_km_kh(knon,klev,yu,yv,yt, &
     1958             ypplay,ypaprs,ytke,ycoefm, ycoefh)
     1959
     1960        ELSE
     1961
     1962
    19531963        CALL coef_diff_turb(dtime, nsrf, knon, ni,  &
    19541964            ypaprs, ypplay, yu, yv, yq, yt, yts, yqsurf, ycdragm, &
    … …  
    19661976       ENDDO
    19671977       ENDIF
     1978
     1979       ENDIF ! iflag_pbl >= 50
     1980
    19681981        IF (prt_level >=10) print *,'coef_diff_turb -> ycoefh ',ycoefh
    19691982!
    … …  
    19811994
    19821995       ENDIF
     1996
     1997        IF (iflag_pbl>=50) THEN
     1998
     1999        CALL atke_compute_km_kh(knon,klev,yu_x,yv_x,yt_x, &
     2000             ypplay,ypaprs,ytke_x,ycoefm_x, ycoefh_x)
     2001
     2002        ELSE
     2003
     2004
    19832005        CALL coef_diff_turb(dtime, nsrf, knon, ni,  &
    19842006            ypaprs, ypplay, yu_x, yv_x, yq_x, yt_x, yts_x, yqsurf_x, ycdragm_x, &
    … …  
    19962018       ENDDO
    19972019       ENDIF
     2020
     2021       ENDIF ! iflag_pbl >= 50
     2022
    19982023        IF (prt_level >=10) print *,'coef_diff_turb -> ycoefh_x ',ycoefh_x
    19992024!
    … …  
    20092034      print *,' args coef_diff_turb: ytke_w ', ytke_w
    20102035       ENDIF
     2036
     2037        IF (iflag_pbl>=50) THEN
     2038
     2039        CALL atke_compute_km_kh(knon,klev,yu_w,yv_w,yt_w, &
     2040             ypplay,ypaprs,ytke_w,ycoefm_w, ycoefh_w)
     2041
     2042        ELSE
     2043
     2044
    20112045        CALL coef_diff_turb(dtime, nsrf, knon, ni,  &
    20122046            ypaprs, ypplay, yu_w, yv_w, yq_w, yt_w, yts_w, yqsurf_w, ycdragm_w, &
    … …  
    20232057       ENDDO
    20242058       ENDIF
     2059
     2060       ENDIF ! iflag_pbl >= 50
     2061
    20252062        IF (prt_level >=10) print *,'coef_diff_turb -> ycoefh_w ',ycoefh_w
    20262063!

  • LMDZ6/trunk/libf/phylmdiso/physiq_mod.F90

    r4389 r4449  
    7979    USE wxios, ONLY: g_ctx, wxios_set_context
    8080#endif
     81    USE atke_turbulence_ini_mod, ONLY : atke_ini
    8182    USE lscp_ini_mod, ONLY : lscp_ini
    8283    USE lscp_mod, ONLY : lscp
    … …  
    18841885!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
    18851886       CALL wake_ini(rg,rd,rv,prt_level)
     1887       CALL atke_ini(prt_level, lunout, RG, RD, RPI)
    18861888       CALL thermcell_ini(iflag_thermals,prt_level,tau_thermals,lunout, &
    18871889   &    RG,RD,RCPD,RKAPPA,RLVTT,RETV)
Note: See TracChangeset for help on using the changeset viewer.