mod_bulk Module Reference

Functions/Subroutines
subroutine	asimple_drag (spdx, spdy, strx, stry)
subroutine	psimple_drag (spdx, spdy, strx, stry)

Function/Subroutine Documentation

asimple_drag()

subroutine mod_bulk::asimple_drag	(	real(sp), dimension(:), intent(in), allocatable, target	spdx,
		real(sp), dimension(:), intent(in), allocatable, target	spdy,
		real(sp), dimension(:), intent(inout), allocatable, target	strx,
		real(sp), dimension(:), intent(inout), allocatable, target	stry
	)

Definition at line 49 of file mod_bulk.f90.

    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)

Here is the call graph for this function:

Here is the caller graph for this function:

psimple_drag()

subroutine mod_bulk::psimple_drag	(	real(sp), dimension(:), intent(in), pointer	spdx,
		real(sp), dimension(:), intent(in), pointer	spdy,
		real(sp), dimension(:), intent(inout), pointer	strx,
		real(sp), dimension(:), intent(inout), pointer	stry
	)

Definition at line 70 of file mod_bulk.f90.

    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
    

Here is the caller graph for this function: