home *** CD-ROM | disk | FTP | other *** search
/ The World of Computer Software / World_Of_Computer_Software-02-387-Vol-3of3.iso / h / house_ii.zip / FOR / TCONC.FOR < prev    next >
Text File  |  1992-04-17  |  3KB  |  124 lines

  1.       SUBROUTINE  TCONC(DEGRAD,DTIM1)
  2. C
  3. C  CALCULATE CONSTANTS USED TO CALCULATE LUMP TEMPERATURES
  4. C  BY FINITE DIFFERENCE EQNS
  5. C
  6. C  LINK THESE CALCULATIONS TO CONSTRUCTS, SINCE NO. OF CONSTANTS ARE
  7. C  A FUNCTION OF NO. OF LUMPS OR NODES WHICH IS SPECIFIED IN CONSTRUCT
  8. C
  9. C - - - PARAMETER CONSTANTS
  10. CMDK NKONST
  11. CMDK NWL
  12. CMDK NZN
  13. CMDK NZW
  14.       PARAMETER (FMODMN=30.)
  15. C
  16. C - - - COMMON BLOCKS
  17. C
  18. CMDK CNSTRK
  19. CMDK ENCBK1
  20. CMDK ENCBLK
  21. CMDK IRDFQ
  22. CMDK IZZQ
  23. CMDK STRUCA
  24. CMDK TCNBLK
  25. CMDK TCONCB
  26. C
  27.       DO 10 I=1,15
  28.    10 KONI(I)=0
  29.       NKON=0
  30. C
  31. C  FIND CONSTRUCT NOS USED IN THIS SIMULATION
  32.       DO 100 J=1,NROOMS
  33.       NRM=NRMA(J)
  34.       N2=NWALLA(NRM)
  35. C
  36.       DO 90 K2=1,N2
  37.       IZW=NENC(NRM,K2)
  38.       KON=KONSTA(IZW)
  39.       IF(NKON.EQ.0)THEN
  40.           KONI(1)=KON
  41.           NKON=1
  42.           GO TO 90
  43.           ENDIF
  44.       IFG=0
  45.       DO 80 I=1,NKON
  46.       IF(KON.EQ.KONI(I))IFG=1
  47.    80 CONTINUE
  48.       IF(IFG.EQ.0)THEN
  49.           NKON=NKON+1
  50.           KONI(NKON)=KON
  51.           ENDIF
  52.    90 CONTINUE
  53.   100 CONTINUE
  54.       NKON=NKON+1
  55.       KONI(NKON)=KONRF(1)
  56.       IF(NATC.EQ.2)THEN
  57.           NKON=NKON+1
  58.           KONI(NKON)=KONRF(2)
  59.           ENDIF
  60. C
  61. C  NOW CALC CONSTANTS, TCON, FOR EACH CONSTRUCT
  62.       DO 150 IK=1,NKON
  63.       KON=KONI(IK)
  64.       NLM=NLMP(KON)
  65. C
  66.       VK1=WK(1,KON)
  67.       VX1=WX(1,KON)
  68.       VRHO1=WRHO(1,KON)
  69.       VCP1=WCP(1,KON)
  70.       IF(NLM.EQ.3)THEN
  71.           VK2=WK(2,KON)
  72.           VX2=WX(2,KON)
  73.           VRHO2=WRHO(2,KON)
  74.           VCP2=WCP(2,KON)
  75.           VK3=WK(3,KON)
  76.           VX3=WX(3,KON)
  77.           VRHO3=WRHO(3,KON)
  78.           VCP3=WCP(3,KON)
  79.           ENDIF
  80. C
  81.       IF(NLM.EQ.1)THEN
  82.           FRCPX=VRHO1*VCP1*VX1
  83.           TCON(1,KON)=8.*DTIM1*(VK1/DEGRAD)/(VX1*FRCPX)
  84.           TCON(2,KON)=DTIM1*4./FRCPX
  85.           TCON(3,KON)=TCON(1,KON)
  86. C  FOURIER MODULUS (MUST BE EQ OR GT 30. FOR STABILITY OF FINITE DIFF CALC
  87. C    THIS IS USED IN IWALLS TO CORRECT POSSIBLE INSTABILITY WHICH MAY
  88. C    BE EXPERIENCED WITH LIGHTWEIGHT INSULATED BARRIER WALLS BETWN ZONES
  89. C    .  THE MINIMUM VALUE OF 30. WAS ESTABLISHED BASED ON EXPERIENCE.
  90.           FMOD=VX1**2*VRHO1*VCP1/(VK1*DTIM1)
  91.           IF(FMOD.LT.FMODMN)THEN
  92.               NKSTAB=NKSTAB+1
  93.               IF(NKSTAB.GT.5)STOP 'TCONC:INCR DIM ON KONSTB AND NRATTM'
  94.               RATTM=FMODMN/FMOD
  95.               NRATTM(NKSTAB)=INT(RATTM)
  96.               IF(FLOAT(NRATTM(NKSTAB)).LT.RATTM)
  97.      +            NRATTM(NKSTAB)=NRATTM(NKSTAB)+1
  98.               KONSTB(NKSTAB)=KON
  99.               WRITE(60,500)KON,NRATTM(NKSTAB),FMOD,FMODMN
  100.               ENDIF
  101.           ENDIF
  102.       IF(NLM.EQ.3)THEN
  103.           FRCPX1=VRHO1*VCP1*VX1
  104.           FRCPX2=VRHO2*VCP2*VX2
  105.           FRCPX3=VRHO3*VCP3*VX3
  106.           TCON(1,KON)=DTIM1*2.*VK1/(VX1*FRCPX1)
  107.           TCON(2,KON)=DTIM1*2./FRCPX1
  108.           FR12=0.5*(FRCPX1+FRCPX2)
  109.           TCON(3,KON)=DTIM1*VK1/(VX1*FR12)
  110.           VKD=VK2/DEGRAD
  111.           TCON(4,KON)=DTIM1*VKD/(VX2*FR12)
  112.           FR23=0.5*(FRCPX2+FRCPX3)
  113.           TCON(5,KON)=DTIM1*VKD/(VX2*FR23)
  114.           TCON(6,KON)=DTIM1*VK3/(VX3*FR23)
  115.           TCON(7,KON)=DTIM1*2.*VK3/(VX3*FRCPX3)
  116.           TCON(8,KON)=DTIM1*2./FRCPX3
  117.           ENDIF
  118.   150 CONTINUE
  119.       RETURN
  120.   500 FORMAT(1H0,10(1H$),'STABILITY CORRECTION IMPLEMENTED FOR',
  121.      +' CONSTRUCT NO. ',I4/11X,'WALL TEMP CALC REPEATED ',I4,
  122.      +' TIMES IN IWALLS.  FOURIER MOD= ',F8.1,' CF MIN OF ',F8.1)
  123.       END
  124.