Changeset 5159 for LMDZ6/branches/Amaury_dev/libf/filtrez/lmdz_filtreg.F90

Timestamp:

Aug 2, 2024, 9:58:25 PM (7 weeks ago)

Author:

abarral

Message:

Put dimensions.h and paramet.h into modules

File:

• LMDZ6/branches/Amaury_dev/libf/filtrez/lmdz_filtreg.F90 (modified) (13 diffs)

LMDZ6/branches/Amaury_dev/libf/filtrez/lmdz_filtreg.F90

@@ -2 +2 @@
 
 MODULE lmdz_filtreg
+  USE lmdz_paramet
   IMPLICIT NONE; PRIVATE
   PUBLIC matriceun, matriceus, matricevn, matricevs, matrinvn, matrinvs, &
@@ -13 +14 @@
   SUBROUTINE filtreg(champ, nlat, nbniv, ifiltre, iaire, &
           griscal, iter)
+  USE lmdz_dimensions, ONLY: iim, jjm, llm, ndm
     USE lmdz_coefils, ONLY: jfiltnu, jfiltnv, jfiltsu, jfiltsv, sddu, sddv, unsddu, unsddv, modfrstv, modfrstu
 
     !=======================================================================
-    !
+
     !   Auteur: P. Le Van        07/10/97
     !   ------
-    !
+
     !   Objet: filtre matriciel longitudinal ,avec les matrices precalculees
     !                 pour l'operateur  Filtre    .
     !   ------
-    !
+
     !   Arguments:
     !   ----------
-    !
+
     !  nblat                 nombre de latitudes a filtrer
     !  nbniv                 nombre de niveaux verticaux a filtrer
@@ -35 +37 @@
     !                        +2  Filtre directe
     !                        -2  Filtre inverse
-    !
+
     !  iaire                 1   si champ intensif
     !                        2   si champ extensif (pondere par les aires)
-    !
+
     !  iter                  1   filtre simple
-    !
+
     !=======================================================================
-    !
-    !
+
+
     !                  Variable Intensive
     !            ifiltre = 1     filtre directe
     !            ifiltre =-1     filtre inverse
-    !
+
     !                  Variable Extensive
     !            ifiltre = 2     filtre directe
     !            ifiltre =-2     filtre inverse
+
     !
-    !
-    INCLUDE "dimensions.h"
-    INCLUDE "paramet.h"
+
+
 
     INTEGER :: nlat, nbniv, ifiltre, iter
@@ -326 +328 @@
 
   SUBROUTINE inifgn(dv)
-    !
+
     !    ...  H.Upadyaya , O.Sharma  ...
-    !
+
     USE lmdz_coefils, ONLY: sddv, sddu, unsddu, unsddv, eignfnv, eignfnu
     USE lmdz_ssum_scopy, ONLY: ssum
     USE lmdz_comgeom
 
+USE lmdz_dimensions, ONLY: iim, jjm, llm, ndm
+  USE lmdz_paramet
     IMPLICIT NONE
 
-    INCLUDE "dimensions.h"
-    INCLUDE "paramet.h"
+
+
 
     REAL :: vec(iim, iim), vec1(iim, iim)
@@ -347 +351 @@
     imm1 = iim - 1
     pi = 2. * ASIN(1.)
-    !
+
     DO i = 1, iim
       dlonu(i) = xprimu(i)
@@ -359 +363 @@
       unsddv(i) = 1. / sddv(i)
     END DO
-    !
+
     DO j = 1, iim
       DO i = 1, iim
@@ -368 +372 @@
       END DO
     END DO
-    !
-    !
+
+
     eignfnv(1, 1) = -1.
     eignfnv(iim, 1) = 1.
@@ -386 +390 @@
       END DO
     END DO
-    !
+
     DO j = 1, iim
       DO i = 1, iim
@@ -398 +402 @@
     ENDDO
 
-    !
+
     CALL jacobi(vec, iim, iim, dv, eignfnv, nrot)
     CALL acc(eignfnv, d, iim)
     CALL eigen_sort(dv, eignfnv, iim, iim)
-    !
+
     CALL jacobi(vec1, iim, iim, du, eignfnu, nrot)
     CALL acc(eignfnu, d, iim)
@@ -408 +412 @@
 
     !c   ancienne version avec appels IMSL
-    !
+
     ! CALL MXM(eignfnu,iim,eignfnv,iim,vec,iim)
     ! CALL MXM(eignfnv,iim,eignfnu,iim,vec1,iim)
@@ -414 +418 @@
     ! CALL acc(eignfnv,d,iim)
     ! CALL eigen(eignfnv,dv)
-    !
+
     ! CALL EVCSF(iim,vec1,iim,du,eignfnu,iim)
     ! CALL acc(eignfnu,d,iim)
@@ -587 +591 @@
     USE lmdz_coefils, ONLY: modfrstv, modfrstu, jfiltnu, jfiltnv, coefilu, coefilv, &
             coefilu2, coefilv2, eignfnv, eignfnu, jfiltsu, jfiltsv
+  USE lmdz_dimensions, ONLY: iim, jjm, llm, ndm
     USE lmdz_comgeom
+    USE lmdz_paramet
 
     !    ... H. Upfiltreg_modadhyaya, O.Sharma   ...
@@ -595 +601 @@
     !     Correction  le 28/10/97    P. Le Van .
     !  -------------------------------------------------------------------
-    INCLUDE "dimensions.h"
-    INCLUDE "paramet.h"
 
     REAL  dlonu(iim), dlatu(jjm)