Index: LMDZ6/trunk/libf/phylmd/alpale_th.F90
===================================================================
--- LMDZ6/trunk/libf/phylmd/alpale_th.F90	(revision 3208)
+++ LMDZ6/trunk/libf/phylmd/alpale_th.F90	(revision 3209)
@@ -3,6 +3,5 @@
                        ale_bl_trig, ale_bl_stat, ale_bl,  &
                        alp_bl, alp_bl_stat, &
-                       proba_notrig, random_notrig, birth_rate,  &
-                       q_alp)
+                       proba_notrig, random_notrig, birth_rate)
 
 ! **************************************************************
@@ -44,5 +43,4 @@
 
   REAL, DIMENSION(klon), INTENT(OUT)                         :: birth_rate
-  REAL, DIMENSION(klon), INTENT(OUT)                         :: q_alp
 
   include "thermcell.h"
@@ -146,5 +144,4 @@
                       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.
@@ -152,5 +149,4 @@
                       random_notrig(i)=0.
                       ale_bl_trig(i)=0.
-                      q_alp(i) = 0.
                    endif
                 enddo
@@ -170,5 +166,4 @@
                       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.
@@ -176,5 +171,4 @@
                       random_notrig(i)=0.
                       ale_bl_trig(i)=0.
-                      q_alp(i) = 0.
                    endif
                 enddo
@@ -272,5 +266,4 @@
 !!!              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))/ &
@@ -282,5 +275,4 @@
               random_notrig(i)=0.
               alp_bl(i)=0.
-              q_alp(i) = 0.
            endif
         enddo
@@ -308,6 +300,6 @@
 
           IF (prt_level .GE. 10) THEN
-             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)
+             print *,'alpale_th: ale_bl_trig, alp_bl_stat, birth_rate ', &
+                      ale_bl_trig(1), alp_bl_stat(1), birth_rate(1)
           ENDIF