C+
C NAME:
C	Probe
C PURPOSE:
C	Make model line of sight brightnesses from a density modelat solar probe distances
C
C CATEGORY:
C	Data processing
C CALLING SEQUENCE:
C	probe
C
C INPUTS:
C
C OUTPUTS:
C
C FUNCTIONS/SUBROUTINES:
C	Get3DTval
C PROCEDURE:
C MODIFICATION HISTORY:
C	NOV, 1995 B. Jackson (STEL,UCSD)
C-
	program probe

	cmpAU = 1.49598E13							! Cm per AU
	AU1 = 1.0								! One AU
	AUD = 10.0								! Density of electrons at 1 AU
	dd  = 0.05								! Distance increment of spacecraft from the Sun in AU
	dde = 0.001								! Distance increment along the line of sight from the spacecraft in AU
	ALS10 = 1.0								! One S10
	AMSUN = -26.73								! V-light magnitude of the Sun
	TEN15 = 1.0E15								! 10^15 factor for print-out							

	do K=1,20
	  distsc = dd*K
	  frac    = (AU1/distsc)**2
	  sbright = (ALS10*(10.0**((10.-AMSUN)/2.5))/TEN15)*frac		! Brightness of the Sun in terms of 10^15 S10
	  write(*,'(1X,A,F6.3,/,A,F13.4,A)') 'Distance of the spacecraft from the Sun =', distsc,
     &		' Solar brightness in S10 =',sbright,'E+15'
	  do I=10,180
	    elong = float(I)							! Elongation viewed
C	    print *, 'Elongation =',elong
	    tbright = 0.0
	    do L=1,3000

	      delect  = (float(L)*dde - dde/2.0)				! Distance to the electron
	      ebright = thom_wtf(distsc,delect,elong,3)*distsc			! Brightness in S10 units

	      esdist  = sqrt(delect*delect + distsc*distsc - 2.0*delect*distsc*cosd(elong))	! Distance of the electron from the Sun in AU (law of cosines)
C
C	R^-2 spherically-symmetric heliosphere
	      electn  = ((AU1/esdist)**2)*AUD					! Numbers of electrons*cm^-3 at that distance from the Sun	      tbright = tbright + ebright*electn*dde


	      tbright = tbright + electn*ebright*dde
C	      if(elong.eq.90.) print *, distsc, elong, delect, ebright, esdist, electn, voAU, volm, tbright 

	    end do	

	    tratio  = tbright/sbright
	    write(*,'(1X,A,F6.3,F6.1,F13.4,F13.4,A)') 'distsc, elong, tbright, tratio',distsc, elong, tbright, tratio,'E-15'
	  end do
	end do
	       
        end