SUBROUTINE CALHTC(TIME,XIN,INP,INFO) C******************************************************************************* C THIS SUBROUTINE CALCULATES FOR ALL SURFACES WITHIN A ZONE C HEAT TRANSFER COEFFICIENTS FOR THE INSIDE SURFACE AND C THE CORRESPONDING OUTSIDE SURFACE OF THAT WALL OR WINDOW C C THIS ROUTINE IS CALLED BY: TYPE56 C THIS ROUTINE CALLS: VALUE C C ADAPTED FROM TYPE80 IN 2004 BY S.HOLST. C******************************************************************************* Cak V16 USE TrnsysData, ONLY:LUR,LUW,IFORM,LUK use set_unitNo Cak V16 USE BPROPS USE VAR USE BGEN USE BREG USE BSURFACE USE BTYPES USE BTRANS USE BMATRIX * SH:15.2.2004 variable declaration -- begin -- implicit none * dummy arguments Cak V16 c real TIME Cak V16 * local variables integer I integer N1 integer N2 integer NS integer N integer NTYPE integer iHTC real HTC real VALUE real dTSA real KFLOORUP real EFLOORUP real KFLOORDOWN real EFLOORDOWN real KCEILUP real ECEILUP real KCEILDOWN real ECEILDOWN real KVERTICAL real EVERTICAL * common block variables * /LUNITS/ Cak V16 c integer LUK c integer IFORM c integer LUW c integer LUR * SH:15.2.2004 variable declaration -- end -- c DOUBLE PRECISION XIN c ak INTEGER*4 INP,INFO INTEGER*4 INP c ak DIMENSION XIN(*),INFO(15) real(8), intent(in) :: XIN(*) integer, intent(inout) :: info(15) real(8) :: time !TW integer :: unitNo, typeNo CHARACTER(len=500) Message c ak COMMON /LUNITS/ LUR,LUW,IFORM,LUK Cak V16 CSH 03/04 set parameters CSH KFLOORUP=2.0; EFLOORUP=0.31; KFLOORDOWN=1.08; EFLOORDOWN=0.31 CSH KCEILUP=1.08; ECEILUP=0.31; KCEILDOWN=2.0; ECEILDOWN=0.31 CSH KVERTICAL=1.6; EVERTICAL=0.3 KFLOORUP=HTC_COEF(1)/3.6; EFLOORUP=HTC_COEF(2) KFLOORDOWN=HTC_COEF(3)/3.6; EFLOORDOWN=HTC_COEF(4) KCEILUP=HTC_COEF(5)/3.6; ECEILUP=HTC_COEF(6) KCEILDOWN=HTC_COEF(7)/3.6; ECEILDOWN=HTC_COEF(8) KVERTICAL=HTC_COEF(9)/3.6; EVERTICAL=HTC_COEF(10) Cak V16 !TW unitNo = info(1) !TW typeNo = info(2) Cak V16 DO I=1,NTZONE !loop all (regular+fictive) zones C *** Walls *** IF(NZWALL(I).GT.0) THEN ! if walls N1=IPWALL(I) N2=N1+NZWALL(I)-1 DO N=N1,N2 ! all walls NS=NSWALL(N) NTYPE=NSTYPE(NS) IF (NFACE(NS).EQ.-1) THEN ! back side of adj/int wall HTC=HCO(NTYPE) ELSE ! front side of adj/int or other wall HTC=HCI(NTYPE) ENDIF ! back or front side of wall IF (HTC.GT.-5500) THEN ! standard HTC calc for inside surface HTCI(NS)=VALUE(TIME,HTC,INP,XIN) CSH 08/04 KD IF (HTC.EQ.-4000) THEN ! gap calculation of active layer kuehldecke HTCI(NS)=UC1(NS) ! use stored 1/R1 value from kuehldecken calc for gap calculation ENDIF ! gap calculation of active layer kuehldecke CSH 08/04 KD ELSE ! automatic HTC calc for inside surface iHTC=INT(HTC-0.1) dTSA=TSI(NS)-TBAR(I) IF (iHTC.EQ.-6000) THEN ! floor surface IF (dTSA.GT.0.) THEN ! floor heat flow up HTCI(NS)=MAX(1.,KHTC(NS)*KFLOORUP*(ABS(dTSA)) . **EFLOORUP)*3.6 ELSE ! floor heat flow down HTCI(NS)=MAX(1.,KHTC(NS)*KFLOORDOWN*(ABS(dTSA)) . **EFLOORDOWN)*3.6 ENDIF ! heat flow direction ENDIF ! floor surface IF (iHTC.EQ.-7000) THEN ! ceiling surface IF (dTSA.GT.0.) THEN ! floor heat flow down HTCI(NS)=MAX(1.,KHTC(NS)*KCEILDOWN*(ABS(dTSA)) . **ECEILDOWN)*3.6 ELSE ! floor heat flow up HTCI(NS)=MAX(1.,KHTC(NS)*KCEILUP*(ABS(dTSA)) . **ECEILUP)*3.6 ENDIF ! heat flow direction ENDIF ! ceiling surface IF (iHTC.EQ.-8000) THEN ! vertical surface HTCI(NS)=MAX(1.,KHTC(NS)*KVERTICAL*(ABS(dTSA)) . **EVERTICAL)*3.6 ENDIF ! vertical surface IF (iHTC.LT.-8000) THEN ! error cak V16 write(Message,'("Error in automatic heat transfer calculation & for surface",i4," or its adjacent surface!")') NS CALL MESSAGES(-1,Message,'FATAL',unitNo,typeNo) RETURN c WRITE(LUW,1600) INFO(1),INFO(2),NS ! error in HTC calc for surface NS c CALL MYSTOP (1001) c RETURN cak V16 ENDIF ! error in HTC calc ENDIF ! standard or automatic HTC calcs for inside surface IF (NFACE(NS).EQ.-1) THEN ! back side of adj/int wall HTC=HCI(NTYPE) ELSE ! front side of adj/int or other wall HTC=HCO(NTYPE) ENDIF ! back or front side of wall IF (HTC.GT.-5500) THEN ! standard HTC calc for outside surface HTCO(NS)=VALUE(TIME,HTC,INP,XIN) ELSE ! automatic HTC calc for outside surface (if adj or int surface) iHTC=INT(HTC-0.1) SELECT CASE (CONTYPE(NS)) ! which wall type CASE (1,4) ! external or boundary cak V16 write(Message,'("Error in automatic heat transfer calculation & for surface",i4," or its adjacent surface!: & Automatic heat transfer calculation & is only allowed on surfaces within a zone. & Please check building input file (.bui)!")') NS CALL MESSAGES(-1,Message,'FATAL',unitNo,typeNo) RETURN c WRITE(LUW,1600) INFO(1),INFO(2),NS ! no automatic HTC calc on outside surf of ext or bound wall c CALL MYSTOP (1001) c RETURN cak V16 CASE (2) ! adjacent dTSA=TSO(NS)-TBAR(INT(BOUCON(NS)/1000.)) CASE (3) ! internal dTSA=TSO(NS)-TBAR(I) CASE (5) ! identical boundaries dTSA=TSO(NS)-TBAR(I) END SELECT ! which wall type IF (iHTC.EQ.-6000) THEN ! floor surface IF (dTSA.GT.0.) THEN ! floor heat flow up HTCO(NS)=MAX(1.,KHTC(NS)*KFLOORUP*(ABS(dTSA)) . **EFLOORUP)*3.6 ELSE ! floor heat flow down HTCO(NS)=MAX(1.,KFLOORDOWN*(ABS(dTSA))**EFLOORDOWN)*3.6 ENDIF ! heat flow direction ENDIF ! floor surface IF (iHTC.EQ.-7000) THEN ! ceiling surface IF (dTSA.GT.0.) THEN ! floor heat flow down HTCO(NS)=MAX(1.,KHTC(NS)*KCEILDOWN*(ABS(dTSA)) . **ECEILDOWN)*3.6 ELSE ! floor heat flow up HTCO(NS)=MAX(1.,KHTC(NS)*KCEILUP*(ABS(dTSA)) . **ECEILUP)*3.6 ENDIF ! heat flow direction ENDIF ! ceiling surface IF (iHTC.EQ.-8000) THEN ! vertical surface HTCO(NS)=MAX(1.,KHTC(NS)*KVERTICAL*(ABS(dTSA)) . **EVERTICAL)*3.6 ENDIF ! vertical surface IF (iHTC.LT.-8000) THEN ! error cak V16 write(Message,'("Error in automatic heat transfer calculation & for surface",i4," or its adjacent surface!: & Automatic heat transfer calculation & is only allowed on surfaces within a zone. & Please check building input file (.bui)!")') NS CALL MESSAGES(-1,Message,'FATAL',unitNo,typeNo) RETURN c WRITE(LUW,1600) INFO(1),INFO(2),NS ! error in HTC calc for surface NS c CALL MYSTOP (1001) c RETURN cak V16 ENDIF ! error in HTC calc ENDIF ! standard or automatic HTC calcs for outside surface ENDDO ! all walls ENDIF ! if walls C *** Windows *** IF(NZWIN(I).GT.0) THEN ! if windows N1=IPWIN(I) N2=N1+NZWIN(I)-1 DO N=N1,N2 ! all windows NS=NSWIN(N) NTYPE=NSTYPE(NS) IF (NFACE(NS).EQ.-1) THEN ! back side of int window HTC=DWIN(3,NTYPE) ELSE ! front side of int or external window HTC=DWIN(2,NTYPE) ENDIF ! back or front side of window IF (HTC.GT.-5500) THEN ! standard HTC calc inside glass surface HTCI(NS)=VALUE(TIME,DWIN(2,NTYPE),INP,XIN) ELSE ! automatic HTC calc iHTC=INT(HTC-0.1) dTSA=TSI(NS)-TBAR(I) IF (iHTC.EQ.-6000) THEN ! floor surface IF (dTSA.GT.0.) THEN ! floor heat flow up HTCI(NS)=MAX(1.,KHTC(NS)*KFLOORUP*(ABS(dTSA)) . **EFLOORUP)*3.6 ELSE ! floor heat flow down HTCI(NS)=MAX(1.,KHTC(NS)*KFLOORDOWN*(ABS(dTSA)) . **EFLOORDOWN)*3.6 ENDIF ! heat flow direction ENDIF ! floor surface IF (iHTC.EQ.-7000) THEN ! ceiling surface IF (dTSA.GT.0.) THEN ! floor heat flow down HTCI(NS)=MAX(1.,KHTC(NS)*KCEILDOWN*(ABS(dTSA)) . **ECEILDOWN)*3.6 ELSE ! floor heat flow up HTCI(NS)=MAX(1.,KHTC(NS)*KCEILUP*(ABS(dTSA)) . **ECEILUP)*3.6 ENDIF ! heat flow direction ENDIF ! ceiling surface IF (iHTC.EQ.-8000) THEN ! vertical surface HTCI(NS)=MAX(1.,KHTC(NS)*KVERTICAL*(ABS(dTSA)) . **EVERTICAL)*3.6 ENDIF ! vertical surface IF (iHTC.LT.-8000) THEN ! error cak V16 write(Message,'("Error in automatic heat transfer calculation & for surface",i4," or its adjacent surface!: & Automatic heat transfer calculation & is only allowed on surfaces within a zone. & Please check building input file (.bui)!")') NS CALL MESSAGES(-1,Message,'FATAL',unitNo,typeNo) RETURN c WRITE(LUW,1600) INFO(1),INFO(2),NS ! error in HTC calc for surface NS c CALL MYSTOP (1001) c RETURN cak V16 ENDIF ! error in HTC calc ENDIF ! standard or automatic HTC calcs for inside glass surface IF (NFACE(NS).EQ.-1) THEN ! back side of int window HTC=DWIN(2,NTYPE) ELSE ! front side of int or external window HTC=DWIN(3,NTYPE) ENDIF ! back or front side of window IF (HTC.GT.-5500) THEN ! standard HTC calc outside glass surface HTCO(NS)=VALUE(TIME,HTC,INP,XIN) ELSE ! automatic HTC calc SELECT CASE (CONTYPE(NS)) ! which window type CASE (1) ! external cak V16 write(Message,'("Error in automatic heat transfer calculation & for surface",i4," or its adjacent surface!: & Automatic heat transfer calculation & is only allowed on surfaces within a zone. & Please check building input file (.bui)!")') NS CALL MESSAGES(-1,Message,'FATAL',unitNo,typeNo) RETURN c WRITE(LUW,1600) INFO(1),INFO(2),NS ! no automatic HTC calc on outside surf of ext window c CALL MYSTOP (1001) c RETURN CASE (2) ! adjacent iHTC=INT(HTC-0.1) dTSA=TSI(NADJ(NS))-TBAR(INT(BOUCON(NS)/1000.)) END SELECT ! which window type IF (iHTC.EQ.-6000) THEN ! floor surface IF (dTSA.GT.0.) THEN ! floor heat flow up HTCO(NS)=MAX(1.,KHTC(NS)*KFLOORUP*(ABS(dTSA)) . **EFLOORUP)*3.6 ELSE ! floor heat flow down HTCO(NS)=MAX(1.,KHTC(NS)*KFLOORDOWN*(ABS(dTSA)) . **EFLOORDOWN)*3.6 ENDIF ! heat flow direction ENDIF ! floor surface IF (iHTC.EQ.-7000) THEN ! ceiling surface IF (dTSA.GT.0.) THEN ! floor heat flow down HTCO(NS)=MAX(1.,KHTC(NS)*KCEILDOWN*(ABS(dTSA)) . **ECEILDOWN)*3.6 ELSE ! floor heat flow up HTCO(NS)=MAX(1.,KHTC(NS)*KCEILUP*(ABS(dTSA)) . **ECEILUP)*3.6 ENDIF ! heat flow direction ENDIF ! ceiling surface IF (iHTC.EQ.-8000) THEN ! vertical surface HTCO(NS)=MAX(1.,KHTC(NS)*KVERTICAL*(ABS(dTSA)) . **EVERTICAL)*3.6 ENDIF ! vertical surface IF (iHTC.LT.-8000) THEN ! error cak V16 write(Message,'("Error in automatic heat transfer calculation & for surface",i4," or its adjacent surface!: & Automatic heat transfer calculation & is only allowed on surfaces within a zone. & Please check building input file (.bui)!")') NS CALL MESSAGES(-1,Message,'FATAL',unitNo,typeNo) RETURN c WRITE(LUW,1600) INFO(1),INFO(2),NS ! error in HTC calc for surface NS c CALL MYSTOP (1001) c RETURN ENDIF ! error in HTC calc ENDIF ! standard or automatic HTC calcs for outside glass surface ENDDO ! all windows ENDIF ! if windows ENDDO !loop all regular+fictive zones RETURN C1600 FORMAT(//,1X,'***** ERROR *****',8X,/1X, C .'UNIT ',I3,' TYPE ',I3,' MULTIZONE BUILDING',/1X, C .'ERROR IN AUTOMATIC HEAT TRANSFER CALCULATION FOR SURFACE ',I4, C .' OR ITS ADJACENT SURFACE') END ! CALHTC