C                                                                               
C ****************************************************************              
C * FORTRAN PROGRAM TO READ AND PRINT THE GRIDDED SURFACE AIR    *              
C * TEMPERATURE ANOMALY FILE FOR THE NORTHERN HEMISPHERE         *              
C ****************************************************************              
C                                                                               
      INTEGER I, J, K, DATA(18), LONG(36), LAT(18), YEAR, MONTH                 
      LONG(1)=0                                                                 
      DO 20 I=2,19                                                              
           LONG(I)=LONG(I-1)+10                                                 
   20 CONTINUE                                                                  
      LONG(20)=170                                                              
      DO 30 I=21,36                                                             
           LONG(I)=LONG(I-1)-10                                                 
   30 CONTINUE                                                                  
      LAT(1)=85                                                                 
      DO 40 I=2,18                                                              
           LAT(I)=LAT(I-1)-5                                                    
   40 CONTINUE                                                                  
   45 DO 70 I=1,3                                                               
           READ (5,400,END=800) YEAR, MONTH                                     
           IF (I.EQ.1) THEN                                                     
                WRITE (6,450) YEAR, MONTH                                       
           END IF                                                               
           IF (I.EQ.2) THEN                                                     
                WRITE (6,460) YEAR, MONTH                                       
           END IF                                                               
           IF (I.EQ.3) THEN                                                     
                WRITE (6,470) YEAR, MONTH                                       
           END IF                                                               
           WRITE (6,500) (LAT(K), K=1,18)                                       
           DO 50 J=1,19                                                         
                READ (5,550) (DATA(K), K=1,18)                                  
                WRITE (6,600) LONG(J), (DATA(K), K=1,18)                        
   50      CONTINUE                                                             
           DO 60 J=20,36                                                        
                READ (5,550) (DATA(K), K=1,18)                                  
                WRITE (6,650) LONG(J), (DATA(K), K=1,18)                        
   60      CONTINUE                                                             
           WRITE (6,700)                                                        
   70 CONTINUE                                                                  
      GOTO 45                                                                   
  400 FORMAT (20X,I4,4X,I2)                                                     
  450 FORMAT (1X,'TEMPERATURE ANOMALY DATA',4X,'YEAR = ',I4,                    
     1     2X,'MONTH = ',I2/)                                                   
  460 FORMAT (1X,'NUMBER OF CONTRIBUTING STATIONS',4X,                          
     1     'YEAR = ',I4,2X,'MONTH = ',I2/)                                      
  470 FORMAT (1X,'MEAN INVERSE DISTANCE OF STATIONS',4X,                        
     1     'YEAR = ',I4,2X,'MONTH = ',I2/)                                      
  500 FORMAT (9X,18(3X,I2,'N')/)                                                
  550 FORMAT (18(1X,I5))                                                        
  600 FORMAT (1X,I3,'E',4X,18(1X,I5))                                           
  650 FORMAT (1X,I3,'W',4X,18(1X,I5))                                           
  700 FORMAT (2(/))                                                             
  800 STOP                                                                      
      END