Changeset 5158 for LMDZ6/branches/Amaury_dev/libf/dyn3d/friction.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/dyn3d/friction.F90 (modified) (7 diffs)

LMDZ6/branches/Amaury_dev/libf/dyn3d/friction.F90

@@ -58 +58 @@
   IF (friction_type==0) THEN
     !   calcul des composantes au carre du vent naturel
-    do j = 1, jjp1
-      do i = 1, iip1
+    DO j = 1, jjp1
+      DO i = 1, iip1
         u2(i, j) = ucov(i, j, 1) * ucov(i, j, 1) * unscu2(i, j)
       enddo
     enddo
-    do j = 1, jjm
-      do i = 1, iip1
+    DO j = 1, jjm
+      DO i = 1, iip1
         v2(i, j) = vcov(i, j, 1) * vcov(i, j, 1) * unscv2(i, j)
       enddo
@@ -70 +70 @@
 
     !   calcul du module de V en dehors des poles
-    do j = 2, jjm
-      do i = 2, iip1
+    DO j = 2, jjm
+      DO i = 2, iip1
         modv(i, j) = sqrt(0.5 * (u2(i - 1, j) + u2(i, j) + v2(i, j - 1) + v2(i, j)))
       enddo
@@ -83 +83 @@
     upols = 0.
     vpols = 0.
-    do i = 2, iip1
+    DO i = 2, iip1
       zco = cos(rlonv(i)) * (rlonu(i) - rlonu(i - 1))
       zsi = sin(rlonv(i)) * (rlonu(i) - rlonu(i - 1))
@@ -95 +95 @@
     vpn = sqrt(upoln * upoln + vpoln * vpoln) / pi
     vps = sqrt(upols * upols + vpols * vpols) / pi
-    do i = 1, iip1
+    DO i = 1, iip1
       ! modv(i,1)=vpn
       ! modv(i,jjp1)=vps
@@ -103 +103 @@
 
     !   calcul du frottement au sol.
-    do j = 2, jjm
-      do i = 1, iim
+    DO j = 2, jjm
+      DO i = 1, iim
         ucov(i, j, 1) = ucov(i, j, 1) &
                 - cfric * pdt * 0.5 * (modv(i + 1, j) + modv(i, j)) * ucov(i, j, 1)
@@ -110 +110 @@
       ucov(iip1, j, 1) = ucov(1, j, 1)
     enddo
-    do j = 1, jjm
-      do i = 1, iip1
+    DO j = 1, jjm
+      DO i = 1, iip1
         vcov(i, j, 1) = vcov(i, j, 1) &
                 - cfric * pdt * 0.5 * (modv(i, j + 1) + modv(i, j)) * vcov(i, j, 1)
@@ -120 +120 @@
 
   IF (friction_type==1) THEN
-    do l = 1, llm
+    DO l = 1, llm
       ucov(:, :, l) = ucov(:, :, l) * (1. - pdt * kfrict(l))
       vcov(:, :, l) = vcov(:, :, l) * (1. - pdt * kfrict(l))