Changeset 4843 for LMDZ6/trunk/libf/phylmd/lmdz_thermcell_alp.F90

Timestamp:

Mar 4, 2024, 6:58:03 PM (3 months ago)

Author:

crio

Message:

Nouvelle formulation du strig et correction thermiques montent trop haut

File:

• LMDZ6/trunk/libf/phylmd/lmdz_thermcell_alp.F90 (modified) (7 diffs)

LMDZ6/trunk/libf/phylmd/lmdz_thermcell_alp.F90

@@ -11 +11 @@
      &                  ,pcon,rhobarz,wth3,wmax_sec,lalim,fm,alim_star,zmax &  ! in
 !
-     &                  ,ale_bl,alp_bl,lalim_conv,wght_th &                    ! out
+     &                  ,zcong,ale_bl,alp_bl,lalim_conv,wght_th &                    ! out
      &                  ,zlcl,fraca0,w0,w_conv,therm_tke_max0,env_tke_max0 &   ! out
-     &                  ,n2,s2,ale_bl_stat &                                   ! out
+     &                  ,n2,s2,strig,ale_bl_stat &                                   ! out
      &                  ,therm_tke_max,env_tke_max &                           ! out
      &                  ,alp_bl_det,alp_bl_fluct_m,alp_bl_fluct_tke &          ! out
@@ -50 +50 @@
       real, intent(in) :: pplay(ngrid,nlay),pplev(ngrid,nlay+1)
       integer, intent(in), dimension(ngrid) ::lmax,lalim
-      real, intent(in), dimension(ngrid) :: zmax
+      real, intent(in), dimension(ngrid) :: zmax, zcong
       real, intent(in), dimension(ngrid,nlay+1) :: zw2
       real, intent(in), dimension(ngrid,nlay+1) :: fraca
@@ -69 +69 @@
       integer, intent(out), dimension(ngrid) :: lalim_conv
       real, intent(out), dimension(ngrid) :: zlcl,fraca0,w0,w_conv
-      real, intent(out), dimension(ngrid) :: therm_tke_max0,env_tke_max0,n2,s2,ale_bl_stat
+      real, intent(out), dimension(ngrid) :: therm_tke_max0,env_tke_max0,n2,s2,ale_bl_stat,strig
       real, intent(out), dimension(ngrid,nlay) :: therm_tke_max,env_tke_max
       real, intent(out), dimension(ngrid) :: alp_bl_det,alp_bl_fluct_m,alp_bl_fluct_tke
@@ -94 +94 @@
       real, parameter :: zmax_moy_coef=0.33
       real, dimension(ngrid) :: depth        ! Epaisseur moyenne du cumulus
+      real, dimension(ngrid) :: zcong_moy 
       real, dimension(ngrid) ::  w_max                 ! Vitesse max statistique 
       real, dimension(ngrid) ::  s_max(ngrid)
@@ -119 +120 @@
      s2(ig) = 0.
      ale_bl_stat(ig) = 0.
+     strig(ig) = 0.
      alp_bl_det(ig) = 0.
      alp_bl_fluct_m(ig) = 0.
@@ -249 +251 @@
    s_max(:)=0.
 
+
+
 !-----Epaisseur du nuage (depth) et détermination de la queue du spectre de panaches (n2,s2) et du panache le plus gros (s_max)
    do ig=1,ngrid
@@ -270 +274 @@
 !   print *,'avant Calcul de Wmax '    !!jyg
 
+!CR: calcul de strig
+   if (iflag_strig.eq.0) then
+      strig(:)=s_trig
+   else if (iflag_strig.eq.1) then
+      do ig=1,ngrid
+!         zcong_moy(ig)=zlcl(ig)+zmax_moy_coef*(zcong(ig)-zlcl(ig))
+!         strig(ig)=(hcoef*(zcong_moy(ig)-zlcl(ig))+hmin(ig))**2  
+         strig(ig)=(zcong(ig)-zlcl(ig))**2  
+      enddo
+   else if (iflag_strig.eq.2) then
+      do ig=1,ngrid 
+         if (h_trig.gt.zlcl(ig)) then
+         strig(ig)=(h_trig-zlcl(ig))**2
+         else
+         strig(ig)=s_trig
+         endif
+      enddo   
+   endif
+
    susqr2pi=su_cst*sqrt(2.*Rpi)
    reuler=exp(1.)