ghostuv.f90

Go to the documentation of this file.

!/===========================================================================/
! Copyright (c) 2007, The University of Massachusetts Dartmouth 
! Produced at the School of Marine Science & Technology 
! Marine Ecosystem Dynamics Modeling group
! All rights reserved.
!
! FVCOM has been developed by the joint UMASSD-WHOI research team. For 
! details of authorship and attribution of credit please see the FVCOM
! technical manual or contact the MEDM group.
!
! 
! This file is part of FVCOM. For details, see http://fvcom.smast.umassd.edu 
! The full copyright notice is contained in the file COPYRIGHT located in the 
! root directory of the FVCOM code. This original header must be maintained
! in all distributed versions.
!
! THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS "AS IS" 
! AND ANY EXPRESS OR  IMPLIED WARRANTIES, INCLUDING,  BUT NOT  LIMITED TO,
! THE IMPLIED WARRANTIES OF MERCHANTABILITY AND  FITNESS FOR A PARTICULAR
! PURPOSE ARE DISCLAIMED.  
!
!/---------------------------------------------------------------------------/
! CVS VERSION INFORMATION
! $Id$
! $Name$
! $Revision$
!/===========================================================================/

!==============================================================================|
!   CALCULATE GHOST VELOCITY FOR EXTERNAL MODE                                 !
!==============================================================================|
   SUBROUTINE ghostuv2(I,JJ,UAKK,VAKK)
!==============================================================================|

   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
   END SUBROUTINE ghostuv2


!==============================================================================|
!   CALCULATE GHOST VELOCITY FOR INTERNAL MODE                                 !
!==============================================================================|
   SUBROUTINE ghostuv3(I,JJ,UAKK,VAKK)
!==============================================================================|

   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
   END SUBROUTINE ghostuv3