mirror of
https://github.com/retro-software/B5500-software.git
synced 2026-03-02 17:44:40 +00:00
2c72f7fd1d
1. Commit library tape images, directories, and extracted text files. 2. Commit additional utilities under Unisys-Emode-Tools.
172 lines
11 KiB
Plaintext
172 lines
11 KiB
Plaintext
INTEGER ARRAY HVC2[0:2] ;
|
|
ARRAY TR1S3, TR1S4[0:2], HVC5, HVC9[0:49], HVC3, HVC4, HVC7,
|
|
HVC8[0:2, 0:6] ;
|
|
ARRAY TC, PC[0:N], HVC1[0:N, 0:4] ;
|
|
FILL HVC2[*] WITH 0, 6, 6 ;
|
|
FILL HVC3[1,*] WITH 0.0, -1.1671875, -1.4894278,
|
|
-1.1559073, -1.3616489, -1.2883117, -1.3995702 ;
|
|
FILL HVC3[2,*] WITH 0.0, -1.2935283, -2.9748991,
|
|
-2.4260566, -1.8061761, -1.5492552, -1.4278796 ;
|
|
FILL HVC4[1,*] WITH 0.0, 0.0, -0.10587647,
|
|
-0.15860289, -0.15248362, -0.11896855, -0.12173253 ;
|
|
FILL HVC4[2,*] WITH 0.0, 0.0, -1.7678550,
|
|
-1.8686137, -1.2165188, -0.67716468, -0.61318349 ;
|
|
FILL HVC5[*] WITH 0.0, -0.68264275, -1.0872225,
|
|
0.047697762, -0.024306088, -0.0090858597, -0.001729015,
|
|
-0.0027521429, -2.2596067, 0.11921812, -0.046204775,
|
|
0.014351991, -0.023163497, 0.025471511, 0.013133272,
|
|
0.78097570, -0.17094139, 0.098209577, -0.044902667,
|
|
0.015565173, -0.010983720, 0.013622787, 0.28893313,
|
|
-0.25370399, 0.0052142570, 0.071616235, 0.046127269,
|
|
-0.11245569, -0.067557734, -2.5307361, -0.039864924,
|
|
-0.12607210, 0.13765330, -0.039550505, 0.0075315843,
|
|
-0.21354763, 2.2361128, -0.36080301, 0.49691013,
|
|
-0.41863117, -0.076349532, 0.43597890, -0.40507180,
|
|
0.26338300, 0.0083735620, -0.016716827, -0.016240225,
|
|
0.10507060, -0.12106818, 0.025741676 ;
|
|
FILL HVC7[1,*] WITH 0.0, -1.2729166, -1.6239677,
|
|
-1.5317194, -1.4788939, -1.4345722, -1.5801600 ;
|
|
FILL HVC7[2,*] WITH 0.0, -6.5718750, -12.149681,
|
|
-11.497100, -11.065928, -10.278189, -10.667091 ;
|
|
FILL HVC8[1,*] WITH 0.0, 0.0, -0.090954100,
|
|
-0.056326597, -0.075512113, -0.071434341, -0.063708461;
|
|
FILL HVC8[2,*] WITH 0.0, 0.0, -41.217490,
|
|
-22.312134, -21.043415, -20.693912, -26.268748 ;
|
|
FILL HVC9[*] WITH 0.0, 39.745103, 6.3744973,
|
|
-0.11410040, -0.0036340168, -8.2198154@-5, -2.1580705@-5,
|
|
2.9001251@-6, -44.671381, -6.7048192, 0.073453717,
|
|
-0.0066641733, 0.0017309322, -9.824585@-5, 1.343254@-5,
|
|
41.351011, 4.5915409, -0.029361153, 0.0018827038,
|
|
-0.0040032069, 7.6001770@-4, -9.750253@-5, -27.356096,
|
|
0.97045128, -0.28080321, 0.032749281, 0.0074142633,
|
|
-0.0017002214, 4.3344771@-4, 39.904305, -4.4069227,
|
|
0.52546214, -0.092967052, 0.0057307519, 0.0034063759,
|
|
-1.2990048@-4, 46.719260, 8.5852725, -0.62558932,
|
|
-0.011829404, -0.041450650, 0.014043180, -9.8770289@-4,
|
|
-3.7232051, -0.30920320, -0.0049008792, 0.0011831959,
|
|
9.1663621@-4, -1.8376908@-4, 3.8674416@-5 ;
|
|
PROCEDURE HVPCOR ( N, T, P, Z, CALCDHV, HV, PHVT, PHVP, PHVY ) ; HVPC 1
|
|
|
|
COMMENT THIS PROCEDURE COMPUTES THE ENTHALPY OF A GASEOUS MIXTURE HVPC 2
|
|
AND, WITH SEPARATE OPTIONS, ITS PARTIAL DERIVATIVES WITH HVPC 3
|
|
RESPECT TO TEMPERATURE, PRESSURE AND MOL FRACTIONS . HVPC 4
|
|
|
|
MARK N. LAYMAN HVPC 5
|
|
PROFESSIONAL SERVICES GROUP, BURROUGHS CORPORATION. HVPC 6
|
|
|
|
CARD SEQUENCE STARTS WITH HVPC-0001. HVPC 7
|
|
FIRST RELEASE 3-20-63. HVPC 8
|
|
|
|
THE INPUT PARAMETERS ARE HVPC 9
|
|
N - NUMBER OF COMPONENTS -INTEGER HVPC 10
|
|
T - TEMPERATURE, F - REAL HVPC 11
|
|
P - PRESSURE, PSIA - REAL HVPC 12
|
|
Z - COMPONENT MOLS OR MOL FRACTIONS - REAL ARRAY [0:N] HVPC 13
|
|
CALCDHV - BOOLEAN ARRAY [0:3]. IF CALCDHV[1], [2] AND [3] HVPC 14
|
|
ARE TRUE, THEN PHVT, PHVP AND PHVY WILL BE HVPC 15
|
|
COMPUTED, RESPECTIVELY . HVPC 16
|
|
IN ADDITION, ARRAYS TC, PC, HVC1 - 5 AND HVC7 - 9 ARE HVPC 17
|
|
NON-LOCAL TO THIS PROCEDURE . HVPC 18
|
|
|
|
|
|
THE OUTPUT PARAMETERS ARE HVPC 19
|
|
HV - VAPOR ENTHALPY, BTU (Z = MOLS) OR BTU/MOL (Z = MOL HVPC 20
|
|
FRACTIONS) - REAL HVPC 21
|
|
PHVT - VAPOR HEAT CAPACITY, BTU/DEG.F (Z = MOLS) OR HVPC 22
|
|
BTU/MOL/DEG.F (Z = MOL FRACTIONS) - REAL HVPC 23
|
|
PHVP - PARTIAL DERIVATIVE OF VAPOR ENTHALPY WITH RESPECT HVPC 24
|
|
TO PRESSURE, BTU/PSIA (Z = MOLS) OR BTU/MOL/PSIA HVPC 25
|
|
(Z = MOL FRACTIONS) - REAL HVPC 26
|
|
PHVY - PARTIAL DERIVATIVE OF VAPOR ENTHALPY WITH RESPECT HVPC 27
|
|
TO MOL FRACTION, BTU/MOL - REAL ARRAY[0:N] ; HVPC 28
|
|
VALUE N, T, P ; HVPC 29
|
|
|
|
INTEGER N ; HVPC 30
|
|
|
|
REAL T, P, HV, PHVT, PHVP ; HVPC 31
|
|
|
|
BOOLEAN ARRAY CALCDHV[0] ; HVPC 32
|
|
|
|
ARRAY Z[0], PHVY[0] ; HVPC 33
|
|
|
|
BEGIN HVPC 34
|
|
|
|
INTEGER I ; HVPC 35
|
|
|
|
BOOLEAN CALCPHVT, CALCPHVP, CALCPHVY ; HVPC 36
|
|
|
|
REAL TPC, PPC, TPR, PPR, Y, SUMZ, SUMZ1, ZI, PHI, TRS1, TRS2, HVPC 37
|
|
B, CT, DT2, HO, THO, PHOT, DHT, T2 ; HVPC 39
|
|
ARRAY DATA, DPHI[0:2] ; HVPC 40
|
|
TPC ~ PPC ~ THO ~ PHOT ~ SUMZ ~ DPHI[0] ~ 0.0 ; HVPC 41
|
|
FOR I ~ 1 STEP 1 UNTIL N DO SUMZ ~ Z[I] + SUMZ ; HVPC 42
|
|
SUMZ1 ~ 1.0 / SUMZ ; CALCPHVY ~ CALCDHV[3]; HVPC 43
|
|
CALCPHVT ~ CALCDHV[1] ; CALCPHVP ~ CALCDHV[2]; HVPC 44
|
|
FOR I ~ 1 STEP 1 UNTIL N DO HVPC 45
|
|
BEGIN HVPC 46
|
|
Y ~ Z[I] | SUMZ1 ; HVPC 47
|
|
TPC ~ TC[I] | Y + TPC ; PPC ~ PC[I] | Y + PPC HVPC 48
|
|
END ; HVPC 49
|
|
DATA[1] ~ TPR ~ (T + 459.6) / TPC ; PPR ~ P / PPC ;HVPC 50
|
|
IF CALCPHVY THEN CALCDHV[1] ~ CALCDHV[2] ~ TRUE ; HVPC 51
|
|
|
|
COMMENT TEST PSEUDO-REDUCED PRESSURE TO CHOOSE PROPER SURFACE FIT;HVPC 52
|
|
|
|
IF PPR { 1.0 THEN HVPC 53
|
|
BEGIN HVPC 54
|
|
DATA[2] ~ - LN(PPR) ; HVPC 55
|
|
|
|
COMMENT COMPUTE THE PRESSURE CORRECTION AND ITS PARTIAL HVPC 56
|
|
DERIVATIVES USING THE LOW RANGE FIT ; HVPC 57
|
|
|
|
PHIPART (2,DATA,HVC3,HVC4,HVC5,HVC2,CALCDHV, PHI,DPHI) ; HVPC 58
|
|
PHI ~ EXP(PHI) ; DHT ~ - TPC | PHI ; HVPC 59
|
|
IF CALCPHVT THEN PHVT ~ - DPHI[1] | PHI ; HVPC 60
|
|
IF CALCPHVP THEN PHVP ~ - DPHI[2] | DHT / P ; HVPC 61
|
|
IF CALCPHVY THEN HVPC 62
|
|
BEGIN HVPC 63
|
|
TRS1 ~ (DPHI[1] | TPR - 1.0) | PHI ; HVPC 64
|
|
TRS2 ~ - DPHI[2] | TPC | PHI / PPC HVPC 65
|
|
END HVPC 66
|
|
END HVPC 67
|
|
ELSE HVPC 68
|
|
BEGIN HVPC 69
|
|
DATA[2] ~ PPR ; HVPC 70
|
|
|
|
COMMENT COMPUTE THE PRESSURE CORRECTION AND ITS PARTIAL HVPC 71
|
|
DERIVATIVES USING THE HIGH RANGE FIT ; HVPC 72
|
|
|
|
PHIPART (2,DATA,HVC7,HVC8,HVC9,HVC2,CALCDHV, PHI,DPHI) ; HVPC 73
|
|
DHT ~ - TPC | PHI / PPR ; HVPC 74
|
|
IF CALCPHVT THEN PHVT ~ - DPHI[1] / PPR ; HVPC 75
|
|
IF CALCPHVP THEN PHVP ~ - (DPHI[2] | TPC + DHT) / P ; HVPC 76
|
|
IF CALCPHVY THEN HVPC 77
|
|
BEGIN HVPC 78
|
|
TRS1 ~ (DPHI[1] | TPR - PHI) / PPR ; HVPC 79
|
|
TRS2 ~ (DPHI[2] / PPC - PHI / P) | TPC HVPC 80
|
|
END HVPC 81
|
|
END ; HVPC 82
|
|
T2 ~ T | T ; HVPC 83
|
|
FOR I ~ 1 STEP 1 UNTIL N DO HVPC 84
|
|
BEGIN HVPC 85
|
|
ZI ~ Z[I] ; B ~ HVC1[I,2] ; HVPC 86
|
|
CT ~ HVC1[I,3] | T ; DT2 ~ HVC1[I,4] | T2 ; HVPC 87
|
|
|
|
COMMENT COMPUTE IDEAL GAS STATE ENTHALPY AND HEAT CAPACITY ; HVPC 88
|
|
|
|
HO ~ HVC1[I,1] + (B + CT + DT2) | T ; HVPC 89
|
|
THO ~ HO | ZI + THO ; HVPC 90
|
|
IF CALCPHVT THEN HVPC 91
|
|
PHOT ~ (DT2 | 3.0 + CT | 2.0 + B) | ZI + PHOT ; HVPC 92
|
|
IF CALCPHVY THEN HVPC 93
|
|
PHVY[I] ~ TC[I] | TRS1 + PC[I] | TRS2 + HO HVPC 94
|
|
END ; HVPC 95
|
|
|
|
COMMENT COMPUTE THE TOTAL PRESSURE CORRECTED ENTHALPY AND ITS HVPC 96
|
|
PARTIAL DERIVATIVES ; HVPC 97
|
|
|
|
HV ~ DHT | SUMZ + THO ; HVPC 98
|
|
IF CALCPHVT THEN PHVT ~ PHVT | SUMZ + PHOT ; HVPC 99
|
|
IF CALCPHVP THEN PHVP ~ PHVP | SUMZ ; HVPC 100
|
|
END HVPCOR ; HVPC 101
|