| 1 | SUBROUTINE PIA(K,TBAR,PNN,PCC,PCN,PHN) |
|---|
| 2 | PARAMETER (NSPECI=46, NTEMPS=14) |
|---|
| 3 | REAL PIANN,PIACC,PIACN,PIAHN |
|---|
| 4 | COMMON /PIAC/ PIANN(NSPECI,NTEMPS),PIACC(NSPECI,NTEMPS) |
|---|
| 5 | & ,PIACN(NSPECI,NTEMPS),PIAHN(NSPECI,NTEMPS),TMIN,TMAX |
|---|
| 6 | TD=(TMAX-TMIN)/(NTEMPS-1) |
|---|
| 7 | IF (TBAR .LT. TMIN) TBAR=TMIN |
|---|
| 8 | IF (TBAR .GT. TMAX) TBAR=TMAX |
|---|
| 9 | DO 100 I=1,NTEMPS |
|---|
| 10 | T=TMIN + (I-1)*TD |
|---|
| 11 | IF (TBAR .GT. T) GO TO 100 |
|---|
| 12 | FACTOR= (T-TBAR)/TD |
|---|
| 13 | PNN= PIANN(K,I) - FACTOR*(PIANN(K,I)-PIANN(K,I-1)) |
|---|
| 14 | PCC= PIACC(K,I) - FACTOR*(PIACC(K,I)-PIACC(K,I-1)) |
|---|
| 15 | PCN= PIACN(K,I) - FACTOR*(PIACN(K,I)-PIACN(K,I-1)) |
|---|
| 16 | PHN= PIAHN(K,I) - FACTOR*(PIAHN(K,I)-PIAHN(K,I-1)) |
|---|
| 17 | RETURN |
|---|
| 18 | 100 CONTINUE |
|---|
| 19 | RETURN |
|---|
| 20 | END |
|---|
