Coker A.K. Fortran Programs for Chemical Process Design, Analysis, and Simulation
Подождите немного. Документ загружается.
22
25
C
C
33
C
C
ACALI=WG*SQRT(TI*ZI)/VAL2
DO 22 I=l, 14
IF (Y(I) .ST. ACALI) THEN
GO TO 25
ELSE
ENDIF
CONTINUE
ADESI=Y(I)
CALCULATE THE MAXIMUM FLOWRATE FROM THE MANUFACTURER'S
ORIFICE AREA
CALL MANUF(ADESI)
WGMAX=ADESI*VAL2/SQRT(TI*ZI)
CALCULATE THE REACTIVE FORCE
CALL REACTF(WGMAX,NI,TI,MWI)
T1 = TI-460.0
RETUR~
END
THIS PROGRAM SIZES THE RELIEF-VALVE FOR STEAM CONDITION
SUBROUTINE STEAM (WS, SP2, OP2, KSH, T2, MW2, N2)
CHARACTER COUNT(14),REGIME(14), REG (14)
CHARACTER*I4 COUNTI(14),COUNLL(14), ACOUNT(14)
COMMON/NAME/I,Y,IJK
COMMON/DATA2/WSMAX,ADES2,ACAL2,C2
COMMON/DATA4/COUNT,REGIME,COUNTI,COUNLL,REG, ACOUNT
DIMENSION Y(14)
REAL KSH,OP2,MW2,N2,KN,KD
KD = THE EFFECTIVE COEFFICIENT OF DISCHARGE (KD=0.975)
KN = CORRECTION FACTOR FOR NAPIER EQUATION (KN=I)
CALCULATE THE GAS COEFFICIENT (C2) AND THE ORIFICE AREA
KN = 1.0
KD = 0.975
T2 =T2+460,0
P2=SP2+(SP2*OP2/100.O)+I4.7
VAL3=N2*(2.0/(N2+I.0))**((N2+I.0)/(N2-1.0))
C2=520.O*SQRT(VAL3)
ACAL2=WS/(51.5*P2*KSH*KD*KN)
DO 33 I=i,14
IF (X(I) .ST. ACAL2) THEN
GO TO 25
ELSE
ENDIF
CONTINUE
CALCULATE THE MAXIMUM FLOWRATE FROM THE MANUFACTURER'S ORIFICE
AREA.
407
25
44
25
ADES2=Y(I)
CALL MANUF(ADES2)
WSMAX=51.5*P2*KSH*KN*KD*ADES2
CALCULATE THE REACTIVE FORCE.
CALL REACTF(WSMAX,N2,T2,MW2)
T2 = T2-460.0
RETURN
END
****************************************************************
THIS PROGRAM SIZES THE RELIEF VALVE FOR LIQUID CONDITION.
SIZING FOR LIQUID RELIEF VALVES REQUIRING LIQUID CAPACITY
CERTIFICATION.
****************************************************************
SUBROUTINE LIQUID (GPM, SP3, PB, SPL, VISL, KV, KW, KD)
DIMENSION Y(14)
CHARACTER COUNT(14),REGIME(14),REG(14),REGI(14)
CHARACTER*I4 COUNTI(14),ACO~T(14),COUNLL(14),ACONT(14)
REAL KD, KW, K~, KVl
COMMON/DATA3/GPIMAX,GP2MAX,ADES3,ACAL3,ACAL31,ADES31,RE
COMMON/NAME/I,Y, IJK
COMMON/DATA4/COUNT,REGIME,COUNTI,COUNLL, REG, ACOUNT
COMMON/DATAS/KVl
COMMON/DATA/REGI,ACONT
CALCITE THE ORIFICE AREA FOR LIQUID CAPACITY CERTIFICATION.
SP3 = SP3+(0.1*SP3)
VAL3 = (SP3-PB)
VAL4=38.0*KD*KW*KV*SQRT(VAL3)
ACAL3=GPM*SQRT(SPL)/VAL4
DO 44 I =i, 14
IF (Y(I) .GT. ACAL3) THEN
GO TO 25
ELSE
ENDIF
CONTINUE
ADES3=Y(I)
IJK=I
CALCULATE THE MANUFACTURER'S ORIFICE AREA AND THE MAXIMUM
LIQUID FLOWRATE.
CALL MANUF(ADES3 )
GPIMAX=( ADES3*VAL4 )/SQRT(SPL)
CALL OUT3 (GPM, SP3, PB, SPL, VISL, KV, KW, KD)
408
C CHECK WHETHER THE VISCOSITY IS EQUAL TO OR GREATER THAN 10.0cP
30
66
77
90
I00
IF (VISL .GE. I0.) THEN
GO TO 3O
ELSE
ENDIF
GO TO 90
CALL AMANU(ADES3)
REGI(I)=REG(I)
ACONT(I)=ACOUNT(I)
GPIMAX=(ADES3*VAL4)/SQRT(SPL)
CALL OUT31(GPM,SP3,PB,SPL,VISL,KV,k~,KP)
RE=(GPM*2800.0*SPL)/(VISL*SQRT(ADES3))
VAL6=(I.2502*ALOG(RE))-(0.17510*(ALOG(RE))**2)
VAL7=(0.01087*(ALOG(RE))**3)-(0.00025*(ALOG(RE))**4)
KVI=-2.38878+VAL6+VAL7
VAL8=38.0*KD*KW*KVI*SQRT(VAL3)
ACAL31=GPM*SQRT(SPL)/VAL8
DO 66 I=i,14
IF (Y(I) .ST. ACAL31) THEN
GO TO 77
ELSE
ENDIF
CONTINUE
ADES31=Y(I)
CALL MANUFL(ADES3,ADES31)
GP2MAX=(ADES31*VAL8)/SQRT(SPL)
RETURN
END
*****************************************************************
THIS PROGRAM PRINTS THE RESULTS OF VAPOR RELIEF SIZING ONTO A
PRINTER
*****************************************************************
SUBROUTINE OUT1 (WG, SPI, OPI, TI, Zl, MWI, Nl)
CHARACTER COUNT( 14 ) ,REGIME( 14 ) ,REG( 14 ) ,REGI ( 14 )
CHARACTER*I4 COUNT1 ( 14 ) ,ACOUNT( 14 ) ,COUNLL( 14 ) ,ACONT(14 )
CHARACTER*I0 CHOICE
COMMON/NAME/I, Y, IJK
COMMON/DATAI/WGMAX,ADESI,ACALI,CI
COMMON/DATA4/COUNT,REGIME,COUNTI,COUNLL,
REG, ACOUNT
COMMON/DATA5/FORCE
COMMON CHOICE
REAL OPI,MWI,N1
DIMENSION Y(14)
WRITE(I, I00)
FORMAT(20X,'VAPOR RELIEF VALVE SIZING')
WRITE(l, ii0)
409
105
ii0
120
130
140
150
160
170
180
190
200
210
220
230
240
250
260
C
20
WRITE (i, 105)CHOICE
FORMAT (3X,'BASIS OF SELECTION:', T65, A)
FORMAT(IH ,76(IH*))
WRITE(l, 120)WG
FORMAT(3X,'REQUIRED VAPOR RATE, ib/hr.:',T65, FI2.0)
WRITE(I, 130) SPI
FORMAT(3X,'VALVE SET PRESSURE, psig:',T65, F9.3)
WRITE(l, 140) OPI
FORMAT(3X,'PERCENTAGE IN OVERPRESSURE, (%):',T65, F9.3)
WRITE(I, 150)TI
FORMAT(3X,'UPSTRE~M FLUID TEMPERATURE, DEG. F:',T65, F9.3)
WRITE(l, 160)Zl
FORMAT(3X,'COMPRESSIBILITY FACTOR:',T65, F9.4)
WRITE(I, 170) HWI
FORMAT(3X,'MOLECULAR WEIGHT OF VAPOR:', T65, F9.3)
WRITE(l, 180)NI
FORMAT(3X,'RATIO OF SPECIFIC HEATS OF FLUID:', T65, F9.3)
WRITE(l, 190) Cl
FORMAT(3X,'GAS CONSTANT DETERMINED BY RATIO OF Cp/Cv:',T65,F9.3)
WRITE(I, 200) ACALI
FORMAT(3X,'CALCULATED VALVE ORIFICE AREA, in^2:',T65, F9.3)
WRITE(I, 210) ADESI
FORMAT(3X,'NEAREST STANDARD ORIFICE AREA, in^2:',T65, F9.3)
WRITE(l, 220) REGIME(I)
FORMAT(3X, 'NEAREST STANDARD ORIFICE SIZE:', T65, A)
WRITE(l, 230) COUNT1 (I)
FORMAT(3X, 'PREFERRED VALVE BODY SIZE: (IN-OUT) inches: ' ,T65,AI4)
WRITE(l, 240) WGMAX
FORMAT(3X, 'MAXIMUM VAPOR FLOWRATE WITH STD. ORIFICE, Ib/hr. :',
T65,FI2.0)
WRITE(I, 250) FORCE
FORMAT(3X,'REACTION FORCE, Ib:',T65, F9.0)
WRITE
( i, 260)
FORMAT ( 76 ( IH- ) )
FORM FEED THE PRINTING PAPER TO THE TOP OF THE NEXT PAGE.
WRITE (i, *)CHAR(12)
RETURN
STOP
END
THIS PROGRAM PRINTS THE RESULTS OF LIQUID RELIEF CONDITION
410
270
280
285
290
300
310
320
330
340
350
360
370
380
390
400
410
ONTO A PRINTER
****************************************************************
SUBROUTINE OUT3 (GPM, SP3, PB, SPL, VISL, KV, KW, KD)
DIMENSION Y(14)
CHARACTER COUNT(14),REGIME(14),REG(14),REGI(14)
CHARACTER*I4 COUNTI(14),ACOUNT(!4),COUNLL(i4),ACONT(14)
CHARACTER*!0 CHOICE
CO~.~ON/N~4E/I, Y, IJK
COMMON/DATA/REG
1, ACONT
CO~HION/DATA]/GPIMAX,
GP2MAX, ADES3, ACAL3, ACAL31, ADES31, RE
COMMON/DATA4/COUNT, REGIME, COUNT1, COUNLL, REG, ACOUNT
COMMON/DATA8/KVI
CO~ON CI{OICE
REAL KV,KW,KD,KVI
~ITE(I, 270)
FOR}IAT(20X,'RELIEF VALVE SIZING FOR LIQUID')
WRITE(l, 280)
FORMAT(1H ,76(IH*))
WRITE (i, 285)CHOICE
FORMAT (3X,'BASIS OF SELECTION:', T65, A)
~{RITE(I, 290) GPM
FORMAT(3X,'LIQUID FLOWRATE, US gpm:',T65, FI2.0)
W~ITE(I, 300) SP3
FORMAT(3X,'VALVE SET PRESSURE, psig:',T65, F9.3)
I~RITE(I,310) PB
FORMAT(3X,'BACK PRESSURE ON VALVE, psig:',T55, F9.3)
WRITE(I, 320) SPL
FORMAT(3X,'FLUID SPECIFIC GRAVITY:', T65, F9.3)
WRITE(I,330) VISL
FO~4AT(3X,'FLUID VISCOSITY, cP:', T65, F9.4)
~RITE(I, 340) KV
FORMAT(3X,'VISCOSITY CORRECTION FACTOR, Kv:',T65, F9.4)
WRITE(l, 350) KW
FORMAT(3X,'BACKPRESSURE CORRECTION FACTOR, Kw:',T65, F9.4)
WRITE(l, 360) KD
FORMAT(3X,'EFFECTIVE COEFFICIENT OF DISCHARGE, Kd:', T65, F9.3)
WRITE(l, 370) ACAL3
FORMAT(3X,'CALCULATED VALVE ORIFICE AREA, in^2:',T65,F9.3)
W~ITE(I, 380) ADES3
FORMAT(3X,'NEAREST STANDARD ORIFICE AREA, in^2: ' , T65,F9.3)
WRITE(l, 390) REGIME (I)
FORMAT(3X,'NEAREST STANDARD ORIFICE SIZE:', T65, A)
Z~ITE(I, 400) COUNLL(I)
FORMAT(3X,'PREFERRED VALVE BODY SIZE, (IN-OUT) inches:',T65,A14)
WRITE(l, 410) GPIMAX
FOR~T(3X,'~XI}KIM LIQUID FLOW~-RATE WITH STD. ORIFICE, gpm:',
T65,FI2.0)
411
420
430
440
450
460
470
480
490
500
510
520
530
540
550
FORM FEED THE PRINTING PAPER TO TOP OF THE NEXT PAGE.
RETURN
END
**********************************************************
SUBROUTINE OUT31 (GPM, SP3, PB, SPL, VISL, KV, KW, KD)
CHARACTER COUNT(14),REGIME(14),REG(14),REGI(!4)
CHARACTER*I4 COUNTI(14),ACOUNT(14),COUNLL(14),ACONT(14)
CHARACTER*I0 CHOICE
DIMENSION Y(14)
COMMON/NAME/I,Y,IJK
COMMON/DATA/REGI,ACONT
COMMON/DATA3/GPIMAX,GP2MAX,ADES3,ACAL3,ACAL31,ADES31,RE
COMMON/DATA4/COUNT,REGIME,COUNTI,COUNLL, REG, ACOUNT
COMMON/DATA8/KVI
COMMON CHOICE
REAL KV,KW,KD,KVI
WRITE(l, 420)
FORMAT(20X,'RELIEF VALVE SIZING FOR LIQUID')
WRITE(l, 430)
FORMAT(IH ,76(IH*))
WRITE (I, 440)CHOICE
FORMAT (3X,'BASIS OF SELECTION:', T65, A)
WRITE(l, 450) SPS
FORMAT(3X,'LIQUID FLOW RATE, U.S gpm:',T65,F12.0)
WRITE(l, 460) SP3
FORMAT(3X,'VALVE SET PRESSURE, psig:',T65,F9.3)
WRITE(I,470) SPL
FO~AT(3X,'FLUID SPECIFIC GRAVITY:',T65, F9.3)
WRITE(l, 480) KV
FORMAT(3X,'VISCOSITY CORRECTION FACTOR Kv=I.0:',T65,F9.3)
WRITE(l, 490) KW
FORMAT(3X,'BACK PRESSURE CORRECTION FACTOR Kw=I.0:',T65,F9.3)
WRITE(l, 500)KD
FORMAT(3X,'COEFFICIENT OF DISCHARGE Kd=0.65:',T65, F9.3)
WRITE(l, 510) ACAL3
FORMAT(3X,'CALCULATED VALVE ORIFICE AREA, in^2:',T65,F9.3)
WRITE(l, 520) ADES3
FORMAT(3X,'NEAREST STANDARD ORIFICE AREA, in^2: ' , T65,F9.3)
WRITE(l, 530) REGI(I)
FORMAT(BX,'NEAREST STANDARD ORIFICE SIZE:',T65, A)
WRITE(l, 540) ACONT(I)
FORMAT(3X,'PREFERRED VALVE BODY SIZE, (IN-OUT) inches:',
T65,AI4)
WRITE(l, 550) GPIMAX
FORMAT(3X,'MAXIMUM LIQUID FLOWRATE WITH STD. ORIFICE, gpm:',
T65,FI2.0)
WRITE(l, 560)
412
560
C
570
580
590
600
610
620
630
640
C
FOR/4AT ( 76 ( IH- ) )
FORM FEED THE PRINTING PAPER TO THE TOP OF THE NEXT PAGE.
WRITE ( i, * ]CHAR( 12 )
RETURN
END
********************************************************
THIS PROGRAM PRINTS THE RESULTS OF LIQUID RELIEF VALVE
SIZING WITH VISCOSITY GREATER THAN 10cP.
SUBROUTINE OUT32 (VISL)
DIMENSION Y(14)
CHARACTER COUNT(!4),REGIME(14),REG(14),REGI(14)
CHARACTER*I4 COUNTI(14),ACOUNT(!4),COUNLL(14),ACONT(!4)
CO~MON/NAME/I,Y, IJK
COMMON/DATA3/GPIMAX,GP2MAX,ADES3,ACAL3,ACAL31,ADES31,RE
COMMON/DATA4/COUNT,REGIME,COUNTI,COUNLL, REG, ACOU~[T
COMMON/DATA8/KVI
REAL KVI
WRITE(I, 570) KVI
FORMAT(3X, 'VISCOSITY CORRECTION FACTOR: ', T65, F9.3)
WRITE(I, 580) RE
FORMAT(3X, 'REYNOLDS NUMBER: ', T65, FI2.0)
WRITE(l, 590) ACAL31
FORMAT(3X,'CALCULATED VALVE ORIFICE AREA, in^2:',T65,F9.3)
WRITE(l, 600) ADES31
FORMAT(3X,'NEAREST STANDARD ORIFICE ~EA, in^2: ' , T65,F9.3)
WRITE(l, 610) REGIME (I)
FORMAT(3X,'NEAREST STANDARD ORIFICE SIZE:', T65, A)
WRITE(I, 620) COUNLL (I)
FORMAT(3X, 'PREFERRED VALVE BODY SIZE, (IN-OUT) inches: ' ,T65,AI4)
WRITE(l, 630) GP2MAX
FORMAT(3X,'}~XIMUM LIQUID FLOWRATE WITH STD. ORIFICE, gpm:',
T65,FI2.0)
WRITE
(I, 640)
FORMAT ( 76 ( IH- ) )
FORM FEED THE PRINTING PAPER TO THE TOP OF THE NEXT PAGE.
WRITE ( i, * )CHAR( 12 ]
RETURN
END
413
650
660
665
670
680
690
700
710
720
730
740
750
760
770
780
790
800
810
C
THIS PROGRAM PRINTS TI{E RESULTS OF STE~! RELIEF SIZING
ONTO A PRI~$TER
***************************************************************
SUBROUTINE OUT2 (WS, SP2, OP2, KSH, T2, MW2, N2)
CHARACTER COUNT(14),REGIME(14), REG(14)
CHARACTE}{*!4 COUNT1 ( 14 ), COUNLL( 14 ), ACOUNT( 14 )
CHARACTER* i0 CHOICE
COMbiON/NAME/I,
Y, IJK
COHHON/DATA2/WSMAX, ADES2, ACAL2, C2
COMMON/DATA4/COUNT,REGIME,COUNTI,COUNLL,REG, ACOUNT
CObN.[ON/DATAS/FORCE
COI.H~ON CHOICE
DIMENSION Y(!4)
REAL KSI[,MW2,N2
WRITE( l, 650)
FOP~T(20X,'RELIEF VALVE SIZING FOR STEAM CONDITION')
WRITE(l, 660)
FORMAT(IN ,76(IH*))
WRITE (i, 665)CHOICE
FORMAT (3X, 'BASIS OF SELECTION:', T65, A)
WRITE(l, 670) WS
FORMAT(3X,'REQUIRED STE~]{ RATE, ib/hr.:',T65,F12.0)
WRITE(l, 680)SP2
FORMAT(3X,'VALVE SET PRESSURE, psig:',T65,F9.3)
WRITE(l, 690) OP2
FORMAT(3X,'PERCENTAGE OF OVERPRESSURE, (%)'',T65,F9.9)
WRITE(I, 700) KSH
FORMAT(3X,'FACTOR FOR SUPERHEATED KSH:0.88/SATL<~t%TED KSH=I:',
T65,F9.3)
WRITE(I, 710) T2
FORMAT(3X,'STEAM TEMPERATURE, DEG. F:',T65,F9.3)
WRITE(l, 720) ~.~42
FO~4AT(3X,'MOLECULAR WEIGHT:', T65, F9.3)
9~ITE(I, 730) N2
FO~-~T(3X,'RATIO OF SPECIFIC HEATS OF STE~{:', T65, F9.3)
WRITE(l, 740) C2
FORMAT(3X,'GAS CONSTANT DETERMINED BY ~hTIO OF Cp/Cv:',T65,F9.3)
WRITE(I, 750) ACAL2
FORMAT(3X,'CALCULATED VALVE ORIFICE AREA, in^2: ' , T65, F9.3)
WRITE(l, 760) ADES2
FOP~AT(3X,'NEAREST STANDARD AR~%, in^2: ' , T65, F9.3)
WRITE(l, 770) REGII4E (I)
FORMAT(3X,'NEAREST STANDARD ORIFICE SIZE:', T65, A)
WRITE(l, 780) COUNTI(I)
FORMAT(3X,'PREFERRED VALVE BODY SIZE, (IN-OUT) inches:',T65,A14)
WRITE(l, 790) WSI,~X
FORMAT(3X,'MAXIMUM STEAM FLOWRATE WITH STD. ORIFICE, ib/hr.:',
T65,FI2.0)
WRITE(I, 800)FORCE
FORMAT(3X,'REACTION FORCE, Ib:',T65,FI2.0)
WRITE (I, 810)
FORY~T (76(1H-))
FORM FEED THE PRINTING PAPER TO THE TOP OF THE NEXT PAGE.
WRITE (i, *)CHAR(12)
RETURN
END
414
44
C
30
****************************************************************
THIS PROGRAM SIZES THE RELIEF VALVE FOR AIR-CONDITION
****************************************************************
SUBROUTINE AIR (VA, T4, SP4, OP4)
REAL OP4,N,MW,KB,KD
CHARACTER COUNT(14),REGIME(14), REG(14)
CHARACTER*I4 COUNT1(14), COUNLL(14),ACOUNT(14)
COMMON/NAME/I, Y, IJK
COMMON/DATA4/COUNT,REGIME,COUNTI,COUNLL, REG, ACOUNT
COMMON/DATA5/FORCE
COMMON/DATA6/ADES4,VAMAX,ACAL4,C4
DIMENSION Y(14)
CALCULATE THE ORIFICE AREA AND THE GAS COEFFICIENT.
KB =I.0
KD=0.953
T = T4+460.0
MW=28.97
N=I.40
P4=SP4+(SP4*OP4/100.0)+I4.7
RHO= (MW*P4)/(I0.72"T)
ACAL4=(VA*SQRT(T))/(41S.0*KB*KD*P4)
VAL4=N*(2.0/(N+I.0))**((N+I.0)/(N-I.0))
C4=520.O*SQRT(VAL4)
DO 44 I=I,14
IF (Y (I) .ST. ACAL4) THEN
GO TO 30
ELSE
ENDIF
CONTINUE
CALCULATE THE MANUFACTURER'S AREA AND MAXIMUM AIR RATE.
ADES4=Y(I)
CALL MANUF(ADES4)
VAMAX=(418.0*ADES4*P4*KB*KD)/SQRT(T)
WAMAX=(60.0*VAMAX*RHO)
CALCULATE THE REACTIVE FORCE.
FORCE = WAMAX/15.9
RETURN
END
THIS PROGRAM PRINTS THE RESULTS OF AIR RELIEF VALVE SIZING ONTO
A PRINTER
****************************************************************
SUBROUTINE OUT4 (VA, T4, SP4, OP4)
CHARACTER*I4 COUNT1(14), COUNLL(14), ACOUNT(14)
CHARACTER*I0 CHOICE
CHARACTER COUNT(14),REGIME(14), REG (147
COMMON/NAME/I,Y,IJK
COMMON/DATA4/COUNT,REGIME,COUNTI,COUNLL,
REG, ACOUNT
COMMON/DATA5/FORCE
COMMON/DATA6/ADES4,VAMAX,ACAL4,C4
COMMON CHOICE
415
820
830
835
84O
85O
860
870
880
890
900
i000
ll00
1200
1300
1400
C
DIMENSION Y(14)
REAL OP4
WRITE(l, 820)
FORMAT(20X,'RELIEF VALVE SIZE FOR AIR CONDITION')
WRITE(I, 830)
FORMAT(IH ,76(IH*))
WRITE (i, 835)CHOICE
FORMAT (3X,'BASIS OF SELECTION:', T65, A)
WRITE(l, 840)VA
FORMAT(3X,'REQUIRED AIR FLOWRATE, S.ft^3/min:',T65,F12.0)
WRITE(l, 850) T4
FORMAT(3X,'AIR TEMPERATURE, DEG. F:',T65,F9.3)
WRITE(l, 860) SP4
FORMAT(3X,'VALVE SET PRESSURE, psig:',T65,F9.3)
WRITE(l, 870) OP4
FORMAT(3X,'PERCENTAGE OF OVERPRESSURE, (%):',T65,Fg.3)
WRITE(l, 880) C4
FORMAT(3X,'GAS CONSTANT DETERMINED BY RATIO OF Cp/Cv:',T65,F9.3)
WRITE(l, 890) ACAL4
FORMAT(3X,'CALCULATED VALVE ORIFICE AREA, in^2:',T65,F9.3)
WRITE(l, 900) ADES4
FORMAT(3X,'NEAREST STANDARD AREA, in^2: ' , T65, F9.3)
WRITE(l, i000) REGIME (I)
FORMAT(3X,'NEAREST STANDARD ORIFICE SIZE:',T65,A)
WRITE(l, ii00) ACOUNT(I)
FORMAT(3X,'PREFERRED VALVE BODY SIZE, (IN-OUT) inches:',T65,A14)
WRITE(l, 1200) VAMAX
FORMAT(3X,'MAXIMUM AIR FLOWRATE WITH STD. ORIFICE, ft^3/min: ' ,
T65,FI2.0)
WRITE(l, 1300) FORCE
FORMAT(3X,'REACTION FORCE, Ib:',T65,FI2.0)
WRITE (i, 1400)
FORMAT (76(iH-))
FORM FEED THE PRINTING PAPER TO THE TOP OF THE NEXT PAGE.
WRITE (i, *)CHAR(12)
RETURN
END
THIS PROGRAM SIZES THE RELIEF VALVE FOR FIRE CONDITION
SUBROUTINE FIRE (F, AREA, LAT, SP5, OP5, T5, Z5, MW5, N5)
REAL LAT,OP5,MW5,N5
CHARACTER COUNT(14),REGIME(14), REG (14)
CHARACTER*I4 COUNTI(14),COUNLL(14), ACOUNT(14)
COMMON/NAME/I,Y,IJK
COMMON/DATAI/WGMAX,ADESI,ACALI,CI
COMMON/DATA4/COUNT,REGIME,COUNTI,COUNLL,REG, ACOUNT
COMMON/DATA5/FORCE
COMMON/DATA7/Q,WF
DIMENSION Y(14)
CALCULATE THE TOTAL HEAT ABSORBED
Q=21000.0*F*(AREA)**0.82
416