Changeset 4843 for LMDZ6/trunk/libf/phylmd/alpale_th.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/alpale_th.F90 (modified) (7 diffs)

LMDZ6/trunk/libf/phylmd/alpale_th.F90

@@ -3 +3 @@
 !
 SUBROUTINE alpale_th ( dtime, lmax_th, t_seri, cell_area,  &
-                       cin, s2, n2,  &
+                       cin, s2, n2, strig,  &
                        ale_bl_trig, ale_bl_stat, ale_bl,  &
                        alp_bl, alp_bl_stat, &
@@ -36 +36 @@
   REAL, DIMENSION(klon), INTENT(IN)                          :: ale_bl_stat
   REAL, DIMENSION(klon), INTENT(IN)                          :: cin
-  REAL, DIMENSION(klon), INTENT(IN)                          :: s2, n2
+  REAL, DIMENSION(klon), INTENT(IN)                          :: s2, n2, strig
                                                                 
   REAL, DIMENSION(klon), INTENT(INOUT)                       :: ale_bl_trig, ale_bl
@@ -161 +161 @@
                    if ( (ale_bl_stat(i) .gt. abs(cin(i))+1.e-10) )  then
                       proba_notrig(i)=proba_notrig(i)* &
-                         (1.-exp(-s_trig/s2(i)))**(n2(i)*dtime/tau_trig(i))
+                         (1.-exp(-strig(i)/s2(i)))**(n2(i)*dtime/tau_trig(i))
                       !        print *, 'proba_notrig(i) ',proba_notrig(i)
                       if (random_notrig(i) .ge. proba_notrig(i)) then 
@@ -168 +168 @@
                          ale_bl_trig(i)=0.
                       endif
-                      birth_rate(i) = n2(i)*exp(-s_trig/s2(i))/(tau_trig(i)*cell_area(i))
+                      birth_rate(i) = n2(i)*exp(-strig(i)/s2(i))/(tau_trig(i)*cell_area(i))
 !!!                      birth_rate(i) = max(birth_rate(i),1.e-18)
                    else
@@ -183 +183 @@
                    if ( (Ale_bl(i) .gt. abs(cin(i))+1.e-10) )  then
                       proba_notrig(i)=proba_notrig(i)* &
-                         (1.-exp(-s_trig/s2(i)))**(n2(i)*dtime/tau_trig(i))
+                         (1.-exp(-strig(i)/s2(i)))**(n2(i)*dtime/tau_trig(i))
                       !        print *, 'proba_notrig(i) ',proba_notrig(i)
                       if (random_notrig(i) .ge. proba_notrig(i)) then 
@@ -190 +190 @@
                          ale_bl_trig(i)=0.
                       endif
-                      birth_rate(i) = n2(i)*exp(-s_trig/s2(i))/(tau_trig(i)*cell_area(i))
+                      birth_rate(i) = n2(i)*exp(-strig(i)/s2(i))/(tau_trig(i)*cell_area(i))
 !!!                      birth_rate(i) = max(birth_rate(i),1.e-18)
                    else
@@ -283 +283 @@
         do i=1,klon
            if ( (ale_bl_ref(i) .gt. abs(cin(i))+1.e-10) )  then
-              birth_number = n2(i)*exp(-s_trig/s2(i))
+              birth_number = n2(i)*exp(-strig(i)/s2(i))
               birth_rate(i) = birth_number/(tau_trig(i)*cell_area(i))
 !!!              birth_rate(i) = max(birth_rate(i),1.e-18)