adjust2d3d.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$
!/===========================================================================/

   SUBROUTINE adjust2d3d(ADJUST_TYPE)
!==============================================================================|
!    ADJUST 3D VELOCITY USING DEFECT BETWEEN UPDATED AND CURRENT VERTICALLY    !
!    AVERAGED VELOCITIES                               !
!                                          !
!    FORMULA IS:                                   !
!                                          !
!      U_adjusted = U_orig + eps*(U_avg_new - U_avg_current)               !
!      eps = 0 : no adjustment                             !
!      eps = 1 : full adjustment                           !
!==============================================================================|
   USE mod_wd
   USE all_vars
   USE mod_utils
   IMPLICIT NONE
   INTEGER, INTENT(IN) :: ADJUST_TYPE
   INTEGER :: I,K
   REAL(SP), PARAMETER :: EPS = 1.0_sp
   REAL(SP) :: UAC,VAC,UTMP,VTMP
!==============================================================================!

   if(dbg_set(dbg_sbr)) write(ipt,*) "Start: adjust2d3d"

   SELECT CASE(adjust_type)

   CASE(1)
   DO i=1,nt
     uac    = sum(u(i,1:kbm1)*dz1(i,1:kbm1))
     vac    = sum(v(i,1:kbm1)*dz1(i,1:kbm1))
     u(i,1:kbm1) = u(i,1:kbm1) + eps*(ua(i) - uac) 
     v(i,1:kbm1) = v(i,1:kbm1) + eps*(va(i) - vac) 
   END DO

   CASE(2)
   uard = uard/float(isplit)
   vard = vard/float(isplit)
!!#  if defined (1)
!!   UARDS = UARDS/FLOAT(ISPLIT)
!!   VARDS = VARDS/FLOAT(ISPLIT)
!!#  endif

                                                                                                                         
   DO i=1,nt
     IF(iswetct(i)*iswetc(i) == 1)THEN
       utmp = 0.0_sp ; vtmp = 0.0_sp
       DO k=1,kbm1
         utmp = utmp + u(i,k)*dz1(i,k)
         vtmp = vtmp + v(i,k)*dz1(i,k)
       END DO
       utmp = utmp*dt1(i)
       vtmp = vtmp*dt1(i)
       DO k=1,kbm1
         u(i,k) = u(i,k) - (utmp-uard(i))/dt1(i)
         v(i,k) = v(i,k) - (vtmp-vard(i))/dt1(i)
       END DO
     END IF
   END DO

!!#  if defined (1)
!!   DO I=1,NT
!!     UTMP = 0.0_SP ; VTMP = 0.0_SP
!!     DO K=1,KBM1
!!       UTMP = UTMP + U(I,K)*DZ1(I,K)
!!       VTMP = VTMP + V(I,K)*DZ1(I,K)
!!     END DO
!!     UTMP = UTMP*DT1(I)
!!     VTMP = VTMP*DT1(I)
!!     DO K=1,KBM1
!!       US(I,K) = U(I,K) - (UTMP-UARDS(I))/DT1(I)
!!       VS(I,K) = V(I,K) - (VTMP-VARDS(I))/DT1(I)
!!     END DO
!!   END DO
!!#  endif

   END SELECT

   if(dbg_set(dbg_sbr)) write(ipt,*) "End: adjust2d3d"

   END SUBROUTINE adjust2d3d
!==============================================================================|