Ignore:
Timestamp:
Nov 17, 2021, 11:23:39 PM (4 years ago)
Author:
evignon
Message:

! Prise en compte de l'orographie sous maille, des heterogeneites de surface
! sur le ratqs + nouvelle version des ratqs interactifs de Louis
! Le tout est dan un module: calcratqs_multi_mod.
! Pour l'instant, les nouvelles contributions peuvent s'activer
! uniquement de facon separee

Location:
LMDZ6/trunk/libf/phylmdiso
Files:
1 added
2 edited

Legend:

Unmodified
Added
Removed
  • LMDZ6/trunk/libf/phylmdiso/phys_local_var_mod.F90

    r4003 r4009  
    1616      REAL, SAVE, ALLOCATABLE :: u_seri(:,:), v_seri(:,:)
    1717      !$OMP THREADPRIVATE(u_seri, v_seri)
    18       REAL, SAVE, ALLOCATABLE :: l_mixmin(:,:,:), l_mix(:,:,:), tke_dissip(:,:,:)
    19       !$OMP THREADPRIVATE(l_mixmin, l_mix, tke_dissip)
     18      REAL, SAVE, ALLOCATABLE :: l_mixmin(:,:,:), l_mix(:,:,:), tke_dissip(:,:,:), wprime(:,:,:)
     19      !$OMP THREADPRIVATE(l_mixmin, l_mix, tke_dissip, wprime)
    2020      REAL, SAVE, ALLOCATABLE :: tr_seri(:,:,:)
    2121      !$OMP THREADPRIVATE(tr_seri)
     
    693693      ALLOCATE(t_seri(klon,klev),q_seri(klon,klev),ql_seri(klon,klev),qs_seri(klon,klev))
    694694      ALLOCATE(u_seri(klon,klev),v_seri(klon,klev))
    695       ALLOCATE(l_mixmin(klon,klev+1,nbsrf), l_mix(klon,klev+1,nbsrf), tke_dissip(klon,klev+1,nbsrf))
    696       l_mix(:,:,:)=0. ; l_mixmin(:,:,:)=0. ; tke_dissip(:,:,:)=0. ! doit etre initialse car pas toujours remplis
     695      ALLOCATE(l_mixmin(klon,klev+1,nbsrf), l_mix(klon,klev+1,nbsrf), tke_dissip(klon,klev+1,nbsrf), wprime(klon,klev+1,nbsrf))
     696      l_mix(:,:,:)=0. ; l_mixmin(:,:,:)=0. ; tke_dissip(:,:,:)=0. ; wprime(:,:,:)=0. ! doit etre initialse car pas toujours remplis
    697697
    698698      ALLOCATE(tr_seri(klon,klev,nbtr))
     
    10901090      DEALLOCATE(t_seri,q_seri,ql_seri,qs_seri)
    10911091      DEALLOCATE(u_seri,v_seri)
    1092       DEALLOCATE(l_mixmin,l_mix, tke_dissip)
     1092      DEALLOCATE(l_mixmin,l_mix, tke_dissip, wprime)
    10931093
    10941094      DEALLOCATE(tr_seri)
  • LMDZ6/trunk/libf/phylmdiso/physiq_mod.F90

    r4003 r4009  
    264264       d_t_vdf_x, d_t_vdf_w, &
    265265       d_q_vdf_x, d_q_vdf_w, &
    266        pbl_tke_input, &
     266       pbl_tke_input, tke_dissip, l_mix, wprime, &
    267267       t_therm, q_therm, u_therm, v_therm, &
    268268       cdragh_x, cdragh_w, &
     
    12881288    integer iostat
    12891289
     1290    REAL, dimension(klon,klev+1) :: tke_dissip_ave, l_mix_ave, wprime_ave
    12901291    REAL zzz
    12911292    !albedo SB >>>
     
    44684469    ! Computation of ratqs, the width (normalized) of the subrid scale
    44694470    ! water distribution
    4470     CALL  calcratqs(klon,klev,prt_level,lunout,        &
     4471
     4472    tke_dissip_ave(:,:)=0.
     4473    l_mix_ave(:,:)=0.
     4474    wprime_ave(:,:)=0.
     4475
     4476    DO nsrf = 1, nbsrf
     4477       DO i = 1, klon
     4478          tke_dissip_ave(i,:) = tke_dissip_ave(i,:) + tke_dissip(i,:,nsrf)*pctsrf(i,nsrf)
     4479          l_mix_ave(i,:) = l_mix_ave(i,:) + l_mix(i,:,nsrf)*pctsrf(i,nsrf)
     4480          wprime_ave(i,:) = wprime_ave(i,:) + wprime(i,:,nsrf)*pctsrf(i,nsrf)
     4481       ENDDO
     4482    ENDDO
     4483
     4484   CALL  calcratqs(klon,klev,prt_level,lunout,        &
    44714485         iflag_ratqs,iflag_con,iflag_cld_th,pdtphys,  &
    44724486         ratqsbas,ratqshaut,ratqsp0, ratqsdp, &
    44734487         tau_ratqs,fact_cldcon,wake_s, wake_deltaq,   &
    44744488         ptconv,ptconvth,clwcon0th, rnebcon0th,     &
    4475          paprs,pplay,q_seri,zqsat,fm_therm, &
     4489         paprs,pplay,t_seri,q_seri, qtc_cv, sigt_cv, zqsat, &
     4490         pbl_tke(:,:,is_ave),tke_dissip_ave,l_mix_ave,wprime_ave,t2m,q2m,fm_therm, &
    44764491         ratqs,ratqsc,ratqs_inter)
     4492
    44774493
    44784494    !
Note: See TracChangeset for help on using the changeset viewer.