1 | SUBROUTINE SETPIA(IPRINT,IFLAG) |
---|
2 | PARAMETER (NSPECI=46,NSPC1I=47, NTEMPS=14) |
---|
3 | REAL PIANN,PIACC,PIACN,PIAHN,T,W |
---|
4 | COMMON /PIAC/ PIANN(NSPECI,NTEMPS),PIACC(NSPECI,NTEMPS) |
---|
5 | & ,PIACN(NSPECI,NTEMPS),PIAHN(NSPECI,NTEMPS),TMIN,TMAX |
---|
6 | COMMON /SPECTI/ BWNI(NSPC1I),WNOI(NSPECI),DWNI(NSPECI) |
---|
7 | & ,WLNI(NSPECI) |
---|
8 | DIMENSION T(NTEMPS),W(NSPECI) |
---|
9 | DATA TMAX,TMIN/190.,60./ |
---|
10 | DO 11 K=1,NSPECI |
---|
11 | DO 11 NT=1,NTEMPS |
---|
12 | PIANN(K,NT)=0.0 |
---|
13 | PIACC(K,NT)=0.0 |
---|
14 | PIACN(K,NT)=0.0 |
---|
15 | PIAHN(K,NT)=0.0 |
---|
16 | 11 CONTINUE |
---|
17 | NWAVES=NSPECI |
---|
18 | DO 99 INU=1,NSPECI |
---|
19 | W(INU)=WNOI(INU) |
---|
20 | IF (WNOI(INU) .LT. 1000.) GOTO 99 |
---|
21 | NWAVES=INU-1 |
---|
22 | GO TO 991 |
---|
23 | 99 CONTINUE |
---|
24 | 991 CONTINUE |
---|
25 | DT=(TMAX-TMIN)/(NTEMPS-1) |
---|
26 | DO 12 I=1,NTEMPS |
---|
27 | T(I)=TMIN + (I-1)*DT |
---|
28 | 12 CONTINUE |
---|
29 | CALL REGIS (W,NSPECI,NWAVES,T,NTEMPS,NTEMPS, |
---|
30 | & PIANN,PIACC,PIACN,PIAHN) |
---|
31 | RETURN |
---|
32 | CC TO BE DELETED... |
---|
33 | 1234 NT=NTEMPS |
---|
34 | NNU=NSPECI |
---|
35 | WRITE(6,101) |
---|
36 | 101 FORMAT(/55X,'N2-N2 ABSORPTION'//) |
---|
37 | WRITE(6,100) (T(IT),IT=1,NT) |
---|
38 | 100 FORMAT(2X,'NU\T',13F9.0/) |
---|
39 | DO 3 INU=1,NNU |
---|
40 | 3 WRITE(6,110) WNOI(INU),(PIANN(INU,IT),IT=1,NT) |
---|
41 | WRITE(6,103) |
---|
42 | 103 FORMAT(/55X,'CH4-CH4 ABSORPTION'//) |
---|
43 | WRITE(6,100) (T(IT),IT=1,NT) |
---|
44 | DO 5 INU=1,NNU |
---|
45 | 5 WRITE(6,110) WNOI(INU),(PIACC(INU,IT),IT=1,NT) |
---|
46 | WRITE(6,105) |
---|
47 | 105 FORMAT(/50X,'N2-CH4 + CH4-N2 ABSORPTION'//) |
---|
48 | DO 7 INU=1,NNU |
---|
49 | 7 WRITE(6,110) WNOI(INU),(PIACN(INU,IT),IT=1,NT) |
---|
50 | WRITE(6,102) |
---|
51 | 102 FORMAT(/55X,'N2-H2 + H2-N2 ABSORPTION'//) |
---|
52 | WRITE(6,100) (T(IT),IT=1,NT) |
---|
53 | DO 4 INU=1,NNU |
---|
54 | 4 WRITE(6,110) WNOI(INU),(PIAHN(INU,IT),IT=1,NT) |
---|
55 | 110 FORMAT(1X,F5.0,2X,13(1PD9.2)) |
---|
56 | RETURN |
---|
57 | END |
---|