Changeset 5136 for LMDZ6/branches/Amaury_dev/libf/dyn3d_common/gr_v_scal.f90

Timestamp:

Jul 28, 2024, 4:17:54 PM (3 months ago)

Author:

abarral

Message:

Put comgeom.h, comgeom2.h into modules

File:

• LMDZ6/branches/Amaury_dev/libf/dyn3d_common/gr_v_scal.f90 (modified) (4 diffs)

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

@@ -1 @@
-
 ! $Header$
 
-SUBROUTINE gr_v_scal(nx,x_v,x_scal)
+SUBROUTINE gr_v_scal(nx, x_v, x_scal)
   !%W%    %G%
   !=======================================================================
@@ -25 +24 @@
   !
   !=======================================================================
+  USE lmdz_comgeom
   IMPLICIT NONE
   !-----------------------------------------------------------------------
@@ -32 +32 @@
   INCLUDE "dimensions.h"
   INCLUDE "paramet.h"
-  INCLUDE "comgeom.h"
 
   !   Arguments:
@@ -38 +37 @@
 
   INTEGER :: nx
-  REAL :: x_v(ip1jm,nx),x_scal(ip1jmp1,nx)
+  REAL :: x_v(ip1jm, nx), x_scal(ip1jmp1, nx)
 
   !   Local:
   !   ------
 
-  INTEGER :: l,ij
+  INTEGER :: l, ij
 
   !-----------------------------------------------------------------------
 
-  DO l=1,nx
-     DO ij=iip2,ip1jm
-        x_scal(ij,l)= &
-              (airev(ij-iip1)*x_v(ij-iip1,l)+airev(ij)*x_v(ij,l)) &
-              /(airev(ij-iip1)+airev(ij))
-     ENDDO
-     DO ij=1,iip1
-        x_scal(ij,l)=0.
-     ENDDO
-     DO ij=ip1jm+1,ip1jmp1
-        x_scal(ij,l)=0.
-     ENDDO
+  DO l = 1, nx
+    DO ij = iip2, ip1jm
+      x_scal(ij, l) = &
+              (airev(ij - iip1) * x_v(ij - iip1, l) + airev(ij) * x_v(ij, l)) &
+                      / (airev(ij - iip1) + airev(ij))
+    ENDDO
+    DO ij = 1, iip1
+      x_scal(ij, l) = 0.
+    ENDDO
+    DO ij = ip1jm + 1, ip1jmp1
+      x_scal(ij, l) = 0.
+    ENDDO
   ENDDO