subroutine dpde2dxdy(mjd,dpsi,deps,dX,dY)

C-------------------------------------------------------------------------C
C									  C
C       Conversion of offsets dpsi, deps to dX, dY                        C
C									  C
C-------------------------------------------------------------------------C
C	 								  C
C	References :                                                      C
C		Lieske et al. 1977, A&A 58, pp. 1-16                      C
C               IERS Conventions 2000, chap. 5                            C
C	 	Herring et al. 2002, JGR 107, B4			  C
C	 								  C
C-------------------------------------------------------------------------C
C	 								  C
C	Input :                                                           C
C		mjd : time in mjd                                         C
C		dpsi, deps : offsets in arcsec ref. to IAU2000A           C
C	Output :                                                          C
C               dX, dY : offsets in arcsec ref. to IAU2000A	          C
C	 								  C
C-------------------------------------------------------------------------C
C	 								  C
C	Subroutine interface : none					  C
C	 								  C
C-------------------------------------------------------------------------C
C	 								  C
C	Author : Sebastien Lambert (Sebastien.Lambert@obspm.fr)		  C
C	Last modified : 25/11/02	 				  C
C	 								  C
C-------------------------------------------------------------------------C

	implicit none

C	Input variables

	double precision               :: mjd,dpsi,deps

C 	Output variables

	double precision               :: dX,dY
	
C	Other variables
	
	double precision               :: jc,eps0,psiA,chiA,epsA
	double precision, parameter    :: a2r=4.84813681109535993D-6
        double precision, parameter    :: r2a=2.06264806247096355D+5
	
C	Julian centuries from J2000.0

	jc=(mjd-51544.5D0)/36525.D0

C	Lieske et al. (1977) expressions in arcsec and converted in radians
C	and improved numerical value of psiA from Herring et al. (2002)

	eps0=84381.448D0

        psiA=5038.47875D0*jc
     .      -1.07259D0*jc*jc
     .      -0.001147D0*jc*jc*jc

        chiA=10.5526D0*jc
     .      -2.38064D0*jc*jc
     .      -0.001125D0*jc*jc*jc

        epsA=eps0
     .      -46.8402D0*jc
     .      -0.00059D0*jc*jc
     .      +0.001813D0*jc*jc*jc

        eps0=eps0*a2r
        psiA=psiA*a2r
        chiA=chiA*a2r
        epsA=epsA*a2r
	
C	Chapter 5, equation (23)

	dX=dpsi*dsin(epsA)
     .    +deps*(psiA*dcos(epsA)-chiA)
	dY=deps
     .    -dpsi*dsin(epsA)*(psiA*dcos(epsA)-chiA)
     
	end