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

MODULE mod_bulk
  USE mod_prec
  implicit none

CONTAINS



  SUBROUTINE asimple_drag(spdx,spdy,strx,stry)
    IMPLICIT NONE
    REAL(SP),ALLOCATABLE, TARGET, INTENT(IN)  :: SPDX(:),SPDY(:)
    REAL(SP),ALLOCATABLE, TARGET, INTENT(INOUT) :: STRX(:),STRY(:)

    REAL(SP), POINTER :: SPDXP(:),SPDYP(:)
    REAL(SP), POINTER :: STRXP(:),STRYP(:)
    
    spdxp => spdx
    spdyp => spdy

    strxp => strx
    stryp => stry
    
    CALL psimple_drag(spdxp,spdyp,strxp,stryp)


  END SUBROUTINE asimple_drag



  SUBROUTINE psimple_drag(spdx,spdy,strx,stry)
    IMPLICIT NONE
    REAL(SP), POINTER,INTENT(IN)  :: SPDX(:),SPDY(:)
    REAL(SP), POINTER, INTENT(INOUT) :: STRX(:),STRY(:)
    INTEGER :: I, N
    REAL(SP) :: CD, WDS, TX, TY


    IF(.not.Associated(spdx)) WRITE(6,*) "SIMPLE DRAG: SPDX is not associated"
    IF(.not.Associated(spdy)) WRITE(6,*) "SIMPLE DRAG: SPDY is not associated"
    IF(.not.Associated(strx)) WRITE(6,*) "SIMPLE DRAG: STRX is not associated"
    IF(.not.Associated(stry)) WRITE(6,*) "SIMPLE DRAG: STRY is not associated"

    n = ubound(spdx,1)


    IF(n /= ubound(spdy,1)) WRITE(6,*) "SIMPLE DRAG: MIS-MATCHED DIMENSIONS"
    IF(n /= ubound(stry,1)) WRITE(6,*) "SIMPLE DRAG: MIS-MATCHED DIMENSIONS"
    IF(n /= ubound(strx,1)) WRITE(6,*) "SIMPLE DRAG: MIS-MATCHED DIMENSIONS"

    DO i=1,n
       tx = spdx(i)
       ty = spdy(i)
       wds=sqrt(tx*tx+ty*ty)
       cd=1.2e-3
       IF (wds >= 11.0_sp) cd=(0.49_sp+0.065_sp*wds)*1.e-3_sp
       IF (wds >= 25.0_sp) cd=(0.49_sp+0.065_sp*25.0_sp)*1.e-3_sp

       strx(i) = 1.2_sp*cd*tx*wds
       stry(i) = 1.2_sp*cd*ty*wds

    END DO
    

  END SUBROUTINE psimple_drag
  


END MODULE mod_bulk