Changeset 5159 for LMDZ6/branches/Amaury_dev/libf/dyn3d_common/diverg_gam.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/dyn3d_common/diverg_gam.f90 (modified) (5 diffs)

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

@@ -4 +4 @@
 SUBROUTINE diverg_gam(klevel,cuvscvgam,cvuscugam,unsairegam , &
         unsapolnga,unsapolsga,  x, y,  div )
-  !
+
   ! P. Le Van
-  !
+
   !  *********************************************************************
   !  ... calcule la divergence a tous les niveaux d'1 vecteur de compos.
@@ -15 +15 @@
   USE lmdz_comgeom
 
+  USE lmdz_dimensions, ONLY: iim, jjm, llm, ndm
+  USE lmdz_paramet
   IMPLICIT NONE
-  !
+
   !  x  et  y  sont des arguments  d'entree pour le s-prog
   !    div      est  un argument  de sortie pour le s-prog
   !
-  INCLUDE "dimensions.h"
-  INCLUDE "paramet.h"
-  !
+
+
+
   !    ..........          variables en arguments    ...................
-  !
+
   INTEGER :: klevel
   REAL :: x( ip1jmp1,klevel ),y( ip1jm,klevel ),div( ip1jmp1,klevel )
   REAL :: cuvscvgam(ip1jm),cvuscugam(ip1jmp1),unsairegam(ip1jmp1)
   REAL :: unsapolnga,unsapolsga
-  !
+
   !    ...............     variables  locales   .........................
 
@@ -37 +39 @@
   !    ...................................................................
 
-  !
+
   DO l = 1,klevel
-  !
+
     DO  ij = iip2, ip1jm - 1
      div( ij + 1, l )     = ( &
@@ -46 +48 @@
            unsairegam( ij+1 )
     ENDDO
-  !
+
   ! ....  correction pour  div( 1,j,l)  ......
   ! ....   div(1,j,l)= div(iip1,j,l) ....
-  !
+
   !DIR$ IVDEP
     DO  ij = iip2,ip1jm,iip1
      div( ij,l ) = div( ij + iim,l )
     ENDDO
-  !
+
   ! ....  calcul  aux poles  .....
-  !
+
     DO  ij  = 1,iim
      aiy1(ij) =    cuvscvgam(    ij       ) * y(     ij     , l )
@@ -63 +65 @@
     sumypn = SSUM ( iim,aiy1,1 ) * unsapolnga
     sumyps = SSUM ( iim,aiy2,1 ) * unsapolsga
-  !
+
     DO  ij = 1,iip1
      div(     ij    , l ) = - sumypn