C --- -----------------------------------------------------------------
      SUBROUTINE EARTH
      PARAMETER ( PI = 3.141593, D2R = PI/180., HPI = PI/2. )
      INCLUDE 'igm_grid.ins'
C --- -record time series at Earth
      CHARACTER FMT*40
      DATA KE,IE / 12,90 /
      DATA FMT / '(F11.0,F7.1,2F7.2,1P,3E11.3,2E10.2)' /

C --- Get Solar "B-angle" & convert it from Carrington to Theta
      CALL BOANG(Time, Bang)
      ThetaE = HPI - Bang*D2R
C --- Get J such that TH(J) .GE. ThetaE .GT. TH(J-1)
      DO 10 J=2,JM
        IF ( TH(J) .GE. ThetaE ) THEN
	  JEP = J
	  JEM = J-1
	  GOTO 20
	ENDIF
   10 CONTINUE
      STOP "ThetaE not found in grid."

   20 Fp = ( ThetaE - TH(JEM) )/DTH
      Fm = 1. - Fp
      VrE = Fp*Vr(KE,JEP,IE) + Fm*Vr(KE,JEM,IE)
      VtE = Fp*Vt(KE,JEP,IE) + Fm*Vt(KE,JEM,IE)
      VpE = Fp*Vp(KE,JEP,IE) + Fm*Vp(KE,JEM,IE)
      BrE = Fp*Br(KE,JEP,IE) + Fm*Br(KE,JEM,IE)
      BtE = Fp*Bt(KE,JEP,IE) + Fm*Bt(KE,JEM,IE)
      BpE = Fp*Bp(KE,JEP,IE) + Fm*Bp(KE,JEM,IE)
      RoE = Fp*Ro(KE,JEP,IE) + Fm*Ro(KE,JEM,IE)
      TpE = Fp*Tp(KE,JEP,IE) + Fm*Tp(KE,JEM,IE)

      OPEN(11,FILE='Earth.ts',ACCESS='APPEND')
      WRITE(11,FMT) Time, VrE,VtE,VpE, BrE,BtE,BpE, RoE,TpE
      CLOSE(11)

      RETURN
      END