ghostuv.f90 File Reference

Go to the source code of this file.

Functions/Subroutines
subroutine	ghostuv2 (I, JJ, UAKK, VAKK)
subroutine	ghostuv3 (I, JJ, UAKK, VAKK)

Function/Subroutine Documentation

ghostuv2()

subroutine ghostuv2	(	integer, intent(in)	I,
		integer, intent(in)	JJ,
		real(sp), intent(out)	UAKK,
		real(sp), intent(out)	VAKK
	)

Definition at line 44 of file ghostuv.f90.

!==============================================================================|

   USE all_vars
   USE bcs
   IMPLICIT NONE
   INTEGER, INTENT(IN)  :: I,JJ
   INTEGER              :: J1,J2
   REAL(SP)             :: DELTX,DELTY,ALPHA1
   REAL(SP)             :: UTMP,VTMP
   REAL(SP), INTENT(OUT):: UAKK,VAKK

   uakk = 0.0_sp; vakk = 0.0_sp

   IF(isbce(i) /= 2)THEN
     j1 = jj+1-int((jj+1)/4)*3
     j2 = jj+2-int((jj+2)/4)*3
     deltx = vx(nv(i,j1))-vx(nv(i,j2))
     delty = vy(nv(i,j1))-vy(nv(i,j2))

     alpha1 = atan2(delty,deltx)

     utmp = ua(i)*cos(alpha1)+va(i)*sin(alpha1)
     vtmp = -ua(i)*sin(alpha1)+va(i)*cos(alpha1)

!     VTMP = -VTMP
     vtmp = 0.0_sp

     uakk = utmp*cos(alpha1)-vtmp*sin(alpha1)
     vakk = utmp*sin(alpha1)+vtmp*cos(alpha1)
   ELSE IF(isbce(i) == 2)THEN
     uakk = ua(i)
     vakk = va(i)
   END IF

   RETURN

Here is the caller graph for this function:

ghostuv3()

subroutine ghostuv3	(	integer, intent(in)	I,
		integer, intent(in)	JJ,
		real(sp), dimension(kb), intent(out)	UAKK,
		real(sp), dimension(kb), intent(out)	VAKK
	)

Definition at line 86 of file ghostuv.f90.

!==============================================================================|

   USE all_vars
   USE bcs
   IMPLICIT NONE
   INTEGER, INTENT(IN)  :: I,JJ
   INTEGER              :: J1,J2,K
   REAL(SP)             :: DELTX,DELTY,ALPHA1
   REAL(SP)             :: UTMP,VTMP
   REAL(SP), INTENT(OUT):: UAKK(KB),VAKK(KB)

   uakk = 0.0_sp; vakk = 0.0_sp

   IF(isbce(i) /= 2)THEN
     j1 = jj+1-int((jj+1)/4)*3
     j2 = jj+2-int((jj+2)/4)*3
     deltx = vx(nv(i,j1))-vx(nv(i,j2))
     delty = vy(nv(i,j1))-vy(nv(i,j2))

     alpha1 = atan2(delty,deltx)

     DO k = 1,kbm1
       utmp = u(i,k)*cos(alpha1)+v(i,k)*sin(alpha1)
       vtmp = -u(i,k)*sin(alpha1)+v(i,k)*cos(alpha1)

!       VTMP = -VTMP
       vtmp = 0.0_sp

       uakk(k) = utmp*cos(alpha1)-vtmp*sin(alpha1)
       vakk(k) = utmp*sin(alpha1)+vtmp*cos(alpha1)
     END DO
     
   ELSE IF(isbce(i) == 2)THEN
     uakk = u(i,:)
     vakk = v(i,:)
   END IF

   RETURN

Here is the caller graph for this function: