Changeset 5158 for LMDZ6/branches/Amaury_dev/libf/phylmd/yamada_c.F90

Timestamp:

Aug 2, 2024, 2:12:03 PM (3 months ago)

Author:

abarral

Message:

Add missing klon on strataer_emiss_mod.F90
Correct various missing explicit declarations
Replace tabs by spaces (tabs are not part of the fortran charset)
Continue cleaning modules
Removed unused arguments and variables

File:

• LMDZ6/branches/Amaury_dev/libf/phylmd/yamada_c.F90 (modified) (9 diffs)

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

@@ -162 +162 @@
   zt(:, :) = pt(:, :) + zalpha * d_t(:, :)
 
-  do k = 1, klev
+  DO k = 1, klev
     exner(:, k) = (play(:, k) / plev(:, 1))**RKAPPA
     masse(:, k) = (plev(:, k) - plev(:, k + 1)) / RG
@@ -170 +170 @@
   ! Atmospheric mass at layer interfaces, where the TKE is computed
   masseb(:, :) = 0.
-  do k = 1, klev
+  DO k = 1, klev
     masseb(:, k) = masseb(:, k) + masse(:, k)
     masseb(:, k + 1) = masseb(:, k + 1) + masse(:, k)
@@ -178 +178 @@
   zlev(:, 1) = 0.
   zlay(:, 1) = RCPD * teta(:, 1) * (1. - exner(:, 1))
-  do k = 1, klev - 1
+  DO k = 1, klev - 1
     zlay(:, k + 1) = zlay(:, k) + 0.5 * RCPD * (teta(:, k) + teta(:, k + 1)) * (exner(:, k) - exner(:, k + 1)) / RG
     zlev(:, k) = 0.5 * (zlay(:, k) + zlay(:, k + 1)) ! PASBO
@@ -187 +187 @@
   fluxt(:, klev + 1) = 0.
 
-  do k = klev, 1, -1
+  DO k = klev, 1, -1
     fluxu(:, k) = fluxu(:, k + 1) + masse(:, k) * d_u(:, k)
     fluxv(:, k) = fluxv(:, k + 1) + masse(:, k) * d_v(:, k)
@@ -197 +197 @@
   dddt(:, 1) = (exner(:, 1) - 1.) * fluxt(:, 1)
 
-  do k = 2, klev
+  DO k = 2, klev
     dddu(:, k) = (zu(:, k) - zu(:, k - 1)) * fluxu(:, k)
     dddv(:, k) = (zv(:, k) - zv(:, k - 1)) * fluxv(:, k)
@@ -263 +263 @@
   !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
 
-  do k = 2, klev
-    do ig = 1, ngrid
+  DO k = 2, klev
+    DO ig = 1, ngrid
       dz(ig, k) = zlay(ig, k) - zlay(ig, k - 1)
       m2(ig, k) = ((zu(ig, k) - zu(ig, k - 1))**2 + (zv(ig, k) - zv(ig, k - 1))**2) / (dz(ig, k) * dz(ig, k))
@@ -301 +301 @@
     !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
 
-    do ig = 1, ngrid
+    DO ig = 1, ngrid
       sq(ig) = 1.e-10
       sqz(ig) = 1.e-10
     enddo
-    do k = 2, klev - 1
-      do ig = 1, ngrid
+    DO k = 2, klev - 1
+      DO ig = 1, ngrid
         zq = sqrt(q2(ig, k))
         sqz(ig) = sqz(ig) + zq * zlev(ig, k) * (zlay(ig, k) - zlay(ig, k - 1))
@@ -312 +312 @@
       enddo
     enddo
-    do ig = 1, ngrid
+    DO ig = 1, ngrid
       l0(ig) = 0.2 * sqz(ig) / sq(ig)
     enddo
-    do k = 2, klev
-      do ig = 1, ngrid
+    DO k = 2, klev
+      DO ig = 1, ngrid
         l(ig, k) = fl(zlev(ig, k), l0(ig), q2(ig, k), n2(ig, k))
       enddo
@@ -329 +329 @@
 
     l0(:) = 150.
-    do k = 2, klev
-      do ig = 1, ngrid
+    DO k = 2, klev
+      DO ig = 1, ngrid
         l(ig, k) = fl(zlev(ig, k), l0(ig), q2(ig, k), n2(ig, k))
       enddo