Changeset 5805 for LMDZ6


Ignore:
Timestamp:
Sep 10, 2025, 3:34:13 AM (2 months ago)
Author:
fhourdin
Message:

Details poches

Location:
LMDZ6/trunk/libf/phylmd
Files:
3 edited

Legend:

Unmodified
Added
Removed

  • LMDZ6/trunk/libf/phylmd/lmdz_wake_dadv.f90

    r5804 r5805  
    6868!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
    6969
    70 print *,'ZZZwake_dadv_IN wk_adv(1) ', wk_adv(1)
    71 print *,'ZZZwake_dadv_IN kupper(1) ', kupper(1)
    72 print *,'ZZZwake_dadv_IN k, thw(1,k), thx(1,k) ', (k, thw(1,k), thx(1,k), k = 1,3)
    73 print *,'ZZZwake_dadv_IN k, deltomg(1,k) ', (k, deltomg(1,k), k = 1,3)
    74 print *,'ZZZwake_dadv_IN k, dp_deltomg(1,k) ', (k, dp_deltomg(1,k), k = 1,3)
    75 print *,'ZZZwake_dadv_IN sigmaw(1) ', sigmaw(1)
    76 print *,'ZZZwake_dadv_IN dsigspread(1) ', dsigspread(1)
     70! print *,'ZZZwake_dadv_IN wk_adv(1) ', wk_adv(1)
     71! print *,'ZZZwake_dadv_IN kupper(1) ', kupper(1)
     72! print *,'ZZZwake_dadv_IN k, thw(1,k), thx(1,k) ', (k, thw(1,k), thx(1,k), k = 1,3)
     73! print *,'ZZZwake_dadv_IN k, deltomg(1,k) ', (k, deltomg(1,k), k = 1,3)
     74! print *,'ZZZwake_dadv_IN k, dp_deltomg(1,k) ', (k, dp_deltomg(1,k), k = 1,3)
     75! print *,'ZZZwake_dadv_IN sigmaw(1) ', sigmaw(1)
     76! print *,'ZZZwake_dadv_IN dsigspread(1) ', dsigspread(1)
    7777
    7878    entr_s(:,:) = 0.
    … …  
    463463  ENDIF! (flag_dadv_implicit)
    464464
    465 print *,'ZZZwake_dadv k, d_deltat_dadv(1,k) ', (k, d_deltat_dadv(1,k), k = 1,3)
     465!print *,'ZZZwake_dadv k, d_deltat_dadv(1,k) ', (k, d_deltat_dadv(1,k), k = 1,3)
    466466
    467467    END SUBROUTINE wake_dadv

  • LMDZ6/trunk/libf/phylmd/lmdz_wake_ini.f90

    r5804 r5805  
    9494
    9595  ! =========================================================================
    96   SUBROUTINE wake_ini(rg_in,rd_in,rv_in,prt_lev)
     96  SUBROUTINE wake_ini(iflag_wake,rg_in,rd_in,rv_in,prt_lev)
    9797  ! =========================================================================
    9898
    … …  
    112112  USE ioipsl_getin_p_mod, ONLY : getin_p
    113113  real eps
    114   integer, intent(in) :: prt_lev
     114  integer, intent(in) :: prt_lev,iflag_wake
    115115  real, intent(in) :: rg_in,rd_in,rv_in
    116116
    … …  
    202202  stark=0.33
    203203  CALL getin_p('stark',stark)
    204   cstart = stark*sqrt(2.*wapecut)
     204  if ( iflag_wake / 10 >= 3 ) then
     205      ! Calcul effectif tenant compte dans l'idee d'une heterogeneite des poches
     206      ! Pour eviter une singularité
     207      cstart = stark*sqrt(2.*wapecut)/2.
     208  else
     209      cstart = stark*sqrt(2.*wapecut)
     210  endif
    205211
    206212  alpk=0.25

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

    r5786 r5805  
    18561856       
    18571857!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
    1858        CALL wake_ini(rg,rd,rv,prt_level)
     1858       CALL wake_ini(iflag_wake,rg,rd,rv,prt_level)
    18591859       CALL yamada_ini(klon,lunout,prt_level)
    18601860       viscom=1.46E-5
Note: See TracChangeset for help on using the changeset viewer.