Changeset 5099 for LMDZ6/branches/Amaury_dev/libf/phylmd/alpale_th.F90

Timestamp:

Jul 22, 2024, 9:29:09 PM (4 months ago)

Author:

abarral

Message:

Replace most uses of CPP_DUST by the corresponding logical defined in lmdz_cppkeys_wrapper.F90
Convert several files from .F to .f90 to allow Dust to compile w/o rrtm/ecrad
Create lmdz_yoerad.f90
(lint) Remove "!" on otherwise empty line

File:

• LMDZ6/branches/Amaury_dev/libf/phylmd/alpale_th.F90 (modified) (14 diffs)

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

@@ -1 @@
-!
+
 ! $Id$
-!
+
 SUBROUTINE alpale_th ( dtime, lmax_th, t_seri, cell_area,  &
                        cin, s2, n2, strig,  &
@@ -20 +20 @@
   USE ioipsl_getin_p_mod, ONLY : getin_p
   USE print_control_mod, ONLY: mydebug=>debug , lunout, prt_level
-!
+
   IMPLICIT NONE
 
@@ -59 +59 @@
   REAL, DIMENSION(klon)                                      :: ale_bl_ref
   REAL, DIMENSION(klon)                                      :: tau_trig
-!
+
     !$OMP THREADPRIVATE(multiply_proba_notrig) 
     !$OMP THREADPRIVATE(random_notrig_max) 
     !$OMP THREADPRIVATE(cv_feed_area) 
     !$OMP THREADPRIVATE(first)
-!
+
  REAL umexp  ! expression of (1.-exp(-x))/x valid for all x, especially when x->0
  REAL x
-!
+
      CHARACTER (LEN=20) :: modname='alpale_th'
      CHARACTER (LEN=80) :: abort_message
@@ -75 +75 @@
 !!!            (1.-max(sign(1.,x-1.e-3),0.))*(-0.5*x*(1.-x/3.*(1.-0.25*x))) !!! bug introduced by mistake  (jyg)
 !!!            (1.-max(sign(1.,x-1.e-3),0.))*(1.-0.5*x*(1.-x/3.*(1.-0.25*x)))  !!! initial correct formula (jyg)
-!
+
 !+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
 !  JYG, 20160513 : Introduction of the Effective Lifting Power (ELP), which 
@@ -83 +83 @@
 ! zero. It is activated by iflag_clos_bl = 3.
 !   The ELP values are stored in the ALP_bl variable.
-!   
+
 !+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
-!
 
     IF (first) THEN
@@ -122 +121 @@
   IF (iflag_clos_bl < 3) THEN
 !---------------------------------------
-!
+
 !      Original code (Nicolas Rochetin)
 !     --------------------------------
@@ -128 +127 @@
           !-----------Stochastic triggering-----------
           if (iflag_trig_bl>=1) then
-             !
+
              IF (prt_level >= 10) THEN
                 WRITE(lunout,*)'cin, ale_bl_stat, alp_bl, alp_bl_stat ', &
@@ -146 +145 @@
                 endif
              enddo
-             !
+
              IF (prt_level >= 10) THEN
                 WRITE(lunout,*)'random_notrig, tau_trig ', &
@@ -202 +201 @@
              ENDIF
 
-             !
              IF (prt_level >= 10) THEN
                 WRITE(lunout,*)'proba_notrig, ale_bl_trig ', &
@@ -233 +231 @@
           endif !(iflag_clos_bl)
 
-!
 !---------------------------------------
   ELSEIF (iflag_clos_bl == 3) THEN  ! (iflag_clos_bl .LT. 3)
 !---------------------------------------
-!
+
 !      New code with Effective Lifting Power
 !     -------------------------------------
@@ -243 +240 @@
           !-----------Stochastic triggering-----------
      if (iflag_trig_bl>=1) then
-        !
+
         IF (prt_level >= 10) THEN
            WRITE(lunout,*)'cin, ale_bl_stat, alp_bl_stat ', &
@@ -272 +269 @@
            endif
         enddo
-        !
+
         IF (prt_level >= 10) THEN
            WRITE(lunout,*)'random_notrig, tau_trig ', &
@@ -308 +305 @@
          enddo
 
-        !
         IF (prt_level >= 10) THEN
            WRITE(lunout,*)'proba_notrig, ale_bl_trig ', &
@@ -326 +322 @@
 
           !cc fin nrlmd le 10/04/2012
-!
+
           !IM/FH: 2011/02/23 
           ! Couplage Thermiques/Emanuel seulement si T<0