adjust2d3d.f90 File Reference

Go to the source code of this file.

Functions/Subroutines
subroutine	adjust2d3d (ADJUST_TYPE)

Function/Subroutine Documentation

adjust2d3d()

subroutine adjust2d3d

(

integer, intent(in)

ADJUST_TYPE

)

Definition at line 41 of file adjust2d3d.f90.

!==============================================================================|
!    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"

Here is the call graph for this function:

Here is the caller graph for this function: