C+
C NAME:
C	AIPS_WTFm
C PURPOSE:
C	freq=IPS Observation Frequency
C	1=327.0MHz: STELab. Nagoya Univ; 2=81.5MHZ: Cambridge Univ.; 3=73.8MHZ: UCSD
C CALLING SEQUENCE:
	function AIPS_WTFmm(fr,ss,zdist) ! Calculate IPS weighting function along los
C INCLUDE:
	include		'sun.h'
	include		'phys.h'
	include		'math.h'
C EXTERNAL TYPE:
	real		IPSFuncQ
C EXTERNAL:
	external	IPSFuncQ		! Spectral function
C COMMON BLOCKS:
	common	/FQIPS/	hachodz,thetadz
!$omp threadprivate(/FQIPS/)
C CALLS:
C	Simpson
C PROCEDURE:
C-
	real    fr 			! frequency in MHz
	real 	ss			! source size in arcseconds
	real	zdist			! distance along the LOS in AU

	hacho = 0.1*PHYS__c/fr		! Radio wavelength (km)
	theta = ss*MATH__RPD/3600	! theta=Source Size(arcsec->radians)
	dz  = zdist*SUN__AU*1.e+8	! Dist along los (AU->km)

	thetadz = theta*dz
	hachodz = hacho*dz

	W = 0
	if (thetadz .ne. 0) W = Simpson(IPSFuncQ,0.0,sqrt(60.0)/thetadz,1.0E-5,13,N)

	AIPS_WTFmm = W

	return
	end