Type | Intent | Optional | Attributes | Name | ||
---|---|---|---|---|---|---|
real(kind=r64), | intent(in) | :: | WindAt10m | |||
real(kind=r64), | intent(in) | :: | WindDir | |||
real(kind=r64), | intent(in) | :: | SurfAzimuth |
Nodes of different colours represent the following:
Solid arrows point from a procedure to one which it calls. Dashed arrows point from an interface to procedures which implement that interface. This could include the module procedures in a generic interface or the implementation in a submodule of an interface in a parent module. Where possible, edges connecting nodes are given different colours to make them easier to distinguish in large graphs.
FUNCTION CalcBlockenWindward(WindAt10m, WindDir, SurfAzimuth) RESULT (Hf)
! FUNCTION INFORMATION:
! AUTHOR Brent Griffith
! DATE WRITTEN Aug 2010
! MODIFIED na
! RE-ENGINEERED na
! PURPOSE OF THIS FUNCTION:
! calculate model equation for forced convection using Blocken correlation
! METHODOLOGY EMPLOYED:
! encapsulate model in function
! REFERENCES:
! Blocken, B., T. Defraeye, D. Derome, J. Carmeliet. 2009.
! High-Resolution CFD Simulations for Forced Convection
! Heat Transfer Coefficients at the Facade of a Low-Rise Building.
! Building and Environment 44 (2009) 2396 - 2412.
! USE STATEMENTS:
! na
IMPLICIT NONE ! Enforce explicit typing of all variables in this routine
! FUNCTION ARGUMENT DEFINITIONS:
REAL(r64) , INTENT(IN) :: WindAt10m
REAL(r64) , INTENT(IN) :: WindDir ! Wind direction measured clockwise from geographhic North
REAL(r64) , INTENT(IN) :: SurfAzimuth ! or Facing, Direction the surface outward normal faces (degrees)
REAL(r64) :: Hf
! FUNCTION PARAMETER DEFINITIONS:
! na
! INTERFACE BLOCK SPECIFICATIONS:
! na
! DERIVED TYPE DEFINITIONS:
! na
! FUNCTION LOCAL VARIABLE DECLARATIONS:
REAL(r64) :: Theta ! angle between wind and surface azimuth
Theta = WindDir - SurfAzimuth - 90.d0 !TODO double check theta
IF (Theta > 180.0d0) Theta = Theta - 360.00
IF (Theta <= 11.25d0) THEN
Hf = 4.6d0 * (WindAt10m**0.89d0)
ELSEIF ((11.25d0 < Theta) .AND. (Theta <= 33.75d0)) THEN
Hf = 5.d0 * (WindAt10m**0.8d0)
ELSEIF ((33.75d0 < Theta) .AND. (Theta <= 56.25d0)) THEN
Hf = 4.6d0 * (WindAt10m**0.84d0)
ELSEIF ((56.25d0 < Theta) .AND. (Theta <= 100.d0)) THEN
Hf = 4.5d0 * (WindAt10m**0.81d0)
ELSE
! should not be used for leeward... check why come here?
Hf = 3.54d0 * (WindAt10m**0.76d0) !emmel model for robustness?
ENDIF
RETURN
END FUNCTION CalcBlockenWindward