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

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

@@ -1 @@
-      SUBROUTINE dudv1_loc ( vorpot, pbaru, pbarv, du, dv )
-      USE parallel_lmdz
-      IMPLICIT NONE
-c
-c-----------------------------------------------------------------------
-c
-c   Auteur:   P. Le Van
-c   -------
-c
-c   Objet:
-c   ------
-c   calcul du terme de  rotation
-c   ce terme est ajoute a  d(ucov)/dt et a d(vcov)/dt  ..
-c   vorpot, pbaru et pbarv sont des arguments d'entree  pour le s-pg ..
-c   du  et dv              sont des arguments de sortie pour le s-pg ..
-c
-c-----------------------------------------------------------------------
+SUBROUTINE dudv1_loc ( vorpot, pbaru, pbarv, du, dv )
+  USE parallel_lmdz
+  IMPLICIT NONE
+  !
+  !-----------------------------------------------------------------------
+  !
+  !   Auteur:   P. Le Van
+  !   -------
+  !
+  !   Objet:
+  !   ------
+  !   calcul du terme de  rotation
+  !   ce terme est ajoute a  d(ucov)/dt et a d(vcov)/dt  ..
+  !   vorpot, pbaru et pbarv 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"
+  INCLUDE "dimensions.h"
+  INCLUDE "paramet.h"
 
-      REAL vorpot( ijb_v:ije_v,llm ) ,pbaru( ijb_u:ije_u,llm ) ,
-     *     pbarv( ijb_v:ije_v,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=ij_begin+iip1
-      if (pole_sud)  ije=ij_end-iip1
-      
-      DO ij = ijb, ije-1
-      du( ij,l ) = 0.125 *(  vorpot(ij-iip1, l) + vorpot( ij, l)  ) *
-     *                    (   pbarv(ij-iip1, l) + pbarv(ij-iim,  l) +
-     *                        pbarv(   ij  , l) + pbarv(ij+ 1 ,  l)   )
-      END DO
-   
- 
-c
-      if (pole_nord) ijb=ij_begin
-      
-      DO ij = ijb, ije-1
-      dv( ij+1,l ) = - 0.125 *(  vorpot(ij, l)  + vorpot(ij+1, l)  ) *
-     *                        (   pbaru(ij, l)  +  pbaru(ij+1   , l) +
-     *                       pbaru(ij+iip1, l)  +  pbaru(ij+iip2, l)  )
-      END DO
-c
-c    .... correction  pour  dv( 1,j,l )  .....
-c    ....   dv(1,j,l)= dv(iip1,j,l) ....
-c
-CDIR$ IVDEP
-      DO ij = ijb, ije, iip1
-      dv( ij,l ) = dv( ij + iim, l )
-      END DO
-c
-      END DO
-c$OMP END DO NOWAIT
-      RETURN
-      END
+  REAL :: vorpot( ijb_v:ije_v,llm ) ,pbaru( ijb_u:ije_u,llm ) , &
+        pbarv( ijb_v:ije_v,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=ij_begin+iip1
+  if (pole_sud)  ije=ij_end-iip1
+
+  DO ij = ijb, ije-1
+  du( ij,l ) = 0.125 *(  vorpot(ij-iip1, l) + vorpot( ij, l)  ) * &
+        (   pbarv(ij-iip1, l) + pbarv(ij-iim,  l) + &
+        pbarv(   ij  , l) + pbarv(ij+ 1 ,  l)   )
+  END DO
+
+
+  !
+  if (pole_nord) ijb=ij_begin
+
+  DO ij = ijb, ije-1
+  dv( ij+1,l ) = - 0.125 *(  vorpot(ij, l)  + vorpot(ij+1, l)  ) * &
+        (   pbaru(ij, l)  +  pbaru(ij+1   , l) + &
+        pbaru(ij+iip1, l)  +  pbaru(ij+iip2, l)  )
+  END DO
+  !
+  !    .... correction  pour  dv( 1,j,l )  .....
+  !    ....   dv(1,j,l)= dv(iip1,j,l) ....
+  !
+  !DIR$ IVDEP
+  DO ij = ijb, ije, iip1
+  dv( ij,l ) = dv( ij + iim, l )
+  END DO
+  !
+  END DO
+!$OMP END DO NOWAIT
+
+END SUBROUTINE dudv1_loc