Changeset 5246 for LMDZ6/trunk/libf/dyn3dmem/bernoui_loc.f90

Timestamp:

Oct 21, 2024, 2:58:45 PM (23 hours ago)

Author:

abarral

Message:

Convert fixed-form to free-form sources .F -> .{f,F}90
(WIP: some .F remain, will be handled in subsequent commits)

File:

• LMDZ6/trunk/libf/dyn3dmem/bernoui_loc.f90 (moved) (moved from LMDZ6/trunk/libf/dyn3dmem/bernoui_loc.F) (1 diff)

LMDZ6/trunk/libf/dyn3dmem/bernoui_loc.f90

@@ -1 @@
-      SUBROUTINE bernoui_loc (ngrid,nlay,pphi,pecin,pbern)
-      USE parallel_lmdz
-      USE mod_filtreg_p
-      IMPLICIT NONE
+SUBROUTINE bernoui_loc (ngrid,nlay,pphi,pecin,pbern)
+  USE parallel_lmdz
+  USE mod_filtreg_p
+  IMPLICIT NONE
 
-c=======================================================================
-c
-c   Auteur:   P. Le Van
-c   -------
-c
-c   Objet:
-c   ------
-c     calcul de la fonction de Bernouilli aux niveaux s  .....
-c     phi  et  ecin  sont des arguments d'entree pour le s-pg .......
-c          bern       est un  argument de sortie pour le s-pg  ......
-c
-c    fonction de Bernouilli = bern = filtre de( geopotentiel + 
-c                              energ.cinet.)
-c
-c=======================================================================
-c
-c-----------------------------------------------------------------------
-c   Decalrations:
-c   -------------
-c
-      include "dimensions.h"
-      include "paramet.h"
-c
-c   Arguments:
-c   ----------
-c
-      INTEGER nlay,ngrid
-      REAL pphi(ijb_u:ije_u,nlay),pecin(ijb_u:ije_u,nlay)
-      REAL pbern(ijb_u:ije_u,nlay)
-c
-c   Local:
-c   ------
-c
-      INTEGER   ij,l,ijb,ije,jjb,jje
-c
-c-----------------------------------------------------------------------
-c   calcul de Bernouilli:
-c   ---------------------
-c
-      ijb=ij_begin
-      ije=ij_end+iip1
-      if (pole_sud) ije=ij_end
+  !=======================================================================
+  !
+  !   Auteur:   P. Le Van
+  !   -------
+  !
+  !   Objet:
+  !   ------
+  ! calcul de la fonction de Bernouilli aux niveaux s  .....
+  ! phi  et  ecin  sont des arguments d'entree pour le s-pg .......
+  !      bern       est un  argument de sortie pour le s-pg  ......
+  !
+  !    fonction de Bernouilli = bern = filtre de( geopotentiel +
+  !                          energ.cinet.)
+  !
+  !=======================================================================
+  !
+  !-----------------------------------------------------------------------
+  !   Decalrations:
+  !   -------------
+  !
+  include "dimensions.h"
+  include "paramet.h"
+  !
+  !   Arguments:
+  !   ----------
+  !
+  INTEGER :: nlay,ngrid
+  REAL :: pphi(ijb_u:ije_u,nlay),pecin(ijb_u:ije_u,nlay)
+  REAL :: pbern(ijb_u:ije_u,nlay)
+  !
+  !   Local:
+  !   ------
+  !
+  INTEGER :: ij,l,ijb,ije,jjb,jje
+  !
+  !-----------------------------------------------------------------------
+  !   calcul de Bernouilli:
+  !   ---------------------
+  !
+  ijb=ij_begin
+  ije=ij_end+iip1
+  if (pole_sud) ije=ij_end
 
-      jjb=jj_begin
-      jje=jj_end+1
-      if (pole_sud) jje=jj_end
+  jjb=jj_begin
+  jje=jj_end+1
+  if (pole_sud) jje=jj_end
 
-c$OMP DO SCHEDULE(STATIC,OMP_CHUNK)                
-      DO l=1,llm
-    
-        DO 4 ij = ijb,ije
-          pbern( ij,l ) =  pphi( ij,l ) + pecin( ij,l )
-   4    CONTINUE
-       
-       ENDDO
-c$OMP END DO NOWAIT
-c
-c-----------------------------------------------------------------------
-c   filtre:
-c   -------
-c
+!$OMP DO SCHEDULE(STATIC,OMP_CHUNK)
+  DO l=1,llm
 
-        
-        CALL filtreg_p( pbern,jjb_u,jje_u,jjb,jje, jjp1, llm, 
-     &                  2,1, .true., 1 )
-c
-c-----------------------------------------------------------------------
-      
-      
-      RETURN
-      END
+    DO ij = ijb,ije
+      pbern( ij,l ) =  pphi( ij,l ) + pecin( ij,l )
+    END DO
+
+   ENDDO
+!$OMP END DO NOWAIT
+  !
+  !-----------------------------------------------------------------------
+  !   filtre:
+  !   -------
+  !
+
+
+    CALL filtreg_p( pbern,jjb_u,jje_u,jjb,jje, jjp1, llm, &
+          2,1, .true., 1 )
+  !
+  !-----------------------------------------------------------------------
+
+
+  RETURN
+END SUBROUTINE bernoui_loc