Changeset 5144 for LMDZ6/branches/Amaury_dev/libf/phylmd/transp_lay.F90

Timestamp:

Jul 29, 2024, 11:01:04 PM (4 months ago)

Author:

abarral

Message:

Put YOMCST.h into modules

File:

• LMDZ6/branches/Amaury_dev/libf/phylmd/transp_lay.F90 (modified) (3 diffs)

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

@@ -1 @@
-
 ! $Header$
 
 SUBROUTINE transp_lay(paprs, tsol, t, q, u, v, geom, vtran_e, vtran_q, &
-    utran_e, utran_q)
+        utran_e, utran_q)
 
   USE dimphy
+  USE lmdz_yomcst
+
   IMPLICIT NONE
   ! ======================================================================
@@ -13 +14 @@
   ! ======================================================================
 
-  include "YOMCST.h"
-
-  REAL paprs(klon, klev+1), tsol(klon)
+  REAL paprs(klon, klev + 1), tsol(klon)
   REAL t(klon, klev), q(klon, klev), u(klon, klev), v(klon, klev)
   REAL utran_e(klon, klev), utran_q(klon, klev)
@@ -35 +34 @@
   DO l = 1, klev
     DO i = 1, klon
-      esh = rcpd*t(i, l) + rlvtt*q(i, l) + geom(i, l)
-      utran_e(i, l) = utran_e(i, l) + u(i, l)*esh*(paprs(i,l)-paprs(i,l+1))/ &
-        rg
-      utran_q(i, l) = utran_q(i, l) + u(i, l)*q(i, l)*(paprs(i,l)-paprs(i,l+1 &
-        ))/rg
-      vtran_e(i, l) = vtran_e(i, l) + v(i, l)*esh*(paprs(i,l)-paprs(i,l+1))/ &
-        rg
-      vtran_q(i, l) = vtran_q(i, l) + v(i, l)*q(i, l)*(paprs(i,l)-paprs(i,l+1 &
-        ))/rg
+      esh = rcpd * t(i, l) + rlvtt * q(i, l) + geom(i, l)
+      utran_e(i, l) = utran_e(i, l) + u(i, l) * esh * (paprs(i, l) - paprs(i, l + 1)) / &
+              rg
+      utran_q(i, l) = utran_q(i, l) + u(i, l) * q(i, l) * (paprs(i, l) - paprs(i, l + 1 &
+              )) / rg
+      vtran_e(i, l) = vtran_e(i, l) + v(i, l) * esh * (paprs(i, l) - paprs(i, l + 1)) / &
+              rg
+      vtran_q(i, l) = vtran_q(i, l) + v(i, l) * q(i, l) * (paprs(i, l) - paprs(i, l + 1 &
+              )) / rg
     END DO
   END DO
 
-
 END SUBROUTINE transp_lay