Nodes of different colours represent the following:
Solid arrows point from a parent (sub)module to the submodule which is descended from it. Dashed arrows point from a module being used to the module or program unit using it. Where possible, edges connecting nodes are given different colours to make them easier to distinguish in large graphs.
Type | Intent | Optional | Attributes | Name | ||
---|---|---|---|---|---|---|
integer, | intent(in) | :: | SurfNum | |||
integer, | intent(in) | :: | BeamDIffFlag | |||
real(kind=r64), | intent(out) | :: | CFSAbs(CFSMAXNL+1,2) |
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.
SUBROUTINE CalcEQLOpticalProperty(SurfNum, BeamDIffFlag, CFSAbs)
! SUBROUTINE INFORMATION:
! AUTHOR Bereket Nigusse
! DATE WRITTEN May 2013
! MODIFIED na
!
! RE-ENGINEERED na
! PURPOSE OF THIS SUBROUTINE:
! Calculates the system optical properties from the individual layers
! properties at each time step. The values returned are the layer-by-layer
! absorptance and system transmittance for both beam and diffuse radiation.
! METHODOLOGY EMPLOYED:
! Uses the net radiation method developed for ASHWAT fenestration
! model (ASHRAE RP-1311) by John Wright, the University of WaterLoo
! REFERENCES:
! na
! USE STATEMENTS:
USE DataEnvironment, ONLY: SOLCOS
USE DaylightingManager, ONLY: ProfileAngle
IMPLICIT NONE ! Enforce explicit typing of all variables in this routine
! SUBROUTINE ARGUMENT DEFINITIONS:
INTEGER, INTENT(IN) :: SurfNum
INTEGER, INTENT(IN) :: BeamDIffFlag ! identifier index of diffuse and beam SW radiation
REAL(r64), INTENT(OUT) :: CFSAbs(CFSMAXNL+1,2) ! absorbed beam solar radiation by layers fraction
! SUBROUTINE PARAMETER DEFINITIONS:
! na
! INTERFACE BLOCK SPECIFICATIONS
! na
! DERIVED TYPE DEFINITIONS
! na
! SUBROUTINE LOCAL VARIABLE DECLARATIONS:
!
REAL(r64) :: ProfAngHor ! Solar profile angle (radians) for horizontal blind
REAL(r64) :: ProfAngVer ! Solar profile angle (radians) for vertical blind
REAL(r64) :: IncAng ! incident angle degree
REAL(r64) :: IncidAngle ! = ACOS(SOLCOS(3))
REAL(r64) :: Abs1( CFSMAXNL+1, 2) !
INTEGER :: Lay ! window layer index
INTEGER :: EQLNum ! equivalent layer window construction index
INTEGER :: ConstrNum ! construction index
INTEGER :: I ! index
INTEGER :: J ! index
! Flow
CFSAbs = 0.0d0
ProfAngHor = 0.0d0
ProfAngVer = 0.0d0
ConstrNum = Surface(SurfNum)%Construction
EQLNum = Construct(Surface(SurfNum)%Construction)%EQLConsPtr
IF ( BeamDIffFlag /= isDIFF) THEN
IF(CosIncAng(SurfNum,HourOfDay,TimeStep) <= 0.0d0) RETURN
DO Lay = 1, CFS(EQLNum)%NL
IF ( ISVBLayer(CFS(EQLNum)%L(Lay)) )THEN
IF (CFS(EQLNum)%L(Lay)%LTYPE == ltyVBHOR) THEN
CALL ProfileAngle(SurfNum,SOLCOS,Horizontal,ProfAngHor)
ELSEIF(CFS(EQLNum)%L(Lay)%LTYPE == ltyVBVER) THEN
CALL ProfileAngle(SurfNum,SOLCOS,Vertical,ProfAngVer)
ENDIF
ENDIF
END DO
! Incident angle
IncAng = ACOS(CosIncAng(SurfNum,HourOfDay,TimeStep))
CALL CalcEQLWindowOpticalProperty(CFS(EQLNum), BeamDIffFlag, Abs1, IncAng, ProfAngVer, ProfAngHor)
CFSAbs(1:CFSMAXNL+1,1)= Abs1(1:CFSMAXNL+1, 1)
CFSAbs(1:CFSMAXNL+1,2)= Abs1(1:CFSMAXNL+1, 2)
ELSE
IF ( EQLDiffPropFlag(EQLNum) ) THEN
DO Lay = 1, CFS(EQLNum)%NL
IF ( ISVBLayer(CFS(EQLNum)%L(Lay)) )THEN
IF (CFS(EQLNum)%L(Lay)%LTYPE == ltyVBHOR) THEN
CALL ProfileAngle(SurfNum,SOLCOS,Horizontal,ProfAngHor)
ELSEIF(CFS(EQLNum)%L(Lay)%LTYPE == ltyVBVER) THEN
CALL ProfileAngle(SurfNum,SOLCOS,Vertical,ProfAngVer)
ENDIF
ENDIF
END DO
CALL CalcEQLWindowOpticalProperty( CFS(EQLNum), BeamDIffFlag, Abs1, IncAng, ProfAngVer, ProfAngHor)
CFSAbs(1:CFSMAXNL+1,:) = Abs1(1:CFSMAXNL+1, :)
CFSDiffAbsTrans(EQLNum,1:CFSMAXNL+1,:) = Abs1(1:CFSMAXNL+1,:)
Construct(ConstrNum)%TransDiff = Abs1(CFS(EQLNum)%NL+1, 1)
Construct(ConstrNum)%AbsDiffFrontEQL(1:CFSMAXNL) = Abs1(1:CFSMAXNL,1)
Construct(ConstrNum)%AbsDiffBackEQL(1:CFSMAXNL) = Abs1(1:CFSMAXNL,2)
Construct(ConstrNum)%ReflectSolDiffFront = CFS(EQLNum)%L(1)%SWP_EL%RHOSFDD
Construct(ConstrNum)%ReflectSolDiffBack = CFS(EQLNum)%L(CFS(EQLNum)%NL)%SWP_EL%RHOSBDD
IF (.NOT. CFS(EQLNum)%ISControlled) EQLDiffPropFlag(EQLNum) = .FALSE.
ELSE
CFSAbs(1:CFSMAXNL+1,:) = CFSDiffAbsTrans(EQLNum,1:CFSMAXNL+1,:)
Construct(ConstrNum)%TransDiff = CFSDiffAbsTrans(EQLNum,CFS(EQLNum)%NL+1,1)
Construct(ConstrNum)%AbsDiffFrontEQL(1:CFSMAXNL) = CFSAbs(1:CFSMAXNL,1)
Construct(ConstrNum)%AbsDiffBackEQL(1:CFSMAXNL) = CFSAbs(1:CFSMAXNL,2)
ENDIF
ENDIF
RETURN
END SUBROUTINE CalcEQLOpticalProperty