Changeset 5158 for LMDZ6/branches/Amaury_dev/libf/dyn3d/interp_horiz.F90

Timestamp:

Aug 2, 2024, 2:12:03 PM (7 weeks 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/interp_horiz.F90 (modified) (6 diffs)

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

@@ -36 +36 @@
   ! Autres variables
   ! """"""""""""""""
-  REAL :: airetest(imn + 1, jmn + 1)
   INTEGER :: ii, jj, l
 
@@ -45 +44 @@
   INTEGER :: iik(kllm), jjk(kllm), jk(kllm), ik(kllm)
   REAL :: intersec(kllm)
-  REAL :: R
   REAL :: totn, tots
-
-  LOGICAL :: firstcall, firsttest, aire_ok
-  save firsttest
-  data firsttest /.TRUE./
-  data aire_ok /.TRUE./
-
-
-
-
 
   ! initialisation
@@ -68 +57 @@
           ktotal, iik, jjk, jk, ik, intersec, airen)
 
-  do l = 1, lm
-    do jj = 1, jmn + 1
-      do ii = 1, imn + 1
+  DO l = 1, lm
+    DO jj = 1, jmn + 1
+      DO ii = 1, imn + 1
         varn(ii, jj, l) = 0.
       END DO
@@ -83 +72 @@
   PRINT *, 'ktotal 1 = ', ktotal
 
-  do k = 1, ktotal
-    do l = 1, lm
+  DO k = 1, ktotal
+    DO l = 1, lm
       varn(iik(k), jjk(k), l) = varn(iik(k), jjk(k), l) &
               + varo(ik(k), jk(k), l) * intersec(k) / airen(iik(k), jjk(k))
@@ -91 +80 @@
 
   ! Une seule valeur au pole pour les variables ! :
-  ! -----------------------------------------------
-  do l = 1, lm
+  DO l = 1, lm
     totn = 0.
     tots = 0.
-    do ii = 1, imn + 1
+    DO ii = 1, imn + 1
       totn = totn + varn(ii, 1, l)
       tots = tots + varn (ii, jmn + 1, l)
     END DO
-    do ii = 1, imn + 1
+    DO ii = 1, imn + 1
       varn(ii, 1, l) = totn / REAL(imn + 1)
       varn(ii, jmn + 1, l) = tots / REAL(imn + 1)
@@ -105 +93 @@
   END DO
 
-
-  !---------------------------------------------------------------
-  !  TEST  TEST  TEST  TEST  TEST  TEST  TEST  TEST  TEST  TEST
-  !!       if (.NOT.(firsttest)) goto 99
-  !!       firsttest = .FALSE.
-  !! !     write (*,*) 'INTERP. HORIZ. : TEST SUR LES AIRES:'
-  !!       do jj =1 , jmn+1
-  !!         do ii=1, imn+1
-  !!           airetest(ii,jj) =0.
-  !!         END DO
-  !!       END DO
-  !!       PRINT *, 'ktotal = ', ktotal
-  !!       PRINT *, 'jmn+1 =', jmn+1, 'imn+1', imn+1
-  !!
-  !!       do k=1,ktotal
-  !!          airetest(iik(k),jjk(k))= airetest(iik(k),jjk(k)) +intersec(k)
-  !!       end DO
-  !!
-  !!
-  !!       PRINT *, 'fin boucle'
-  !!       do jj =1 , jmn+1
-  !!        do ii=1, imn+1
-  !!          r = airen(ii,jj)/airetest(ii,jj)
-  !!          if ((r.gt.1.001).OR.(r.lt.0.999)) THEN
-  !! !             write (*,*) '********** PROBLEME D'' AIRES !!!',
-  !! !     &                   ' DANS L''INTERPOLATION HORIZONTALE'
-  !! !             WRITE(*,*)'ii,jj,airen,airetest',
-  !! !     &          ii,jj,airen(ii,jj),airetest(ii,jj)
-  !!              aire_ok = .FALSE.
-  !!          end if
-  !!        END DO
-  !!       END DO
-  !! !      if (aire_ok) WRITE(*,*) 'INTERP. HORIZ. : AIRES OK'
-  !!  99   continue
-
-  ! FIN TEST  FIN TEST  FIN TEST  FIN TEST  FIN TEST  FIN TEST  FIN TEST
-  !---------------------------------------------------------------
-
-
 END SUBROUTINE  interp_horiz