wreal.f90 File Reference

Go to the source code of this file.

Functions/Subroutines
subroutine	wreal

Function/Subroutine Documentation

wreal()

subroutine wreal

(

)

Definition at line 44 of file wreal.f90.

!==============================================================================|
   USE all_vars
   USE mod_utils
   USE mod_wd
   USE mod_nesting
   IMPLICIT NONE
   REAL(SP) :: DDDX,DDDY,DEDX,DEDY,ETF1AA,WW1,WW2 
   INTEGER  :: I,K,J1,J2,J3
   INTEGER  :: J, II
   REAL(SP) :: U_TMP, V_TMP
!==============================================================================|

   IF(dbg_set(dbg_sbr)) WRITE(ipt,*) "Start: wreal"

!------------------------------------------------------------------------------!
!  SAVE OMEGA VELOCITY FROM PREVIOUS TIMESTEP (USED FOR LAGRANGIAN TRACKING)   !
!------------------------------------------------------------------------------!

   wtts = wts

!------------------------------------------------------------------------------!
!  CALCULATE A NEW OMEGA VELOCITY                                              !
!------------------------------------------------------------------------------!

!!===========yding==================


   DO i=1,n
    IF(iswetc(i) == 1)THEN
     j1=nv(i,1)
     j2=nv(i,2)
     j3=nv(i,3)
     dddx=awx(i,1) * d(j1)+awx(i,2) * d(j2)+awx(i,3)*d(j3)
     dddy=awy(i,1) * d(j1)+awy(i,2) * d(j2)+awy(i,3)*d(j3)
     dedx=awx(i,1)*elf(j1)+awx(i,2)*elf(j2)+awx(i,3)*elf(j3)
     dedy=awy(i,1)*elf(j1)+awy(i,2)*elf(j2)+awy(i,3)*elf(j3)
     etf1aa=one_third*(el(nv(i,1))+el(nv(i,2))+el(nv(i,3)))
     DO k=1,kbm1
      ww1=0.5_sp*(w(i,k)+w(i,k+1))+u(i,k)*(zz1(i,k)*dddx+dedx)+ &
                                v(i,k)*(zz1(i,k)*dddy+dedy)
      ww2=(zz1(i,k)+1.)*(etf1aa-et1(i))/dti
      ww(i,k)=ww1+ww2
     END DO
    ELSE
     DO k=1,kbm1
      ww(i,k)=0.0_sp
     END DO
    END IF
   END DO



   IF(dbg_set(dbg_sbr)) WRITE(ipt,*) "End: wreal"

Here is the call graph for this function:

Here is the caller graph for this function: