alew1.F File Reference

#include "implicit_f.inc"
#include "com01_c.inc"
#include "com04_c.inc"
#include "com08_c.inc"
#include "tabsiz_c.inc"

Go to the source code of this file.

Functions/Subroutines
subroutine	alew1 (d, v, w, ale_nn_connect, nale, nodft, nodlt, nbrcvois, nbsdvois, lnrcvois, lnsdvois)

Function/Subroutine Documentation

alew1()

subroutine alew1	(		d,
			v,
			w,
		type(t_connectivity), intent(in)	ale_nn_connect,
		integer, dimension(numnod)	nale,
		integer	nodft,
		integer	nodlt,
		integer, dimension(*)	nbrcvois,
		integer, dimension(*)	nbsdvois,
		integer, dimension(*)	lnrcvois,
		integer, dimension(*)	lnsdvois )

Definition at line 33 of file alew1.F.

C-----------------------------------------------
C   M o d u l e s
C-----------------------------------------------      
      USE ale_connectivity_mod
      USE ale_mod
C-----------------------------------------------
C   D e s c r i p t i o n
C-----------------------------------------------
C Compute Grid for /ALE/GRID/DISP
C
C     X,D,V are allocated to SX,SD,DV=3*(NUMNOD_L+NUMVVOIS_L)
C      in grid subroutine it may needed to access nodes which
C      are connected to a remote elem. They are sored in X(1:3,NUMNOD+1:)
C      Consequently X is defined here X(3,SX/3) instead of X(3,NUMNOD) as usually
C-----------------------------------------------
C   I m p l i c i t   T y p e s
C-----------------------------------------------
#include      "implicit_f.inc"
C-----------------------------------------------
C   C o m m o n   B l o c k s
C-----------------------------------------------
#include      "com01_c.inc"
#include      "com04_c.inc"
#include      "com08_c.inc"
#include      "tabsiz_c.inc"
C-----------------------------------------------
C   D u m m y   A r g u m e n t s
C-----------------------------------------------
! SPMD CASE : SX >= 3*NUMNOD    (SX = 3*(NUMNOD_L+NRCVVOIS_L))
! X(1:3,1:NUMNOD) : local nodes
!  (1:3, NUMNOD+1:) additional nodes (also on adjacent domains but connected to the boundary of the current domain)
!      idem with D(SD), and V(SV)
C-----------------------------------------------
      INTEGER NALE(NUMNOD), NODFT, NODLT,
     .        NBRCVOIS(*),NBSDVOIS(*),
     .        LNRCVOIS(*),LNSDVOIS(*)
      my_real d(3,sd/3), v(3,sv/3), w(3,sw/3)
      TYPE(t_connectivity), INTENT(IN) :: ALE_NN_CONNECT
C-----------------------------------------------
C   L o c a l   V a r i a b l e s
C-----------------------------------------------
      INTEGER I, LENCOM
      my_real dx, dy, dz,fac
      INTEGER :: IAD, IAD1, IAD2, NODE_ID
      my_real :: dx_max, dx_min
      my_real :: dy_max, dy_min
      my_real :: dz_max, dz_min
C-----------------------------------------------
C   S o u r c e   L i n e s
C-----------------------------------------------
 
C  SPMD : COMM D
      IF(nspmd > 1)THEN
!$OMP SINGLE
         lencom = nbrcvois(nspmd+1)+nbsdvois(nspmd+1)
         CALL spmd_xvois(d,nbrcvois,nbsdvois,lnrcvois,lnsdvois,lencom)
!$OMP END SINGLE
      END IF
      fac=one/dt2/three
      DO i = nodft, nodlt
         IF(iabs(nale(i)) == 1) THEN
            iad1 = ale_nn_connect%IAD_CONNECT(i)
            iad2 = ale_nn_connect%IAD_CONNECT(i + 1) - 1
            dx_max = -huge(zero)
            dx_min = huge(zero)
            dy_max = -huge(zero)
            dy_min = huge(zero)
            dz_max = -huge(zero)
            dz_min = huge(zero)
            DO iad = iad1, iad2
               node_id = ale_nn_connect%CONNECTED(iad)
               dx_max = max(dx_max, d(1, node_id))
               dx_min = min(dx_min, d(1, node_id))
               dy_max = max(dy_max, d(2, node_id))
               dy_min = min(dy_min, d(2, node_id))
               dz_max = max(dz_max, d(3, node_id))
               dz_min = min(dz_min, d(3, node_id))
            ENDDO
            dx = fac * (dx_max + dx_min - two * d(1, i))
            dy = fac * (dy_max + dy_min - two * d(2, i))
            dz = fac * (dz_max + dz_min - two * d(3, i))
            w(1,i)=sign(one,dx)*min(abs(dx),ale%GRID%ALPHA)
            w(2,i)=sign(one,dy)*min(abs(dy),ale%GRID%ALPHA)
            w(3,i)=sign(one,dz)*min(abs(dz),ale%GRID%ALPHA)
         ELSEIF(nale(i) == 0)THEN
            w(1,i)=v(1,i)
            w(2,i)=v(2,i)
            w(3,i)=v(3,i)
         ELSE
            w(1,i)=zero
            w(2,i)=zero
            w(3,i)=zero     
         ENDIF
C     
      ENDDO
C     
      RETURN