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

Timestamp:

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

Author:

abarral

Message:

Put YOMCST.h into modules

File:

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

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

@@ -1 @@
-
 ! $Id$
 
-SUBROUTINE transp(paprs, tsol, t, q, ql, qs, u, v, geom, & 
-                  utran_e, vtran_e, utran_q, vtran_q, utran_w, vtran_w)
+SUBROUTINE transp(paprs, tsol, t, q, ql, qs, u, v, geom, &
+        utran_e, vtran_e, utran_q, vtran_q, utran_w, vtran_w)
 
   USE dimphy
+  USE lmdz_yomcst
+
   IMPLICIT NONE
   ! ======================================================================
@@ -13 +14 @@
   ! ======================================================================
 
-  include "YOMCST.h"
-
   !--inputs
-  REAL, INTENT(IN)  :: paprs(klon, klev+1), tsol(klon), geom(klon, klev)
-  REAL, INTENT(IN)  :: t(klon, klev), q(klon, klev), ql(klon, klev), qs(klon, klev)
-  REAL, INTENT(IN)  :: u(klon, klev), v(klon, klev)
+  REAL, INTENT(IN) :: paprs(klon, klev + 1), tsol(klon), geom(klon, klev)
+  REAL, INTENT(IN) :: t(klon, klev), q(klon, klev), ql(klon, klev), qs(klon, klev)
+  REAL, INTENT(IN) :: u(klon, klev), v(klon, klev)
   !--outputs
   REAL, INTENT(OUT) :: utran_e(klon), vtran_e(klon) !--lateral flux of dry static energy (J m-1 s-1)
@@ -39 +38 @@
   DO l = 1, klev
     DO i = 1, klon
-      dm= (paprs(i,l)-paprs(i,l+1))/RG  !--mass of layer kg m-2
+      dm = (paprs(i, l) - paprs(i, l + 1)) / RG  !--mass of layer kg m-2
       !--moist static energy
-!      e = rcpd*t(i, l) + rlvtt*q(i, l) + geom(i, l)
+      !      e = rcpd*t(i, l) + rlvtt*q(i, l) + geom(i, l)
       !--dry static energy
-      e = rcpd*t(i, l) + geom(i, l)
-      utran_e(i) = utran_e(i) + u(i, l)*e*dm
-      vtran_e(i) = vtran_e(i) + v(i, l)*e*dm
+      e = rcpd * t(i, l) + geom(i, l)
+      utran_e(i) = utran_e(i) + u(i, l) * e * dm
+      vtran_e(i) = vtran_e(i) + v(i, l) * e * dm
       !--water vapour
-      utran_q(i) = utran_q(i) + u(i, l)*q(i,l)*dm 
-      vtran_q(i) = vtran_q(i) + v(i, l)*q(i,l)*dm
+      utran_q(i) = utran_q(i) + u(i, l) * q(i, l) * dm
+      vtran_q(i) = vtran_q(i) + v(i, l) * q(i, l) * dm
       !--total water
-      utran_w(i) = utran_w(i) + u(i, l)*(q(i,l)+ql(i,l)+qs(i,l))*dm 
-      vtran_w(i) = vtran_w(i) + v(i, l)*(q(i,l)+ql(i,l)+qs(i,l))*dm
+      utran_w(i) = utran_w(i) + u(i, l) * (q(i, l) + ql(i, l) + qs(i, l)) * dm
+      vtran_w(i) = vtran_w(i) + v(i, l) * (q(i, l) + ql(i, l) + qs(i, l)) * dm
     ENDDO
   ENDDO
 
-
 END SUBROUTINE transp