! beam_eff.astro ! Calculate CSO beam efficiency ! S. Radford 2012-06-11, 2013-05-09 ! based on Class macro by D. C. Lis ! CSO receivers are DSB; Data are calibrated for lines with signal in only one SB. ! Planets have signal in both SB, so must halve observed temperature. say "Input parameters:" say "&1: observing frequency (GHz)" say "&2: planet name" say "&3: (ssb) planet antenna temperature (K)" ! say "&4: planet brightness temperature (K)" ! say "&5: beam size (arcs) (<0=diffraction beam)" let obs_freq &1 /new double define character*12 pla_name let pla_name &2 let ssb_temp &3 /new double ! let pla_temp &4 /new double ! let fwhm &5 /new double ! ! Kooi 24.4 @ 345 GHz; diffraction limited 21 => times 1.16 say "Observatory: "'observatory' let frequency obs_freq planet 'pla_name' def double fwhm def double pla_temp def double pla_diam let fwhm astro%planet%beam let pla_temp astro%planet%tb let pla_diam sqrt(astro%planet%maj*astro%planet%min) ! if fwhm.lt.0 then ! let fwhm 1.22*3e10|10.4e2|1e9|obs_freq*206265.!*1.16 ! endif let round int(10*fwhm) /new double say "Assumed beam size: "'round|10'" arcsec" let cfac 1-exp(-log(2.)*(pla_diam|fwhm)**2) /new double let mbef ssb_temp|2.|pla_temp|cfac*100. /new double let round int(mbef*10) say "Beam efficiency : "'round|10'" %"