Changeset 5159 for LMDZ6/branches/Amaury_dev/libf/dyn3dmem/diverg_p.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/dyn3dmem/diverg_p.f90 (modified) (8 diffs)

LMDZ6/branches/Amaury_dev/libf/dyn3dmem/diverg_p.f90

@@ -1 +1 @@
 SUBROUTINE diverg_p(klevel,x,y,div)
-  !
+
   ! P. Le Van
-  !
+
   !  *********************************************************************
   !  ... calcule la divergence a tous les niveaux d'1 vecteur de compos.
@@ -12 +12 @@
   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 )
   INTEGER :: l,ij
-  !
+
   !    ...............     variables  locales   .........................
 
@@ -33 +35 @@
   INTEGER :: ijb,ije
   !    ...................................................................
-  !
-  !
+
+
   ijb=ij_begin
   ije=ij_end
@@ -42 +44 @@
 !$OMP DO SCHEDULE(STATIC,OMP_CHUNK)
   DO l = 1,klevel
-  !
+
     DO  ij = ijb, ije - 1
      div( ij + 1, l )     = &
@@ -48 +50 @@
            cuvsurcv(ij-iim) * y(ij-iim,l) - cuvsurcv(ij+1) * y(ij+1,l)
     ENDDO
-  !
+
   ! ....  correction pour  div( 1,j,l)  ......
   ! ....   div(1,j,l)= div(iip1,j,l) ....
-  !
+
   !DIR$ IVDEP
     DO  ij = ijb,ije,iip1
      div( ij,l ) = div( ij + iim,l )
     ENDDO
-  !
+
   ! ....  calcul  aux poles  .....
-  !
+
     IF (pole_nord) THEN
       DO  ij  = 1,iim
@@ -64 +66 @@
       ENDDO
       sumypn = SSUM ( iim,aiy1,1 ) / apoln
-  !
+
       DO  ij = 1,iip1
        div(     ij    , l ) = - sumypn
@@ -75 +77 @@
       ENDDO
       sumyps = SSUM ( iim,aiy2,1 ) / apols
-  !
+
       DO  ij = 1,iip1
        div( ij + ip1jm, l ) =   sumyps
@@ -88 +90 @@
   !cc        CALL filtreg( div, jjp1, klevel, 2, 2, .TRUE., 1 )
 
-  !
+
 !$OMP DO SCHEDULE(STATIC,OMP_CHUNK)
     DO l = 1, klevel