Oppdif Revised
Oppdif Revised
131 WRITE (LOUT, *) ' FATAL ERROR, NOT ENOUGH WORK SPACE PROVIDED' 196 C
132 STOP 197 C
133 ENDIF 198 C SET THE POINTERS INTO THE SOLUTION VECTOR
134 C 199 C
135 C 200 NT = 1
136 CALL CKSYMS (C(ICC), LOUT, KSYM, IERR) 201 NG = 2
137 C 202 NF = 3
138 CALL CKCRAY (LINE, KK, KSYM, LOUT, 4, NRAD, NFD, IERR) 203 NH = 4
139 C 204 NYS = 4
140 C 205 NY = 5
141 C JUST FOR TEST 206 C
142 C 207 C CALL FOR CHEMKIN WORK ARRAY LENGTHS
143 WRITE (LOUT,'(1X,A)') 'RADIATING SPECIES:' 208 C
144 209 CALL CKLEN (LINKCK, LOUT, LENICK, LENRCK, LENCCK, IFLAG1)
145 WRITE (LOUT,'(8(1X,A))') 210 CALL MCLEN (LINKMC, LOUT, LENIMC, LENRMC, IFLAG2)
146 211 IF (IFLAG1.GT.0 .OR. IFLAG2.GT.0) STOP
147 1 ( KSYM( NRAD(K) ), K = 1, 4) 212 C
148 213 C real chemkin work space
149 214 NCKW = 1
150 C END TEST 215 C real transport work space
151 216 NMCW = NCKW + LENRCK
152 CALL FLDRIV (LENITW, LENTWP, LIN, LOUT, LREST, 217 NTOT = NMCW + LENRMC
153 1 LSAVE, LRCRVR, KK, II, MM, NATJ, NMAX, LENICK, LENRCK, 218 C
154 2 R(NCKW), R(NMCW), R(NYOX), R(NYFL), R(NWT), R(NFL), 219 C integer chemkin space
155 3 R(NPD), R(NOX), R(NSCH), R(NX), R(NVIS), R(NCON), 220 ICKW = 1
156 4 R(NTGV), R(NXGV), R(ND), R(NDKJ), R(NTDR), R(NYV), 221 C integer transport space
157 5 R(NABV), R(NBLW), R(NBUF), R(NTWP), R(NS), R(NSN), 222 IMCW = ICKW + LENICK
158 6 R(NFF), R(NFN), R(NDS), R(NSSV), R(NA), 223 ITOT = IMCW + LENIMC
159 7 I(ICKW), I(IMCW), I(ITPW), C(IKS), C(ICC), 224 C
160 8 I(IIP), L(LAC), L(LMK), 225 C characcter chemkin space
161 C EQUIL ARRAYS... 226 ICC = 1
162 9 I(IEQW), LENIEQ, R(NEQW), LENREQ, R(NREAC), 227 ICTOT = ICC + LENCCK
163 1 R(NDCS), C(IATOM) ) 228 C
164 229 IF (ITOT.LT.LENIWK .AND. NTOT.LT.LENRWK .AND. ICTOT.LT.LENCWK)
165 C 230 1 THEN
166 RETURN 231 CALL CKINIT (LENICK, LENRCK, LENCCK, LINKCK, LOUT, I, R, C,
167 END 232 1 IFLAG)
168 C 233 IF (IFLAG .GT. 0) STOP
169 C---------------------------------------------------------------------- 234 ILINK = 'CKLINK'
170 C 235 WRITE (LSAVE) ILINK
171 SUBROUTINE POINT (LINKCK, LINKMC, NMAX, LOUT, LSAVE, KK, II, MM, 236 CALL CKSAVE (LOUT, LSAVE, I, R, C)
172 1 NATJ, LENIWK, LENRWK, LENCWK, LENITW, LENTWP, 237 CALL CKINDX (I, R, MM, KK, II, NFIT)
173 2 LENICK, LENRCK, LENIEQ, LENREQ, 238 C
174 3 LTOT, ITOT, NTOT, ICTOT, I, R, C) 239 CALL MCINIT (LINKMC, LOUT, LENIMC, LENRMC, I(IMCW), R(NMCW),
175 C 240 1 IFLAG)
176 C*****precision > double 241 IF (IFLAG .GT. 0) STOP
177 IMPLICIT DOUBLE PRECISION (A-H, O-Z), INTEGER (I-N) 242 ILINK = 'MCLINK'
178 C*****END precision > double 243 WRITE (LSAVE) ILINK
179 C 244 CALL MCSAVE (LOUT, LSAVE, I(IMCW), R(NMCW))
180 C*****precision > single 245 ENDIF
181 C IMPLICIT REAL (A-H, O-Z), INTEGER (I-N) 246 C
182 C*****END precision > single 247 C TWOPOINT WORK SPACE LENGTHS
183 C 248 C
184 DIMENSION I(*), R(*) 249 NATJ = KK+4
185 CHARACTER C(*)*(*), ILINK*16 250 LENITW = 3 * NMAX
186 C 251 LENTWP = (7*NATJ+2) * NMAX
187 COMMON /FLFLFL/ NCKW, NMCW, NYOX, NYFL, NWT, NFL, NPD, NOX, 252 C
188 1 NSCH, NX, NVIS, NCON, NTGV, NXGV, ND, NDKJ, 253 C STANJAN WORK SPACE LENGTHS
189 2 NTDR, NYV, NABV, NBLW, NBUF, NTWP, NS, NSN, 254 C
190 3 NFF, NFN, NDS, NSSV, NA, 255 NCON = 0
191 4 ICKW, IMCW, ITPW, IKS, ICC, IIP, LAC, LMK, 256 CALL EQLEN (MM, KK, NCON, LENIEQ, LENREQ)
192 C ADDED FOR EQUIL... 257 C
193 5 IEQW, IATOM, NREAC, NDCS, NEQW 258 C APPORTION THE FLOATING POINT SPACE
194 C 259 C
195 COMMON /LOCS/ NT, NG, NF, NH, NYS, NY 260 NYOX = NTOT
261 NYFL = NYOX + KK 326 7 IEQWRK, LENIEQ, REQWRK, LENREQ, REAC,
262 NWT = NYFL + KK 327 8 DCS, ATOM )
263 NFL = NWT + KK 328 C
264 NPD = NFL + KK 329 C*****precision > double
265 NOX = NPD + KK 330 IMPLICIT DOUBLE PRECISION (A-H, O-Z), INTEGER (I-N)
266 NSCH = NOX + KK 331 C*****END precision > double
267 NX = NSCH + KK*5 332 C
268 NVIS = NX + NMAX 333 C*****precision > single
269 NCON = NVIS + NMAX 334 C IMPLICIT REAL (A-H, O-Z), INTEGER (I-N)
270 NTGV = NCON + NMAX 335 C*****END precision > single
271 NXGV = NTGV + NMAX 336 C
272 ND = NXGV + NMAX 337 DIMENSION PMAC(7,4), NRAD(4), NNARROW(6)
273 NDKJ = ND + KK*NMAX 338 C
274 NTDR = NDKJ + KK*KK*NMAX 339 COMMON /LOCS/ NT, NG, NF, NH, NYS, NY
275 NYV = NTDR + KK*NMAX 340 COMMON /RAD/ PMAC, NRAD ,NNARROW
276 NABV = NYV + KK*NMAX 341 COMMON/CLAYER/ QL(100),QLS(100,600)
277 NBLW = NABV + NATJ*NMAX 342
278 NBUF = NBLW + NATJ*NMAX 343 C
279 NTWP = NBUF + NATJ*NMAX 344 C
280 NS = NTWP + LENTWP 345 DIMENSION ICKWRK(*), RCKWRK(*), IMCWRK(*), RMCWRK(*), YOXID(*),
281 NSN = NS + NATJ*NMAX 346 1 YFUEL(*), WT(*), FUEL(*), PROD(*), OXID(*), X(*),
282 NFF = NSN + NATJ*NMAX 347 2 VISC(*), COND(*), D(KK,*), TDR(KK,*), YV(KK,*),
283 NFN = NFF + NATJ*NMAX 348 3 DKJ(KK,KK,*), S(NATJ,*), SN(NATJ,*), F(NATJ,*),
284 NDS = NFN + NATJ*NMAX 349 4 XGIVEN(*), TGIVEN(*), SCRCHK(KK,5), KW(6), XMF(KK)
285 NSSV = NDS + NMAX 350 C
286 NA = NSSV + NMAX 351 C EQUIL SPACE
287 C NTOT = NA + (6*NATJ-2) * (NATJ*NMAX) - 1 352 C
288 NEQW = NA + (6*NATJ-2) * (NATJ*NMAX) - 1 353 DIMENSION IEQWRK(LENIEQ), REQWRK(LENREQ), REAC(KK,2), DCS(KK)
289 NREAC = NEQW + LENREQ 354 C
290 NDCS = NREAC + 2*KK 355 C DIMENSION NEWTON SPACE
291 NTOT = NDCS + KK 356 C
292 C 357 DIMENSION ABOVE(NATJ,*), BELOW(NATJ,*), BUFFER(NATJ,*),
293 C APPORTION THE INTEGER SPACE 358 1 TWPWK(*), ITWPWK(*), A(6 * NATJ - 2, *), IP(NATJ, *),
294 C 359 3 FN(NATJ,*), DS(*), SSAVE(*), LEVEL(2)
295 ITPW = ITOT 360 C
296 IIP = ITPW + LENITW 361 CHARACTER KSYM(*)*(*), CCKWRK(*)*(*), ATOM(*)*(*), LINE*80
297 C ITOT = IIP + NATJ*NMAX - 1 362 C
298 IEQW = IIP + NATJ*NMAX - 1 363 LOGICAL LTIME, LTIME2, LENRGY, LMULTI, LTDIF, LUMESH, LRSTRT,
299 ITOT = IEQW + LENIEQ 364 1 LCNTUE, LSEN, LVARMC, RSTCNT, ERROR, FUNCTN, JACOBN,
300 C 365 2 REENTR, SOLVE, STORE, SUCCES, ADAPT, SHOW, SAVE, UPDATE,
301 C APPORTION THE LOGICAL SPACE 366 3 ACTIVE(*), MARK(*), ENERGY, LUSTGV, LREGRD, LRADI,LNARR
302 C 367 C
303 LAC = 1 368 INTEGER CALL, CALLS
304 LMK = LAC + NATJ 369 C
305 LTOT = LMK + NMAX - 1 370 DOUBLE PRECISION XFRACT(KK)
306 C 371
307 C APPORTION THE CHARACTER SPACE 372 CHARACTER*16 ISOLUT, ISENSI, ICHR, ICKLNK, IMCLNK
308 C 373 LOGICAL IERR, LEQUIL
309 IKS = ICTOT 374 DATA ISOLUT/'SOLUTION '/, ISENSI/'SENSITIVITY '/,
310 IATOM = IKS + KK 375 1 ICKLNK/'CKLINK '/, IMCLNK/'MCLINK '/
311 ICTOT = IATOM + MM 376 DATA ADAFLR /1.0E-08/, LCNTUE/.FALSE./, RFAC/1.0/
312 RETURN 377 DATA LINE/'CO2 H2O CH4 CO O2 N2'/
313 END 378 C
314 C 379 C
315 C---------------------------------------------------------------------- 380 CALL CKSYMS (CCKWRK, LOUT, KSYM, IERR)
316 C 381 C
317 SUBROUTINE FLDRIV (LENITW, LENTWP, LIN, LOUT, 382 CALL CKCRAY (LINE, KK, KSYM, LOUT, 6, NNARROW, NFD, IERR)
318 1 LREST, LSAVE, LRCRVR, KK, II, MM, NATJ, NMAX, 383
319 2 LENICK, LENRCK, RCKWRK, RMCWRK, YOXID, YFUEL, 384 C COMPUTE THE UNIT ROUNDOFF, AND THE RELATIVE AND ABSOLUTE
320 3 WT, FUEL, PROD, OXID, SCRCHK, X, VISC, COND, 385 C PERTURBATIONS FOR THE JACOBIAN EVALUATION.
321 4 TGIVEN, XGIVEN, D, DKJ, TDR, YV, ABOVE, BELOW, 386 C
322 5 BUFFER, TWPWK, S, SN, F, FN, DS, SSAVE, A, 387 U = 1.0
323 6 ICKWRK, IMCWRK, ITWPWK, KSYM, CCKWRK, IP, 388 10 CONTINUE
324 7 ACTIVE, MARK, 389 U = U*0.5
325 C EQUIL ARRAYS... 390 COMP = 1.0 + U
911 CALL EQSOL (KK, REQWRK, REAC(1,2), DCS, SOLN(NT,J), PDUM, 976 ETAJ = (J-1)*DETA
912 1 HDUM, VDUM, SDUM, WMDUM, CDUM, CDET) 977 DO 70 I = ISTART, JJOLD
913 C 978 IF (ETAJ .LE. WORK(I)) THEN
914 CALL CKXTY (REAC(1,2), ICKWRK, RCKWRK, SOLN(NY,J)) 979 DEL = (ETAJ-WORK(I-1))/(WORK(I)-WORK(I-1))
915 100 CONTINUE 980 XNEW(J) = XOLD(I-1)+(XOLD(I)-XOLD(I-1))*DEL
916 C 981 GO TO 80
917 RETURN 982 ELSE
918 END 983 ISTART = I
919 C 984 ENDIF
920 C--------------------------------------------------------------------- 985 70 CONTINUE
921 C 986 WRITE (6, *) ' *** VALUE OF ETA NOT FOUND ***'
922 SUBROUTINE REGRID (SNEW, SOLD, NVAR, JJNEW, JJOLD, XNEW, XOLD, 987 80 CONTINUE
923 1 WORK, PR, R, N ) 988 C
924 C 989 C interpolate solution...
925 C*****precision > double 990 C
926 IMPLICIT DOUBLE PRECISION (A-H, O-Z), INTEGER (I-N) 991 CALL INTPL8 (SOLD, SNEW, XOLD, XNEW, NVAR, JJOLD, JJNEW)
927 C*****END precision > double 992 C
928 C 993 RETURN
929 C*****precision > single 994 END
930 C IMPLICIT REAL (A-H, O-Z), INTEGER (I-N) 995 C
931 C*****END precision > single 996 C------------------------------------------------------------------
932 C 997 SUBROUTINE INTPL8 (F1, F2, X1, X2, MVAR, N1, N2)
933 DIMENSION SNEW(NVAR,*), SOLD(NVAR,*), XNEW(*), XOLD(*), WORK(*) 998 C------------------------------------------------------------------
934 C 999 C
935 C interpolate the solution onto an equidistributing mesh 1000 C interpolate to get f2(x2) from f1(x1)
936 C 1001 C
937 C 1002 C
938 C compute weight coefficients 1003 C*****precision > double
939 C 1004 IMPLICIT DOUBLE PRECISION (A-H, O-Z), INTEGER (I-N)
940 R0 = 1. - PR 1005 C*****END precision > double
941 R1 = PR *R /(R+1.) 1006 C
942 R2 = PR - R1 1007 C*****precision > single
943 TV1 = 0. 1008 C IMPLICIT REAL (A-H, O-Z), INTEGER (I-N)
944 DO 10 I = 2, JJOLD 1009 C*****END precision > single
945 TV1 = TV1 + ABS( SOLD(N,I) - SOLD(N,I-1) ) 1010 C
946 10 CONTINUE 1011 DIMENSION F1(MVAR,*), F2(MVAR,*), X1(*), X2(*)
947 TV2 = 0. 1012 DO 100 J = 2, N2-1
948 DO 20 I = 2, JJOLD-1 1013 XVAL = X2(J)
949 TV2 = TV2 + ABS( (SOLD(N,I+1)-SOLD(N,I))/(XOLD(I+1)-XOLD(I)) 1014 DO 50 I = 2, N1
950 1 - (SOLD(N,I)-SOLD(N,I-1))/(XOLD(I)-XOLD(I-1)) ) 1015 IF(XVAL.LE.X1(I)) THEN
951 20 CONTINUE 1016 DEL = (XVAL-X1(I-1))/(X1(I)-X1(I-1))
952 XLEN = XOLD(JJOLD) - XOLD(1) 1017 DO 40 L = 1, MVAR
953 B1 = R1 * XLEN /(R0 * TV1) 1018 F2(L,J)=F1(L,I-1)+(F1(L,I)-F1(L,I-1))*DEL
954 B2 = R2 * XLEN /(R0 * TV2) 1019 40 CONTINUE
955 C 1020 GOTO 80
956 C compute partial sums of weight function 1021 ENDIF
957 C 1022 50 CONTINUE
958 WORK(1) = 0. 1023 WRITE(7,*) ' *** STOP...INTERPOLATION ERROR.'
959 DO 50 I = 2, JJOLD 1024 80 CONTINUE
960 DX = XOLD(I) -XOLD(I-1) 1025 100 CONTINUE
961 WORK(I) = DX + B1*ABS( SOLD(N,I) - SOLD(N,I-1) ) + WORK(I-1) 1026 C
962 1 + B2*ABS( (SOLD(N,I+1)-SOLD(N,I))/(XOLD(I+1)-XOLD(I)) 1027 C endpoints..
963 2 - (SOLD(N,I)-SOLD(N,I-1))/(XOLD(I)-XOLD(I-1)) ) 1028 C
964 50 CONTINUE 1029 DO 110 L = 1, MVAR
965 DO 65 I = 2, JJOLD 1030 F2(L,1) = F1(L,1)
966 WORK(I) = WORK(I)/WORK(JJOLD) 1031 F2(L,N2) = F1(L,N1)
967 65 CONTINUE 1032 110 CONTINUE
968 C 1033 END
969 C interpolate onto uniform eta grid to find new x 1034 C
970 C 1035 C-------------------------------------------------------------------
971 XNEW(1) = XOLD(1) 1036 C
972 XNEW(JJNEW) = XOLD(JJOLD) 1037 SUBROUTINE JACOB(KK,POINTS,COMPS,LENRGY,LRADI,LNARR,LMULTI,LTDIF,
973 ISTART = 2 1038 1 LTIME, PR, WT, YOXID, YFUEL, DT,
974 DETA = 1./FLOAT(JJNEW-1) 1039 2 NTEMP, XGIVEN, TGIVEN, MESH, SN, X0, WNDFAC,
975 DO 80 J = 2, JJNEW-1 1040 3 ABSOL, RELAT, SCRTCH, YV,
1041 4 TFUEL, TOXID, VFUEL, VOXID, AFUEL, AOXID, 1106 PERTRB(INDEX) = ABS(X0(COMP, POINT)) * RELAT + ABSOL
1042 4 VISC, COND, D, DKJ, 1107 X0(COMP, POINT) = X0(COMP, POINT) + PERTRB(INDEX)
1043 5 TDR, ICKWRK, RCKWRK, IMCWRK, RMCWRK, Y1, Y0, 1108 0100 CONTINUE
1044 6 A, PERTRB, SAVE, RFAC) 1109 C
1045 C 1110 C/// CALL THE FUNCTION AT THE PERTRBED X0 AND STORE THE RESULT IN Y1.
1046 C*****precision > double 1111 C
1047 IMPLICIT DOUBLE PRECISION (A-H, O-Z), INTEGER (I-N) 1112 C WRITE(*,*) 'Third call to fun!,Line1049'
1048 DOUBLE PRECISION MESH 1113 CALL FUN(KK,POINTS,COMPS,LENRGY,LRADI,LNARR,LMULTI,LTDIF,.FALSE.,
1049 C*****END precision > double 1114 1 LTIME, PR, WT, YOXID, YFUEL, DT, NTEMP, XGIVEN,
1050 C 1115 2 TGIVEN, MESH, SN, X0, WNDFAC, SCRTCH(1, 1),
1051 C*****precision > single 1116 3 SCRTCH(1,5), YV, SCRTCH(1, 2), SCRTCH(1, 3),
1052 C IMPLICIT REAL (A-H, O-Z), INTEGER (I-N) 1117 4 SCRTCH(1, 4), TFUEL, TOXID, VFUEL, VOXID,
1053 C REAL MESH 1118 4 AFUEL, AOXID,
1054 C*****END precision > single 1119 5 VISC, COND, D, DKJ, TDR, ICKWRK, RCKWRK,
1055 C 1120 6 IMCWRK, RMCWRK, Y1, RFAC)
1056 INTEGER CLASS, COMP, COMPS, INDEX, LOWER, POINT, POINTS, ROWS, 1121 C
1057 1 UPPER 1122 C/// RESTORE X0 TO ITS ORIGINAL VALUE.
1058 LOGICAL LENRGY, LRADI, LNARR, LMULTI, LTDIF, LTIME 1123 C
1059 C 1124 INDEX = 0
1060 PARAMETER (LOWER = 1, UPPER = 1, MODULS = LOWER + 1 + UPPER) 1125 DO 0200 POINT = CLASS, POINTS, MODULS
1061 C 1126 INDEX = INDEX + 1
1062 DIMENSION A((2 * LOWER + UPPER + 3) * COMPS - 2, *), 1127 X0(COMP, POINT) = SAVE(INDEX)
1063 1 MESH(*), PERTRB(*), SAVE(*), SCRTCH(KK,*), 1128 0200 CONTINUE
1064 2 X0(COMPS, *), Y0(COMPS, *), Y1(COMPS, *) 1129 C
1065 C 1130 C/// DIFFERENCE TO GET THE COLUMNS OF THE JACOBIAN.
1066 KOFSET = (LOWER + UPPER + 2) * COMPS - 1 1131 C
1067 ROWS = (2 * LOWER + UPPER + 3) * COMPS - 2 1132 INDEX = 0
1068 C 1133 DO 0300 POINT = CLASS, POINTS, MODULS
1069 C/// ZERO THE MATRIX STORAGE SPACE. 1134 INDEX = INDEX + 1
1070 C 1135 TEMP = 1.0 / PERTRB(INDEX)
1071 ZERO = 0.0 1136 K = COMP + (POINT - 1) * COMPS
1072 C*****precision > double 1137 DO 0300 J2 = MAX (POINT - LOWER, 1),
1073 CALL DCOPY (ROWS * COMPS * POINTS, ZERO, 0, A, 1) 1138 + MIN (POINT + UPPER, POINTS)
1074 C*****END precision > double 1139 JOFSET = (J2 - 1) * COMPS - K + KOFSET
1075 C 1140 DO 0300 J1 = 1, COMPS
1076 C*****precision > single 1141 A(J1 + JOFSET, K) = (Y1(J1, J2) - Y0(J1, J2)) * TEMP
1077 C CALL SCOPY (ROWS * COMPS * POINTS, ZERO, 0, A, 1) 1142 0300 CONTINUE
1078 C*****END precision > single 1143 C
1079 C 1144 C/// BOTTOM OF THE LOOPS OVER THE RESIDUE CLASSES AND SOLUTION
1080 C 1145 C/// COMPONENTS.
1081 C/// CALL THE FUNCTION AT X0 AND STORE IN Y0. 1146 C
1082 C 1147 0400 CONTINUE
1083 C WRITE(*,*) 'Second call to fun!,Line1020' 1148 C
1084 CALL FUN (KK,POINTS,COMPS,LENRGY,LRADI,LNARR,LMULTI,LTDIF,.TRUE., 1149 RETURN
1085 1 LTIME, PR, WT, YOXID, YFUEL, DT, NTEMP, XGIVEN, 1150 END
1086 2 TGIVEN, MESH, SN, X0, WNDFAC, SCRTCH(1, 1), 1151 C
1087 3 SCRTCH(1,5), YV, SCRTCH(1, 2), SCRTCH(1, 3), 1152 C---------------------------------------------------------------------
1088 4 SCRTCH(1, 4), TFUEL, TOXID, VFUEL, VOXID, 1153 C
1089 4 AFUEL, AOXID, 1154 SUBROUTINE FUN (KK, JJ, NATJ, LENRGY,LRADI,LNARR,LMULTI, LTDIF,
1090 5 VISC, COND, D, DKJ, TDR, ICKWRK, RCKWRK, IMCWRK, 1155 1 LVARMC, LTIME, PR, WT, YOXID, YFUEL, DT, NTEMP,
1091 6 RMCWRK, Y0, RFAC) 1156 2 XGIVEN, TGIVEN, X, SN, S, WNDFAC, YAV, XAV, YV,
1092 C 1157 3 WDOT, CP, H, TFUEL, TOXID, VFUEL, VOXID,
1093 C 1158 4 AFUEL, AOXID, VISC,
1094 C/// TOP OF THE LOOPS OVER THE RESIDUE CLASSES AND SOLUTION COMPONENTS. 1159 5 COND, D, DKJ, TDR, ICKWRK, RCKWRK, IMCWRK, RMCWRK,
1095 C 1160 6 F, RFAC)
1096 DO 0400 CLASS = 1, MODULS 1161 C
1097 DO 0400 COMP = 1, COMPS 1162 C This subroutine forms the function for the differential equations
1098 C 1163 C Indicies:
1099 C/// FOR A GIVEN RESIDUE CLASS AND A GIVEN SOLUTION COMPONENT, 1164 C NT : Temperature
1100 C/// PERTRB THE X0 VECTOR AT POINTS IN THE SAME RESIDUE CLASS. 1165 C NG : Velocity gradient
1101 C 1166 C NF : Velocity
1102 INDEX = 0 1167 C NH : Eigenvalue, pressure gradient
1103 DO 0100 POINT = CLASS, POINTS, MODULS 1168 C KK : no. of species
1104 INDEX = INDEX + 1 1169 C JJ : no. of gridpoints
1105 SAVE(INDEX) = X0(COMP, POINT) 1170 C NATJ : no. variables at each point
1171 C LENERGY = .TRUE. if solve energy equation 1236 CALL CKRHOY (PR, TAV, YAV, ICKWRK, RCKWRK, RHOP)
1172 C LMULTI . 1237 C
1173 C LTDIFF 1238 DO 30 K = 1, KK
1174 C LVARMC 1239 C?? F(NYS+K, 1) = RHOP*YFUEL(K)*VFUEL
1175 C LTIME 1240 C?? 1 - RHOP*S(NYS+K,1)*VFUEL - RHOP*YV(K,1)
1176 C VFUEL, VOXID, TFUEL, TOXID, YFUEL, YOXID... 1241 F(NYS+K, 1) = S(NYS+K, 1) - YFUEL(K)
1177 C 1242 30 CONTINUE
1178 C--------------------------------------------------------------------- 1243 IF (LENRGY) THEN
1179 C 1244 F(NT,1) = TFUEL - S(NT,1)
1180 C*****precision > double 1245 ELSE
1181 IMPLICIT DOUBLE PRECISION (A-H, O-Z), INTEGER (I-N) 1246 CALL TEMP (NTEMP, X(1), XGIVEN, TGIVEN, TI)
1182 C*****END precision > double 1247 F(NT, 1) = S(NT,1) - TI
1183 C 1248 ENDIF
1184 C*****precision > single 1249 F(NF,1) = S(NF,1) - VFUEL*RHOP/2.
1185 C IMPLICIT REAL (A-H, O-Z), INTEGER (I-N) 1250 F(NG,1) = AFUEL*RHOP - S(NG,1)
1186 C*****END precision > single 1251 F(NH,1) = S(NH,2) - S(NH,1)
1187 C 1252 C
1188 DIMENSION PMAC(7,4),NRAD(4), NNARROW(6) 1253 C INTERIOR MESH POINTS
1189 COMMON /LOCS/ NT, NG, NF, NH, NYS, NY 1254 C
1190 C 1255 HASNARR=0.0
1191 COMMON/CLAYER/ QL(100),QLS(100,600) 1256 C
1192 C 1257 IF (LNARR) THEN
1193 COMMON/RCDATA/NOL,XLL(100),TL(100),RSC(7,100),WALLT,WAVMAX,WAVMIN 1258 C
1194 1259 C OPEN(7,FILE='RC.DAT')
1195 COMMON /RAD/ PMAC, NRAD, NNARROW 1260 C
1196 1261 C WRITE(7,9999) JJ-1
1197 DIMENSION WT(*), YOXID(*), YFUEL(*), XGIVEN(*), TGIVEN(*), X(*), 1262 C9999 FORMAT(I2)
1198 1 S(NATJ,*), SN(NATJ,*), YAV(*), XAV(*), YV(KK,*), 1263 C
1199 2 WDOT(*), CP(*), H(*), VISC(*), COND(*), D(KK,*), 1264 NOL = JJ-1
1200 3 TDR(KK,*), DKJ(KK,KK,*), ICKWRK(*), RCKWRK(*), 1265
1201 4 IMCWRK(*), RMCWRK(*), F(NATJ,*) 1266 DO 9998 I=1,JJ-1
1202 C 1267
1203 LOGICAL LENRGY,LRADI,LNARR,LMULTI, LTDIF, LVARMC, LTIME 1268 CALL CKYTX(S(NY,I),ICKWRK, RCKWRK,XFRACT)
1204 DOUBLE PRECISION XFRACT(KK) 1269
1205 C 1270 DELTAX=(X(I+1)-X(I))/100
1206 C 1271 C
1207 C timelimit 1272 C WRITE(7,9997) DELTAX,S(NT,I),
1208 C 1273 C 1 (XFRACT(NNARROW(K))*101.,K=1,6),0.0
1209 C*****unicos timelimit 1274 C
1210 C call tremain (secs) 1275 C9997 FORMAT(F8.5,1X,F7.2,1X,7(D9.3,1X))
1211 C if (secs .lt. 60.0) stop 1276 C
1212 C*****END unicos timelimit 1277 XLL(I) = DELTAX
1213 C 1278
1214 C EVALUATE AND STORE THE TRANSPORT COEFFICIENTS 1279 TL(I) = S(NT,I)
1215 C 1280
1216 C WRITE(*,*) 'Now in fun()!' 1281 DO 9997,K=1,6
1217 IF (LVARMC .AND. .NOT.LTIME) THEN 1282
1218 CALL MTRNPR (KK, JJ, NATJ, LENRGY, LMULTI, LTDIF, PR, X, S, 1283 RSC(K,I) = XFRACT(NNARROW(K))*101.
1219 1 YAV, ICKWRK, RCKWRK, IMCWRK, RMCWRK, 1284
1220 2 H, CP, XAV, VISC, COND, D, TDR, DKJ) 1285 9997 CONTINUE
1221 ENDIF 1286
1222 C 1287 RSC(7,I) = 0.
1223 C EVALUATE AND STORE THE DIFFUSION VELOCITIES 1288
1224 C (in the following call h(*) and cp(*) are used 1289 9998 CONTINUE
1225 C temporarily for scratch space) 1290 C
1226 C 1291 C
1227 CALL MDIFV (KK, JJ, NATJ, LMULTI, LTDIF, X, S, WT, YAV, 1292 C WRITE(7,9996) 298.0,500.0,10000.0
1228 1 H, CP, PR, D, TDR, ICKWRK, RCKWRK, YV, DKJ) 1293 C
1229 C 1294 C9996 FORMAT(F5.1,1X,F5.1,1X,F7.1)
1230 C FUEL BOUNDARY (X=0) 1295 C
1231 C 1296 C WRITE(7,9995) 0
1232 TAV = 0.5 * (S(NT,1) + S(NT,2)) 1297 C
1233 DO 10 K = 1, KK 1298 C9995 FORMAT(I1)
1234 YAV(K) = 0.5 * (S(NYS+K,1) + S(NYS+K,2)) 1299 C
1235 10 CONTINUE 1300 C
1301 WALLT = 298. 1366 WINDA = 0.
1302 1367 WINDB = 1.
1303 WAVMAX = 10000. 1368 WINDC = -1.
1304 1369 ENDIF
1305 WAVMIN = 500. 1370 DXWIND = WINDA*X(J+1) + WINDB*X(J) + WINDC*X(J-1)
1306 1371 ENDIF
1307 C CLOSE(7) 1372 C
1308 C 1373 C SPECIES CONSERVATION EQUATION
1309 CCCC CALL RADCAL 1374 C
1310 C 1375 IF (WNDFAC .EQ. 1.0) THEN
1311 1376 DO 200 K = 1, KK
1312 HASNARR = 1.0 1377 DYK = WINDA*S(NYS+K,J+1) + WINDB*S(NYS+K,J)
1313 C 1378 1 + WINDC*S(NYS+K,J-1)
1314 LENRGY = .TRUE. 1379 F(NYS+K,J) = 2.0*S(NF,J) * DYK / DXWIND
1315 1380 1 - WDOT(K) * WT(K)
1316 ENDIF 1381 2 + (RHOP*YV(K,J) - RHOM*YV(K,J-1)) / DXAV
1317 1382 200 CONTINUE
1318 DO 1000 J = 2, JJ-1 1383 ELSE
1319 C 1384 DO 210 K = 1, KK
1320 TAV = 0.5 * (S(NT,J) + S(NT,J+1)) 1385 F(NYS+K,J) = 2.0*S(NF,J) *
1321 DO 100 K = 1, KK 1386 1 (CENDFP*S(NYS+K,J+1) + CENDFC*S(NYS+K,J) +
1322 YAV(K) = 0.5 * (S(NYS+K,J) + S(NYS+K,J+1)) 1387 1 CENDFM*S(NYS+K,J-1) ) -
1323 100 CONTINUE 1388 2 WDOT(K) * WT(K) +
1324 C 1389 3 (RHOP*YV(K,J) - RHOM*YV(K,J-1)) / DXAV
1325 RHOM = RHOP 1390 210 CONTINUE
1326 CALL CKRHOY (PR, TAV, YAV, ICKWRK, RCKWRK, RHOP) 1391 ENDIF
1327 CALL CKRHOY (PR, S(NT,J), S(NY,J), ICKWRK, RCKWRK, RHOJ) 1392 C
1328 CALL CKRHOY (PR, S(NT,J+1), S(NY,J+1), ICKWRK, RCKWRK, RHOJP1) 1393 C STREAMFUNCTION EQUATION
1329 CALL CKRHOY (PR, S(NT,J-1), S(NY,J-1), ICKWRK, RCKWRK, RHOJM1) 1394 C
1330 C 1395 F(NF,J) = (S(NG,J)+S(NG,J-1))/2. - (S(NF,J) - S(NF,J-1))/DXM
1331 C FORM THE CHEMICAL RATE TERMS 1396 C
1332 C 1397 C MOMENTUM EQUATION
1333 CALL CKWYP (PR, S(NT,J), S(NY,J), ICKWRK, RCKWRK, WDOT) 1398 C
1334 CALL CKHML (S(NT,J), ICKWRK, RCKWRK, H) 1399 BIGCP = VISC(J)
1335 CALL CKCPBS (S(NT,J), S(NY,J), ICKWRK, RCKWRK, CPB) 1400 BIGCM = VISC(J-1)
1336 CALL CKCPMS (S(NT,J), ICKWRK, RCKWRK, CP) 1401 IF (WNDFAC .EQ. 1.0) THEN
1337 1402 DUGR = WINDA*( S(NF,J+1)*S(NG,J+1)/RHOJP1 )
1338 C 1403 1 + WINDB*( S(NF,J )*S(NG,J )/RHOJ )
1339 C ACCELERATE RATES TO FACILITATE SOLUTION? 1404 2 + WINDC*( S(NF,J-1)*S(NG,J-1)/RHOJM1 )
1340 C 1405 F(NG,J) = - 3.0*S(NG,J)**2/RHOJ +
1341 DO 120 K = 1, KK 1406 2 2.0 * DUGR / DXWIND
1342 WDOT(K) = WDOT(K) *RFAC 1407 3 - ( BIGCP * (S(NG,J+1)/RHOJP1 - S(NG,J)/RHOJ) / DXP -
1343 120 CONTINUE 1408 3 BIGCM * (S(NG,J)/RHOJ -S(NG,J-1)/RHOJM1)/DXM ) /
1344 C 1409 3 DXAV - S(NH,J)
1345 C FORM THE MESH DIFFERENCES 1410 ELSE
1346 C 1411 F(NG,J) = - 3.0*S(NG,J)**2/RHOJ +
1347 DXP = (X(J+1) - X(J) ) 1412 2 2.0 * (CENDFP*S(NG,J+1)*S(NF,J+1)/RHOJP1 +
1348 DXM = (X(J) - X(J-1)) 1413 2 CENDFC*S(NG,J)*S(NF,J)/RHOJ +
1349 DXAV = 0.5 * (X(J+1) - X(J-1)) 1414 2 CENDFM*S(NG,J-1)*S(NF,J-1)/RHOJM1 ) -
1350 DXPM = (X(J+1) - X(J-1)) 1415 3 ( BIGCP * (S(NG,J+1)/RHOJP1 - S(NG,J)/RHOJ) / DXP -
1351 C 1416 3 BIGCM * (S(NG,J)/RHOJ -S(NG,J-1)/RHOJM1)/DXM ) /
1352 C FORM THE COEFFICIENTS FOR CENTRAL DIFFERENCES 1417 3 DXAV -
1353 C 1418 4 S(NH,J)
1354 CENDFM = - DXP / (DXM*DXPM) 1419 ENDIF
1355 CENDFC = (DXP-DXM) / (DXP*DXM) 1420 C
1356 CENDFP = DXM / (DXP*DXPM) 1421 C ENERGY EQUATION
1357 C 1422 C
1358 C WINDWARD DIFFERENCING FACTORS 1423
1359 C 1424
1360 IF (WNDFAC .EQ. 1.0) THEN 1425 IF (LENRGY) THEN
1361 IF (S(NF,J) .LE. 0.) THEN 1426
1362 WINDA = 1. 1427 C Form the radiation coefficient of flame
1363 WINDB = -1. 1428
1364 WINDC = 0. 1429 CALL RADEMIS(KK, ICKWRK, RCKWRK, PR, S(NT,J), S(NY,J), PMAC,
1365 ELSE 1430 1 NRAD, 4, GASEMI)
1431 C 1496 CALL CKRHOY (PR, S(NT,JJ), S(NY,JJ), ICKWRK, RCKWRK, RHOJ)
1432 HASRAD = 0.0 1497 C
1433 1498 DO 1200 K = 1, KK
1434 IF (LRADI) THEN 1499 C.. F(NYS+K, JJ) = S(NYS+K,JJ) - YOXID(K)
1435 HASRAD = 1.0 1500 F(NYS+K, JJ) = RHOJ*YOXID(K)*VOXID
1436 ENDIF 1501 1 - RHOJ*S(NYS+K,JJ)*VOXID - RHOJ*YV(K,JJ-1)
1437 1502 1200 CONTINUE
1438 C*****precision > double 1503 C
1439 TDOT = DDOT (KK, WDOT, 1, H, 1) 1504 IF (LENRGY) THEN
1440 C*****END precision > double 1505 F(NT, JJ) = S(NT,JJ) - TOXID
1441 C*****precision > single 1506 ELSE
1442 C TDOT = SDOT (KK, WDOT, 1, H, 1) 1507 CALL TEMP (NTEMP, X(JJ), XGIVEN, TGIVEN, TI)
1443 C*****END precision > single 1508 F(NT, JJ) = S(NT,JJ) - TI
1444 C 1509 ENDIF
1445 SUM = 0.0 1510 C
1446 DO 400 K = 1, KK 1511 F(NH,JJ) = RHOJ*VOXID/2.0 - S(NF,JJ)
1447 SUM = SUM + 0.5 * (YV(K,J) + YV(K,J-1)) * 1512 F(NG,JJ) = S(NG,JJ) - AOXID*RHOJ
1448 1 CP(K) * (CENDFP*S(NT,J+1) + CENDFC*S(NT,J) + 1513 F(NF,JJ) = (S(NG,JJ)+S(NG,JJ-1))/2. - (S(NF,JJ) - S(NF,JJ-1))/DXP
1449 2 CENDFM*S(NT,J-1) ) 1514 C
1450 400 CONTINUE 1515 C ADD THE TIME STEP, IF NEEDED
1451 C 1516 C
1452 C WRITE(*,2999) HASNARR*QL(J)*3.14159*2/(X(J+1)-X(J))*1000 1517 IF (.NOT. LTIME) RETURN
1453 2999 FORMAT(8X,'Pr=',E12.6,'erg/cm3') 1518 C
1454 1519 DO 2500 J = 2, JJ-1
1455 IF (WNDFAC .EQ. 1.0) THEN 1520 C
1456 DTEMP = WINDA*S(NT,J+1) + WINDB*S(NT,J) + WINDC*S(NT,J-1) 1521 CALL CKRHOY (PR, S(NT,J), S(NY,J), ICKWRK, RCKWRK, RHO)
1457 F(NT,J) = 1522 C
1458 1 2.0*S(NF,J) * DTEMP / DXWIND - 1523 DO 2200 K = 1, KK-1
1459 1 ( COND(J)*(S(NT,J+1)-S(NT,J))/DXP - 1524 DYDT = (S(NYS+K,J) - SN(NYS+K,J)) / DT
1460 2 COND(J-1)*(S(NT,J)-S(NT,J-1))/DXM ) / 1525 F(NYS+K,J) = F(NYS+K,J) + RHO*DYDT
1461 3 (CPB*DXAV) + 1526 2200 CONTINUE
1462 4 RHOJ * SUM / CPB + 1527
1463 5 TDOT / CPB 1528 DGDT = (S(NG,J) - SN(NG,J)) / DT
1464 6 + HASRAD*(S(NT,J)**4-298.**4)*(5.67E-8)*4*10*GASEMI/CPB 1529 F(NG,J) = F(NG,J) + DGDT
1465 7 + HASNARR*QL(J)*3.14159*2/(X(J+1)-X(J))*1000./CPB 1530 C
1466 1531 IF (LENRGY) THEN
1467 ELSE 1532 DTDT = (S(NT,J) - SN(NT,J)) / DT
1468 F(NT,J) = 2.0*S(NF,J) * 1533 F(NT,J) = F(NT,J) + RHO*DTDT
1469 1 (CENDFP*S(NT,J+1) + CENDFC*S(NT,J) + 1534 ENDIF
1470 1 CENDFM*S(NT,J-1) ) - 1535 2500 CONTINUE
1471 1 ( COND(J)*(S(NT,J+1)-S(NT,J))/DXP - 1536 C
1472 2 COND(J-1)*(S(NT,J)-S(NT,J-1))/DXM ) / 1537 RETURN
1473 3 (CPB*DXAV) + 1538 END
1474 4 RHOJ * SUM / CPB + 1539 C--------------------------------------------------------------------
1475 5 TDOT / CPB 1540 SUBROUTINE MTRNPR (KK, JJ, NATJ, LENRGY, LMULTI, LTDIF, PR, X, S,
1476 6 + HASRAD*(S(NT,J)**4-298.**4)*(5.67E-8)*4*10*GASEMI/CPB 1541 1 YAV, ICKWRK, RCKWRK, IMCWRK, RMCWRK,
1477 7 + HASNARR*QL(J)*3.14159*2/(X(J+1)-X(J))*1000./CPB 1542 2 XMF, XMFP, XAV, VISC, COND, D, TDR, DKJ)
1478 1543 C
1479 ENDIF 1544 C*****precision > double
1480 1545 IMPLICIT DOUBLE PRECISION (A-H, O-Z), INTEGER (I-N)
1481 1546 C*****END precision > double
1482 C 1547 C
1483 ELSE 1548 C*****precision > single
1484 CALL TEMP (NTEMP, X(J), XGIVEN, TGIVEN, TI) 1549 C IMPLICIT REAL (A-H, O-Z), INTEGER (I-N)
1485 F(NT, J) = S(NT,J) - TI 1550 C*****END precision > single
1486 ENDIF 1551 C
1487 C 1552 COMMON /LOCS/ NT, NG, NF, NH, NYS, NY
1488 C EIGENVALUE EQUATION, H 1553 C
1489 C 1554 DIMENSION X(*), S(NATJ,*), YAV(*), XAV(*), XMF(*),
1490 F(NH,J) = S(NH,J+1) - S(NH,J) 1555 1 XMFP(*), ICKWRK(*), RCKWRK(*), IMCWRK(*), RMCWRK(*),
1491 C 1556 2 VISC(*), COND(*), DKJ(KK,KK,*), D(KK,*), TDR(KK,*)
1492 1000 CONTINUE 1557 LOGICAL LENRGY, LTDIF, LMULTI
1493 C 1558 C
1494 C OXIDIZER BOUNDARY (X=L) 1559 DATA EPS/1.0E-30/
1495 C 1560 C
1561 IF (LMULTI) THEN 1626 1 XMF, XMFP, PR, D, TDR, ICKWRK, RCKWRK, YV, DKJ)
1562 C-----------------MULTICOMPONENT FORMULAS--------------------- 1627 C
1563 CALL CKYTX (S(NY,1), ICKWRK, RCKWRK, XMFP) 1628 C*****precision > double
1564 DO 200 J = 1, JJ-1 1629 IMPLICIT DOUBLE PRECISION (A-H, O-Z), INTEGER (I-N)
1565 TAV = 0.5 * (S(NT,J) + S(NT,J+1)) 1630 C*****END precision > double
1566 C 1631 C
1567 C DIMENSIONAL TEMPERATURE AT THE GRID POINTS 1632 C*****precision > single
1568 C 1633 C IMPLICIT REAL (A-H, O-Z), INTEGER (I-N)
1569 DO 100 K = 1, KK 1634 C*****END precision > single
1570 YAV(K) = 0.5 * (S(NYS+K,J) + S(NYS+K,J+1)) 1635 C
1571 XMF(K) = XMFP(K) 1636 COMMON /LOCS/ NT, NG, NF, NH, NYS, NY
1572 100 CONTINUE 1637 C
1573 CALL CKYTX (YAV, ICKWRK, RCKWRK, XAV) 1638 DIMENSION X(*), S(NATJ,*), WT(*), YAV(*), XMF(*), XMFP(*),
1574 CALL CKYTX (S(NY,J+1), ICKWRK, RCKWRK, XMFP) 1639 1 D(KK,*), TDR(KK,*), DKJ(KK,KK,*), YV(KK,*),
1575 CALL CKMMWX (XAV, ICKWRK, RCKWRK, WTMAV) 1640 2 ICKWRK(*), RCKWRK(*)
1576 CALL MCMDIF (PR, TAV, XAV, KK, IMCWRK, RMCWRK, DKJ(1,1,J)) 1641 LOGICAL LTDIF, LDIRCT, LMULTI
1577 CALL MCAVIS (TAV, XAV, RMCWRK, VISC(J)) 1642 C
1578 DO 75 K = 1, KK 1643 LDIRCT = .TRUE.
1579 SUMN = 0.0 1644 C
1580 DO 50 L = 1, KK 1645 C
1581 CC SUMN = SUMN + WT(L) * DKJ(K,L,J) * 1646 CALL CKYTX (S(NY,1), ICKWRK, RCKWRK, XMFP)
1582 CC 1 (XMFP(L) - XMF(L)) / (X(J+1) - X(J)) 1647 C
1583 SUMN = SUMN + DKJ(K,L,J) * 1648 C LOOP OVER ALL MESH POINTS, COMPUTING THE DIFFUSION
1584 1 (S(NYS+L,J+1) - S(NYS+L,J)) / (X(J+1) - X(J)) 1649 C VELOCITY AT THE MID POINTS. THE INDEXING IS SUCH THAT
1585 50 CONTINUE 1650 C YV(K,J) IS THE DIFFUSION VELOCITY OF THE KTH SPECIES
1586 CC DENOM = - (XMFP(K) - XMF(K)) / (X(J+1) - X(J)) 1651 C MIDWAY BETWEEN NODES J AND J+1.
1587 DENOM = - (S(NYS+K,J+1) - S(NYS+K,J)) / (X(J+1) - X(J)) 1652 C
1588 D(K,J) = (SUMN + EPS) / ( WTMAV * (DENOM + EPS)) 1653 DO 1000 J = 1, JJ-1
1589 75 CONTINUE 1654 C
1590 C 1655 TMF = S(NT, J)
1591 C DETERMINE THE MIXTURE CONDUCTIVITY AND 1656 TMFP = S(NT, J+1)
1592 C THERMAL DIFFUSION COEFFICIENT AT J 1657 TAV = 0.5 * (TMF + TMFP)
1593 C 1658 DO 300 K = 1, KK
1594 IF (LENRGY .OR. LTDIF) 1659 YAV(K) = 0.5 * (S(NYS+K,J) + S(NYS+K,J+1))
1595 1 CALL MCMCDT (PR, TAV, XAV, IMCWRK, RMCWRK, 1660 XMF(K) = XMFP(K)
1596 2 ICKWRK, RCKWRK, TDR(1,J), COND(J)) 1661 300 CONTINUE
1597 C 1662 C
1598 200 CONTINUE 1663 CALL CKMMWY (YAV, ICKWRK, RCKWRK, WTM)
1599 RETURN 1664 CALL CKRHOY (PR, TAV, YAV, ICKWRK, RCKWRK, RHOAV)
1600 ELSE 1665 CALL CKYTX (S(NY,J+1), ICKWRK, RCKWRK, XMFP)
1601 C-----------------MIXTURE-AVERAGED FORMULAS--------------------- 1666 C
1602 C 1667 IF (LMULTI. AND. LDIRCT) THEN
1603 DO 400 J = 1, JJ-1 1668 C EVALUATE THE MULTICOMPONENT DIFFUSION VELOCITY DIRECTLY,
1604 TAV = 0.5 * (S(NT,J) + S(NT,J+1)) 1669 C RATHER THAN USE THE MIXTURE-AVERAGED FORM FOR D(K,J)
1605 C 1670 DO 475 K = 1, KK
1606 C DIMENSIONAL TEMPERATURE AT THE GRID POINTS 1671 SUM = 0.0
1607 C 1672 DO 450 L = 1, KK
1608 DO 300 K = 1, KK 1673 SUM = SUM + WT(L) * DKJ(K,L,J) *
1609 YAV(K) = 0.5 * (S(NYS+K,J) + S(NYS+K,J+1)) 1674 1 (XMFP(L)-XMF(L)) / (X(J+1)-X(J))
1610 300 CONTINUE 1675 450 CONTINUE
1611 CALL CKYTX (YAV, ICKWRK, RCKWRK, XAV) 1676 YV(K,J) = (WT(K)/WTM**2) * SUM
1612 CALL MCADIF (PR, TAV, XAV, RMCWRK, D(1,J)) 1677 475 CONTINUE
1613 CALL MCAVIS (TAV, XAV, RMCWRK, VISC(J)) 1678 ELSE
1614 C 1679 C USE MIXTURE-AVERAGED FORM FOR FICKIAN DIFFUSION,
1615 C DETERMINE THE MIXTURE CONDUCTIVITY AT J 1680 C WHETHER WE ARE USING THE MULTICOMPONENT FORMALISM
1616 C 1681 C OR MIXTURE-AVERAGED
1617 IF (LENRGY) CALL MCACON (TAV, XAV, RMCWRK, COND(J)) 1682 DO 500 K = 1, KK
1618 IF (LTDIF) CALL MCATDR (TAV, XAV, IMCWRK, RMCWRK, TDR(1,J)) 1683 YV(K,J) = - D(K,J) * (WT(K)/WTM) *
1619 C 1684 1 (XMFP(K)-XMF(K)) / (X(J+1)-X(J))
1620 400 CONTINUE 1685 500 CONTINUE
1621 ENDIF 1686 ENDIF
1622 RETURN 1687 C
1623 END 1688 C ADD THE THERMAL DIFFUSION, IF REQUESTED
1624 C-------------------------------------------------------------------- 1689 C
1625 SUBROUTINE MDIFV (KK, JJ, NATJ, LMULTI, LTDIF, X, S, WT, YAV, 1690 IF (LTDIF) THEN
2211 C 2276 C
2212 C-------------------------------------------------------------- 2277 C ABSOLUTE NEWTON ITERATION CONVERGENCE CRITERIA
2213 C 2278 C
2214 C READ NEXT INPUT LINE 2279 ELSE IF (KEYWRD .EQ. 'ATOL') THEN
2215 C 2280 CALL CKXNUM (LINE, 1, LOUT, NVAL, ATOL, IERR)
2216 WRITE (LOUT,'(/A/)') ' KEYWORD INPUT ' 2281 KERR = KERR.OR.IERR
2217 C 2282 C
2218 90 CONTINUE 2283 C RELATIVE NEWTON ITERATION CONVERGENCE CRITERIA
2219 KEYWRD = ' ' 2284 C
2220 LINE = ' ' 2285 ELSE IF (KEYWRD .EQ. 'RTOL') THEN
2221 READ (LIN, 7000) KEYWRD, LINE 2286 CALL CKXNUM (LINE, 1, LOUT, NVAL, RTOL, IERR)
2222 WRITE (LOUT, 8000) KEYWRD, LINE 2287 KERR = KERR.OR.IERR
2223 C 2288 C
2224 C IS THIS A KEYWORD COMMENT? 2289 C ABSOLUTE NEWTON CONVERGENCE CRITERIA FOR TIMESTEPS
2225 C 2290 C
2226 IF (KEYWRD(1:1) .EQ. '.' .OR. KEYWRD(1:1) .EQ. '/' 2291 ELSE IF (KEYWRD .EQ. 'ATIM') THEN
2227 1 .OR. KEYWRD(1:1) .EQ. '!') GO TO 90 2292 CALL CKXNUM (LINE, 1, LOUT, NVAL, ATIM, IERR)
2228 CALL UPCASE (KEYWRD) 2293 KERR = KERR.OR.IERR
2229 C 2294 C
2230 C--------------PROBLEM TYPE KEYWORDS-------------------- 2295 C RELATIVE NEWTON CONVERGENCE CRITERIA FOR TIMESTEPS
2231 C 2296 C
2232 C ENERGY EQUATION IS NOT INCLUDED 2297 ELSE IF (KEYWRD .EQ. 'RTIM') THEN
2233 C 2298 CALL CKXNUM (LINE, 1, LOUT, NVAL, RTIM, IERR)
2234 IF (KEYWRD .EQ. 'TGIV') THEN 2299 KERR = KERR.OR.IERR
2235 LENRGY = .FALSE. 2300 C
2236 C 2301 C TIME STEP STARTING PROCEDURE
2237 C ENERGY EQUATION IS INCLUDED 2302 C
2238 C 2303 ELSE IF (KEYWRD .EQ. 'TIME') THEN
2239 ELSE IF (KEYWRD .EQ. 'ENRG') THEN 2304 LTIME = .TRUE.
2240 LENRGY = .TRUE. 2305 CALL CKXNUM (LINE, 2, LOUT, NVAL, VALUE, IERR)
2241 C 2306 KERR = KERR.OR.IERR
2242 ELSE IF (KEYWRD .EQ. 'RADI') THEN 2307 NUMDT = INT(VALUE(1))
2243 LRADI = .TRUE. 2308 DT = VALUE(2)
2244 2309 C
2245 ELSE IF (KEYWRD .EQ. 'NARR') THEN 2310 C TIME STEPPING, AFTER ADDING THE ENERGY EQUATION
2246 LNARR = .TRUE. 2311 C
2247 2312 ELSE IF (KEYWRD .EQ. 'TIM2') THEN
2248 C FACTOR ON REACTION RATES 2313 LTIME2 = .TRUE.
2249 C 2314 CALL CKXNUM (LINE, 2, LOUT, NVAL, VALUE, IERR)
2250 ELSE IF (KEYWRD .EQ. 'GFAC') THEN 2315 KERR = KERR.OR.IERR
2251 CALL CKXNUM (LINE, 1, LOUT, NVAL, RATEF, IERR) 2316 NUMDT2 = INT(VALUE(1))
2252 KERR = KERR.OR.IERR 2317 DT2 = VALUE(2)
2253 C 2318 C
2254 C USE EQUILIBRIUM SOLUTION ESTIMATE 2319 C TIMESTEP INCREASE WHEN TIMESTEP DOES NOT CHANGE SOLUTION
2255 C 2320 C
2256 ELSE IF (KEYWRD .EQ. 'EQUL') THEN 2321 ELSE IF (KEYWRD .EQ. 'UFAC') THEN
2257 LEQUIL = .TRUE. 2322 CALL CKXNUM (LINE, 1, LOUT, NVAL, UFAC, IERR)
2258 C 2323 KERR = KERR.OR.IERR
2259 C FIX TEMP IN EQUIL SOLUTION ESTIMATE 2324 C
2260 C 2325 C TIMESTEP DECREASE WHEN NEWTON FAILS CONVERGENCE ON TIMESTEP
2261 ELSE IF (KEYWRD .EQ. 'FXEQ') THEN 2326 C
2262 LFXTEQ = .TRUE. 2327 ELSE IF (KEYWRD .EQ. 'DFAC') THEN
2263 C 2328 CALL CKXNUM (LINE, 1, LOUT, NVAL, DFAC, IERR)
2264 C LINEAR PROFILE ON INITIAL SOLUTION 2329 KERR = KERR.OR.IERR
2265 C 2330 C
2266 ELSE IF (KEYWRD .EQ. 'LINE') THEN 2331 C MINIMUM TIMESTEP
2267 IPROFL = 2 2332 C
2268 C 2333 ELSE IF (KEYWRD .EQ. 'DTMN') THEN
2269 C PROFILE ON INITIAL SOLUTION IS PLATEAU 2334 CALL CKXNUM (LINE, 1, LOUT, NVAL, DTMIN, IERR)
2270 C 2335 KERR = KERR.OR.IERR
2271 ELSE IF (KEYWRD .EQ. 'PLAT') THEN 2336 C
2272 IPROFL = 3 2337 C MAXIMUM TIMESTEP
2273 C 2338 C
2274 C--------------METHOD OPTIONS KEYWORDS-------------------- 2339 ELSE IF (KEYWRD .EQ. 'DTMX') THEN
2275 C 2340 CALL CKXNUM (LINE, 1, LOUT, NVAL, DTMAX, IERR)
2341 KERR = KERR.OR.IERR 2406 CALL CKXNUM (LINE, 1, LOUT, NVAL, GRAD, IERR)
2342 C 2407 KERR = KERR.OR.IERR
2343 C DO NOT DO THE FIXED TEMPERATURE PROBLEM 2408 C
2344 C 2409 C CURVATURE MESH ADAPTION PARAMETER
2345 ELSE IF (KEYWRD .EQ. 'NOFT') THEN 2410 C
2346 N1CALL = 2 2411 ELSE IF (KEYWRD .EQ. 'CURV') THEN
2347 C 2412 CALL CKXNUM (LINE, 1, LOUT, NVAL, CURV, IERR)
2348 C WINDWARD DIFFERENCING 2413 KERR = KERR.OR.IERR
2349 C 2414 C
2350 ELSE IF (KEYWRD .EQ. 'WDIF') THEN 2415 C NUMBER OF MESH POINTS TO ADD AT ONCE
2351 WNDFAC = 1.0 2416 C
2352 C 2417 ELSE IF (KEYWRD .EQ. 'NADP') THEN
2353 C CENTRAL DIFFERENCING 2418 CALL CKXNUM (LINE, 1, LOUT, NVAL, VALUE, IERR)
2354 C 2419 KERR = KERR.OR.IERR
2355 ELSE IF (KEYWRD .EQ. 'CDIF') THEN 2420 NADP = INT(VALUE(1))
2356 WNDFAC = 0.0 2421 C
2357 C 2422 C RETIREMENT PERIOD BEFORE INCREASING THE TIMESTEP
2358 C FLOOR VALUE FOR THE SPECIES BOUNDS 2423 C
2359 C 2424 ELSE IF (KEYWRD .EQ. 'IRET') THEN
2360 ELSE IF (KEYWRD .EQ. 'SFLR') THEN 2425 CALL CKXNUM (LINE, 1, LOUT, NVAL, VALUE, IERR)
2361 CALL CKXNUM (LINE, 1, LOUT, NVAL, SFLR, IERR) 2426 KERR = KERR.OR.IERR
2362 KERR = KERR.OR.IERR 2427 IRETIR = INT(VALUE(1))
2363 C 2428 C
2364 C OXIDIZER 2429 C INITIAL TIME STEPS BEFORE NEWTON
2365 C 2430 C
2366 ELSE IF (KEYWRD .EQ. 'KOUT') THEN 2431 ELSE IF (KEYWRD .EQ. 'ISTP') THEN
2367 CALL CKCRAY (LINE, KK, KSYM, LOUT, 6, KW, NFD, IERR) 2432 CALL CKXNUM (LINE, 1, LOUT, NVAL, VALUE, IERR)
2368 IF (IERR) THEN 2433 KERR = KERR.OR.IERR
2369 WRITE (LOUT,'(A)') 2434 NINIT = INT(VALUE(1))
2370 1 ' ERROR READING KEYWRD '//KEYWRD 2435 C
2371 KERR = .TRUE. 2436 C RETIREMENT AGE OF JACOBIAN DURING NEWTON
2372 IF (NFD .GT. 6) THEN 2437 C
2373 WRITE (LOUT,'(A,I4,A)') 'FOUND ', NFD,' SPECIES' 2438 ELSE IF (KEYWRD .EQ. 'NJAC') THEN
2374 NW = 6 2439 CALL CKXNUM (LINE, 1, LOUT, NVAL, VALUE, IERR)
2375 ENDIF 2440 KERR = KERR.OR.IERR
2376 ELSE 2441 NJAC = INT(VALUE(1))
2377 NW = NFD 2442 C
2378 ENDIF 2443 C RETIREMENT AGE OF JACOBIAN DURING TIME STEPPING
2379 C 2444 C
2380 C--------------GRID PARAMETER KEYWORDS-------------------- 2445 ELSE IF (KEYWRD .EQ. 'TJAC') THEN
2381 C 2446 CALL CKXNUM (LINE, 1, LOUT, NVAL, VALUE, IERR)
2382 C 2447 KERR = KERR.OR.IERR
2383 C NUMBER OF INITIAL MESH POINTS 2448 ITJAC = INT(VALUE(1))
2384 C (THIS IS OVERWRITTEN 'GRID' INPUT) 2449 C
2385 C 2450 C CENTER OF MIXING REGION
2386 ELSE IF (KEYWRD .EQ. 'NPTS') THEN 2451 C
2387 CALL CKXNUM (LINE, 1, LOUT, NVAL, VALUE, IERR) 2452 ELSE IF (KEYWRD .EQ. 'XCEN') THEN
2388 KERR = KERR.OR.IERR 2453 NOPT(2) = .TRUE.
2389 NPTS = INT(VALUE(1)) 2454 CALL CKXNUM (LINE, 1, LOUT, NVAL, XCEN, IERR)
2390 C 2455 KERR = KERR.OR.IERR
2391 C INITIAL MESH 2456 C
2392 C 2457 C DISTANCE AT WHICH END BOUNDARY CONDITION IS APPLIED
2393 ELSE IF (KEYWRD .EQ. 'GRID') THEN 2458 C
2394 LUMESH = .FALSE. 2459 ELSE IF (KEYWRD .EQ. 'XEND') THEN
2395 CALL CKXNUM (LINE, 1, LOUT, NVAL, VALUE, IERR) 2460 NEC(1) = .TRUE.
2396 IF (IERR .OR. NP+1.GT.NMAX) THEN 2461 CALL CKXNUM (LINE, 1, LOUT, NVAL, XEND, IERR)
2397 KERR = .TRUE. 2462 KERR = KERR.OR.IERR
2398 ELSE 2463 C
2399 NP = NP + 1 2464 C WIDTH OF MIXING ZONE
2400 X(NP) = VALUE(1) 2465 C
2401 ENDIF 2466 ELSE IF (KEYWRD .EQ. 'WMIX') THEN
2402 C 2467 NOPT(1) = .TRUE.
2403 C GRADIENT MESH ADAPTION PARAMETER 2468 CALL CKXNUM (LINE, 1, LOUT, NVAL, WMIX, IERR)
2404 C 2469 KERR = KERR.OR.IERR
2405 ELSE IF (KEYWRD .EQ. 'GRAD') THEN 2470 C
2731 C 2796 C
2732 C MAKE SURE THE INITIAL GRID POINTS ARE IN ORDER 2797 C THIS SUBROUTINE USES BISECTION TO LINEARLY INTERPOLATE
2733 C 2798 C AN ARRAY OF XX,TT PAIRS. GIVEN AN XX,TT PAIR THIS ROUTINE
2734 IF ((.NOT. CNTNUD) .AND. (.NOT. LUMESH)) THEN 2799 C RETURNS THE INTERPOLATED VALUE OF THE T AT THE POINT X.
2735 NPTS = NP 2800 C
2736 DO 5550 N = 2, NPTS 2801 C INPUT-
2737 IF (X(N-1) .GE. X(N)) THEN 2802 C NPTS - NUMBER OF XX,TT PAIRS.
2738 WRITE (LOUT, *) 2803 C X - LOCATION AT WHICH INTERPOLATED T IS DESIRED.
2739 1 ' ERROR...INITIAL GRID IS OUT OF ORDER' 2804 C XX - ARRAY OF X POINTS AT WHICH TT ARE GIVEN.
2740 STOP 2805 C TT - ARRAY OF T VALUES AT THE XX LOCATIONS.
2741 ENDIF 2806 C
2742 5550 CONTINUE 2807 C OUTPUT-
2743 ENDIF 2808 C T - INTERPOLATED T AT POINT X
2744 C 2809 C
2745 IF (NTEMP.GT.0) THEN 2810 CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC
2746 C 2811 C
2747 C MAKE SURE THE (X,T) PAIRS ARE IN ORDER 2812 C check for x outside (1,npts)
2748 C 2813 C
2749 DO 5500 N = 2, NTEMP 2814 IF (X .LE. XX(2)) THEN
2750 IF (XX(N-1) .GE. XX(N)) THEN 2815 N = 2
2751 WRITE (LOUT, *) 2816 S = (TT(N) - TT(N-1)) / (XX(N) - XX(N-1))
2752 1 ' ERROR...SPECIFIED TEMPERATURES ARE OUT OF ORDER' 2817 ELSEIF (X .GE. XX(NPTS-1)) THEN
2753 STOP 2818 N = NPTS-1
2754 ENDIF 2819 S = (TT(N+1) - TT(N)) / (XX(N+1) - XX(N))
2755 5500 CONTINUE 2820 ELSE
2756 C 2821 NLO = 1
2757 C MAKE SURE THE GIVEN TEMPERATURES SPAN THE XEND-XSTR DOMAIN 2822 NHI = NPTS
2758 C 2823 S = 0.0
2759 IF (.NOT.LRSTRT .OR. .NOT.CNTNUD .OR. LUSTGV) THEN 2824 C
2760 IF (XX(1).GT.0.0 .OR. XX(NTEMP).LT.XEND) THEN 2825 C bisect interval
2761 WRITE (LOUT, *) 2826 C
2762 1 ' ERROR...GIVEN TEMPERATURE PROFILE DOES NOT SPAN XEND-XSTR' 2827 50 CONTINUE
2763 STOP 2828 N = (NLO+NHI)/2
2764 ENDIF 2829 IF (X .LT. XX(N)) THEN
2765 ENDIF 2830 IF (X .LT. XX(N-1)) THEN
2766 ENDIF 2831 NHI = N
2767 C 2832 GO TO 50
2768 C SET OPTIONAL INPUT IF NEEDED 2833 ELSEIF (X .EQ. XX(N-1)) THEN
2769 C 2834 N = N-1
2770 IF (.NOT. NOPT(1)) WMIX = (XEND-0.0)*0.50 2835 ELSE
2771 C 2836 S = (TT(N) - TT(N-1)) / (XX(N) - XX(N-1))
2772 IF (.NOT. NOPT(2)) XCEN = (XEND-0.0)*0.35 2837 ENDIF
2773 C 2838 ELSEIF (X .GT. XX(N)) THEN
2774 C formats 2839 IF (X .GT. XX(N+1)) THEN
2775 C 2840 NLO = N
2776 RETURN 2841 GO TO 50
2777 7000 FORMAT(A4, A) 2842 ELSEIF (X .EQ. XX(N+1)) THEN
2778 8000 FORMAT(10X, A4, A76) 2843 N = N + 1
2779 END 2844 ELSE
2780 C 2845 S = (TT(N+1) - TT(N)) / (XX(N+1) - XX(N))
2781 C---------------------------------------------------------------------- 2846 ENDIF
2782 C 2847 ENDIF
2783 SUBROUTINE TEMP(NPTS, X, XX, TT, T) 2848 ENDIF
2784 C 2849 C
2785 C*****precision > double 2850 100 CONTINUE
2786 IMPLICIT DOUBLE PRECISION (A-H, O-Z), INTEGER (I-N) 2851 T = TT(N) + S * (X - XX(N))
2787 C*****END precision > double 2852 RETURN
2788 C 2853 END
2789 C*****precision > single 2854 C
2790 C IMPLICIT REAL (A-H, O-Z), INTEGER (I-N) 2855 C--------------------------------------------------------------------
2791 C*****END precision > single 2856 C
2792 C 2857 SUBROUTINE REASEN (II,KK,JJ,NATJ,LDA,LENRGY,LRADI,LNARR,LMULTI,
2793 DIMENSION XX(*), TT(*) 2858 1 LTDIF, LSAVE, LOUT, LVARMC, LTIME, PR,
2794 C 2859 2 WT, YOXID, YFUEL, DT, NTEMP,
2795 CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 2860 3 XGIVEN, TGIVEN, X, SN, S, WNDFAC, ABS, REL,
2861 4 SCRCHK, YV, VISC, COND, D, DKJ, TDR, ICKWRK, 2926 C CALL SGBFA (A, LDA, NATJ*JJ, ML, MU, IP, INFO)
2862 5 RCKWRK, IMCWRK, RMCWRK, F, FN, A, DS, SSAVE, 2927 C*****END precision > single
2863 6 IP, RFAC) 2928 C
2864 C 2929 IF (INFO .NE. 0) THEN
2865 C*****precision > double 2930 WRITE(LOUT,*) 'ERROR IN FACTORING THE JACOBIAN, INFO = ',INFO
2866 IMPLICIT DOUBLE PRECISION (A-H, O-Z), INTEGER (I-N) 2931 STOP
2867 C*****END precision > double 2932 ENDIF
2868 C 2933 C
2869 C*****precision > single 2934 WRITE (LSAVE) ISENSI
2870 C IMPLICIT REAL (A-H, O-Z), INTEGER (I-N) 2935 C
2871 C*****END precision > single 2936 C COMPUTE THE RAW SENSITIVITY COEFFICIENTS WITH RESPECT TO
2872 C 2937 C THE RATE CONSTANTS, D(MASS FRACTION)/D(RATE CONSTANT).
2873 DIMENSION WT(*), YOXID(*), YFUEL(*), XGIVEN(*), TGIVEN(*), X(*), 2938 C
2874 1 SN(NATJ,*), S(NATJ, *), F(NATJ,*), FN(NATJ,*), 2939 DO 1000 I = 1, II
2875 2 SCRCHK(KK,*), YV(KK,*), VISC(*), COND(*), D(KK,*), 2940 C
2876 3 TDR(KK,*), DKJ(KK,KK,*), DS(*), SSAVE(*), IP(*), 2941 CALL CKRDEX (I, RCKWRK, SAVEP)
2877 4 ICKWRK(*), RCKWRK(*), IMCWRK(*), RMCWRK(*), A(LDA, *) 2942 DP = REL*SAVEP + ABS
2878 C 2943 CALL CKRDEX (-I, RCKWRK, SAVEP+DP)
2879 COMMON /LOCS/ NT, NG, NF, NH, NYS, NY 2944 C
2880 C 2945 CALL FUN (KK,JJ,NATJ,LENRGY,LRADI,LNARR,LMULTI,LTDIF, LVARMC,
2881 LOGICAL LENRGY, LRADI,LNARR,LMULTI, LTDIF, LVARMC, LTIME 2946 1 LTIME,PR,WT,YOXID, YFUEL, DT,NTEMP,XGIVEN,TGIVEN,
2882 C 2947 2 X, SN, S, WNDFAC, SCRCHK(1,1), SCRCHK(1,5), YV,
2883 CHARACTER*18 ISENSI 2948 3 SCRCHK(1,2), SCRCHK(1,3), SCRCHK(1,4),
2884 DATA ISENSI /'SENSITIVITY '/ 2949 4 TFUEL, TOXID, VFUEL, VOXID, AFUEL, AOXID,
2885 C 2950 4 VISC, COND, D, DKJ,
2886 C*****precision > single 2951 5 TDR, ICKWRK, RCKWRK, IMCWRK, RMCWRK, F, RFAC)
2887 C DATA ONE/1.0E0/ 2952 C
2888 C*****END precision > single 2953 CALL CKRDEX (-I, RCKWRK, SAVEP)
2889 C*****precision > double 2954 C
2890 DATA ONE/1.0D0/ 2955 DO 100 J = 1, JJ
2891 C*****END precision > double 2956 DO 100 N = 1, NATJ
2892 2957 SN(N,J) = - (F(N,J)-FN(N,J)) / DP
2893 LTIME = .FALSE. 2958 100 CONTINUE
2894 LVARMC = .TRUE. 2959 C
2895 ML = 2*NATJ - 1 2960 C
2896 MU = 2*NATJ - 1 2961 C*****precision > double
2897 C 2962 CALL DGBSL (A, LDA, NATJ*JJ, ML, MU, IP, SN, 0)
2898 C FORM THE JACOBIAN, AND AT THE SAME TIME EVALUATE THE 2963 C*****END precision > double
2899 C UNPERTURBED FUNCTION FN 2964 C
2900 C 2965 C*****precision > single
2901 CALL JACOB (KK, JJ, NATJ,LENRGY ,LRADI,LNARR, LMULTI,LTDIF,LTIME, 2966 C CALL SGBSL (A, LDA, NATJ*JJ, ML, MU, IP, SN, 0)
2902 1 PR, WT, YOXID, YFUEL, DT, NTEMP, XGIVEN, TGIVEN, 2967 C*****END precision > single
2903 2 X, SN, S, WNDFAC, ABS, REL, SCRCHK, YV, 2968 C
2904 3 TFUEL, TOXID, VFUEL, VOXID, AFUEL, AOXID, 2969 C
2905 3 VISC, COND, D, DKJ, 2970 C SN(N,J) NOW CONTAINS THE RAW SENSIVITY MATRIX DS(N,J)/DAi
2906 4 TDR, ICKWRK, RCKWRK, IMCWRK, RMCWRK, F, FN, A, DS, 2971 C
2907 5 SSAVE, RFAC) 2972 C
2908 2973 C NORMALIZE THE SENSIVITY COEFFICIENTS
2909 ! SUBROUTINE JACOB(KK,POINTS,COMPS,LENRGY,LRADI,LNARR,LMULTI,LTDIF, 2974 C
2910 ! 1 LTIME, PR, WT, YOXID, YFUEL, DT, 2975 DO 200 J = 1, JJ
2911 ! 2 NTEMP, XGIVEN, TGIVEN, MESH, SN, X0, WNDFAC, 2976 SN(NT,J) = SN(NT,J) * MAX(ONE,SAVEP) / S(NT,J)
2912 ! 3 ABSOL, RELAT, SCRTCH, YV, 2977 200 CONTINUE
2913 ! 4 TFUEL, TOXID, VFUEL, VOXID, AFUEL, AOXID, 2978 C
2914 ! 4 VISC, COND, D, DKJ, 2979 DO 500 J = 1, JJ
2915 ! 5 TDR, ICKWRK, RCKWRK, IMCWRK, RMCWRK, Y1, Y0, 2980 SUM = 0.E0
2916 ! 6 A, PERTRB, SAVE, RFAC) 2981 DO 300 L = 1, KK
2917 ! 2982 SUM = SUM + SN(NYS+L,J) / WT(L)
2918 C 2983 300 CONTINUE
2919 C FACTOR THE JACOBIAN 2984 CALL CKMMWY (S(NY,J), ICKWRK, RCKWRK, WTM)
2920 C 2985 DO 400 K = 1, KK
2921 C*****precision > double 2986 SN(NYS+K,J) = SN(NYS+K,J) *
2922 CALL DGBFA (A, LDA, NATJ*JJ, ML, MU, IP, INFO) 2987 1 MAX(ONE,SAVEP) / S(NYS+K,J)
2923 C*****END precision > double 2988 2 - MAX(ONE,SAVEP) * WTM * SUM
2924 C 2989 400 CONTINUE
2925 C*****precision > single 2990 500 CONTINUE
2991 C 3056 C
2992 WRITE (LSAVE) I, ((SN(N,J),N=1,NATJ), J=1,JJ) 3057 GASEMI = GASEMI * PRESS/1013250
2993 C 3058 C
2994 1000 CONTINUE 3059 ELSE
2995 C 3060 GASEMI = 0.50
2996 RETURN 3061 ENDIF
2997 END 3062 C
2998 3063 C WRITE (*,8101) PCO,PCH4,PCO2,PH2O
2999 C 3064 C8101 format (1x,'PCO=',f7.3,' PCH4=',f7.3
3000 C------------------------------------------------------------------------------------ 3065 C 1 ,' PCO2=',f7.3,' PH2O=',f7.3)
3001 SUBROUTINE RADEMIS(KK,ICKWRK, RCKWRK,PRESS,TEMP,MASSF,POLYCO, 3066 C
3002 1 GASNUM,NGAS,GASEMI) 3067 C WRITE (*,8102) MOLEF(GASNUM(1)),MOLEF(GASNUM(2)),
3003 C Parameters Illustration: 3068 C 1 MOLEF(GASNUM(3)),MOLEF(GASNUM(4))
3004 C 3069 C8102 format (1x,'XCO=',f8.3,' XCH4=',f8.3
3005 C PR: Pressure (Double) 3070 C 1 ,' XCO2=',f8.3,' XH2O=',f8.3)
3006 C TEMP: Temperature (Double) 3071 C
3007 C MASSF: Double Array of Species Mass Fractions 3072 C WRITE (*,8100) GASEMI
3008 C Must be at leaet KK of Dimension 3073 C8100 format (1x,'GASEMI=',f10.4)
3009 C POLYCO: A 4*7 Double Array, Polymial Coefficients 3074 C
3010 C GASNUM: Index of Radiative Gases in Species List 3075 C WRITE (*,8107) TEMP1
3011 C Integer array that has Dimension of NGAS 3076 C8107 format (1x,'TEMP1=',f10.4)
3012 C NGAS: Number of radiative gases 3077 C
3013 C GASEMI: OUTPUT---The Mixture Emissivity 3078 RETURN
3014 IMPLICIT DOUBLE PRECISION (A-H, O-Z), INTEGER (I-N) 3079 C
3015 C 3080 END
3016 DOUBLE PRECISION MASSF(*), MOLEF(KK) 3081 C
3017 INTEGER GASNUM(*) 3082 C
3018 DIMENSION POLYCO(7,4) 3083
3019 C 3084 C**************************************************************************
3020 C 3085
3021 CALL CKYTX(MASSF,ICKWRK, RCKWRK,MOLEF) 3086 C
3022 C 3087
3023 C 3088 CCC SUBROUTINE RADCAL
3024 TEMP1 = (TEMP/1000)*9/5 3089 CCC
3025 IF( (TEMP1 .GE. 0.0) .AND. (TEMP1 .LE. 10.0)) THEN 3090 CCC IMPLICIT DOUBLE PRECISION (A-H, O-Z), INTEGER (I-N)
3026 C 3091 CCC
3027 C 3092 CCC real :: ommin
3028 PCO=0.0 3093 CCC real :: ommax
3029 DO 1996 I=1,7 3094 CCC
3030 PCO = PCO + POLYCO(I,1)*TEMP1**(I-1) 3095 CCC
3031 1996 CONTINUE 3096 CCC DOUBLE PRECISION SDWEAK,GDINV,GDDINV,XC,AOM,Q,QW(600),TTAU(600),
3032 3097 CCC
3033 C 3098 CCC 2 XTOT(100),XT(600),XSTAR(100),X(4,100),UK,TAU(100)
3034 PCH4=0.0 3099 CCC
3035 DO 1997 I=1,7 3100 CCC COMMON/CMAIN/OMMIN,OMMAX,NOM,TWALL,P(6,100),SPECIE(5),NPRINT
3036 PCH4 = PCH4 + POLYCO(I,2)*TEMP1**(I-1) 3101 CCC
3037 1997 CONTINUE 3102 CCC COMMON/CPARAM/GAMMA(4,7)
3038 C 3103 CCC
3039 3104 CCC COMMON/CPART/W(100),XPART(100),T(100),DD(100),NPT
3040 PCO2=0.0 3105 CCC
3041 DO 1998 I=1,7 3106 CCC COMMON/CLAYER/ QL(100),QLS(100,600)
3042 PCO2 = PCO2 + POLYCO(I,3)*TEMP1**(I-1) 3107 CCC
3043 1998 CONTINUE 3108 CCC DIMENSION U(4,100),AC(100),AD(100),GC(4,100),AMBDA(600),TAUS(100),
3044 C 3109 CCC
3045 3110 CCC 2 AB(600),PKPA(6)
3046 PH2O=0.0 3111 CCC
3047 DO 1999 I=1,7 3112 CCC COMMON/RCDATA/NOL,XLL(100),TL(100),RSC(7,100),WALLT,WAVMAX,WAVMIN
3048 PH2O = PH2O + POLYCO(I,4)*TEMP1**(I-1) 3113 CCC
3049 1999 CONTINUE 3114 CCCC
3050 C 3115 CCCC [NOTE: THE TOTAL INTENSITY CALCULATED IS THAT WHICH LEAVES INTERVAL J=1.
3051 3116 CCCC P(I,J) IS PARTIAL PRESSURE, ATM, OF SPECIES I IN INTERVAL J.
3052 GASEMI = PCO*MOLEF(GASNUM(1))+ 3117 CCCC I=1,2,3,4,5, OR 6 IMPLIES SPECIES IS CO2, H2O, CH4, CO, O2, OR N2, RESP.]
3053 1 PCH4*MOLEF(GASNUM(2))+ 3118 CCCC
3054 2 PCO2*MOLEF(GASNUM(3))+ 3119 CCCC OPEN (7,FILE='RC.DAT')
3055 3 PH2O*MOLEF(GASNUM(4)) 3120 CCCC READ(7,*)NPT
3121 CCC 3186 CCCC REPLACE THE PREVIOUS WITH THE FOLLOWING
3122 CCC NPT = NOL 3187 CCC
3123 CCC 3188 CCC NOM=IFIX((OMMAX-OMMIN)/600.)
3124 CCCC IF(NPT.EQ.0)GO TO 3000 3189 CCC
3125 CCC 3190 CCC GO TO 9906
3126 CCC SPECIE(1)=0. 3191 CCC
3127 CCC SPECIE(2)=0. 3192 CCC9901 NOM=IFIX((OMMAX-OMMIN)/5.)
3128 CCC SPECIE(3)=0. 3193 CCC
3129 CCC SPECIE(4)=0. 3194 CCC GO TO 9906
3130 CCC SPECIE(5)=0. 3195 CCC
3131 CCC 3196 CCC9902 NOM=IFIX((OMMAX-OMMIN)/50.)
3132 CCC DO 9950 J=1,NPT 3197 CCC
3133 CCC 3198 CCC GO TO 9906
3134 CCCC READ(7,*)DD(J),T(J),(P(I,J),I=1,6),W(J) 3199 CCC
3135 CCC 3200 CCC9903 NOM=IFIX((1100.-OMMIN)/5.)+IFIX((5000.-1100.)/25.)
3136 CCC DD(J) = XLL(J) 3201 CCC
3137 CCC 3202 CCC 2 +IFIX((OMMAX-5000.)/50.)
3138 CCC T(J) = TL(J) 3203 CCC
3139 CCC 3204 CCC GO TO 9906
3140 CCC DO 9949 I=1,6 3205 CCC
3141 CCC 3206 CCC9904 NOM=IFIX((1100.-OMMIN)/5.)+IFIX((OMMAX-1100.)/25.)
3142 CCC P(I,J) = RSC(I,J) 3207 CCC
3143 CCC 3208 CCC GO TO 9906
3144 CCC9949 CONTINUE 3209 CCC
3145 CCC 3210 CCC9905 NOM=IFIX((5000.-OMMIN)/25.)+IFIX((OMMAX-5000.)/50.)
3146 CCC W(J) = RSC(7,J) 3211 CCC
3147 CCC 3212 CCC9906 NPRINT=1
3148 CCC DO 9948 I=1,6 3213 CCC
3149 CCC 3214 CCCC
3150 CCC9948 P(I,J)=P(I,J)/101. 3215 CCC DOM=5.0
3151 CCC 3216 CCC
3152 CCC SPECIE(1)=P(1,J)+SPECIE(1) 3217 CCC OMEGA=OMMIN-DOM
3153 CCC 3218 CCC
3154 CCC SPECIE(2)=P(2,J)+SPECIE(2) 3219 CCC NM=NOM-1
3155 CCC 3220 CCC
3156 CCC SPECIE(3)=P(3,J)+SPECIE(3) 3221 CCCC WRITE (*,2999) NOM,OMMIN,OMMAX
3157 CCC 3222 CCC
3158 CCC SPECIE(4)=P(4,J)+SPECIE(4) 3223 CCCC2999 FORMAT ('NOM=',I5,'OMMIN=',I5,'OMMAX=',I10)
3159 CCC 3224 CCCC
3160 CCC SPECIE(5)=W(J)+SPECIE(5) 3225 CCC
3161 CCC 3226 CCCC LOOP 1000 COMPUTES EACH SPECTRAL CONTRIBUTION
3162 CCC9950 CONTINUE 3227 CCC
3163 CCC 3228 CCCC *********************************************
3164 CCCC READ(7,*)TWALL,OMMIN,OMMAX 3229 CCC
3165 CCC TWALL = WALLT 3230 CCC
3166 CCC 3231 CCC
3167 CCC OMMIN = WAVMIN 3232 CCC DO 1000 KK=1,NOM
3168 CCC 3233 CCC
3169 CCC OMMAX = WAVMAX 3234 CCC OMEGA=OMEGA+DOM
3170 CCC 3235 CCC
3171 CCCC CLOSE(7) 3236 CCC IF(OMEGA.LE.1100.)GO TO 109
3172 CCCC THESE STATEMENTS ARE COMMENTED OUT. GET A FINE WAVELENGTH INCREMENT 3237 CCC
3173 CCC 3238 CCC OMEGA=OMEGA+20.
3174 CCC IF(OMMAX.LT.1100.)GO TO 9901 3239 CCC
3175 CCC 3240 CCC IF(OMEGA.LE.5000.)GO TO 109
3176 CCC IF(OMMIN.GT.5000.)GO TO 9902 3241 CCC
3177 CCC 3242 CCC OMEGA=OMEGA+25.
3178 CCC IF(OMMIN.LT.1100..AND.OMMAX.GT.5000.)GO TO 9903 3243 CCC
3179 CCC 3244 CCC109 AMBDA(KK)=10000./OMEGA
3180 CCC IF(OMMIN.LT.1100.)GO TO 9904 3245 CCC
3181 CCC 3246 CCC ABGAS=0.
3182 CCC IF(OMMAX.GT.5000.)GO TO 9905 3247 CCC
3183 CCC 3248 CCCC
3184 CCC NOM=IFIX((OMMAX-OMMIN)/25.) 3249 CCC
3185 CCC 3250 CCCC LOOP 200 COMPUTES THE CONTRIBUTION OF EACH SPECIES TO TAU
4811 CCC 4876 CCC 1 .680E-07, .620E-03, .590E-01, .443E+00, .103E+01, .156E+01, 750
4812 CCC DATA A1/ 4877 CCC
4813 CCC 4878 CCC 1 .385E-07, .275E-03, .450E-01, .330E+00, .845E+00, .136E+01, 775
4814 CCC 1 .950E+00, .103E+00, .420E-01, .114E-01, .450E-02, .300E-02, 50 4879 CCC
4815 CCC 4880 CCC 1 .670E-07, .113E-03, .355E-01, .242E+00, .695E+00, .117E+01, 800
4816 CCC 1 .208E+01, .365E+00, .113E+00, .375E-01, .195E-01, .134E-01, 75 4881 CCC
4817 CCC 4882 CCC 1 .113E-06, .500E-04, .289E-01, .174E+00, .560E+00, .100E+01/ 825
4818 CCC 1 .368E+01, .990E+00, .300E+00, .104E+00, .577E-01, .365E-01, 100 4883 CCC
4819 CCC 4884 CCC DATA A5/
4820 CCC 1 .650E+01, .201E+01, .650E+00, .214E+00, .128E+00, .845E-01, 125 4885 CCC
4821 CCC 4886 CCC 1 .195E-06, .230E-04, .245E-01, .123E+00, .450E+00, .855E+00, 850
4822 CCC 1 .825E+01, .325E+01, .121E+01, .415E+00, .260E+00, .168E+00, 150 4887 CCC
4823 CCC 4888 CCC 1 .328E-06, .103E-04, .214E-01, .100E+00, .357E+00, .718E+00, 875
4824 CCC 1 .870E+01, .452E+01, .189E+01, .765E+00, .450E+00, .289E+00, 175 4889 CCC
4825 CCC 4890 CCC 1 .560E-06, .460E-05, .189E-01, .830E-01, .278E+00, .595E+00, 900
4826 CCC 1 .810E+01, .540E+01, .261E+01, .126E+01, .695E+00, .460E+00, 200 4891 CCC
4827 CCC 4892 CCC 1 .950E-06, .205E-05, .174E-01, .730E-01, .239E+00, .492E+00, 925
4828 CCC 1 .682E+01, .600E+01, .337E+01, .179E+01, .101E+01, .679E+00/ 225 4893 CCC
4829 CCC 4894 CCC 1 .160E-05, .140E-05, .166E-01, .665E-01, .211E+00, .405E+00, 950
4830 CCC DATA A2/ 4895 CCC
4831 CCC 4896 CCC 1 .275E-05, .350E-05, .165E-01, .630E-01, .195E+00, .352E+00, 975
4832 CCC 1 .493E+01, .622E+01, .407E+01, .230E+01, .135E+01, .935E+00, 250 4897 CCC
4833 CCC 4898 CCC 1 .470E-05, .850E-05, .167E-01, .620E-01, .190E+00, .312E+00, 1000
4834 CCC 1 .316E+01, .592E+01, .456E+01, .281E+01, .172E+01, .122E+01, 275 4899 CCC
4835 CCC 4900 CCC 1 .810E-05, .215E-04, .175E-01, .630E-01, .191E+00, .289E+00/ 1025
4836 CCC 1 .199E+01, .528E+01, .479E+01, .328E+01, .213E+01, .149E+01, 300 4901 CCC
4837 CCC 4902 CCC DATA A6/
4838 CCC 1 .113E+01, .450E+01, .484E+01, .361E+01, .249E+01, .179E+01, 325 4903 CCC
4839 CCC 4904 CCC 1 .136E-04, .570E-04, .188E-01, .675E-01, .194E+00, .281E+00, 1050
4840 CCC 1 .585E+00, .370E+01, .471E+01, .383E+01, .284E+01, .208E+01, 350 4905 CCC
4841 CCC 4906 CCC 1 .235E-04, .150E-03, .208E-01, .745E-01, .202E+00, .283E+00, 1075
4842 CCC 1 .293E+00, .289E+01, .443E+01, .394E+01, .312E+01, .237E+01, 375 4907 CCC
4843 CCC 4908 CCC 1 .400E-04, .380E-03, .233E-01, .865E-01, .223E+00, .314E+00, 1100
4844 CCC 1 .138E+00, .205E+01, .400E+01, .396E+01, .330E+01, .260E+01, 400 4909 CCC
4845 CCC 4910 CCC 1 .680E-04, .950E-03, .268E-01, .122E+00, .260E+00, .380E+00, 1125
4846 CCC 1 .620E-01, .143E+01, .347E+01, .388E+01, .341E+01, .280E+01/ 425 4911 CCC
4847 CCC 4912 CCC 1 .120E-03, .245E-02, .343E-01, .176E+00, .328E+00, .461E+00, 1150
4848 CCC DATA A3/ 4913 CCC
4849 CCC 4914 CCC 1 .200E-03, .620E-02, .638E-01, .251E+00, .411E+00, .511E+00, 1175
4850 CCC 1 .255E-01, .950E+00, .292E+01, .370E+01, .345E+01, .295E+01, 450 4915 CCC
4851 CCC 4916 CCC 1 .365E-03, .140E-01, .107E+00, .330E+00, .458E+00, .542E+00, 1200
4852 CCC 1 .940E-02, .610E+00, .236E+01, .343E+01, .342E+01, .304E+01, 475 4917 CCC
4853 CCC 4918 CCC 1 .680E-03, .330E-01, .166E+00, .405E+00, .487E+00, .571E+00/ 1225
4854 CCC 1 .340E-02, .386E+00, .188E+01, .310E+01, .334E+01, .309E+01, 500 4919 CCC
4855 CCC 4920 CCC DATA A7/
4856 CCC 1 .105E-02, .236E+00, .145E+01, .274E+01, .319E+01, .307E+01, 525 4921 CCC
4857 CCC 4922 CCC 1 .130E-02, .635E-01, .244E+00, .459E+00, .535E+00, .557E+00, 1250
4858 CCC 1 .350E-03, .144E+00, .110E+01, .238E+01, .300E+01, .301E+01, 550 4923 CCC
4859 CCC 4924 CCC 1 .250E-02, .123E+00, .341E+00, .477E+00, .502E+00, .562E+00, 1275
4860 CCC 1 .126E-03, .820E-01, .818E+00, .204E+01, .276E+01, .289E+01, 575 4925 CCC
4861 CCC 4926 CCC 1 .500E-02, .212E+00, .407E+00, .547E+00, .531E+00, .514E+00, 1300
4862 CCC 1 .430E-04, .445E-01, .598E+00, .174E+01, .248E+01, .275E+01, 600 4927 CCC
4863 CCC 4928 CCC 1 .103E-01, .285E+00, .489E+00, .592E+00, .497E+00, .486E+00, 1325
4864 CCC 1 .150E-04, .242E-01, .427E+00, .145E+01, .222E+01, .260E+01/ 625 4929 CCC
4865 CCC 4930 CCC 1 .219E-01, .328E+00, .491E+00, .558E+00, .489E+00, .485E+00, 1350
4866 CCC DATA A4/ 4931 CCC
4867 CCC 4932 CCC 1 .485E-01, .345E+00, .505E+00, .521E+00, .477E+00, .484E+00, 1375
4868 CCC 1 .510E-05, .127E-01, .294E+00, .118E+01, .195E+01, .241E+01, 650 4933 CCC
4869 CCC 4934 CCC 1 .114E+00, .361E+00, .538E+00, .563E+00, .503E+00, .502E+00, 1400
4870 CCC 1 .170E-05, .630E-02, .200E+00, .950E+00, .169E+01, .221E+01, 675 4935 CCC
4871 CCC 4936 CCC 1 .249E+00, .460E+00, .621E+00, .624E+00, .538E+00, .538E+00/ 1425
4872 CCC 1 .570E-06, .300E-02, .134E+00, .748E+00, .146E+01, .200E+01, 700 4937 CCC
4873 CCC 4938 CCC DATA A8/
4874 CCC 1 .195E-06, .140E-02, .902E-01, .580E+00, .124E+01, .178E+01, 725 4939 CCC
4875 CCC 4940 CCC 1 .397E+00, .569E+00, .749E+00, .768E+00, .581E+00, .565E+00, 1450
4941 CCC 5006 CCC 1 .480E-03, .298E-02, .931E-02, .189E-01, .329E-01, .403E-01, 2200
4942 CCC 1 .418E+00, .627E+00, .824E+00, .849E+00, .640E+00, .594E+00, 1475 5007 CCC
4943 CCC 5008 CCC 1 .405E-03, .175E-02, .696E-02, .152E-01, .295E-01, .365E-01/ 2225
4944 CCC 1 .108E+01, .125E+01, .113E+01, .940E+00, .807E+00, .663E+00, 1500 5009 CCC
4945 CCC 5010 CCC DATA A12/
4946 CCC 1 .165E+01, .155E+01, .118E+01, .670E+00, .562E+00, .483E+00, 1525 5011 CCC
4947 CCC 5012 CCC 1 .321E-03, .120E-02, .452E-02, .101E-01, .252E-01, .331E-01, 2250
4948 CCC 1 .142E+01, .675E+00, .557E+00, .349E+00, .276E+00, .263E+00, 1550 5013 CCC
4949 CCC 5014 CCC 1 .229E-03, .721E-03, .364E-02, .930E-02, .225E-01, .305E-01, 2275
4950 CCC 1 .451E+00, .202E+00, .132E+00, .118E+00, .134E+00, .156E+00, 1575 5015 CCC
4951 CCC 5016 CCC 1 .195E-03, .544E-03, .318E-02, .750E-02, .202E-01, .284E-01, 2300
4952 CCC 1 .603E-01, .538E-01, .863E-01, .112E+00, .120E+00, .125E+00, 1600 5017 CCC
4953 CCC 5018 CCC 1 .154E-03, .375E-03, .185E-02, .603E-02, .175E-01, .269E-01, 2325
4954 CCC 1 .501E+00, .252E+00, .118E+00, .112E+00, .131E+00, .140E+00/ 1625 5019 CCC
4955 CCC 5020 CCC 1 .101E-03, .263E-03, .119E-02, .480E-02, .156E-01, .253E-01, 2350
4956 CCC DATA A9/ 5021 CCC
4957 CCC 5022 CCC 1 .852E-04, .185E-03, .909E-03, .360E-02, .133E-01, .241E-01, 2375
4958 CCC 1 .730E+00, .430E+00, .237E+00, .191E+00, .171E+00, .170E+00, 1650 5023 CCC
4959 CCC 5024 CCC 1 .763E-04, .137E-03, .711E-03, .316E-02, .122E-01, .237E-01, 2400
4960 CCC 1 .149E+01, .506E+00, .294E+00, .238E+00, .210E+00, .201E+00, 1675 5025 CCC
4961 CCC 5026 CCC 1 .615E-04, .126E-03, .610E-03, .257E-02, .101E-01, .218E-01/ 2425
4962 CCC 1 .100E+01, .553E+00, .434E+00, .340E+00, .260E+00, .220E+00, 1700 5027 CCC
4963 CCC 5028 CCC DATA A13/
4964 CCC 1 .802E+00, .658E+00, .528E+00, .411E+00, .300E+00, .240E+00, 1725 5029 CCC
4965 CCC 5030 CCC 1 .480E-04, .113E-03, .518E-03, .201E-02, .920E-02, .200E-01, 2450
4966 CCC 1 .580E+00, .527E+00, .460E+00, .378E+00, .322E+00, .283E+00, 1750 5031 CCC
4967 CCC 5032 CCC 1 .372E-04, .106E-03, .435E-03, .168E-02, .785E-02, .183E-01, 2475
4968 CCC 1 .330E+00, .403E+00, .430E+00, .356E+00, .318E+00, .270E+00, 1775 5033 CCC
4969 CCC 5034 CCC 1 .355E-04, .101E-03, .376E-03, .168E-02, .669E-02, .166E-01, 2500
4970 CCC 1 .250E+00, .393E+00, .405E+00, .342E+00, .301E+00, .275E+00, 1800 5035 CCC
4971 CCC 5036 CCC 1 .358E-04, .990E-04, .366E-03, .167E-02, .651E-02, .156E-01, 2525
4972 CCC 1 .147E+00, .249E+00, .313E+00, .318E+00, .291E+00, .268E+00/ 1825 5037 CCC
4973 CCC 5038 CCC 1 .389E-04, .102E-03, .376E-03, .167E-02, .641E-02, .152E-01, 2550
4974 CCC DATA A10/ 5039 CCC
4975 CCC 5040 CCC 1 .422E-04, .106E-03, .373E-03, .168E-02, .656E-02, .150E-01, 2575
4976 CCC 1 .910E-01, .252E+00, .298E+00, .295E+00, .269E+00, .253E+00, 1850 5041 CCC
4977 CCC 5042 CCC 1 .521E-04, .111E-03, .371E-03, .170E-02, .673E-02, .152E-01, 2600
4978 CCC 1 .580E-01, .158E+00, .214E+00, .244E+00, .244E+00, .245E+00, 1875 5043 CCC
4979 CCC 5044 CCC 1 .646E-04, .121E-03, .384E-03, .179E-02, .798E-02, .179E-01/ 2625
4980 CCC 1 .370E-01, .113E+00, .184E+00, .218E+00, .214E+00, .218E+00, 1900 5045 CCC
4981 CCC 5046 CCC DATA A14/
4982 CCC 1 .244E-01, .118E+00, .156E+00, .188E+00, .195E+00, .200E+00, 1925 5047 CCC
4983 CCC 5048 CCC 1 .742E-04, .129E-03, .479E-03, .201E-02, .788E-02, .175E-01, 2650
4984 CCC 1 .162E-01, .606E-01, .976E-01, .141E+00, .166E+00, .179E+00, 1950 5049 CCC
4985 CCC 5050 CCC 1 .953E-04, .165E-03, .544E-03, .249E-02, .945E-02, .204E-01, 2675
4986 CCC 1 .112E-01, .425E-01, .903E-01, .133E+00, .148E+00, .156E+00, 1975 5051 CCC
4987 CCC 5052 CCC 1 .101E-03, .190E-03, .761E-03, .324E-02, .106E-01, .231E-01, 2700
4988 CCC 1 .780E-02, .400E-01, .765E-01, .112E+00, .129E+00, .137E+00, 2000 5053 CCC
4989 CCC 5054 CCC 1 .147E-03, .272E-03, .892E-03, .441E-02, .125E-01, .257E-01, 2725
4990 CCC 1 .540E-02, .352E-01, .647E-01, .876E-01, .110E+00, .118E+00/ 2025 5055 CCC
4991 CCC 5056 CCC 1 .195E-03, .326E-03, .100E-02, .499E-02, .147E-01, .295E-01, 2750
4992 CCC DATA A11/ 5057 CCC
4993 CCC 5058 CCC 1 .261E-03, .421E-03, .145E-02, .568E-02, .161E-01, .306E-01, 2775
4994 CCC 1 .380E-02, .252E-01, .507E-01, .705E-01, .888E-01, .100E+00, 2050 5059 CCC
4995 CCC 5060 CCC 1 .305E-03, .515E-03, .195E-02, .754E-02, .185E-01, .363E-01, 2800
4996 CCC 1 .260E-02, .179E-01, .377E-01, .546E-01, .724E-01, .828E-01, 2075 5061 CCC
4997 CCC 5062 CCC 1 .362E-03, .645E-03, .237E-02, .830E-02, .205E-01, .373E-01/ 2825
4998 CCC 1 .180E-02, .123E-01, .294E-01, .443E-01, .608E-01, .686E-01, 2100 5063 CCC
4999 CCC 5064 CCC DATA A15/
5000 CCC 1 .127E-02, .850E-02, .212E-01, .378E-01, .579E-01, .640E-01, 2125 5065 CCC
5001 CCC 5066 CCC 1 .507E-03, .850E-03, .274E-02, .888E-02, .234E-01, .431E-01, 2850
5002 CCC 1 .880E-03, .680E-02, .152E-01, .275E-01, .449E-01, .521E-01, 2150 5067 CCC
5003 CCC 5068 CCC 1 .799E-03, .118E-02, .322E-02, .110E-01, .262E-01, .451E-01, 2875
5004 CCC 1 .620E-02, .400E-02, .107E-01, .214E-01, .374E-01, .453E-01, 2175 5069 CCC
5005 CCC 5070 CCC 1 .935E-03, .160E-02, .386E-02, .126E-01, .292E-01, .530E-01, 2900
5331 CCC 1 .139E+00, .105E+00, .844E-01, .687E-01, .593E-01, .560E-01, 5350 5396 CCC
5332 CCC 5397 CCC 1 .988E-05, .420E-04, .210E-03, .832E-03, .317E-02, .842E-02, 6075
5333 CCC 1 .774E-01, .710E-01, .683E-01, .618E-01, .556E-01, .534E-01, 5375 5398 CCC
5334 CCC 5399 CCC 1 .991E-05, .425E-04, .219E-03, .877E-03, .340E-02, .888E-02, 6100
5335 CCC 1 .858E-01, .483E-01, .579E-01, .547E-01, .503E-01, .495E-01, 5400 5400 CCC
5336 CCC 5401 CCC 1 .102E-04, .435E-04, .231E-03, .937E-03, .361E-02, .929E-02, 6125
5337 CCC 1 .985E-01, .575E-01, .589E-01, .510E-01, .451E-01, .449E-01/ 5425 5402 CCC
5338 CCC 5403 CCC 1 .110E-04, .486E-04, .244E-03, .971E-03, .402E-02, .994E-02, 6150
5339 CCC DATA A28/ 5404 CCC
5340 CCC 5405 CCC 1 .127E-04, .579E-04, .257E-03, .111E-02, .437E-02, .104E-01, 6175
5341 CCC 1 .996E-01, .682E-01, .539E-01, .489E-01, .454E-01, .446E-01, 5450 5406 CCC
5342 CCC 5407 CCC 1 .131E-04, .612E-04, .277E-03, .113E-02, .465E-02, .110E-01, 6200
5343 CCC 1 .680E-01, .680E-01, .548E-01, .495E-01, .460E-01, .458E-01, 5475 5408 CCC
5344 CCC 5409 CCC 1 .150E-04, .783E-04, .353E-03, .116E-02, .510E-02, .116E-01/ 6225
5345 CCC 1 .325E-01, .520E-01, .515E-01, .483E-01, .449E-01, .454E-01, 5500 5410 CCC
5346 CCC 5411 CCC DATA A32/
5347 CCC 1 .150E-01, .350E-01, .451E-01, .464E-01, .452E-01, .449E-01, 5525 5412 CCC
5348 CCC 5413 CCC 1 .178E-04, .922E-04, .394E-03, .157E-02, .555E-02, .123E-01, 6250
5349 CCC 1 .620E-02, .238E-01, .369E-01, .408E-01, .414E-01, .417E-01, 5550 5414 CCC
5350 CCC 5415 CCC 1 .203E-04, .115E-03, .481E-03, .188E-02, .601E-02, .131E-01, 6275
5351 CCC 1 .270E-02, .158E-01, .282E-01, .339E-01, .366E-01, .384E-01, 5575 5416 CCC
5352 CCC 5417 CCC 1 .230E-04, .145E-03, .617E-03, .183E-02, .644E-02, .139E-01, 6300
5353 CCC 1 .113E-02, .101E-01, .203E-01, .263E-01, .303E-01, .333E-01, 5600 5418 CCC
5354 CCC 5419 CCC 1 .280E-04, .187E-03, .723E-03, .202E-02, .686E-02, .146E-01, 6325
5355 CCC 1 .829E-03, .590E-02, .148E-01, .206E-01, .247E-01, .295E-01/ 5625 5420 CCC
5356 CCC 5421 CCC 1 .305E-04, .209E-03, .811E-03, .243E-02, .779E-02, .157E-01, 6350
5357 CCC DATA A29/ 5422 CCC
5358 CCC 5423 CCC 1 .455E-04, .244E-03, .935E-03, .243E-02, .844E-02, .166E-01, 6375
5359 CCC 1 .365E-03, .310E-02, .969E-02, .154E-01, .203E-01, .258E-01, 5650 5424 CCC
5360 CCC 5425 CCC 1 .661E-04, .320E-03, .989E-03, .288E-02, .902E-02, .173E-01, 6400
5361 CCC 1 .240E-03, .130E-02, .589E-02, .112E-01, .164E-01, .222E-01, 5675 5426 CCC
5362 CCC 5427 CCC 1 .723E-04, .397E-03, .122E-02, .359E-02, .100E-01, .184E-01/ 6425
5363 CCC 1 .158E-03, .400E-03, .417E-02, .850E-02, .134E-01, .190E-01, 5700 5428 CCC
5364 CCC 5429 CCC DATA A33/
5365 CCC 1 .103E-03, .262E-03, .208E-02, .594E-02, .109E-01, .162E-01, 5725 5430 CCC
5366 CCC 5431 CCC 1 .847E-04, .481E-03, .143E-02, .429E-02, .108E-01, .192E-01, 6450
5367 CCC 1 .741E-04, .181E-03, .142E-02, .455E-02, .907E-02, .141E-01, 5750 5432 CCC
5368 CCC 5433 CCC 1 .103E-03, .591E-03, .174E-02, .488E-02, .116E-01, .200E-01, 6475
5369 CCC 1 .625E-04, .135E-03, .816E-03, .316E-02, .698E-02, .121E-01, 5775 5434 CCC
5370 CCC 5435 CCC 1 .131E-03, .703E-03, .247E-02, .549E-02, .124E-01, .205E-01, 6500
5371 CCC 1 .499E-04, .111E-03, .624E-03, .230E-02, .551E-02, .102E-01, 5800 5436 CCC
5372 CCC 5437 CCC 1 .165E-03, .872E-03, .265E-02, .641E-02, .131E-01, .211E-01, 6525
5373 CCC 1 .325E-04, .677E-04, .425E-03, .124E-02, .385E-02, .818E-02/ 5825 5438 CCC
5374 CCC 5439 CCC 1 .205E-03, .110E-02, .298E-02, .749E-02, .140E-01, .218E-01, 6550
5375 CCC DATA A30/ 5440 CCC
5376 CCC 5441 CCC 1 .253E-03, .130E-02, .346E-02, .811E-02, .150E-01, .230E-01, 6575
5377 CCC 1 .231E-04, .563E-04, .278E-03, .986E-03, .290E-02, .672E-02, 5850 5442 CCC
5378 CCC 5443 CCC 1 .338E-03, .150E-02, .445E-02, .890E-02, .159E-01, .237E-01, 6600
5379 CCC 1 .165E-04, .481E-04, .247E-03, .944E-03, .253E-02, .612E-02, 5875 5444 CCC
5380 CCC 5445 CCC 1 .437E-03, .170E-02, .491E-02, .107E-01, .170E-01, .245E-01/ 6625
5381 CCC 1 .126E-04, .432E-04, .241E-03, .886E-03, .220E-02, .582E-02, 5900 5446 CCC
5382 CCC 5447 CCC DATA A34/
5383 CCC 1 .118E-04, .420E-04, .235E-03, .847E-03, .209E-02, .571E-02, 5925 5448 CCC
5384 CCC 5449 CCC 1 .581E-03, .190E-02, .537E-02, .116E-01, .179E-01, .254E-01, 6650
5385 CCC 1 .110E-04, .408E-04, .226E-03, .812E-03, .221E-02, .604E-02, 5950 5450 CCC
5386 CCC 5451 CCC 1 .685E-03, .220E-02, .578E-02, .128E-01, .189E-01, .263E-01, 6675
5387 CCC 1 .101E-04, .400E-04, .213E-03, .805E-03, .239E-02, .641E-02, 5975 5452 CCC
5388 CCC 5453 CCC 1 .900E-03, .250E-02, .649E-02, .134E-01, .195E-01, .275E-01, 6700
5389 CCC 1 .983E-05, .395E-04, .186E-03, .801E-03, .247E-02, .691E-02, 6000 5454 CCC
5390 CCC 5455 CCC 1 .121E-02, .280E-02, .722E-02, .142E-01, .202E-01, .281E-01, 6725
5391 CCC 1 .979E-05, .401E-04, .193E-03, .805E-03, .260E-02, .732E-02/ 6025 5456 CCC
5392 CCC 5457 CCC 1 .152E-02, .330E-02, .813E-02, .161E-01, .212E-01, .288E-01, 6750
5393 CCC DATA A31/ 5458 CCC
5394 CCC 5459 CCC 1 .185E-02, .370E-02, .907E-02, .168E-01, .222E-01, .292E-01, 6775
5395 CCC 1 .976E-05, .410E-04, .201E-03, .814E-03, .285E-02, .776E-02, 6050 5460 CCC
5461 CCC 1 .220E-02, .430E-02, .929E-02, .183E-01, .233E-01, .294E-01, 6800 5526 CCC
5462 CCC 5527 CCC 1 .772E-03, .751E-03, .384E-02, .575E-02, .537E-02, .594E-02, 7525
5463 CCC 1 .255E-02, .500E-02, .114E-01, .195E-01, .245E-01, .289E-01/ 6825 5528 CCC
5464 CCC 5529 CCC 1 .491E-03, .600E-03, .301E-02, .453E-02, .380E-02, .434E-02, 7550
5465 CCC DATA A35/ 5530 CCC
5466 CCC 5531 CCC 1 .275E-03, .410E-03, .193E-02, .366E-02, .319E-02, .332E-02, 7575
5467 CCC 1 .290E-02, .580E-02, .167E-01, .215E-01, .260E-01, .291E-01, 6850 5532 CCC
5468 CCC 5533 CCC 1 .185E-01, .280E-03, .131E-02, .232E-02, .247E-02, .256E-02, 7600
5469 CCC 1 .320E-02, .670E-02, .208E-01, .237E-01, .274E-01, .293E-01, 6875 5534 CCC
5470 CCC 5535 CCC 1 .101E-03, .160E-03, .915E-03, .150E-02, .186E-02, .197E-02/ 7625
5471 CCC 1 .360E-02, .880E-02, .220E-01, .253E-01, .282E-01, .300E-01, 6900 5536 CCC
5472 CCC 5537 CCC DATA A39/
5473 CCC 1 .400E-02, .920E-02, .238E-01, .273E-01, .290E-01, .304E-01, 6925 5538 CCC
5474 CCC 5539 CCC 1 .691E-04, .110E-03, .565E-03, .114E-02, .205E-02, .192E-02, 7650
5475 CCC 1 .460E-02, .108E-01, .272E-01, .279E-01, .298E-01, .310E-01, 6950 5540 CCC
5476 CCC 5541 CCC 1 .476E-04, .750E-04, .114E-02, .124E-02, .175E-02, .187E-02, 7675
5477 CCC 1 .530E-02, .128E-01, .304E-01, .292E-01, .297E-01, .312E-01, 6975 5542 CCC
5478 CCC 5543 CCC 1 .305E-04, .590E-04, .529E-03, .114E-02, .160E-02, .185E-02, 7700
5479 CCC 1 .620E-02, .152E-01, .344E-01, .303E-01, .293E-01, .310E-01, 7000 5544 CCC
5480 CCC 5545 CCC 1 .240E-04, .480E-04, .293E-03, .842E-03, .141E-02, .184E-02, 7725
5481 CCC 1 .760E-02, .182E-01, .341E-01, .297E-01, .290E-01, .300E-01/ 7025 5546 CCC
5482 CCC 5547 CCC 1 .170E-04, .360E-04, .122E-03, .435E-03, .124E-02, .182E-02, 7750
5483 CCC DATA A36/ 5548 CCC
5484 CCC 5549 CCC 1 .120E-04, .240E-04, .121E-03, .435E-03, .118E-02, .187E-02, 7775
5485 CCC 1 .980E-02, .222E-01, .398E-01, .318E-01, .291E-01, .294E-01, 7050 5550 CCC
5486 CCC 5551 CCC 1 .810E-05, .170E-04, .103E-03, .439E-03, .126E-02, .192E-02, 7800
5487 CCC 1 .132E-01, .271E-01, .402E-01, .294E-01, .274E-01, .282E-01, 7075 5552 CCC
5488 CCC 5553 CCC 1 .550E-05, .120E-04, .866E-04, .367E-03, .119E-02, .193E-02/ 7825
5489 CCC 1 .190E-01, .335E-01, .421E-01, .286E-01, .262E-01, .269E-01, 7100 5554 CCC
5490 CCC 5555 CCC DATA A40/
5491 CCC 1 .240E-01, .432E-01, .431E-01, .276E-01, .245E-01, .257E-01, 7125 5556 CCC
5492 CCC 5557 CCC 1 .390E-05, .900E-05, .716E-04, .351E-03, .116E-02, .194E-02, 7850
5493 CCC 1 .288E-01, .570E-01, .458E-01, .270E-01, .228E-01, .243E-01, 7150 5558 CCC
5494 CCC 5559 CCC 1 .295E-05, .830E-05, .373E-04, .254E-03, .114E-02, .196E-02, 7875
5495 CCC 1 .323E-01, .740E-01, .449E-01, .261E-01, .214E-01, .221E-01, 7175 5560 CCC
5496 CCC 5561 CCC 1 .230E-05, .800E-05, .465E-04, .298E-03, .117E-02, .201E-02, 7900
5497 CCC 1 .570E-01, .890E-01, .435E-01, .225E-01, .199E-01, .196E-01, 7200 5562 CCC
5498 CCC 5563 CCC 1 .225E-05, .820E-05, .367E-04, .252E-03, .116E-02, .205E-02, 7925
5499 CCC 1 .216E-01, .680E-01, .378E-01, .239E-01, .195E-01, .192E-01/ 7225 5564 CCC
5500 CCC 5565 CCC 1 .220E-05, .840E-05, .371E-04, .268E-03, .127E-02, .211E-02, 7950
5501 CCC DATA A37/ 5566 CCC
5502 CCC 5567 CCC 1 .223E-05, .920E-05, .396E-04, .273E-03, .128E-02, .216E-02, 7975
5503 CCC 1 .126E-01, .475E-01, .364E-01, .238E-01, .197E-01, .192E-01, 7250 5568 CCC
5504 CCC 5569 CCC 1 .235E-05, .103E-04, .415E-04, .263E-03, .121E-02, .221E-02, 8000
5505 CCC 1 .117E-01, .369E-01, .385E-01, .249E-01, .212E-01, .204E-01, 7275 5570 CCC
5506 CCC 5571 CCC 1 .280E-05, .125E-04, .633E-04, .363E-03, .136E-02, .231E-02/ 8025
5507 CCC 1 .140E-01, .370E-01, .419E-01, .272E-01, .228E-01, .213E-01, 7300 5572 CCC
5508 CCC 5573 CCC DATA A41/
5509 CCC 1 .425E-01, .418E-01, .440E-01, .280E-01, .248E-01, .229E-01, 7325 5574 CCC
5510 CCC 5575 CCC 1 .310E-05, .150E-04, .979E-04, .492E-03, .150E-02, .241E-02, 8050
5511 CCC 1 .640E-01, .460E-01, .427E-01, .290E-01, .263E-01, .238E-01, 7350 5576 CCC
5512 CCC 5577 CCC 1 .370E-05, .180E-04, .120E-03, .580E-03, .167E-02, .251E-02, 8075
5513 CCC 1 .385E-01, .385E-01, .374E-01, .259E-01, .235E-01, .224E-01, 7375 5578 CCC
5514 CCC 5579 CCC 1 .420E-05, .200E-04, .987E-04, .509E-03, .171E-02, .257E-02, 8100
5515 CCC 1 .182E-01, .179E-01, .282E-01, .231E-01, .211E-01, .214E-01, 7400 5580 CCC
5516 CCC 5581 CCC 1 .510E-05, .240E-04, .134E-03, .547E-03, .173E-02, .267E-02, 8125
5517 CCC 1 .170E-01, .810E-02, .191E-01, .175E-01, .181E-01, .194E-01/ 7425 5582 CCC
5518 CCC 5583 CCC 1 .600E-05, .270E-04, .121E-03, .534E-03, .172E-02, .274E-02, 8150
5519 CCC DATA A38/ 5584 CCC
5520 CCC 5585 CCC 1 .720E-05, .300E-04, .204E-03, .684E-03, .184E-02, .285E-02, 8175
5521 CCC 1 .161E-01, .370E-02, .105E-01, .127E-01, .152E-01, .171E-01, 7450 5586 CCC
5522 CCC 5587 CCC 1 .820E-05, .330E-04, .276E-03, .819E-03, .199E-02, .297E-02, 8200
5523 CCC 1 .145E-01, .170E-02, .554E-02, .855E-02, .113E-01, .131E-01, 7475 5588 CCC
5524 CCC 5589 CCC 1 .100E-04, .380E-04, .317E-03, .859E-03, .214E-02, .308E-02/ 8225
5525 CCC 1 .175E-02, .140E-02, .385E-02, .595E-02, .803E-02, .945E-02, 7500 5590 CCC
5851 CCC 1 .210E+00, .780E+00, .127E+01, .133E+01, .137E+01, .132E+01, 710 5916 CCC
5852 CCC 5917 CCC DATA B10/
5853 CCC 1 .190E+00, .780E+00, .140E+01, .146E+01, .147E+01, .142E+01, 715 5918 CCC
5854 CCC 5919 CCC 1 .000E+00, .000E+00, .000E+00, .330E-01, .560E-01, .750E-01, 860
5855 CCC 1 .900E+00, .106E+01, .140E+01, .150E+01, .155E+01, .134E+01, 720 5920 CCC
5856 CCC 5921 CCC 1 .000E+00, .000E+00, .000E+00, .300E-01, .530E-01, .750E-01, 865
5857 CCC 1 .720E-01, .300E+00, .800E+00, .100E+01, .115E+01, .126E+01, 725 5922 CCC
5858 CCC 5923 CCC 1 .000E+00, .000E+00, .000E+00, .290E-01, .530E-01, .850E-01, 870
5859 CCC 1 .640E-01, .210E+00, .560E+00, .720E+00, .860E+00, .102E+01, 730 5924 CCC
5860 CCC 5925 CCC 1 .000E+00, .000E+00, .000E+00, .240E-01, .470E-01, .900E-01, 875
5861 CCC 1 .680E-01, .210E+00, .530E+00, .670E+00, .790E+00, .101E+01/ 735 5926 CCC
5862 CCC 5927 CCC 1 .000E+00, .000E+00, .000E+00, .220E-01, .450E-01, .860E-01, 880
5863 CCC DATA B7/ 5928 CCC
5864 CCC 5929 CCC 1 .000E+00, .000E+00, .000E+00, .000E+00, .000E+00, .000E+00,
5865 CCC 1 .690E-01, .210E+00, .540E+00, .690E+00, .820E+00, .910E+00, 740 5930 CCC
5866 CCC 5931 CCC 1 .000E+00, .000E+00, .000E+00, .000E+00, .000E+00, .000E+00,
5867 CCC 1 .330E-01, .140E+00, .390E+00, .530E+00, .690E+00, .770E+00, 745 5932 CCC
5868 CCC 5933 CCC 1 .000E+00, .000E+00, .000E+00, .000E+00, .000E+00, .000E+00/
5869 CCC 1 .230E-01, .780E-01, .270E+00, .410E+00, .560E+00, .890E+00, 750 5934 CCC
5870 CCC 5935 CCCC THE FOLLOWING DATA ARE FOR THE 7.7 MICRON BAND OF CH4
5871 CCC 1 .300E-01, .860E-01, .280E+00, .400E+00, .520E+00, .710E+00, 755 5936 CCC
5872 CCC 5937 CCCC TEMP,K= 290 600 850 WAVE NO.
5873 CCC 1 .175E-01, .620E-01, .225E+00, .335E+00, .450E+00, .660E+00, 760 5938 CCC
5874 CCC 5939 CCC DATA C1/
5875 CCC 1 .105E-01, .450E-01, .180E+00, .280E+00, .380E+00, .600E+00, 765 5940 CCC
5876 CCC 5941 CCC 1 0., 0., 0.,
5877 CCC 1 .450E-02, .300E-01, .148E+00, .240E+00, .345E+00, .570E+00, 770 5942 CCC
5878 CCC 5943 CCC 1 0., 0., 0.03,
5879 CCC 1 .000E+00, .140E-01, .124E+00, .205E+00, .285E+00, .430E+00/ 775 5944 CCC
5880 CCC 5945 CCC 1 0., 0., 0.22,
5881 CCC DATA B8/ 5946 CCC
5882 CCC 5947 CCC 1 0.16, 0.20, 0.47,
5883 CCC 1 .000E+00, .115E-01, .110E+00, .185E+00, .260E+00, .375E+00, 780 5948 CCC
5884 CCC 5949 CCC 1 0.34, 0.34, 0.62,
5885 CCC 1 .000E+00, .135E-01, .840E-01, .140E+00, .205E+00, .335E+00, 785 5950 CCC
5886 CCC 5951 CCC 1 0.69, 0.53, 0.65,
5887 CCC 1 .000E+00, .430E-02, .650E-01, .120E+00, .185E+00, .325E+00, 790 5952 CCC
5888 CCC 5953 CCC 1 1.27, 0.88, 1.09,
5889 CCC 1 .000E+00, .000E+00, .540E-01, .115E+00, .180E+00, .315E+00, 795 5954 CCC
5890 CCC 5955 CCC 1 1.68, 1.38, 0.87/
5891 CCC 1 .000E+00, .000E+00, .440E-01, .950E-01, .150E+00, .270E+00, 800 5956 CCC
5892 CCC 5957 CCC DATA C2/
5893 CCC 1 .000E+00, .000E+00, .360E-01, .790E-01, .125E+00, .205E+00, 805 5958 CCC
5894 CCC 5959 CCC 1 0.55, 0.28, 0.40,
5895 CCC 1 .000E+00, .000E+00, .250E-01, .650E-01, .110E+00, .178E+00, 810 5960 CCC
5896 CCC 5961 CCC 1 1.25, 0.86, 0.93,
5897 CCC 1 .000E+00, .000E+00, .180E-01, .620E-01, .103E+00, .153E+00/ 815 5962 CCC
5898 CCC 5963 CCC 1 0.34, 0.59, 0.75,
5899 CCC DATA B9/ 5964 CCC
5900 CCC 5965 CCC 1 0., 0.13, 0.25,
5901 CCC 1 .000E+00, .000E+00, .320E-01, .580E-01, .860E-01, .147E+00, 820 5966 CCC
5902 CCC 5967 CCC 1 0., 0., 0.06,
5903 CCC 1 .000E+00, .000E+00, .800E-02, .510E-01, .870E-01, .134E+00, 825 5968 CCC
5904 CCC 5969 CCC 1 0., 0., 0.,
5905 CCC 1 .000E+00, .000E+00, .600E-02, .480E-01, .830E-01, .133E+00, 830 5970 CCC
5906 CCC 5971 CCC 1 0., 0., 0.,
5907 CCC 1 .000E+00, .000E+00, .000E+00, .430E-01, .780E-01, .118E+00, 835 5972 CCC
5908 CCC 5973 CCC 1 0., 0., 0./
5909 CCC 1 .000E+00, .000E+00, .000E+00, .420E-01, .700E-01, .108E+00, 840 5974 CCC
5910 CCC 5975 CCCC THE FOLLOWING DATA ARE FOR THE 3.3 MICRON BAND OF CH4
5911 CCC 1 .000E+00, .000E+00, .000E+00, .360E-01, .640E-01, .980E-01, 845 5976 CCC
5912 CCC 5977 CCCC TEMP, K= 290 600 850
5913 CCC 1 .000E+00, .000E+00, .000E+00, .350E-01, .610E-01, .870E-01, 850 5978 CCC
5914 CCC 5979 CCC DATA C3/
5915 CCC 1 .000E+00, .000E+00, .000E+00, .320E-01, .580E-01, .860E-01/ 855 5980 CCC
5981 CCC 1 0., 0., 0.03, 6046 CCC
5982 CCC 6047 CCC 1 0., 0., 0.03,
5983 CCC 1 0., 0., 0.03, 6048 CCC
5984 CCC 6049 CCC 1 0., 0., 0./
5985 CCC 1 0., 0., 0.03, 6050 CCC
5986 CCC 6051 CCC END
5987 CCC 1 0., 0., 0.06, 6052 CCC
5988 CCC 6053
5989 CCC 1 0.03, 0.03, 0.09,
5990 CCC
5991 CCC 1 0.07, 0.07, 0.12,
5992 CCC
5993 CCC 1 0.09, 0.09, 0.12,
5994 CCC
5995 CCC 1 0.14, 0.15, 0.22/
5996 CCC
5997 CCC DATA C4/
5998 CCC
5999 CCC 1 0.18, 0.22, 0.28,
6000 CCC
6001 CCC 1 0.24, 0.31, 0.37,
6002 CCC
6003 CCC 1 0.33, 0.44, 0.47,
6004 CCC
6005 CCC 1 0.45, 0.50, 0.53,
6006 CCC
6007 CCC 1 0.59, 0.62, 0.62,
6008 CCC
6009 CCC 1 0.74, 0.70, 0.68,
6010 CCC
6011 CCC 1 0.91, 0.77, 0.72,
6012 CCC
6013 CCC 1 1.00, 0.81, 0.75/
6014 CCC
6015 CCC DATA C5/
6016 CCC
6017 CCC 1 1.03, 0.84, 0.78,
6018 CCC
6019 CCC 1 1.03, 0.84, 0.78,
6020 CCC
6021 CCC 1 1.00, 0.81, 0.75,
6022 CCC
6023 CCC 1 0.94, 0.77, 0.72,
6024 CCC
6025 CCC 1 0.72, 0.68, 0.68,
6026 CCC
6027 CCC 1 0.52, 0.63, 0.63,
6028 CCC
6029 CCC 1 0.33, 0.50, 0.56,
6030 CCC
6031 CCC 1 0.25, 0.42, 0.50/
6032 CCC
6033 CCC DATA C6/
6034 CCC
6035 CCC 1 0.17, 0.26, 0.37,
6036 CCC
6037 CCC 1 0.08, 0.18, 0.31,
6038 CCC
6039 CCC 1 0.04, 0.11, 0.22,
6040 CCC
6041 CCC 1 0., 0.06, 0.16,
6042 CCC
6043 CCC 1 0., 0.02, 0.12,
6044 CCC
6045 CCC 1 0., 0., 0.06,