grvdef Subroutine

subroutine grvdef(tjd, loc, pos1, pobs, pos2)

SUBROUTINE GRVDEF VERSION 2. THIS SUBROUTINE COMPUTES THE TOTAL GRAVITATIONAL DEFLECTION OF LIGHT FOR THE OBSERVED OBJECT DUE TO THE MAJOR GRAVITATING BODIES IN THE SOLAR SYSTEM. THIS VERSION IS A DUMMY. NO CORRECTION IS APPLIED.

 TJD    = (NOT USED)
 LOC    = (NOT USED)
 POS1   = POSITION VECTOR OF OBSERVED OBJECT, WITH RESPECT TO
          ORIGIN AT OBSERVER (OR THE GEOCENTER), REFERRED
          TO ICRS AXES, COMPONENTS IN AU (IN)
 POBS   = (NOT USED)
 POS2   = POSITION VECTOR OF OBSERVED OBJECT, WITH RESPECT TO
          ORIGIN AT OBSERVER (OR THE GEOCENTER), REFERRED
          TO ICRS AXES, CORRECTED FOR GRAVITATIONAL
          DEFLECTION, COMPONENTS IN AU (OUT)

NOTE: IN THIS VERSION, POS2 = POS1. THE GRAVITATIONAL DEFLECTION THAT IS NEGLECTED HERE CAN REACH 1.8 ARCSECONDS AT THE LIMB OF THE SUN, BUT IS LESS THAN 10 MILLIARCSECONDS OVER THE AREA OF SKY MORE THAN 45 DEGREES FROM THE SUN. SEE TABLE 3.26.1 ON PAGE 138 OF THE EXPLANATORY SUPPLEMENT TO THE ASTRONOMICAL ALMANAC (1992).

Arguments

Type IntentOptional Attributes Name
double precision :: tjd
integer :: loc
double precision, dimension(3) :: pos1
double precision, dimension(3) :: pobs
double precision, dimension(3) :: pos2

Source Code

subroutine grvdef (tjd,loc,pos1,pobs,pos2)

double precision tjd,pos1,pobs,pos2
dimension pos1(3), pobs(3), pos2(3)

do j = 1, 3
    pos2(j) = pos1(j)
end do

end subroutine grvdef