Changeset 3209 for LMDZ6/trunk/libf

Timestamp:

Feb 16, 2018, 12:50:44 PM (7 years ago)

Author:

jyg

Message:

error in the prvious commit

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, birth_rate,  &
-                       q_alp)
+                       proba_notrig, random_notrig, birth_rate)
 
 ! **************************************************************
@@ -44 +43 @@
 
   REAL, DIMENSION(klon), INTENT(OUT)                         :: birth_rate
-  REAL, DIMENSION(klon), INTENT(OUT)                         :: q_alp
 
   include "thermcell.h"
@@ -146 +144 @@
                       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 +149 @@
                       random_notrig(i)=0.
                       ale_bl_trig(i)=0.
-                      q_alp(i) = 0.
                    endif
                 enddo
@@ -170 +166 @@
                       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 +171 @@
                       random_notrig(i)=0.
                       ale_bl_trig(i)=0.
-                      q_alp(i) = 0.
                    endif
                 enddo
@@ -272 +266 @@
 !!!              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 +275 @@
               random_notrig(i)=0.
               alp_bl(i)=0.
-              q_alp(i) = 0.
            endif
         enddo
@@ -308 +300 @@
 
           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