DIMENSION IYEAR(138),NHEM(138),SHEM(138),GLOBE(138),
     1           NENSO(138),SENSO(138),GENSO(138)
       OPEN(UNIT=5,FILE='ndp022r2.dat')
       WRITE(6,5)
  5     FORMAT(1X,'HEMISPHERIC AND GLOBAL SURFACE TEMPERATURE ANOMALIES',
     1   /,6X,'(RELATIVE TO A 1950-1979 REFERENCE PERIOD)',/,/,38X,
     2   'ENSO-CORRECTED',/,/,6X,'NORTHERN',1X,'SOUTHERN',7X,'NORTHERN',
     3   1X,'SOUTHERN',/,1X,'YEAR',4X,'HEM.',4X,'HEM.',2X,'GLOBAL',4X,
     4   'HEM.',4X,'HEM.',2X,'GLOBAL')
        DO 20 I=1,138
         READ(5,10,END=99) IYEAR(I),NHEM(I),SHEM(I),GLOBE(I),NENSO(I),
     1  SENSO(I),GENSO(I)
 10      FORMAT(1X,I4,6(F6.2))
         WRITE(6,15) IYEAR(I),NHEM(I),SHEM(I),GLOBE(I),NENSO(I),
     1  SENSO(I),GENSO(I)
 15      FORMAT(1X,I4,2X,6(F6.2,2X))
 20     CONTINUE
 99   CONTINUE
      WRITE(6,25)
 25    FORMAT(/,1X,'ALL ANOMALY ESTIMATES EXPRESSED IN DEGREES CELSIUS.
     1MISSING',/,' VALUES DENOTED BY -99.99.')
      CLOSE(UNIT=5)
      CLOSE(UNIT=6)
        STOP
        END