Changeset 3241 for LMDZ6/trunk/libf/phylmd/cosp/mod_cosp_utils.F90

Timestamp:

Mar 6, 2018, 10:46:13 AM (7 years ago)

Author:

idelkadi

Message:

• Corrections pour tourner en mode debug, controle par la cles ok_debug_cosp dans mod_cosp_constants.F90 la valeur par defaut etant .false.
• Rajouter "implicit none" pour les routines ne le comportant pas
• Nettoyage du code (suppression des commentaires unitils)

File:

• LMDZ6/trunk/libf/phylmd/cosp/mod_cosp_utils.F90 (modified) (2 diffs)

LMDZ6/trunk/libf/phylmd/cosp/mod_cosp_utils.F90

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