      SUBROUTINE PIA(K,TBAR,PNN,PCC,PCN,PHN)
      PARAMETER (NSPECI=46, NTEMPS=14)
      REAL PIANN,PIACC,PIACN,PIAHN
      COMMON /PIAC/ PIANN(NSPECI,NTEMPS),PIACC(NSPECI,NTEMPS)
     & ,PIACN(NSPECI,NTEMPS),PIAHN(NSPECI,NTEMPS),TMIN,TMAX
      TD=(TMAX-TMIN)/(NTEMPS-1)
      IF (TBAR .LT. TMIN) TBAR=TMIN
      IF (TBAR .GT. TMAX) TBAR=TMAX
      DO 100 I=1,NTEMPS
      T=TMIN + (I-1)*TD
      IF (TBAR .GT. T) GO TO 100
      FACTOR= (T-TBAR)/TD
      PNN= PIANN(K,I) - FACTOR*(PIANN(K,I)-PIANN(K,I-1))
      PCC= PIACC(K,I) - FACTOR*(PIACC(K,I)-PIACC(K,I-1))
      PCN= PIACN(K,I) - FACTOR*(PIACN(K,I)-PIACN(K,I-1))
      PHN= PIAHN(K,I) - FACTOR*(PIAHN(K,I)-PIAHN(K,I-1))
      RETURN
 100  CONTINUE
      RETURN
      END
