Changeset 5144 for LMDZ6/branches/Amaury_dev/libf/phylmd/qcheck.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/qcheck.F90 (modified) (3 diffs)

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

@@ -1 @@
-FUNCTION qcheck(klon,klev,paprs,q,ql,aire)
+FUNCTION qcheck(klon, klev, paprs, q, ql, aire)
+  USE lmdz_yomcst
+
   IMPLICIT NONE
 
@@ -5 +7 @@
   ! la conservation de l'eau
 
-  include "YOMCST.h"
-  INTEGER,INTENT(IN) :: klon,klev
-  REAL,INTENT(IN) :: paprs(klon,klev+1), q(klon,klev), ql(klon,klev)
-  REAL,INTENT(IN) :: aire(klon)
+  INTEGER, INTENT(IN) :: klon, klev
+  REAL, INTENT(IN) :: paprs(klon, klev + 1), q(klon, klev), ql(klon, klev)
+  REAL, INTENT(IN) :: aire(klon)
   REAL qtotal, zx, qcheck
   INTEGER i, k
@@ -14 +15 @@
   zx = 0.0
   DO i = 1, klon
-     zx = zx + aire(i)
+    zx = zx + aire(i)
   ENDDO
   qtotal = 0.0
   DO k = 1, klev
-     DO i = 1, klon
-        qtotal = qtotal + (q(i,k)+ql(i,k)) * aire(i) &
-             *(paprs(i,k)-paprs(i,k+1))/RG
-     ENDDO
+    DO i = 1, klon
+      qtotal = qtotal + (q(i, k) + ql(i, k)) * aire(i) &
+              * (paprs(i, k) - paprs(i, k + 1)) / RG
+    ENDDO
   ENDDO
 
-  qcheck = qtotal/zx
+  qcheck = qtotal / zx
 
 END FUNCTION qcheck