	function tver2000(mjd,choice)
	
C-------------------------------------------------------------------------C
C									  C
C       Tidal variations in the Earth rotation                            C
C       from IERS Convention 2000 model                                   C
C									  C
C-------------------------------------------------------------------------C
C                                                                         C
C	Input :  time in mjd                                              C
C	 	 choice = 1 --> correction on UT1 in ms			  C
C	 	 choice = 2 --> correction on LOD in ms			  C
C	 	 choice = 3 --> correction on Omega in rad/s              C
C	 	 choice = 4 --> correction on dC/C in rad                 C
C                                                                         C
C       Omega is the mean Earth velocity (7.292115D-5 rad/s)              C
C       C is the axial moment of inertia of the whole Earth               C
C       (8.037D37 kg.m2.s-1)                                              C
C                                                                         C
C	Output : correction                                               C
C                                                                         C
C	NB : to obtain variations free from tidal effects, the            C
C	correction should be substracted from the observed value          C
C	 								  C
C-------------------------------------------------------------------------C
C	 								  C
C	Subroutine interface : none			 		  C
C	 								  C
C-------------------------------------------------------------------------C
C	 								  C
C	Reference : IERS Conv. 2000, chapter 8                            C
C	 								  C
C-------------------------------------------------------------------------C
C	 								  C
C	Author : Sebastien Lambert					  C
C	Last modified : 02/05/03	 				  C
C	 								  C
C-------------------------------------------------------------------------C

	implicit none
	
C	Input variables

	double precision                    :: mjd
	integer                             :: choice

C 	Output variables

	double precision tver2000
	
C	Other variables
	
	integer  i,j
	double precision  jc,phi,a2r,r2a,corr,Om
	double precision  TAB1 (7,62) 
	double precision  TAB2(7,62)
	double precision  TAB3(7,62)
	double precision  ARG(5)
	
C	Mean angular velocity of the Earth in rad/s

	Om=7.292115D-5
	
C	Conversion radians <--> arcseconds

        a2r=4.84813681109535993D-6
        r2a=2.06264806247096355D+5
	
C	Julian centuries from J2000.0

	jc=(mjd-51544.5D0)/36525.D0
	
C	Delaunay arguments from IERS Conventions (1996)

	ARG(1)=134.96340251D0*3600.D0
     .        +1717915923.2178D0*jc
     .        +31.8792D0*jc*jc
     .        +0.051635D0*jc*jc*jc
     .        -0.00024470D0*jc*jc*jc*jc
	ARG(2)=357.52910918D0*3600.D0
     .        +129596581.0481D0*jc
     .        -0.5532D0*jc*jc
     .        +0.000136D0*jc*jc*jc
     .        -0.00001149D0*jc*jc*jc*jc
	ARG(3)=93.27209062D0*3600.D0
     .        +1739527262.8478D0*jc
     .        -12.7512D0*jc*jc
     .        -0.001037D0*jc*jc*jc
     .        +0.00000417D0*jc*jc*jc*jc
	ARG(4)=297.85019547D0*3600.D0
     .        +1602961601.2090D0*jc
     .        -6.3706D0*jc*jc
     .        +0.006593D0*jc*jc*jc
     .        -0.00003169D0*jc*jc*jc*jc
	ARG(5)=125.04455501D0*3600.D0
     .        -6962890.5431D0*jc
     .        +7.4722D0*jc*jc
     .        +0.007702D0*jc*jc*jc
     .        -0.00005939D0*jc*jc*jc*jc
	
	do i=1,5
	   ARG(i)=ARG(i)*a2r
	end do

C-------------------------------------------------------------------------C
C-------------------------------------------------------------------------C
C-------------------------------------------------------------------------C

C	Defraigne & Smits (1999) in 1.D-4 second for UT1

C-------------------------------------------------------------------------C
C-------------------------------------------------------------------------C
C-------------------------------------------------------------------------C

        data ((TAB1(i,j),i=1,7),j=   1,  10) /
     .  1.D0,  0.D0,  2.D0,  2.D0,  2.D0,   -0.020000D0,    0.000000D0,
     .  2.D0,  0.D0,  2.D0,  0.D0,  1.D0,   -0.040000D0,    0.000000D0,
     .  2.D0,  0.D0,  2.D0,  0.D0,  2.D0,   -0.100000D0,    0.000000D0,
     .  0.D0,  0.D0,  2.D0,  2.D0,  1.D0,   -0.050000D0,    0.000000D0,
     .  0.D0,  0.D0,  2.D0,  2.D0,  2.D0,   -0.120000D0,    0.000000D0,
     .  1.D0,  0.D0,  2.D0,  0.D0,  0.D0,   -0.040000D0,    0.000000D0,
     .  1.D0,  0.D0,  2.D0,  0.D0,  1.D0,   -0.410000D0,    0.000000D0,
     .  1.D0,  0.D0,  2.D0,  0.D0,  2.D0,   -1.000000D0,    0.010000D0,
     .  3.D0,  0.D0,  0.D0,  0.D0,  0.D0,   -0.020000D0,    0.000000D0,
     . -1.D0,  0.D0,  2.D0,  2.D0,  1.D0,   -0.080000D0,    0.000000D0 /
        data ((TAB1(i,j),i=1,7),j=  11,  20) /
     . -1.D0,  0.D0,  2.D0,  2.D0,  2.D0,   -0.200000D0,    0.000000D0,
     .  1.D0,  0.D0,  0.D0,  2.D0,  0.D0,   -0.080000D0,    0.000000D0,
     .  2.D0,  0.D0,  2.D0, -2.D0,  2.D0,    0.020000D0,    0.000000D0,
     .  0.D0,  1.D0,  2.D0,  0.D0,  2.D0,    0.030000D0,    0.000000D0,
     .  0.D0,  0.D0,  2.D0,  0.D0,  0.D0,   -0.300000D0,    0.000000D0,
     .  0.D0,  0.D0,  2.D0,  0.D0,  1.D0,   -3.220000D0,    0.020000D0,
     .  0.D0,  0.D0,  2.D0,  0.D0,  2.D0,   -7.790000D0,    0.050000D0,
     .  2.D0,  0.D0,  0.D0,  0.D0, -1.D0,    0.020000D0,    0.000000D0,
     .  2.D0,  0.D0,  0.D0,  0.D0,  0.D0,   -0.340000D0,    0.000000D0,
     .  2.D0,  0.D0,  0.D0,  0.D0,  1.D0,    0.020000D0,    0.000000D0 /
        data ((TAB1(i,j),i=1,7),j=  21,  30) /
     .  0.D0, -1.D0,  2.D0,  0.D0,  2.D0,   -0.020000D0,    0.000000D0,
     .  0.D0,  0.D0,  0.D0,  2.D0, -1.D0,    0.050000D0,    0.000000D0,
     .  0.D0,  0.D0,  0.D0,  2.D0,  0.D0,   -0.740000D0,    0.000000D0,
     .  0.D0,  0.D0,  0.D0,  2.D0,  1.D0,   -0.050000D0,    0.000000D0,
     .  0.D0, -1.D0,  0.D0,  2.D0,  0.D0,   -0.050000D0,    0.000000D0,
     .  1.D0,  0.D0,  2.D0, -2.D0,  1.D0,    0.050000D0,    0.000000D0,
     .  1.D0,  0.D0,  2.D0, -2.D0,  2.D0,    0.100000D0,    0.000000D0,
     .  1.D0,  1.D0,  0.D0,  0.D0,  0.D0,    0.040000D0,    0.000000D0,
     . -1.D0,  0.D0,  2.D0,  0.D0,  0.D0,    0.050000D0,    0.000000D0,
     . -1.D0,  0.D0,  2.D0,  0.D0,  1.D0,    0.180000D0,    0.000000D0 /
        data ((TAB1(i,j),i=1,7),j=  31,  40) /
     . -1.D0,  0.D0,  2.D0,  0.D0,  2.D0,    0.440000D0,    0.000000D0,
     .  1.D0,  0.D0,  0.D0,  0.D0, -1.D0,    0.540000D0,    0.000000D0,
     .  1.D0,  0.D0,  0.D0,  0.D0,  0.D0,   -8.330000D0,    0.060000D0,
     .  1.D0,  0.D0,  0.D0,  0.D0,  1.D0,    0.550000D0,    0.000000D0,
     .  0.D0,  0.D0,  0.D0,  1.D0,  0.D0,    0.050000D0,    0.000000D0,
     .  1.D0, -1.D0,  0.D0,  0.D0,  0.D0,   -0.060000D0,    0.000000D0,
     . -1.D0,  0.D0,  0.D0,  2.D0, -1.D0,    0.120000D0,    0.000000D0,
     . -1.D0,  0.D0,  0.D0,  2.D0,  0.D0,   -1.840000D0,    0.010000D0,
     . -1.D0,  0.D0,  0.D0,  2.D0,  1.D0,    0.130000D0,    0.000000D0,
     .  1.D0,  0.D0, -2.D0,  2.D0, -1.D0,    0.020000D0,    0.000000D0 /
        data ((TAB1(i,j),i=1,7),j=  41,  50) /
     . -1.D0, -1.D0,  0.D0,  2.D0,  0.D0,   -0.090000D0,    0.000000D0,
     .  0.D0,  2.D0,  2.D0, -2.D0,  2.D0,   -0.060000D0,    0.000000D0,
     .  0.D0,  1.D0,  2.D0, -2.D0,  1.D0,    0.030000D0,    0.000000D0,
     .  0.D0,  1.D0,  2.D0, -2.D0,  2.D0,   -1.910000D0,    0.020000D0,
     .  0.D0,  0.D0,  2.D0, -2.D0,  0.D0,    0.260000D0,    0.000000D0,
     .  0.D0,  0.D0,  2.D0, -2.D0,  1.D0,    1.180000D0,   -0.010000D0,
     .  0.D0,  0.D0,  2.D0, -2.D0,  2.D0,  -49.060000D0,    0.430000D0,
     .  0.D0,  2.D0,  0.D0,  0.D0,  0.D0,   -0.200000D0,    0.000000D0,
     .  2.D0,  0.D0,  0.D0, -2.D0, -1.D0,    0.050000D0,    0.000000D0,
     .  2.D0,  0.D0,  0.D0, -2.D0,  0.D0,   -0.560000D0,    0.010000D0 /
        data ((TAB1(i,j),i=1,7),j=  51,  60) /
     .  2.D0,  0.D0,  0.D0, -2.D0,  1.D0,    0.040000D0,    0.000000D0,
     .  0.D0, -1.D0,  2.D0, -2.D0,  1.D0,   -0.050000D0,    0.000000D0,
     .  0.D0,  1.D0,  0.D0,  0.D0, -1.D0,    0.090000D0,    0.000000D0,
     .  0.D0, -1.D0,  2.D0, -2.D0,  2.D0,    0.820000D0,   -0.010000D0,
     .  0.D0,  1.D0,  0.D0,  0.D0,  0.D0,  -15.650000D0,    0.150000D0,
     .  0.D0,  1.D0,  0.D0,  0.D0,  1.D0,   -0.140000D0,    0.000000D0,
     .  1.D0,  0.D0,  0.D0, -1.D0,  0.D0,    0.030000D0,    0.000000D0,
     .  2.D0,  0.D0, -2.D0,  0.D0,  0.D0,   -0.140000D0,    0.000000D0,
     . -2.D0,  0.D0,  2.D0,  0.D0,  1.D0,    0.430000D0,   -0.010000D0,
     . -1.D0,  1.D0,  0.D0,  1.D0,  0.D0,   -0.040000D0,    0.000000D0 /
        data ((TAB1(i,j),i=1,7),j=  61,  62) /
     .  0.D0,  0.D0,  0.D0,  0.D0,  2.D0,    8.200000D0,    0.110000D0,
     .  0.D0,  0.D0,  0.D0,  0.D0,  1.D0,-1689.540000D0,  -25.040000D0 /
	
C-------------------------------------------------------------------------C
C-------------------------------------------------------------------------C
C-------------------------------------------------------------------------C

C-------------------------------------------------------------------------C
C-------------------------------------------------------------------------C
C-------------------------------------------------------------------------C

C	Defraigne & Smits (1999) in 1.D-5 second for LOD

C-------------------------------------------------------------------------C
C-------------------------------------------------------------------------C
C-------------------------------------------------------------------------C

        data ((TAB2(i,j),i=1,7),j=   1,  10) /
     .  1.D0,  0.D0,  2.D0,  2.D0,  2.D0,    0.000000D0,    0.260000D0,
     .  2.D0,  0.D0,  2.D0,  0.D0,  1.D0,    0.000000D0,    0.380000D0,
     .  2.D0,  0.D0,  2.D0,  0.D0,  2.D0,    0.000000D0,    0.910000D0,
     .  0.D0,  0.D0,  2.D0,  2.D0,  1.D0,    0.000000D0,    0.450000D0,
     .  0.D0,  0.D0,  2.D0,  2.D0,  2.D0,    0.010000D0,    1.090000D0,
     .  1.D0,  0.D0,  2.D0,  0.D0,  0.D0,    0.000000D0,    0.270000D0,
     .  1.D0,  0.D0,  2.D0,  0.D0,  1.D0,    0.020000D0,    2.840000D0,
     .  1.D0,  0.D0,  2.D0,  0.D0,  2.D0,    0.040000D0,    6.850000D0,
     .  3.D0,  0.D0,  0.D0,  0.D0,  0.D0,    0.000000D0,    0.120000D0,
     . -1.D0,  0.D0,  2.D0,  2.D0,  1.D0,    0.000000D0,    0.540000D0 /
        data ((TAB2(i,j),i=1,7),j=  11,  20) /
     . -1.D0,  0.D0,  2.D0,  2.D0,  2.D0,    0.010000D0,    1.300000D0,
     .  1.D0,  0.D0,  0.D0,  2.D0,  0.D0,    0.000000D0,    0.500000D0,
     .  2.D0,  0.D0,  2.D0, -2.D0,  2.D0,    0.000000D0,   -0.110000D0,
     .  0.D0,  1.D0,  2.D0,  0.D0,  2.D0,    0.000000D0,   -0.120000D0,
     .  0.D0,  0.D0,  2.D0,  0.D0,  0.D0,    0.010000D0,    1.390000D0,
     .  0.D0,  0.D0,  2.D0,  0.D0,  1.D0,    0.090000D0,   14.860000D0,
     .  0.D0,  0.D0,  2.D0,  0.D0,  2.D0,    0.220000D0,   35.840000D0,
     .  2.D0,  0.D0,  0.D0,  0.D0, -1.D0,    0.000000D0,   -0.100000D0,
     .  2.D0,  0.D0,  0.D0,  0.D0,  0.D0,    0.010000D0,    1.550000D0,
     .  2.D0,  0.D0,  0.D0,  0.D0,  1.D0,    0.000000D0,   -0.080000D0 /
        data ((TAB2(i,j),i=1,7),j=  21,  30) /
     .  0.D0, -1.D0,  2.D0,  0.D0,  2.D0,    0.000000D0,    0.110000D0,
     .  0.D0,  0.D0,  0.D0,  2.D0, -1.D0,    0.000000D0,   -0.200000D0,
     .  0.D0,  0.D0,  0.D0,  2.D0,  0.D0,    0.020000D0,    3.140000D0,
     .  0.D0,  0.D0,  0.D0,  2.D0,  1.D0,    0.000000D0,    0.220000D0,
     .  0.D0, -1.D0,  0.D0,  2.D0,  0.D0,    0.000000D0,    0.210000D0,
     .  1.D0,  0.D0,  2.D0, -2.D0,  1.D0,    0.000000D0,   -0.130000D0,
     .  1.D0,  0.D0,  2.D0, -2.D0,  2.D0,    0.000000D0,   -0.260000D0,
     .  1.D0,  1.D0,  0.D0,  0.D0,  0.D0,    0.000000D0,   -0.100000D0,
     . -1.D0,  0.D0,  2.D0,  0.D0,  0.D0,    0.000000D0,   -0.110000D0,
     . -1.D0,  0.D0,  2.D0,  0.D0,  1.D0,    0.000000D0,   -0.410000D0 /
        data ((TAB2(i,j),i=1,7),j=  31,  40) /
     . -1.D0,  0.D0,  2.D0,  0.D0,  2.D0,   -0.010000D0,   -1.020000D0,
     .  1.D0,  0.D0,  0.D0,  0.D0, -1.D0,   -0.010000D0,   -1.230000D0,
     .  1.D0,  0.D0,  0.D0,  0.D0,  0.D0,    0.130000D0,   18.990000D0,
     .  1.D0,  0.D0,  0.D0,  0.D0,  1.D0,   -0.010000D0,   -1.250000D0,
     .  0.D0,  0.D0,  0.D0,  1.D0,  0.D0,    0.000000D0,   -0.110000D0,
     .  1.D0, -1.D0,  0.D0,  0.D0,  0.D0,    0.000000D0,    0.120000D0,
     . -1.D0,  0.D0,  0.D0,  2.D0, -1.D0,    0.000000D0,   -0.240000D0,
     . -1.D0,  0.D0,  0.D0,  2.D0,  0.D0,    0.020000D0,    3.630000D0,
     . -1.D0,  0.D0,  0.D0,  2.D0,  1.D0,    0.000000D0,   -0.260000D0,
     .  1.D0,  0.D0, -2.D0,  2.D0, -1.D0,    0.000000D0,   -0.040000D0 /
        data ((TAB2(i,j),i=1,7),j=  41,  50) /
     . -1.D0, -1.D0,  0.D0,  2.D0,  0.D0,    0.000000D0,    0.160000D0,
     .  0.D0,  2.D0,  2.D0, -2.D0,  2.D0,    0.000000D0,    0.040000D0,
     .  0.D0,  1.D0,  2.D0, -2.D0,  1.D0,    0.000000D0,   -0.020000D0,
     .  0.D0,  1.D0,  2.D0, -2.D0,  2.D0,    0.010000D0,    0.980000D0,
     .  0.D0,  0.D0,  2.D0, -2.D0,  0.D0,    0.000000D0,   -0.090000D0,
     .  0.D0,  0.D0,  2.D0, -2.D0,  1.D0,    0.000000D0,   -0.420000D0,
     .  0.D0,  0.D0,  2.D0, -2.D0,  2.D0,    0.150000D0,   16.880000D0,
     .  0.D0,  2.D0,  0.D0,  0.D0,  0.D0,    0.000000D0,    0.070000D0,
     .  2.D0,  0.D0,  0.D0, -2.D0, -1.D0,    0.000000D0,   -0.020000D0,
     .  2.D0,  0.D0,  0.D0, -2.D0,  0.D0,    0.000000D0,    0.170000D0 /
        data ((TAB2(i,j),i=1,7),j=  51,  60) /
     .  2.D0,  0.D0,  0.D0, -2.D0,  1.D0,    0.000000D0,   -0.010000D0,
     .  0.D0, -1.D0,  2.D0, -2.D0,  1.D0,    0.000000D0,    0.010000D0,
     .  0.D0,  1.D0,  0.D0,  0.D0, -1.D0,    0.000000D0,   -0.020000D0,
     .  0.D0, -1.D0,  2.D0, -2.D0,  2.D0,    0.000000D0,   -0.140000D0,
     .  0.D0,  1.D0,  0.D0,  0.D0,  0.D0,    0.030000D0,    2.690000D0,
     .  0.D0,  1.D0,  0.D0,  0.D0,  1.D0,    0.000000D0,    0.020000D0,
     .  1.D0,  0.D0,  0.D0, -1.D0,  0.D0,    0.000000D0,    0.000000D0,
     .  2.D0,  0.D0, -2.D0,  0.D0,  0.D0,    0.000000D0,   -0.010000D0,
     . -2.D0,  0.D0,  2.D0,  0.D0,  1.D0,    0.000000D0,   -0.020000D0,
     . -1.D0,  1.D0,  0.D0,  1.D0,  0.D0,    0.000000D0,    0.000000D0 /
        data ((TAB2(i,j),i=1,7),j=  61,  62) /
     .  0.D0,  0.D0,  0.D0,  0.D0,  2.D0,    0.000000D0,    0.150000D0,
     .  0.D0,  0.D0,  0.D0,  0.D0,  1.D0,    0.230000D0,  -15.620000D0 /
	
C-------------------------------------------------------------------------C
C-------------------------------------------------------------------------C
C-------------------------------------------------------------------------C

C-------------------------------------------------------------------------C
C-------------------------------------------------------------------------C
C-------------------------------------------------------------------------C

C	Defraigne & Smits (1999) in 1.D-14 rad/s for Omega

C-------------------------------------------------------------------------C
C-------------------------------------------------------------------------C
C-------------------------------------------------------------------------C

        data ((TAB3(i,j),i=1,7),j=   1,  10) /
     .  1.D0,  0.D0,  2.D0,  2.D0,  2.D0,   -0.220000D0,    0.000000D0,
     .  2.D0,  0.D0,  2.D0,  0.D0,  1.D0,   -0.320000D0,    0.000000D0,
     .  2.D0,  0.D0,  2.D0,  0.D0,  2.D0,   -0.760000D0,    0.000000D0,
     .  0.D0,  0.D0,  2.D0,  2.D0,  1.D0,   -0.380000D0,    0.000000D0,
     .  0.D0,  0.D0,  2.D0,  2.D0,  2.D0,   -0.920000D0,   -0.010000D0,
     .  1.D0,  0.D0,  2.D0,  0.D0,  0.D0,   -0.220000D0,    0.000000D0,
     .  1.D0,  0.D0,  2.D0,  0.D0,  1.D0,   -2.400000D0,   -0.010000D0,
     .  1.D0,  0.D0,  2.D0,  0.D0,  2.D0,   -5.780000D0,   -0.030000D0,
     .  3.D0,  0.D0,  0.D0,  0.D0,  0.D0,   -0.110000D0,    0.000000D0,
     . -1.D0,  0.D0,  2.D0,  2.D0,  1.D0,   -0.460000D0,    0.000000D0 /
        data ((TAB3(i,j),i=1,7),j=  11,  20) /
     . -1.D0,  0.D0,  2.D0,  2.D0,  2.D0,   -1.100000D0,   -0.010000D0,
     .  1.D0,  0.D0,  0.D0,  2.D0,  0.D0,   -0.420000D0,    0.000000D0,
     .  2.D0,  0.D0,  2.D0, -2.D0,  2.D0,    0.090000D0,    0.000000D0,
     .  0.D0,  1.D0,  2.D0,  0.D0,  2.D0,    0.100000D0,    0.000000D0,
     .  0.D0,  0.D0,  2.D0,  0.D0,  0.D0,   -1.170000D0,   -0.010000D0,
     .  0.D0,  0.D0,  2.D0,  0.D0,  1.D0,  -12.540000D0,   -0.080000D0,
     .  0.D0,  0.D0,  2.D0,  0.D0,  2.D0,  -30.250000D0,   -0.180000D0,
     .  2.D0,  0.D0,  0.D0,  0.D0, -1.D0,    0.080000D0,    0.000000D0,
     .  2.D0,  0.D0,  0.D0,  0.D0,  0.D0,   -1.310000D0,   -0.010000D0,
     .  2.D0,  0.D0,  0.D0,  0.D0,  1.D0,    0.070000D0,    0.000000D0 /
        data ((TAB3(i,j),i=1,7),j=  21,  30) /
     .  0.D0, -1.D0,  2.D0,  0.D0,  2.D0,   -0.090000D0,    0.000000D0,
     .  0.D0,  0.D0,  0.D0,  2.D0, -1.D0,    0.170000D0,    0.000000D0,
     .  0.D0,  0.D0,  0.D0,  2.D0,  0.D0,   -2.650000D0,   -0.020000D0,
     .  0.D0,  0.D0,  0.D0,  2.D0,  1.D0,   -0.190000D0,    0.000000D0,
     .  0.D0, -1.D0,  0.D0,  2.D0,  0.D0,   -0.170000D0,    0.000000D0,
     .  1.D0,  0.D0,  2.D0, -2.D0,  1.D0,    0.110000D0,    0.000000D0,
     .  1.D0,  0.D0,  2.D0, -2.D0,  2.D0,    0.220000D0,    0.000000D0,
     .  1.D0,  1.D0,  0.D0,  0.D0,  0.D0,    0.080000D0,    0.000000D0,
     . -1.D0,  0.D0,  2.D0,  0.D0,  0.D0,    0.090000D0,    0.000000D0,
     . -1.D0,  0.D0,  2.D0,  0.D0,  1.D0,    0.350000D0,    0.000000D0 /
        data ((TAB3(i,j),i=1,7),j=  31,  40) /
     . -1.D0,  0.D0,  2.D0,  0.D0,  2.D0,    0.860000D0,    0.010000D0,
     .  1.D0,  0.D0,  0.D0,  0.D0, -1.D0,    1.040000D0,    0.010000D0,
     .  1.D0,  0.D0,  0.D0,  0.D0,  0.D0,  -16.030000D0,   -0.110000D0,
     .  1.D0,  0.D0,  0.D0,  0.D0,  1.D0,    1.050000D0,    0.010000D0,
     .  0.D0,  0.D0,  0.D0,  1.D0,  0.D0,    0.090000D0,    0.000000D0,
     .  1.D0, -1.D0,  0.D0,  0.D0,  0.D0,   -0.100000D0,    0.000000D0,
     . -1.D0,  0.D0,  0.D0,  2.D0, -1.D0,    0.200000D0,    0.000000D0,
     . -1.D0,  0.D0,  0.D0,  2.D0,  0.D0,   -3.070000D0,   -0.020000D0,
     . -1.D0,  0.D0,  0.D0,  2.D0,  1.D0,    0.220000D0,    0.000000D0,
     .  1.D0,  0.D0, -2.D0,  2.D0, -1.D0,    0.030000D0,    0.000000D0 /
        data ((TAB3(i,j),i=1,7),j=  41,  50) /
     . -1.D0, -1.D0,  0.D0,  2.D0,  0.D0,   -0.130000D0,    0.000000D0,
     .  0.D0,  2.D0,  2.D0, -2.D0,  2.D0,   -0.030000D0,    0.000000D0,
     .  0.D0,  1.D0,  2.D0, -2.D0,  1.D0,    0.010000D0,    0.000000D0,
     .  0.D0,  1.D0,  2.D0, -2.D0,  2.D0,   -0.830000D0,   -0.010000D0,
     .  0.D0,  0.D0,  2.D0, -2.D0,  0.D0,    0.080000D0,    0.000000D0,
     .  0.D0,  0.D0,  2.D0, -2.D0,  1.D0,    0.350000D0,    0.000000D0,
     .  0.D0,  0.D0,  2.D0, -2.D0,  2.D0,  -14.250000D0,   -0.130000D0,
     .  0.D0,  2.D0,  0.D0,  0.D0,  0.D0,   -0.060000D0,    0.000000D0,
     .  2.D0,  0.D0,  0.D0, -2.D0, -1.D0,    0.010000D0,    0.000000D0,
     .  2.D0,  0.D0,  0.D0, -2.D0,  0.D0,   -0.140000D0,    0.000000D0 /
        data ((TAB3(i,j),i=1,7),j=  51,  60) /
     .  2.D0,  0.D0,  0.D0, -2.D0,  1.D0,    0.010000D0,    0.000000D0,
     .  0.D0, -1.D0,  2.D0, -2.D0,  1.D0,   -0.010000D0,    0.000000D0,
     .  0.D0,  1.D0,  0.D0,  0.D0, -1.D0,    0.010000D0,    0.000000D0,
     .  0.D0, -1.D0,  2.D0, -2.D0,  2.D0,    0.120000D0,    0.000000D0,
     .  0.D0,  1.D0,  0.D0,  0.D0,  0.D0,   -2.270000D0,   -0.020000D0,
     .  0.D0,  1.D0,  0.D0,  0.D0,  1.D0,   -0.020000D0,    0.000000D0,
     .  1.D0,  0.D0,  0.D0, -1.D0,  0.D0,    0.000000D0,    0.000000D0,
     .  2.D0,  0.D0, -2.D0,  0.D0,  0.D0,    0.010000D0,    0.000000D0,
     . -2.D0,  0.D0,  2.D0,  0.D0,  1.D0,    0.020000D0,    0.000000D0,
     . -1.D0,  1.D0,  0.D0,  1.D0,  0.D0,    0.000000D0,    0.000000D0 /
        data ((TAB3(i,j),i=1,7),j=  61,  62) /
     .  0.D0,  0.D0,  0.D0,  0.D0,  2.D0,   -0.130000D0,    0.000000D0,
     .  0.D0,  0.D0,  0.D0,  0.D0,  1.D0,   13.180000D0,   -0.200000D0 /

C-------------------------------------------------------------------------C
C-------------------------------------------------------------------------C
C-------------------------------------------------------------------------C

C	Implementation for time jc in second

	corr=0.D0
		
	select case (choice)
	
c	Correction to remove from UT1 in ms

	case default
	
	do j=1,62
	   phi=0.D0
	   do i=1,5
	      phi=phi+ARG(i)*TAB1(i,j)
	   end do
	   corr=corr+TAB1(6,j)*dsin(phi)+TAB1(7,j)*dcos(phi)
	end do
	corr=corr*1.D-1
	
c	Correction to remove from LOD in ms

	case (2)
	
	do j=1,62
	   phi=0.D0
	   do i=1,5
	      phi=phi+ARG(i)*TAB2(i,j)
	   end do
	   corr=corr+TAB2(6,j)*dsin(phi)+TAB2(7,j)*dcos(phi)
	end do
	corr=corr*1.D-2
	
	case (3)
	
c	Correction to remove from Om in rad/s

	do j=1,62
	   phi=0.D0
	   do i=1,5
	      phi=phi+ARG(i)*TAB3(i,j)
	   end do
	   corr=corr+TAB3(6,j)*dsin(phi)+TAB3(7,j)*dcos(phi)
	end do
	corr=corr*1.D-14
	
c	Correction to remove from dC/C in rad

	do j=1,62
	   phi=0.D0
	   do i=1,5
	      phi=phi+ARG(i)*TAB3(i,j)
	   end do
	   corr=corr+TAB3(6,j)*dsin(phi)+TAB3(7,j)*dcos(phi)
	end do
	corr=corr*1.D-14
	corr=corr/Om
	
	end select
	
	tver2000=corr
	
	end
