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

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

@@ -1 @@
-
 ! $Header$
 
@@ -35 +34 @@
 
         ! aire totale
-        airetot = airetot + aire(i)*msk(i)
+        airetot = airetot + aire(i) * msk(i)
 
         ! ponderation par la masse
-        moyglo = moyglo + champ(i)*aire(i)*msk(i)
+        moyglo = moyglo + champ(i) * aire(i) * msk(i)
       END IF
     END DO
@@ -49 +48 @@
 
       ! ponderation par la masse
-      moyglo = moyglo + champ(i)*aire(i)
+      moyglo = moyglo + champ(i) * aire(i)
     END DO
 
@@ -55 +54 @@
 
   ! moyenne ponderee par l'aire
-  moyglo = moyglo/airetot
-
+  moyglo = moyglo / airetot
 
 END SUBROUTINE moyglo_pondaire
@@ -62 +60 @@
 SUBROUTINE moyglo_pondaima(nhori, nvert, champ, aire, pbord, moyglo)
   USE dimphy
+  USE lmdz_yomcst
+
   IMPLICIT NONE
   ! ==================================================================
@@ -72 +72 @@
 
   ! ==================================================================
-  include "YOMCST.h"
   INTEGER i, k, nhori, nvert
   REAL champ(klon, klev), aire(klon)
-  REAL pbord(klon, klev+1)
+  REAL pbord(klon, klev + 1)
   REAL moyglo
 
@@ -87 +86 @@
   DO k = 1, nvert
     DO i = 1, nhori
-      moyglo = moyglo + champ(i, k)*(pbord(i,k)-pbord(i,k+1))/rg*aire(i)
+      moyglo = moyglo + champ(i, k) * (pbord(i, k) - pbord(i, k + 1)) / rg * aire(i)
     END DO
   END DO
@@ -98 +97 @@
 
   ! moyenne par mettre carre avec ponderation par la masse
-  moyglo = moyglo/airetot
-
+  moyglo = moyglo / airetot
 
 END SUBROUTINE moyglo_pondaima
@@ -105 +103 @@
 SUBROUTINE moyglo_pondmass(nhori, nvert, champ, aire, pbord, moyglo)
   USE dimphy
+  USE lmdz_yomcst
+
   IMPLICIT NONE
   ! ==================================================================
@@ -115 +115 @@
 
   ! ==================================================================
-  include "YOMCST.h"
   INTEGER i, k, nhori, nvert
   REAL champ(klon, klev), aire(klon)
-  REAL pbord(klon, klev+1)
+  REAL pbord(klon, klev + 1)
   REAL moyglo
 
@@ -130 +129 @@
   DO k = 1, nvert
     DO i = 1, nhori
-      moyglo = moyglo + champ(i, k)*(pbord(i,k)-pbord(i,k+1))/rg*aire(i)
+      moyglo = moyglo + champ(i, k) * (pbord(i, k) - pbord(i, k + 1)) / rg * aire(i)
     END DO
   END DO
@@ -138 +137 @@
   DO k = 1, nvert
     DO i = 1, nhori
-      massetot = massetot + (pbord(i,k)-pbord(i,k+1))/rg*aire(i)
+      massetot = massetot + (pbord(i, k) - pbord(i, k + 1)) / rg * aire(i)
     END DO
   END DO
 
   ! moyenne par mettre carre avec ponderation par la masse
-  moyglo = moyglo/massetot
-
+  moyglo = moyglo / massetot
 
 END SUBROUTINE moyglo_pondmass