Changeset 5116 for LMDZ6/branches/Amaury_dev/libf/phylmd/lmdz_ratqs_multi.F90

Timestamp:

Jul 24, 2024, 2:54:37 PM (4 months ago)

Author:

abarral

Message:

rename modules properly lmdz_*
move ismin, ismax, minmax into new lmdz_libmath.f90
(lint) uppercase fortran keywords

File:

• LMDZ6/branches/Amaury_dev/libf/phylmd/lmdz_ratqs_multi.F90 (modified) (6 diffs)

LMDZ6/branches/Amaury_dev/libf/phylmd/lmdz_ratqs_multi.F90

@@ -61 +61 @@
 REAL :: taup0 = 70000
 REAL :: taudp = 500
-integer :: lev_out=10
+INTEGER :: lev_out=10
 REAL, DIMENSION (klon,klev) :: zmasse,entr0,detr0,detraincv,dqp,detrain_p,q0,qd0,tau_diss
 REAL, DIMENSION (klon,klev+1) :: fm0
@@ -127 +127 @@
             var_var(i,k) = detr0(i,k)*(varq_therm(i,k)-povariance(i,k))/zmasse(i,k)
             var_det_conv(i,k) =  a_ratqs_cv*(detrain_cv(i,k)/zmasse(i,k))
-            if (sigd(i)/=0) then
+            if (sigd(i)/=0) THEN
                var_desc_prec(i,k) = sigd(i)*(1-sigd(i))*(fqd(i,k)*tau_cumul/sigd(i))**2/tau_cumul
             else
@@ -145 +145 @@
       var_desc_conv(:,klev) = var_desc_conv(:,klev-1)
       
-      if (klein) then
+      if (klein) THEN
          do k=1,klev-1
             do i=1,klon
               qd0(:,:) = 0.0
-              if (sigd(i)/=0) then
+              if (sigd(i)/=0) THEN
                 qd0(i,k) = fqd(i,k)*tau_cumul/sigd(i) 
               endif
@@ -172 +172 @@
                  qd0(i,k) = fqd(i,k)*tau_cumul/sigd(i) 
                endif
-               if (sigt_cv(i,k)/=0) then
+               if (sigt_cv(i,k)/=0) THEN
                  q0(i,k) = fqcomp(i,k)*tau_cumul/sigt_cv(i,k)
                endif
@@ -195 +195 @@
       do k=1,klev
         do i=1,klon
-           if(q_seri(i,k)>=1E-7) then
+           IF(q_seri(i,k)>=1E-7) THEN
                ratqs_inter_(i,k) = abs(povariance(i,k))**0.5/q_seri(i,k)    
                sigma_therm(i,k) = abs(varq_therm(i,k))**0.5     ! sigma dans les thermiques
@@ -205 +205 @@
       enddo
       
-return
+RETURN
 end