c_ ---------------------------------------------------------------------
c_ RCS lines preceded by "c_ "
c_ ---------------------------------------------------------------------
c_
c_ $Source: /www/html/ipsl/OCMIP/phase2/simulations/CFC/boundcond/RCS/read_cfcatm.f,v $
c_ $Revision: 1.2 $
c_ $Date: 1998/07/20 15:33:32 $   ;  $State: Exp $
c_ $Author: jomce $ ;  $Locker:  $
c_
c_ ---------------------------------------------------------------------
c_ $Log: read_cfcatm.f,v $
c_ Revision 1.2  1998/07/20 15:33:32  jomce
c_ Used new data set from Walker et al. (sent by J. Bullister 17/7/1998)
c_ from mid 1931 to mid 1997.  Format of cfc1112.atm has changed.  This
c_ file modified accordingly, and old atm record moved to old/cfc1112.atm.old
c_
c_ Revision 1.1  1998/07/07 15:22:00  orr
c_ Initial revision
c_
c_ ---------------------------------------------------------------------
c_
      PROGRAM read_cfcatm

C     --------------------------------------------------------------------
C     Reads temporal history of atmospheric CFC-11 and CFC-12 (both in pptv)
C
C     year         = year of measurement (mid-year)
C     p11(year,nt) = pcfc11  in northern (1) and southern (2) hemisphere 
C     p12(year,nt) = pcfc12  in northern (1) and southern (2) hemisphere 
C     --------------------------------------------------------------------

c      Jean-Claude Dutay, LSCE/CEA-CNRS, Saclay, France, 10 May 1998
c      - modfied by J. Orr (9 June 1998)

      IMPLICIT NONE 

      INTEGER i, n, nyrbeg, nyrend, nt, iu
      PARAMETER (nyrbeg=31, nyrend=97, nt=2)
      PARAMETER (iu=10)

      REAL year(nyrbeg:nyrend)
      REAL p11(nyrbeg:nyrend,nt), p12(nyrbeg:nyrend,nt)
      CHARACTER*80 filen
C
C     OPEN FILE
C     ---------
      filen='cfc1112.atm'
      OPEN(UNIT=iu,FILE=filen,STATUS='old')
c
c     Skip over 1st six descriptor lines
c     ----------------------------------
      DO i=1,6
        READ(iu,99)
      END do    

C
C     READ FILE
C     ---------
      DO n = nyrbeg, nyrend
        READ(iu,*)
     +      year(n), p11(n,1), p12(n,1)
     +     ,         p11(n,2), p12(n,2)
        WRITE(6,100)
     +      year(n), p11(n,1), p12(n,1)
     +     ,         p11(n,2), p12(n,2)
      END DO 

 99   FORMAT(1x)
 100  FORMAT(f7.2, 4f8.2)

      STOP 
      END