Changeset 5099 for LMDZ6/branches/Amaury_dev/libf/phylmdiso/climb_hq_mod.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/phylmdiso/climb_hq_mod.F90 (modified) (23 diffs)

LMDZ6/branches/Amaury_dev/libf/phylmdiso/climb_hq_mod.F90

@@ -1 +1 @@
 MODULE climb_hq_mod
-!
+
 ! Module to solve the verctical diffusion of "q" and "H"; 
 ! specific humidity and potential energi.
-!
+
   USE dimphy
 #ifdef ISO
@@ -42 +42 @@
 
 CONTAINS
-!
-!****************************************************************************************
-!
+
+!****************************************************************************************
+
   SUBROUTINE climb_hq_down(knon, coefhq, paprs, pplay, &
        delp, temp, q, dtime, &
@@ -67 +67 @@
 ! for the quantity X=[Q,H] in equation X(k) = C(k) + D(k)*X(k-1), where k is
 ! the index of the vertical layer.
-!
+
 ! Input arguments
 !****************************************************************************************
@@ -147 +147 @@
 ! 1)
 ! Allocation at first time step only
-!   
+
 !****************************************************************************************
 
@@ -203 +203 @@
 ! 2)
 ! Definition of the coeficient K 
-!
+
 !****************************************************************************************
     Kcoefhq(:,:) = 0.0
@@ -217 +217 @@
 ! 3)
 ! Calculation of gama for "Q" and "H"
-!
+
 !****************************************************************************************
 !   surface pressure is used as reference
@@ -278 +278 @@
 ! 4)
 ! Calculte the coefficients C and D for specific humidity, q
-!
+
 !****************************************************************************************
     
@@ -340 +340 @@
 ! 5)
 ! Calculte the coefficients C and D for potentiel entalpie, H 
-!
+
 !****************************************************************************************
     h_old(:,:) = 0.0
@@ -358 +358 @@
 ! 6)
 ! Return the first layer in output variables
-!
+
 !****************************************************************************************
     Acoef_H_out = Acoef_H
@@ -372 +372 @@
 ! 7)
 ! If Pbl is split, return also the other layers in output variables
-!
+
 !****************************************************************************************
 !!! jyg le 07/02/2012
@@ -415 +415 @@
 
   END SUBROUTINE climb_hq_down
-!
-!****************************************************************************************
-!
+
+!****************************************************************************************
+
   SUBROUTINE calc_coef(knon, Kcoef, gama, delp, X, Ccoef, Dcoef, Acoef, Bcoef)
-!
+
 ! Calculate the coefficients C and D in : X(k) = C(k) + D(k)*X(k-1)
 ! where X is H or Q, and k the vertical level k=1,klev
-!
+
     INCLUDE "YOMCST.h"
 ! Input arguments
@@ -443 +443 @@
 !****************************************************************************************
 ! Niveau au sommet, k=klev
-!
+
 !****************************************************************************************
     Ccoef(:,:) = 0.0
@@ -458 +458 @@
 !****************************************************************************************
 ! Niveau  (klev-1) <= k <= 2
-!
+
 !****************************************************************************************
 
@@ -472 +472 @@
 !****************************************************************************************
 ! Niveau k=1
-!
+
 !****************************************************************************************
 
@@ -482 +482 @@
 
   END SUBROUTINE calc_coef
-!
-!****************************************************************************************
-!
+
+!****************************************************************************************
+
   SUBROUTINE climb_hq_up(knon, dtime, t_old, q_old, &
        flx_q1, flx_h1, paprs, pplay, &
@@ -505 +505 @@
 USE isotopes_verif_mod
 #endif
-! 
+
 ! This routine calculates the flux and tendency of the specific humidity q and 
 ! the potential engergi H. 
@@ -511 +511 @@
 ! X(k) = C(k) + D(k)*X(k-1) for X=[q,H], where the coefficients 
 ! C and D are known from before and k is index of the vertical layer.
-!   
 
 ! Input arguments
@@ -564 +563 @@
 ! Definition of some variables
     REAL, DIMENSION(klon,klev)               :: d_h, zairm
-!
+
 !****************************************************************************************
 
@@ -645 +644 @@
 ! 2)
 ! Calculation of Q and H
-!
+
 !****************************************************************************************
 
@@ -715 +714 @@
 ! 3)
 ! Calculation of the flux for Q and H
-!
+
 !****************************************************************************************
 
@@ -748 +747 @@
 ! 4)
 ! Calculation of tendency for Q and H
-!
+
 !****************************************************************************************
     d_h_col_vdf(:) = 0.0
@@ -798 +797 @@
 !****************************************************************************************
 ! Some deallocations
-!
+
 !****************************************************************************************
     IF (last) THEN
@@ -822 +821 @@
     END IF
   END SUBROUTINE climb_hq_up
-!
-!****************************************************************************************
-!
+
+!****************************************************************************************
+
 END MODULE climb_hq_mod