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

!==============================================================================|
!  Compute Cartesian Vertical Velocity                                         |
!==============================================================================|
   SUBROUTINE wreal               
!==============================================================================|
   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"

   END SUBROUTINE wreal
!==============================================================================|