SUBROUTINE MPDIR(X1,Y1,GDIST,GAZ,MC,LSK,TT,MPDIST,MPAZ,X2,Y2)             
C                                                                               
C        THIS SUBROUTINE WILL SOLVE THE DIRECT PROBLEM OF GEODETIC POSITIONING  
C      ON THE MAPPING PLANE.                                                    
C                                                                               
C        INPUT:                                                                 
C                                                                               
C              X1,Y1 - X,Y COORDINATES OF THE INITIAL POINT.                    
C              GDIST - GEODESIC DISTANCE.                                       
C              GEODETIC AZIMUTH (IN RADIANS).                                   
C              MC - MERIDIAN CONVERENCE (IN RADIANS).                           
C              LSK - LINE SCALE FACTOR.                                         
C              TT - T-T CORRECTION (IN RADIANS).                                
C                                                                               
C        OUTPUT:                                                                
C                                                                               
C              MPDIST - DISTANCE FROM POINT 1 TO 2 ON THE MAPPING PLANE.        
C              MPAZ - AZIMUTH ON THE MAPPING PLANE (IN RADIANS).                
C              X2 Y2 - X,Y,COORDINATES  OF THE OBSERVED POINT.                  
C                                                                               
C                                           WRITTEN BY G.BOWIE, DEC. 1977.      
C                                                                               
      IMPLICIT REAL*8 (A-Z)                                                     
      MPAZ=GAZ-MC-TT                                                            
      MPDIST=GDIST*LSK                                                          
      X2=X1+DSIN(MPAZ)*MPDIST                                                   
      Y2=Y1+DCOS(MPAZ)*MPDIST                                                   
      RETURN                                                                    
      END