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

LMDZ6/branches/Amaury_dev/libf/phylmd/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
 
@@ -30 +30 @@
 
 CONTAINS
-!
-!****************************************************************************************
-!
+
+!****************************************************************************************
+
   SUBROUTINE climb_hq_down(knon, coefhq, paprs, pplay, &
        delp, temp, q, dtime, &
@@ -44 +44 @@
 ! 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
 !****************************************************************************************
@@ -90 +90 @@
 ! 1)
 ! Allocation at first time step only
-!   
+
 !****************************************************************************************
 
@@ -132 +132 @@
 ! 2)
 ! Definition of the coeficient K 
-!
+
 !****************************************************************************************
     Kcoefhq(:,:) = 0.0
@@ -146 +146 @@
 ! 3)
 ! Calculation of gama for "Q" and "H"
-!
+
 !****************************************************************************************
 !   surface pressure is used as reference
@@ -180 +180 @@
 ! 4)
 ! Calculte the coefficients C and D for specific humidity, q
-!
+
 !****************************************************************************************
     
@@ -189 +189 @@
 ! 5)
 ! Calculte the coefficients C and D for potentiel entalpie, H 
-!
+
 !****************************************************************************************
     h_old(:,:) = 0.0
@@ -207 +207 @@
 ! 6)
 ! Return the first layer in output variables
-!
+
 !****************************************************************************************
     Acoef_H_out = Acoef_H
@@ -217 +217 @@
 ! 7)
 ! If Pbl is split, return also the other layers in output variables
-!
+
 !****************************************************************************************
 !!! jyg le 07/02/2012
@@ -244 +244 @@
 
   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
@@ -272 +272 @@
 !****************************************************************************************
 ! Niveau au sommet, k=klev
-!
+
 !****************************************************************************************
     Ccoef(:,:) = 0.0
@@ -287 +287 @@
 !****************************************************************************************
 ! Niveau  (klev-1) <= k <= 2
-!
+
 !****************************************************************************************
 
@@ -301 +301 @@
 !****************************************************************************************
 ! Niveau k=1
-!
+
 !****************************************************************************************
 
@@ -311 +311 @@
 
   END SUBROUTINE calc_coef
-!
-!****************************************************************************************
-!
+
+!****************************************************************************************
+
   SUBROUTINE climb_hq_up(knon, dtime, t_old, q_old, &
        flx_q1, flx_h1, paprs, pplay, &
@@ -322 +322 @@
 !!!
        flux_q, flux_h, d_q, d_t)
-! 
+
 ! This routine calculates the flux and tendency of the specific humidity q and 
 ! the potential engergi H. 
@@ -328 +328 @@
 ! 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
@@ -365 +364 @@
 ! Definition of some variables
     REAL, DIMENSION(klon,klev)               :: d_h, zairm
-!
+
 !****************************************************************************************
     flux_q(:,:) = 0.0
@@ -406 +405 @@
 ! 2)
 ! Calculation of Q and H
-!
+
 !****************************************************************************************
 
@@ -423 +422 @@
 ! 3)
 ! Calculation of the flux for Q and H
-!
+
 !****************************************************************************************
 
@@ -444 +443 @@
 ! 4)
 ! Calculation of tendency for Q and H
-!
+
 !****************************************************************************************
     d_h_col_vdf(:) = 0.0
@@ -461 +460 @@
 !****************************************************************************************
 ! Some deallocations
-!
+
 !****************************************************************************************
     IF (last) THEN
@@ -476 +475 @@
     END IF
   END SUBROUTINE climb_hq_up
-!
-!****************************************************************************************
-!
+
+!****************************************************************************************
+
 END MODULE climb_hq_mod