Changeset 5105 for LMDZ6/branches/Amaury_dev/libf/dyn3dmem/dudv2_loc.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/dyn3dmem/dudv2_loc.f90 (moved) (moved from LMDZ6/branches/Amaury_dev/libf/dyn3dmem/dudv2_loc.F) (1 diff)

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

@@ -1 @@
-      SUBROUTINE dudv2_loc ( teta, pkf, bern, du, dv  )
-      USE parallel_lmdz
-      IMPLICIT NONE
-c
-c=======================================================================
-c
-c   Auteur:  P. Le Van
-c   -------
-c
-c   Objet:
-c   ------
-c
-c   *****************************************************************
-c   ..... calcul du terme de pression (gradient de p/densite )   et
-c          du terme de ( -gradient de la fonction de Bernouilli ) ...
-c   *****************************************************************
-c          Ces termes sont ajoutes a  d(ucov)/dt et a d(vcov)/dt  ..
-c
-c
-c    teta , pkf, bern  sont des arguments d'entree  pour le s-pg  ....
-c    du et dv          sont des arguments de sortie pour le s-pg  ....
-c
-c=======================================================================
-c
-      include "dimensions.h"
-      include "paramet.h"
+SUBROUTINE dudv2_loc ( teta, pkf, bern, du, dv  )
+  USE parallel_lmdz
+  IMPLICIT NONE
+  !
+  !=======================================================================
+  !
+  !   Auteur:  P. Le Van
+  !   -------
+  !
+  !   Objet:
+  !   ------
+  !
+  !   *****************************************************************
+  !   ..... calcul du terme de pression (gradient de p/densite )   et
+  !      du terme de ( -gradient de la fonction de Bernouilli ) ...
+  !   *****************************************************************
+  !      Ces termes sont ajoutes a  d(ucov)/dt et a d(vcov)/dt  ..
+  !
+  !
+  !    teta , pkf, bern  sont des arguments d'entree  pour le s-pg  ....
+  !    du et dv          sont des arguments de sortie pour le s-pg  ....
+  !
+  !=======================================================================
+  !
+  include "dimensions.h"
+  include "paramet.h"
 
-      REAL teta( ijb_u:ije_u,llm ),pkf( ijb_u:ije_u,llm )
-      REAL bern( ijb_u:ije_u,llm )
-      REAL du( ijb_u:ije_u,llm ),  dv( ijb_v:ije_v,llm )
-      INTEGER  l,ij,ijb,ije
-c
-c
-c$OMP DO SCHEDULE(STATIC,OMP_CHUNK)
-      DO l = 1,llm
-c
-      ijb=ij_begin
-      ije=ij_end
-      if (pole_nord) ijb=ijb+iip1
-      if (pole_sud)  ije=ije-iip1
+  REAL :: teta( ijb_u:ije_u,llm ),pkf( ijb_u:ije_u,llm )
+  REAL :: bern( ijb_u:ije_u,llm )
+  REAL :: du( ijb_u:ije_u,llm ),  dv( ijb_v:ije_v,llm )
+  INTEGER :: l,ij,ijb,ije
+  !
+  !
+!$OMP DO SCHEDULE(STATIC,OMP_CHUNK)
+  DO l = 1,llm
+  !
+  ijb=ij_begin
+  ije=ij_end
+  if (pole_nord) ijb=ijb+iip1
+  if (pole_sud)  ije=ije-iip1
 
-      DO ij  = ijb, ije - 1
-       du(ij,l) = du(ij,l) + 0.5* ( teta( ij,l ) + teta( ij+1,l ) ) *
-     * ( pkf( ij,l ) - pkf(ij+1,l) )  + bern(ij,l) - bern(ij+1,l)
-      END DO
-c
-c
-c    .....  correction  pour du(iip1,j,l),  j=2,jjm   ......
-c    ...          du(iip1,j,l) = du(1,j,l)                 ...
-c
-CDIR$ IVDEP
-      DO ij = ijb+iip1-1, ije, iip1
-      du( ij,l ) = du( ij - iim,l )
-      END DO
-c
-c
-      if (pole_nord) ijb=ijb-iip1
+  DO ij  = ijb, ije - 1
+   du(ij,l) = du(ij,l) + 0.5* ( teta( ij,l ) + teta( ij+1,l ) ) * &
+         ( pkf( ij,l ) - pkf(ij+1,l) )  + bern(ij,l) - bern(ij+1,l)
+  END DO
+  !
+  !
+  !    .....  correction  pour du(iip1,j,l),  j=2,jjm   ......
+  !    ...          du(iip1,j,l) = du(1,j,l)                 ...
+  !
+  !DIR$ IVDEP
+  DO ij = ijb+iip1-1, ije, iip1
+  du( ij,l ) = du( ij - iim,l )
+  END DO
+  !
+  !
+  if (pole_nord) ijb=ijb-iip1
 
-      DO ij  = ijb,ije
-      dv( ij,l) = dv(ij,l) + 0.5 * ( teta(ij,l) + teta( ij+iip1,l ) ) *
-     *                             ( pkf(ij+iip1,l) - pkf(  ij,l  ) )
-     *                           +   bern( ij+iip1,l ) - bern( ij  ,l )
-      END DO
-c
-      END DO
-c$OMP END DO NOWAIT 
-c
-      RETURN
-      END
+  DO ij  = ijb,ije
+  dv( ij,l) = dv(ij,l) + 0.5 * ( teta(ij,l) + teta( ij+iip1,l ) ) * &
+        ( pkf(ij+iip1,l) - pkf(  ij,l  ) ) &
+        +   bern( ij+iip1,l ) - bern( ij  ,l )
+  END DO
+  !
+  END DO
+!$OMP END DO NOWAIT
+  !
+
+END SUBROUTINE dudv2_loc