Changeset 3723 for LMDZ6/trunk/libf/phylmd/cospv2/cosp_utils.F90

Timestamp:

Jun 18, 2020, 10:44:07 AM (4 years ago)

Author:

idelkadi

Message:

Debugging COSP v2 for simulators Calipso, Parasol, Cloudsat

File:

• LMDZ6/trunk/libf/phylmd/cospv2/cosp_utils.F90 (modified) (2 diffs)

LMDZ6/trunk/libf/phylmd/cospv2/cosp_utils.F90

@@ -56 +56 @@
     integer :: i,j,k
     real(wp) :: sigma,one_over_xip1,xi,rho0,rho,lambda_x,gamma_4_3_2,delta
+
+    real(wp) :: seuil
+
+    if (ok_debug_cosp) then
+       seuil=1.e-15
+    else
+       seuil=0.0
+    endif
+
+
     
     mxratio = 0.0
@@ -75 +85 @@
                         mxratio(i,j,k)=mxratio(i,j,k)/rho
                         ! Compute effective radius
-                        if ((reff(i,j,k) <= 0._wp).and.(flux(i,k) /= 0._wp)) then
+!                        if ((reff(i,j,k) <= 0._wp).and.(flux(i,k) /= 0._wp)) then
+                         if ((reff(i,j,k) <= 0._wp).and.(flux(i,k) > seuil)) then
                            lambda_x = (a_x*c_x*((rho0/rho)**g_x)*n_ax*gamma1/flux(i,k))**(1._wp/delta)
                            reff(i,j,k) = gamma_4_3_2/lambda_x