Changeset 5105 for LMDZ6/branches/Amaury_dev/libf/dyn3dmem/laplacien_loc.f90

Timestamp:

Jul 23, 2024, 7:14:34 PM (3 months 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/dyn3dmem/laplacien_loc.f90 (moved) (moved from LMDZ6/branches/Amaury_dev/libf/dyn3dmem/laplacien_loc.F) (1 diff)

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

@@ -1 @@
-      SUBROUTINE laplacien_loc ( klevel, teta, divgra )
-c
-c     P. Le Van
-c
-c   ************************************************************
-c    ....     calcul de  (div( grad ))   de   teta  .....
-c   ************************************************************
-c     klevel et teta  sont des arguments  d'entree pour le s-prog
-c      divgra     est  un argument  de sortie pour le s-prog
-c
-      USE parallel_lmdz
-      USE mod_filtreg_p
-      IMPLICIT NONE
-c
-      INCLUDE "dimensions.h"
-      INCLUDE "paramet.h"
-      INCLUDE "comgeom.h"
+SUBROUTINE laplacien_loc ( klevel, teta, divgra )
+  !
+  ! P. Le Van
+  !
+  !   ************************************************************
+  !    ....     calcul de  (div( grad ))   de   teta  .....
+  !   ************************************************************
+  ! klevel et teta  sont des arguments  d'entree pour le s-prog
+  !  divgra     est  un argument  de sortie pour le s-prog
+  !
+  USE parallel_lmdz
+  USE mod_filtreg_p
+  IMPLICIT NONE
+  !
+  INCLUDE "dimensions.h"
+  INCLUDE "paramet.h"
+  INCLUDE "comgeom.h"
 
-c
-c    .........      variables  en arguments   ..............
-c
-      INTEGER klevel
-      REAL teta( ijb_u:ije_u,klevel ), divgra( ijb_u:ije_u,klevel )
-      INTEGER :: l
-c
-c    ............     variables  locales      ..............
-c
-      REAL ghy(ijb_v:ije_v,llm), ghx(ijb_u:ije_u,llm)
-c    .......................................................
+  !
+  !    .........      variables  en arguments   ..............
+  !
+  INTEGER :: klevel
+  REAL :: teta( ijb_u:ije_u,klevel ), divgra( ijb_u:ije_u,klevel )
+  INTEGER :: l
+  !
+  !    ............     variables  locales      ..............
+  !
+  REAL :: ghy(ijb_v:ije_v,llm), ghx(ijb_u:ije_u,llm)
+  !    .......................................................
 
-      
-      INTEGER :: ijb,ije,jjb,jje
-c
-c      CALL SCOPY ( ip1jmp1 * klevel, teta, 1, divgra, 1 )
 
-      ijb=ij_begin-iip1
-      ije=ij_end+iip1
-      if (pole_nord) ijb=ij_begin
-      if (pole_sud ) ije=ij_end
+  INTEGER :: ijb,ije,jjb,jje
+  !
+  !  CALL SCOPY ( ip1jmp1 * klevel, teta, 1, divgra, 1 )
 
-c$OMP DO SCHEDULE(STATIC,OMP_CHUNK)
-      DO l=1,klevel      
-        divgra(ijb:ije,l)=teta(ijb:ije,l)
-      ENDDO
-c$OMP END DO NOWAIT
-      
-      jjb=jj_begin-1
-      jje=jj_end+1
-      if (pole_nord) jjb=jj_begin
-      if (pole_sud ) jje=jj_end
-      
-      CALL filtreg_p( divgra,jjb_u,jje_u,jjb,jje,jjp1, 
-     &                klevel,  2, 1, .TRUE., 1 )
-      CALL   grad_loc ( klevel,divgra,   ghx , ghy              )
-      CALL  divergf_loc ( klevel, ghx , ghy  , divgra           )
+  ijb=ij_begin-iip1
+  ije=ij_end+iip1
+  if (pole_nord) ijb=ij_begin
+  if (pole_sud ) ije=ij_end
 
-      RETURN
-      END
+!$OMP DO SCHEDULE(STATIC,OMP_CHUNK)
+  DO l=1,klevel
+    divgra(ijb:ije,l)=teta(ijb:ije,l)
+  ENDDO
+!$OMP END DO NOWAIT
+
+  jjb=jj_begin-1
+  jje=jj_end+1
+  if (pole_nord) jjb=jj_begin
+  if (pole_sud ) jje=jj_end
+
+  CALL filtreg_p( divgra,jjb_u,jje_u,jjb,jje,jjp1, &
+        klevel,  2, 1, .TRUE., 1 )
+  CALL   grad_loc ( klevel,divgra,   ghx , ghy              )
+  CALL  divergf_loc ( klevel, ghx , ghy  , divgra           )
+
+
+END SUBROUTINE laplacien_loc