Changeset 5105 for LMDZ6/branches/Amaury_dev/libf/dyn3d_common/divgrad.f90

Timestamp:

Jul 23, 2024, 7:14:34 PM (8 weeks ago)

Author:

abarral

Message:

Replace 1DUTILS.h by module lmdz_1dutils.f90
Replace 1DConv.h by module lmdz_old_1dconv.f90 (it's only used by old_* files)
Convert *.F to *.f90
Fix gradsdef.h formatting
Remove unnecessary "RETURN" at the end of functions/subroutines

File:

• LMDZ6/branches/Amaury_dev/libf/dyn3d_common/divgrad.f90 (moved) (moved from LMDZ6/branches/Amaury_dev/libf/dyn3d_common/divgrad.F) (1 diff)

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

@@ -2 +2 @@
 ! $Header$
 
-      SUBROUTINE divgrad (klevel,h, lh, divgra )
-      IMPLICIT NONE
-c
-c=======================================================================
-c
-c  Auteur :   P. Le Van
-c  ----------
-c
-c                              lh
-c      calcul de  (div( grad ))   de h  .....
-c      h  et lh  sont des arguments  d'entree pour le s-prog
-c      divgra     est  un argument  de sortie pour le s-prog
-c
-c=======================================================================
-c
-c   declarations:
-c   -------------
-c
-      include "dimensions.h"
-      include "paramet.h"
-      include "comgeom.h"
-      include "comdissipn.h"
-c
-      INTEGER klevel
-      REAL h( ip1jmp1,klevel ), divgra( ip1jmp1,klevel )
-c
-      REAL ghy(ip1jm,llm), ghx(ip1jmp1,llm)
+SUBROUTINE divgrad (klevel,h, lh, divgra )
+  IMPLICIT NONE
+  !
+  !=======================================================================
+  !
+  !  Auteur :   P. Le Van
+  !  ----------
+  !
+  !                          lh
+  !  calcul de  (div( grad ))   de h  .....
+  !  h  et lh  sont des arguments  d'entree pour le s-prog
+  !  divgra     est  un argument  de sortie pour le s-prog
+  !
+  !=======================================================================
+  !
+  !   declarations:
+  !   -------------
+  !
+  include "dimensions.h"
+  include "paramet.h"
+  include "comgeom.h"
+  include "comdissipn.h"
+  !
+  INTEGER :: klevel
+  REAL :: h( ip1jmp1,klevel ), divgra( ip1jmp1,klevel )
+  !
+  REAL :: ghy(ip1jm,llm), ghx(ip1jmp1,llm)
 
-      INTEGER  l,ij,iter,lh
-c
-c
-c
-      CALL SCOPY ( ip1jmp1*klevel,h,1,divgra,1 )
-c
-      DO iter = 1,lh
+  INTEGER :: l,ij,iter,lh
+  !
+  !
+  !
+  CALL SCOPY ( ip1jmp1*klevel,h,1,divgra,1 )
+  !
+  DO iter = 1,lh
 
-      CALL filtreg ( divgra,jjp1,klevel,2,1,.TRUE.,1  )
+  CALL filtreg ( divgra,jjp1,klevel,2,1,.TRUE.,1  )
 
-      CALL    grad (klevel,divgra, ghx  , ghy          )
-      CALL  diverg (klevel,  ghx , ghy  , divgra       )
+  CALL    grad (klevel,divgra, ghx  , ghy          )
+  CALL  diverg (klevel,  ghx , ghy  , divgra       )
 
-      CALL filtreg ( divgra,jjp1,klevel,2,1,.TRUE.,1)
+  CALL filtreg ( divgra,jjp1,klevel,2,1,.TRUE.,1)
 
-      DO l = 1,klevel
-      DO ij = 1, ip1jmp1
-      divgra( ij,l ) = - cdivh * divgra( ij,l )
-      END DO
-      END DO
-c
-      END DO
-      RETURN
-      END
+  DO l = 1,klevel
+  DO ij = 1, ip1jmp1
+  divgra( ij,l ) = - cdivh * divgra( ij,l )
+  END DO
+  END DO
+  !
+  END DO
+  RETURN
+END SUBROUTINE divgrad