Changeset 5159 for LMDZ6/branches/Amaury_dev/libf/dyn3d_common/divgrad2.f90

Timestamp:

Aug 2, 2024, 9:58:25 PM (3 months ago)

Author:

abarral

Message:

Put dimensions.h and paramet.h into modules

File:

• LMDZ6/branches/Amaury_dev/libf/dyn3d_common/divgrad2.f90 (modified) (5 diffs)

LMDZ6/branches/Amaury_dev/libf/dyn3d_common/divgrad2.f90

@@ -2 +2 @@
 
 SUBROUTINE divgrad2(klevel, h, deltapres, lh, divgra)
-  !
+
   ! P. Le Van
-  !
+
   !   ***************************************************************
-  !
+
   ! .....   calcul de  (div( grad ))   de (  pext * h ) .....
   !   ****************************************************************
   !   h ,klevel,lh et pext  sont des arguments  d'entree pour le s-prg
   !     divgra     est  un argument  de sortie pour le s-prg
-  !
+
   USE lmdz_ssum_scopy, ONLY: scopy
   USE lmdz_comdissipn, ONLY: tetaudiv, tetaurot, tetah, cdivu, crot, cdivh
   USE lmdz_comgeom2
 
+USE lmdz_dimensions, ONLY: iim, jjm, llm, ndm
+  USE lmdz_paramet
   IMPLICIT NONE
   !
-  INCLUDE "dimensions.h"
-  INCLUDE "paramet.h"
+
+
 
   !    .......    variables en arguments   .......
-  !
+
   INTEGER :: klevel
   REAL :: h(ip1jmp1, klevel), deltapres(ip1jmp1, klevel)
   REAL :: divgra(ip1jmp1, klevel)
-  !
+
   !    .......    variables  locales    ..........
-  !
+
   REAL :: signe, nudivgrs, sqrtps(ip1jmp1, llm)
   INTEGER :: l, ij, iter, lh
   !    ...................................................................
 
-  !
+
   signe = (-1.)**lh
   nudivgrs = signe * cdivh
@@ -39 +41 @@
   CALL SCOPY (ip1jmp1 * klevel, h, 1, divgra, 1)
 
-  !
+
   CALL laplacien(klevel, divgra, divgra)
 
@@ -47 +49 @@
     ENDDO
   ENDDO
-  !
+
   DO l = 1, klevel
     DO ij = 1, ip1jmp1
@@ -55 +57 @@
 
   !    ........    Iteration de l'operateur  laplacien_gam    ........
-  !
+
   DO  iter = 1, lh - 2
     CALL laplacien_gam (klevel, cuvscvgam2, cvuscugam2, unsair_gam2, &
             unsapolnga2, unsapolsga2, divgra, divgra)
   ENDDO
-  !
+
   !    ...............................................................
 
@@ -68 +70 @@
     ENDDO
   ENDDO
-  !
+
   CALL laplacien (klevel, divgra, divgra)
-  !
+
   DO l = 1, klevel
     DO ij = 1, ip1jmp1