Changeset 4843


Ignore:
Timestamp:
Mar 4, 2024, 6:58:03 PM (10 months ago)
Author:
crio
Message:

Nouvelle formulation du strig et correction thermiques montent trop haut

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

Legend:

Unmodified
Added
Removed

  • LMDZ6/trunk/libf/phylmd/alpale.h

    r4090 r4843  
    55
    66! dans alealp_th, thermcell_alp, physiq_mod, conf_phys
    7       integer            :: iflag_trig_bl,iflag_clos_bl
     7      integer            :: iflag_trig_bl,iflag_clos_bl,iflag_strig
    88      integer            :: tau_trig_shallow,tau_trig_deep
    9       real               :: s_trig
     9      real               :: s_trig, h_trig
    1010! thermcell_alp et convection ...
    1111      integer            :: iflag_coupl,iflag_clos,iflag_wake
    … …  
    1313      real               :: alp_bl_k
    1414
    15       common/calpale1/iflag_trig_bl,iflag_clos_bl,tau_trig_shallow,tau_trig_deep
    16       common/calpale2/s_trig,iflag_coupl,iflag_clos,iflag_wake,alp_bl_k
     15      common/calpale1/iflag_trig_bl,iflag_clos_bl,tau_trig_shallow,tau_trig_deep,iflag_strig
     16      common/calpale2/s_trig,iflag_coupl,iflag_clos,iflag_wake,alp_bl_k,h_trig
    1717
    1818!$OMP THREADPRIVATE(/calpale1/,/calpale2/)

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

    r4827 r4843  
    33!
    44SUBROUTINE alpale_th ( dtime, lmax_th, t_seri, cell_area,  &
    5                        cin, s2, n2, &
     5                       cin, s2, n2, strig, &
    66                       ale_bl_trig, ale_bl_stat, ale_bl,  &
    77                       alp_bl, alp_bl_stat, &
    … …  
    3636  REAL, DIMENSION(klon), INTENT(IN)                          :: ale_bl_stat
    3737  REAL, DIMENSION(klon), INTENT(IN)                          :: cin
    38   REAL, DIMENSION(klon), INTENT(IN)                          :: s2, n2
     38  REAL, DIMENSION(klon), INTENT(IN)                          :: s2, n2, strig
    3939                                                               
    4040  REAL, DIMENSION(klon), INTENT(INOUT)                       :: ale_bl_trig, ale_bl
    … …  
    161161                   if ( (ale_bl_stat(i) .gt. abs(cin(i))+1.e-10) )  then
    162162                      proba_notrig(i)=proba_notrig(i)* &
    163                          (1.-exp(-s_trig/s2(i)))**(n2(i)*dtime/tau_trig(i))
     163                         (1.-exp(-strig(i)/s2(i)))**(n2(i)*dtime/tau_trig(i))
    164164                      !        print *, 'proba_notrig(i) ',proba_notrig(i)
    165165                      if (random_notrig(i) .ge. proba_notrig(i)) then
    … …  
    168168                         ale_bl_trig(i)=0.
    169169                      endif
    170                       birth_rate(i) = n2(i)*exp(-s_trig/s2(i))/(tau_trig(i)*cell_area(i))
     170                      birth_rate(i) = n2(i)*exp(-strig(i)/s2(i))/(tau_trig(i)*cell_area(i))
    171171!!!                      birth_rate(i) = max(birth_rate(i),1.e-18)
    172172                   else
    … …  
    183183                   if ( (Ale_bl(i) .gt. abs(cin(i))+1.e-10) )  then
    184184                      proba_notrig(i)=proba_notrig(i)* &
    185                          (1.-exp(-s_trig/s2(i)))**(n2(i)*dtime/tau_trig(i))
     185                         (1.-exp(-strig(i)/s2(i)))**(n2(i)*dtime/tau_trig(i))
    186186                      !        print *, 'proba_notrig(i) ',proba_notrig(i)
    187187                      if (random_notrig(i) .ge. proba_notrig(i)) then
    … …  
    190190                         ale_bl_trig(i)=0.
    191191                      endif
    192                       birth_rate(i) = n2(i)*exp(-s_trig/s2(i))/(tau_trig(i)*cell_area(i))
     192                      birth_rate(i) = n2(i)*exp(-strig(i)/s2(i))/(tau_trig(i)*cell_area(i))
    193193!!!                      birth_rate(i) = max(birth_rate(i),1.e-18)
    194194                   else
    … …  
    283283        do i=1,klon
    284284           if ( (ale_bl_ref(i) .gt. abs(cin(i))+1.e-10) )  then
    285               birth_number = n2(i)*exp(-s_trig/s2(i))
     285              birth_number = n2(i)*exp(-strig(i)/s2(i))
    286286              birth_rate(i) = birth_number/(tau_trig(i)*cell_area(i))
    287287!!!              birth_rate(i) = max(birth_rate(i),1.e-18)

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

    r4837 r4843  
    1212     &      ,pbl_tke,pctsrf,omega,airephy &
    1313     &      ,zlcl_th,fraca0,w0,w_conv,therm_tke_max0,env_tke_max0 &
    14      &      ,n2,s2,ale_bl_stat &
     14     &      ,n2,s2,strig,zcong,ale_bl_stat &
    1515     &      ,therm_tke_max,env_tke_max &
    1616     &      ,alp_bl_det,alp_bl_fluct_m,alp_bl_fluct_tke &
    … …  
    8787      real zthl(klon,klev)
    8888      real wmax_sec(klon)
     89      real zcong(klon)
    8990      real zmax_sec(klon)
    9091      real f_sec(klon)
    … …  
    116117      real zlcl_th(klon),fraca0(klon),w0(klon),w_conv(klon)
    117118      real therm_tke_max0(klon),env_tke_max0(klon)
    118       real n2(klon),s2(klon)
     119      real n2(klon),s2(klon),strig(klon)
    119120      real ale_bl_stat(klon)
    120121      real therm_tke_max(klon,klev),env_tke_max(klon,klev)
    … …  
    302303     &      ,zmax0,f0,zw2,fraca,ztv,zpspsk &
    303304     &      ,ztla,zthl,ztva &
    304      &      ,pcon,rhobarz,wth3,wmax_sec,lalim,fm,alim_star,zmax &
     305     &      ,pcon,rhobarz,wth3,wmax_sec,lalim,fm,alim_star,zmax,zcong &
    305306#ifdef ISO         
    306307     &      ,xt_seri,d_xt_the &
    … …  
    314315     &        ,zw2,fraca &                                           ! in
    315316     &        ,pcon,rhobarz,wth3,wmax_sec,lalim,fm,alim_star,zmax &  ! in
    316      &        ,ale,alp,lalim_conv,wght_th &                          ! out
     317     &        ,zcong,ale,alp,lalim_conv,wght_th &                          ! out
    317318     &        ,zlcl_th,fraca0,w0,w_conv,therm_tke_max0,env_tke_max0 &! out
    318      &        ,n2,s2,ale_bl_stat &                                   ! out
     319     &        ,n2,s2,strig,ale_bl_stat &                                   ! out
    319320     &        ,therm_tke_max,env_tke_max &                           ! out
    320321     &        ,alp_bl_det,alp_bl_fluct_m,alp_bl_fluct_tke &          ! out

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

    r4754 r4843  
    144144    REAL,SAVE :: tau_thermals_omp,alp_bl_k_omp
    145145    ! nrlmd le 10/04/2012
    146     INTEGER,SAVE :: iflag_trig_bl_omp,iflag_clos_bl_omp
     146    INTEGER,SAVE :: iflag_trig_bl_omp,iflag_clos_bl_omp,iflag_strig_omp
    147147    INTEGER,SAVE :: tau_trig_shallow_omp,tau_trig_deep_omp
    148     REAL,SAVE    :: s_trig_omp
     148    REAL,SAVE    :: s_trig_omp,h_trig_omp
    149149    ! fin nrlmd le 10/04/2012
    150150    REAL :: alp_offset
    … …  
    15591559    CALL getin('iflag_trig_bl',iflag_trig_bl_omp)
    15601560
     1561    !Config Key  = iflag_strig
     1562    !Config Desc = 
     1563    !Config Def  = 0
     1564    !Config Help =
     1565    !
     1566    iflag_strig_omp = 0
     1567    CALL getin('iflag_strig',iflag_strig_omp)
     1568
    15611569    !Config Key  = s_trig_bl
    15621570    !Config Desc = 
    … …  
    15661574    s_trig_omp = 2e7
    15671575    CALL getin('s_trig',s_trig_omp)
     1576
     1577    !Config Key  = h_trig
     1578    !Config Desc = 
     1579    !Config Def  = 0
     1580    !Config Help =
     1581    !
     1582    h_trig_omp = 6000.
     1583    CALL getin('h_trig',h_trig_omp)
    15681584
    15691585    !Config Key  = tau_trig_shallow
    … …  
    24392455    ! nrlmd le 10/04/2012
    24402456    iflag_trig_bl = iflag_trig_bl_omp
     2457    iflag_strig = iflag_strig_omp
    24412458    s_trig = s_trig_omp
     2459    h_trig = h_trig_omp
    24422460    tau_trig_shallow = tau_trig_shallow_omp
    24432461    tau_trig_deep = tau_trig_deep_omp
    … …  
    28352853    ! nrlmd le 10/04/2012
    28362854    WRITE(lunout,*) ' iflag_trig_bl = ', iflag_trig_bl
     2855    WRITE(lunout,*) ' iflag_strig = ', iflag_strig
    28372856    WRITE(lunout,*) ' s_trig = ', s_trig
     2857    WRITE(lunout,*) ' h_trig = ', h_trig
    28382858    WRITE(lunout,*) ' tau_trig_shallow = ', tau_trig_shallow
    28392859    WRITE(lunout,*) ' tau_trig_deep = ', tau_trig_deep

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

    r4593 r4843  
    1111     &                  ,pcon,rhobarz,wth3,wmax_sec,lalim,fm,alim_star,zmax &  ! in
    1212!
    13      &                  ,ale_bl,alp_bl,lalim_conv,wght_th &                    ! out
     13     &                  ,zcong,ale_bl,alp_bl,lalim_conv,wght_th &                    ! out
    1414     &                  ,zlcl,fraca0,w0,w_conv,therm_tke_max0,env_tke_max0 &   ! out
    15      &                  ,n2,s2,ale_bl_stat &                                   ! out
     15     &                  ,n2,s2,strig,ale_bl_stat &                                   ! out
    1616     &                  ,therm_tke_max,env_tke_max &                           ! out
    1717     &                  ,alp_bl_det,alp_bl_fluct_m,alp_bl_fluct_tke &          ! out
    … …  
    5050      real, intent(in) :: pplay(ngrid,nlay),pplev(ngrid,nlay+1)
    5151      integer, intent(in), dimension(ngrid) ::lmax,lalim
    52       real, intent(in), dimension(ngrid) :: zmax
     52      real, intent(in), dimension(ngrid) :: zmax, zcong
    5353      real, intent(in), dimension(ngrid,nlay+1) :: zw2
    5454      real, intent(in), dimension(ngrid,nlay+1) :: fraca
    … …  
    6969      integer, intent(out), dimension(ngrid) :: lalim_conv
    7070      real, intent(out), dimension(ngrid) :: zlcl,fraca0,w0,w_conv
    71       real, intent(out), dimension(ngrid) :: therm_tke_max0,env_tke_max0,n2,s2,ale_bl_stat
     71      real, intent(out), dimension(ngrid) :: therm_tke_max0,env_tke_max0,n2,s2,ale_bl_stat,strig
    7272      real, intent(out), dimension(ngrid,nlay) :: therm_tke_max,env_tke_max
    7373      real, intent(out), dimension(ngrid) :: alp_bl_det,alp_bl_fluct_m,alp_bl_fluct_tke
    … …  
    9494      real, parameter :: zmax_moy_coef=0.33
    9595      real, dimension(ngrid) :: depth        ! Epaisseur moyenne du cumulus
     96      real, dimension(ngrid) :: zcong_moy
    9697      real, dimension(ngrid) ::  w_max                 ! Vitesse max statistique
    9798      real, dimension(ngrid) ::  s_max(ngrid)
    … …  
    119120     s2(ig) = 0.
    120121     ale_bl_stat(ig) = 0.
     122     strig(ig) = 0.
    121123     alp_bl_det(ig) = 0.
    122124     alp_bl_fluct_m(ig) = 0.
    … …  
    249251   s_max(:)=0.
    250252
     253
     254
    251255!-----Epaisseur du nuage (depth) et détermination de la queue du spectre de panaches (n2,s2) et du panache le plus gros (s_max)
    252256   do ig=1,ngrid
    … …  
    270274!   print *,'avant Calcul de Wmax '    !!jyg
    271275
     276!CR: calcul de strig
     277   if (iflag_strig.eq.0) then
     278      strig(:)=s_trig
     279   else if (iflag_strig.eq.1) then
     280      do ig=1,ngrid
     281!         zcong_moy(ig)=zlcl(ig)+zmax_moy_coef*(zcong(ig)-zlcl(ig))
     282!         strig(ig)=(hcoef*(zcong_moy(ig)-zlcl(ig))+hmin(ig))**2 
     283         strig(ig)=(zcong(ig)-zlcl(ig))**2 
     284      enddo
     285   else if (iflag_strig.eq.2) then
     286      do ig=1,ngrid
     287         if (h_trig.gt.zlcl(ig)) then
     288         strig(ig)=(h_trig-zlcl(ig))**2
     289         else
     290         strig(ig)=s_trig
     291         endif
     292      enddo   
     293   endif
     294
    272295   susqr2pi=su_cst*sqrt(2.*Rpi)
    273296   reuler=exp(1.)

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

    r4590 r4843  
    33
    44   SUBROUTINE thermcell_env(ngrid,nlay,po,pt,pu,pv,pplay,  &
    5      &           pplev,zo,zh,zl,ztv,zthl,zu,zv,zpspsk,pqsat,lev_out)
     5     &           pplev,zo,zh,zl,ztv,zthl,zu,zv,zpspsk,pqsat,lcong,lintercong,lev_out)
    66
    77!--------------------------------------------------------------
    … …  
    2222   real,intent(out), dimension(ngrid,nlay) :: zo,zl,zh,ztv,zthl
    2323   real,intent(out), dimension(ngrid,nlay) :: zpspsk,zu,zv,pqsat
    24    
     24   real, intent(out), dimension(ngrid) :: lintercong
     25   integer, intent(out), dimension(ngrid) :: lcong   
    2526! Local
    2627
    … …  
    7879      enddo
    7980   enddo
     81!CR: Calcul du niveau de congelation
     82   do ig=1,ngrid
     83      lcong(ig)=1
     84   enddo
     85   do ig=1,ngrid
     86      do ll=1,nlay-1
     87         if ((zh(ig,ll).gt.273.15).and.(zh(ig,ll+1).le.273.15)) then
     88            lcong(ig)=ll+1
     89            lintercong(ig)=(ll*(zh(ig,ll+1)-zh(ig,ll))  &
     90     &               -zh(ig,ll)+273.15)/(zh(ig,ll+1)-zh(ig,ll))
     91         endif
     92      enddo
     93   enddo
    8094 
    8195 RETURN

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

    r4590 r4843  
    22CONTAINS
    33
    4       SUBROUTINE thermcell_height(ngrid,nlay,lalim,lmin,linter,lmix,  &
    5      &           zw2,zlev,lmax,zmax,zmax0,zmix,wmax)
     4      SUBROUTINE thermcell_height(ngrid,nlay,lalim,lmin,linter,lcong,lintercong,lmix,  &
     5     &           zw2,zlev,lmax,zmax,zmax0,zmix,wmax,zcong)
    66      IMPLICIT NONE
    77
    … …  
    1414! Entree
    1515      integer, intent(in) :: ngrid,nlay
    16       real, intent(in), dimension(ngrid) :: linter
     16      real, intent(in), dimension(ngrid) :: linter,lintercong
    1717      real, intent(in), dimension(ngrid,nlay+1) :: zlev
    1818! Sortie
    19       real, intent(out), dimension(ngrid) :: wmax,zmax,zmax0,zmix
    20       integer, intent(out), dimension(ngrid) :: lmax
     19      real, intent(out), dimension(ngrid) :: wmax,zmax,zmax0,zmix,zcong
     20      integer, intent(out), dimension(ngrid) :: lmax,lcong
    2121! Les deux
    2222     integer, intent(inout), dimension(ngrid) :: lmix,lalim,lmin
    … …  
    2424
    2525! local
    26      real, dimension(ngrid) :: num,denom,zlevinter
     26     real, dimension(ngrid) :: num,denom,zlevinter,zlevintercong
    2727     integer ig,l
    2828
    … …  
    113113      enddo
    114114
     115!CR:calcul du zcong
     116      do  ig=1,ngrid
     117! calcul de zlevintercong
     118         zlevintercong(ig)=(zlev(ig,lcong(ig)+1)-zlev(ig,lcong(ig)))*  &
     119     &    lintercong(ig)+zlev(ig,lcong(ig))-lcong(ig)*(zlev(ig,lcong(ig)+1)  &
     120     &    -zlev(ig,lcong(ig)))
     121         zcong(ig)=zlevintercong(ig)-zlev(ig,1)
     122!         print*,"calcul zcong",lcong(ig),lintercong(ig),zlevintercong(ig),zcong(ig)
     123      enddo
    115124
    116125      endif

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

    r4692 r4843  
    1515     &                  ,zmax0, f0,zw2,fraca,ztv &
    1616     &                  ,zpspsk,ztla,zthl,ztva &
    17      &                  ,pcon,rhobarz,wth3,wmax_sec,lalim,fm,alim_star,zmax &
     17     &                  ,pcon,rhobarz,wth3,wmax_sec,lalim,fm,alim_star,zmax,zcong &
    1818#ifdef ISO         
    1919     &      ,xtpo,xtpdoadj &
    … …  
    134134      real, intent(out), dimension(ngrid,nlay+1) :: fm
    135135      real, intent(out), dimension(ngrid,nlay) :: alim_star
    136       real, intent(out), dimension(ngrid) :: zmax
     136      real, intent(out), dimension(ngrid) :: zmax,zcong
    137137
    138138!   local:
    … …  
    150150      integer ig,k,l,ierr,ll
    151151      logical sorties
    152       real, dimension(ngrid) :: linter,zmix, zmax_sec
    153       integer,dimension(ngrid) :: lmin,lmix,lmix_bis,nivcon
     152      real, dimension(ngrid) :: linter,zmix, zmax_sec,lintercong
     153      integer,dimension(ngrid) :: lmin,lmix,lmix_bis,nivcon, lcong
    154154      real, dimension(ngrid,nlay) :: ztva_est
    155155      real, dimension(ngrid,nlay) :: deltaz,zlay,zdthladj,zu,zv,z_o,zl,zva,zua,z_oa
    … …  
    227227
    228228          CALL thermcell_env(ngrid,nlay,p_o,ptemp_env,puwind,pvwind,pplay,  &
    229          &           pplev,z_o,ztemp_env,zl,ztv,zthl,zu,zv,zpspsk,zqsat,lev_out)
     229         &           pplev,z_o,ztemp_env,zl,ztv,zthl,zu,zv,zpspsk,zqsat,lcong,lintercong,lev_out)
    230230
    231231      else
    … …  
    456456!-------------------------------------------------------------------------------
    457457!
    458       CALL thermcell_height(ngrid,nlay,lalim,lmin,linter,lmix,zw2,  &
    459      &           zlev,lmax,zmax,zmax0,zmix,wmax)
     458      CALL thermcell_height(ngrid,nlay,lalim,lmin,linter,lcong,lintercong,lmix,zw2,  &
     459     &           zlev,lmax,zmax,zmax0,zmix,wmax,zcong)
    460460! Attention, w2 est transforme en sa racine carree dans cette routine
    461461! Le probleme vient du fait que linter et lmix sont souvent egaux a 1.
    … …  
    729729        endif
    730730      enddo
    731       if (ierr==1) then
    732            abort_message = 'thermcellV0_main: les thermiques vont trop haut '
    733            CALL abort_physic (modname,abort_message,1)
    734       endif
     731!      if (ierr==1) then
     732!           abort_message = 'thermcellV0_main: les thermiques vont trop haut '
     733!           CALL abort_physic (modname,abort_message,1)
     734!      endif
    735735
    736736      if (prt_level.ge.1) print*,'14b OK convect8'

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

    r4835 r4843  
    402402      REAL,ALLOCATABLE,SAVE,DIMENSION(:) :: alp_bl_fluct_m,alp_bl_fluct_tke
    403403!$OMP THREADPRIVATE(alp_bl_fluct_m,alp_bl_fluct_tke)
    404       REAL,ALLOCATABLE,SAVE,DIMENSION(:) :: alp_bl_stat, n2, s2
    405 !$OMP THREADPRIVATE(alp_bl_stat, n2, s2)
     404      REAL,ALLOCATABLE,SAVE,DIMENSION(:) :: alp_bl_stat, n2, s2, strig, zcong, zlcl_th
     405!$OMP THREADPRIVATE(alp_bl_stat, n2, s2, strig, zcong, zlcl_th)
    406406      REAL,ALLOCATABLE,SAVE,DIMENSION(:) :: proba_notrig, random_notrig
    407407!$OMP THREADPRIVATE(proba_notrig, random_notrig)
    … …  
    877877      ALLOCATE(alp_bl_fluct_m(klon), alp_bl_fluct_tke(klon))
    878878      alp_bl_fluct_m(:)=0 ; alp_bl_fluct_tke(:)= 0.
    879       ALLOCATE(alp_bl_stat(klon), n2(klon), s2(klon))
     879      ALLOCATE(alp_bl_stat(klon), n2(klon), s2(klon), strig(klon), zcong(klon), zlcl_th(klon))
    880880      alp_bl_stat(:)=0
    881881      ALLOCATE(proba_notrig(klon), random_notrig(klon))
    … …  
    12051205      DEALLOCATE(alp_bl_conv,alp_bl_det)
    12061206      DEALLOCATE(alp_bl_fluct_m,alp_bl_fluct_tke)
    1207       DEALLOCATE(alp_bl_stat, n2, s2)
     1207      DEALLOCATE(alp_bl_stat, n2, s2, strig, zcong, zlcl_th)
    12081208      DEALLOCATE(proba_notrig, random_notrig)
    12091209!FC

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

    r4835 r4843  
    936936  TYPE(ctrl_out), SAVE :: o_s2 = ctrl_out((/ 1, 6, 6, 6, 10, 10, 11, 11, 11, 11/), &
    937937    's2', 'Surface moyenne des panaches de type 2', 'm2', (/ ('', i=1, 10) /))
    938              
     938  TYPE(ctrl_out), SAVE :: o_strig = ctrl_out((/ 1, 6, 6, 6, 10, 10, 11, 11, 11, 11/), &
     939    'strig', 'Surface moyenne pour atteindre niveau de congelation', 'm2', (/ ('', i=1, 10) /))
     940  TYPE(ctrl_out), SAVE :: o_zcong = ctrl_out((/ 1, 6, 6, 6, 10, 10, 11, 11, 11, 11/), &
     941    'zcong', 'niveau de congelation', 'm', (/ ('', i=1, 10) /)) 
     942  TYPE(ctrl_out), SAVE :: o_zlcl_th = ctrl_out((/ 1, 6, 6, 6, 10, 10, 11, 11, 11, 11/), &
     943    'zlcl_th', 'niveau de condensation', 'm', (/ ('', i=1, 10) /))                   
    939944!-------Déclenchement stochastique
    940945  TYPE(ctrl_out), SAVE :: o_proba_notrig = ctrl_out((/ 1, 6, 6, 6, 10, 10, 11, 11, 11, 11/), &

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

    r4835 r4843  
    9898         o_qtaa, o_Clwaa, &
    9999         o_ftd, o_fqd, o_wdtrainA, o_wdtrainS, o_wdtrainM, &
    100          o_n2, o_s2, o_proba_notrig, &
     100         o_n2, o_s2, o_strig, o_zcong, o_zlcl_th, o_proba_notrig, &
    101101         o_random_notrig, o_ale_bl_stat, &
    102102         o_ale_bl_trig, o_alp_bl_det, &
    … …  
    335335         wake_h, &
    336336         wake_omg, d_t_wake, d_q_wake, Vprecip, qtaa, Clw, &
    337          wdtrainA, wdtrainS, wdtrainM, n2, s2, proba_notrig, &
     337         wdtrainA, wdtrainS, wdtrainM, n2, s2, strig, zcong, zlcl_th, proba_notrig, &
    338338         random_notrig, &
    339339         qclr, qcld, qss, qvc, rnebclr, rnebss, gamma_ss, &
    … …  
    15481548               CALL histwrite_phy(o_kh         ,kh         )
    15491549               CALL histwrite_phy(o_kh_x       ,kh_x       )
     1550          CALL histwrite_phy(o_strig, strig)
     1551          CALL histwrite_phy(o_zcong, zcong)
     1552          CALL histwrite_phy(o_zlcl_th, zlcl_th)
    15501553               CALL histwrite_phy(o_kh_w       ,kh_w       )
    15511554             ENDIF   ! (iflag_pbl_split>=1)

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

    r4835 r4843  
    307307       alp_bl_conv,alp_bl_det,  &
    308308       alp_bl_fluct_m,alp_bl_fluct_tke,  &
    309        alp_bl_stat, n2, s2,  &
     309       alp_bl_stat, n2, s2,  strig, zcong, zlcl_th, &
    310310       proba_notrig, random_notrig,  &
    311311!!       cv_gen,  &  !moved to phys_state_var_mod
    … …  
    750750    !--------Stochastic Boundary Layer Triggering: ALE_BL--------
    751751    !---Propri\'et\'es du thermiques au LCL
    752     real zlcl_th(klon)          ! Altitude du LCL calcul\'e
     752!    real zlcl_th(klon)          ! Altitude du LCL calcul\'e
    753753    ! continument (pcon dans
    754754    ! thermcell_main.F90)
    … …  
    37253725               ,pbl_tke_input,pctsrf,omega,cell_area &
    37263726               ,zlcl_th,fraca0,w0,w_conv,therm_tke_max0,env_tke_max0 &
    3727                ,n2,s2,ale_bl_stat &
     3727               ,n2,s2,strig,zcong,ale_bl_stat &
    37283728               ,therm_tke_max,env_tke_max &
    37293729               ,alp_bl_det,alp_bl_fluct_m,alp_bl_fluct_tke &
    … …  
    37693769!
    37703770          CALL alpale_th( phys_tstep, lmax_th, t_seri, cell_area,  &
    3771                           cin, s2, n2, &
     3771                          cin, s2, n2, strig, &
    37723772                          ale_bl_trig, ale_bl_stat, ale_bl,  &
    37733773                          alp_bl, alp_bl_stat, &
Note: See TracChangeset for help on using the changeset viewer.