calc_vort.f90 File Reference

Go to the source code of this file.

Functions/Subroutines
subroutine	calc_vort (VTYPE)

Function/Subroutine Documentation

calc_vort()

subroutine calc_vort

(

integer, intent(in), optional

VTYPE

)

Definition at line 50 of file calc_vort.f90.

!==============================================================================|
   USE all_vars
   IMPLICIT NONE
   INTEGER, OPTIONAL, INTENT(IN) :: VTYPE
   REAL(SP) :: XFLUX(0:M)
   REAL(SP) :: UIJ,VIJ,EXFLUX,AVE
   INTEGER  :: I,J,IA,IB,I1,CNT,JNODE,MY_VTYPE
   INTEGER, ALLOCATABLE  :: LIST(:)
!==============================================================================|

  !---------------------------------------------
  !Process input arguments
  !---------------------------------------------
  my_vtype = 1
  IF(present(vtype))Then
    my_vtype = vtype
  ENDIF

  !---------------------------------------------
  !Initialize Fields
  !---------------------------------------------
  ALLOCATE(list(m)) ; list = 0
  vort = 0.0

  !===================================================================
  ! Calculate vorticity at all nodes using c.v. half edges
  !   Note minus sign (nodes are order clockwise [against convention])
  !===================================================================

  SELECT CASE(my_vtype)

  !----------------------------------------------
  CASE(1) !Curl of Vert-Avged Velocity
  !----------------------------------------------
  DO i=1,ncv
    i1  = ntrg(i)
    ia  = niec(i,1)
    ib  = niec(i,2)
    uij = ua(i1)
    vij = va(i1) 
    exflux = -(uij*dltxe(i) + vij*dltye(i))  
    vort(ia) = vort(ia)-exflux
    vort(ib) = vort(ib)+exflux
  END DO
  !----------------------------------------------
  CASE(2) !Curl of Flux
  !----------------------------------------------
  DO i=1,ncv
    i1  = ntrg(i)
    ia  = niec(i,1)
    ib  = niec(i,2)
    uij = ua(i1)
    vij = va(i1) 
    exflux = -d1(i1)*(uij*dltxe(i) + vij*dltye(i))  
    vort(ia) = vort(ia)-exflux
    vort(ib) = vort(ib)+exflux
  END DO
  !----------------------------------------------
  CASE DEFAULT
    !ERROR AND HALT 
  END SELECT

  vort = vort/art1

  !===================================================================
  !Correction at Boundaries (May no longer be necessary for 2.5+)
  !===================================================================

   do i=1,m
     if(isonb(i) > 0)then
     ave = 0.
     cnt = 0
     do j=1,ntsn(i)
       jnode = nbsn(i,j)
       if(jnode /= i .and. isonb(jnode) == 0)then
         ave = ave + vort(jnode)
         cnt = cnt + 1
       end if
     end do
     vort(i) = ave/float(cnt)
     if(cnt == 0) list(i) = 1
     end if
   end do

   do i=1,m
     if(list(i) > 0)then
     ave = 0.
     cnt = 0
     do j=1,ntsn(i)
       jnode = nbsn(i,j)
       if(jnode /= i .and. list(jnode) == 0)then
         ave = ave + vort(jnode)
         cnt = cnt + 1
       end if
     end do
     vort(i) = ave/float(cnt)
     end if
   end do

       
   RETURN