Changeset 3208 for LMDZ6/trunk/libf/phylmd/alpale_th.F90

Timestamp:

Feb 16, 2018, 12:42:18 PM (7 years ago)

Author:

jyg

Message:

Implementation of a first crude model of the
dynamic of wake population.

File:

• LMDZ6/trunk/libf/phylmd/alpale_th.F90 (modified) (9 diffs)

LMDZ6/trunk/libf/phylmd/alpale_th.F90

@@ -3 +3 @@
                        ale_bl_trig, ale_bl_stat, ale_bl,  &
                        alp_bl, alp_bl_stat, &
-                       proba_notrig, random_notrig)
+                       proba_notrig, random_notrig, birth_rate,  &
+                       q_alp)
 
 ! **************************************************************
@@ -42 +43 @@
   REAL, DIMENSION(klon), INTENT(OUT)                         :: random_notrig
 
+  REAL, DIMENSION(klon), INTENT(OUT)                         :: birth_rate
+  REAL, DIMENSION(klon), INTENT(OUT)                         :: q_alp
+
   include "thermcell.h"
 
@@ -53 +57 @@
   REAL, DIMENSION(klon)                                      :: ale_bl_ref
   REAL, DIMENSION(klon)                                      :: tau_trig
-  REAL, DIMENSION(klon)                                      :: birth_rate
 !
     !$OMP THREADPRIVATE(random_notrig_max) 
@@ -103 +106 @@
              !
              IF (prt_level .GE. 10) THEN
-                print *,'cin, ale_bl_stat, alp_bl_stat ', &
-                     cin, ale_bl_stat, alp_bl_stat
+                print *,'cin, ale_bl_stat, alp_bl, alp_bl_stat ', &
+                     cin, ale_bl_stat, alp_bl, alp_bl_stat
              ENDIF
 
@@ -141 +144 @@
                          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) = max(birth_rate(i),1.e-18)
+                      q_alp(i) = alp_bl(i)/max(birth_rate(i),1.e-18)
                    else
 !!jyg                      proba_notrig(i)=1.
+                      birth_rate(i) = 0.
                       random_notrig(i)=0.
                       ale_bl_trig(i)=0.
+                      q_alp(i) = 0.
                    endif
                 enddo
@@ -160 +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) = max(birth_rate(i),1.e-18)
+                      q_alp(i) = alp_bl(i)/max(birth_rate(i),1.e-18)
                    else
 !!jyg                      proba_notrig(i)=1.
+                      birth_rate(i) = 0.
                       random_notrig(i)=0.
                       ale_bl_trig(i)=0.
+                      q_alp(i) = 0.
                    endif
                 enddo
@@ -257 +270 @@
               birth_number = n2(i)*exp(-s_trig/s2(i))
               birth_rate(i) = birth_number/(tau_trig(i)*cell_area(i))
+!!!              birth_rate(i) = max(birth_rate(i),1.e-18)
               proba_notrig(i)=proba_notrig(i)*exp(-birth_number*dtime/tau_trig(i))
+              q_alp(i) = alp_bl(i)/max(birth_rate(i),1.e-18)
               Alp_bl(i) = Alp_bl(i)* &
                           umexp(-birth_number*cv_feed_area/cell_area(i))/ &
@@ -264 +279 @@
           else 
 !!jyg              proba_notrig(i)=1.
+              birth_rate(i)=0.
               random_notrig(i)=0.
               alp_bl(i)=0.
+              q_alp(i) = 0.
            endif
         enddo
@@ -291 +308 @@
 
           IF (prt_level .GE. 10) THEN
-             print *,'ale_bl_trig, alp_bl_stat ',ale_bl_trig, alp_bl_stat
+             print *,'alpale_th: ale_bl_trig, alp_bl_stat, birth_rate, q_alp ', &
+                      ale_bl_trig(1), alp_bl_stat(1), birth_rate(1), q_alp(1)
           ENDIF