observing a source

import%20marimo%0A%0A__generated_with%20%3D%20%220.24.0%22%0Aapp%20%3D%20marimo.App(width%3D%22full%22)%0A%0A%0A%40app.cell(hide_code%3DTrue)%0Aasync%20def%20_()%3A%0A%20%20%20%20import%20sys%0A%0A%20%20%20%20%23%20WASM%20has%20no%20site-packages%2C%20so%20divtel%20has%20to%20come%20from%20a%20wheel.%0A%20%20%20%20if%20sys.platform%20%3D%3D%20%22emscripten%22%3A%0A%20%20%20%20%20%20%20%20import%20json%0A%0A%20%20%20%20%20%20%20%20import%20micropip%0A%20%20%20%20%20%20%20%20from%20pyodide.http%20import%20pyfetch%0A%0A%20%20%20%20%20%20%20%20%23%20Pyodide%20runs%20from%20a%20worker%20under%20assets%2F%2C%20so%20step%20up%20one%20level%0A%20%20%20%20%20%20%20%20%23%20to%20reach%20the%20site%20root%20where%20the%20wheel%20is%20published.%0A%20%20%20%20%20%20%20%20base%20%3D%20%22..%2Fpypi%2F%22%0A%20%20%20%20%20%20%20%20%23%20The%20wheel's%20filename%20carries%20a%20version%2C%20so%20read%20it%20from%20a%0A%20%20%20%20%20%20%20%20%23%20manifest%20instead%20of%20hardcoding%20it.%0A%20%20%20%20%20%20%20%20response%20%3D%20await%20pyfetch(base%20%2B%20%22manifest.json%22)%0A%20%20%20%20%20%20%20%20if%20response.status%20!%3D%20200%3A%0A%20%20%20%20%20%20%20%20%20%20%20%20raise%20RuntimeError(%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20f%22could%20not%20fetch%20%7Bbase%7Dmanifest.json%20(HTTP%20%7Bresponse.status%7D).%20%22%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%22The%20divtel%20wheel%20is%20missing%20from%20the%20published%20site.%22%0A%20%20%20%20%20%20%20%20%20%20%20%20)%0A%20%20%20%20%20%20%20%20manifest%20%3D%20json.loads(await%20response.string())%0A%20%20%20%20%20%20%20%20await%20micropip.install(base%20%2B%20manifest%5B%22divtel%22%5D)%0A%20%20%20%20return%0A%0A%0A%40app.cell(hide_code%3DTrue)%0Adef%20_()%3A%0A%20%20%20%20from%20importlib.resources%20import%20files%0A%0A%20%20%20%20import%20astropy.units%20as%20u%0A%20%20%20%20import%20marimo%20as%20mo%0A%20%20%20%20import%20matplotlib.pyplot%20as%20plt%0A%20%20%20%20import%20numpy%20as%20np%0A%20%20%20%20from%20astropy.coordinates%20import%20AltAz%2C%20SkyCoord%2C%20get_body%0A%20%20%20%20from%20astropy.utils%20import%20iers%0A%0A%20%20%20%20from%20divtel.layout%20import%20load_array%0A%20%20%20%20from%20divtel.observation%20import%20Observation%2C%20pointing_coord%0A%20%20%20%20from%20divtel.visualization%20import%20display_hyper_fov%0A%0A%20%20%20%20%23%20Pyodide%20has%20no%20network%2C%20so%20skip%20astropy's%20default%20fetch%20of%20the%20current%0A%20%20%20%20%23%20Earth-orientation%20table.%20The%20bundled%20table%20is%20accurate%20to%20well%20under%0A%20%20%20%20%23%20an%20arcsecond%20near%20the%20release%20date%2C%20plenty%20for%20pointing%20an%20array.%20With%0A%20%20%20%20%23%20no%20download%2C%20astropy%20also%20refuses%20to%20extrapolate%20more%20than%2030%20days%20past%0A%20%20%20%20%23%20the%20bundled%20table's%20predictive%20range%20and%20raises%20instead%3B%20disable%20that%0A%20%20%20%20%23%20check%20too%3B%20a%20dropdown%20offering%20nights%20a%20year%20out%20would%20otherwise%20crash.%0A%20%20%20%20iers.conf.auto_download%20%3D%20False%0A%20%20%20%20iers.conf.auto_max_age%20%3D%20None%0A%0A%20%20%20%20plt.rcParams%5B%22savefig.format%22%5D%20%3D%20%22svg%22%0A%20%20%20%20return%20(%0A%20%20%20%20%20%20%20%20AltAz%2C%0A%20%20%20%20%20%20%20%20Observation%2C%0A%20%20%20%20%20%20%20%20SkyCoord%2C%0A%20%20%20%20%20%20%20%20display_hyper_fov%2C%0A%20%20%20%20%20%20%20%20files%2C%0A%20%20%20%20%20%20%20%20get_body%2C%0A%20%20%20%20%20%20%20%20load_array%2C%0A%20%20%20%20%20%20%20%20mo%2C%0A%20%20%20%20%20%20%20%20np%2C%0A%20%20%20%20%20%20%20%20plt%2C%0A%20%20%20%20%20%20%20%20u%2C%0A%20%20%20%20)%0A%0A%0A%40app.cell(hide_code%3DTrue)%0Adef%20_(mo)%3A%0A%20%20%20%20mo.Html(%0A%20%20%20%20%20%20%20%20%22%22%22%3Cstyle%3E%0A%20%20%20%20%20%20%20%20%20%20img%2C%20svg%20%7B%20max-width%3A%20100%25%3B%20height%3A%20auto%3B%20%7D%0A%20%20%20%20%20%20%20%20%3C%2Fstyle%3E%22%22%22%0A%20%20%20%20)%0A%20%20%20%20return%0A%0A%0A%40app.cell(hide_code%3DTrue)%0Adef%20_(mo)%3A%0A%20%20%20%20mo.md(r%22%22%22%0A%20%20%20%20%23%20Observing%20a%20real%20source%0A%0A%20%20%20%20The%20geometry%20so%20far%20lives%20in%20a%20frame%20fixed%20to%20the%20ground%3A%20an%20array%0A%20%20%20%20pointed%20at%2070%20degrees%20altitude%20stays%20there%20whatever%20the%20hour.%20A%0A%20%20%20%20source%20doesn't.%20It%20rises%2C%20crosses%20the%20sky%2C%20and%20sets%2C%20so%20the%20direction%0A%20%20%20%20you%20point%20in%20changes%20all%20night.%0A%0A%20%20%20%20%60Observation%60%20bridges%20the%20two%3A%20give%20it%20a%20place%20and%20a%20time%2C%20and%20it%0A%20%20%20%20turns%20a%20sky%20position%20into%20the%20altitude%20and%20azimuth%0A%20%20%20%20%60divergent_pointing%60%20wants.%0A%0A%20%20%20%20Pick%20a%20target%20and%20a%20night%3A%0A%20%20%20%20%22%22%22)%0A%20%20%20%20return%0A%0A%0A%40app.cell(hide_code%3DTrue)%0Adef%20_(mo)%3A%0A%20%20%20%20%23%20Hardcoded%20rather%20than%20resolved%20by%20name%2C%20since%20this%20runs%20in%20a%20browser%0A%20%20%20%20%23%20sandbox%20with%20no%20network.%20All%20ICRS%2C%20in%20degrees.%0A%20%20%20%20SOURCES%20%3D%20%7B%0A%20%20%20%20%20%20%20%20%22Crab%20Nebula%22%3A%20(83.633%2C%2022.015)%2C%0A%20%20%20%20%20%20%20%20%22Markarian%20421%22%3A%20(166.114%2C%2038.209)%2C%0A%20%20%20%20%20%20%20%20%22Markarian%20501%22%3A%20(253.468%2C%2039.760)%2C%0A%20%20%20%20%20%20%20%20%22Galactic%20centre%22%3A%20(266.417%2C%20-29.008)%2C%0A%20%20%20%20%20%20%20%20%22PKS%202155-304%22%3A%20(329.717%2C%20-30.226)%2C%0A%20%20%20%20%20%20%20%20%22Centaurus%20A%22%3A%20(201.365%2C%20-43.019)%2C%0A%20%20%20%20%7D%0A%0A%20%20%20%20source%20%3D%20mo.ui.dropdown(%0A%20%20%20%20%20%20%20%20list(SOURCES)%2C%20value%3D%22Crab%20Nebula%22%2C%20label%3D%22source%22%2C%20full_width%3DTrue%0A%20%20%20%20)%0A%20%20%20%20site%20%3D%20mo.ui.dropdown(%0A%20%20%20%20%20%20%20%20%7B%22CTAO-North%20(La%20Palma)%22%3A%20%22north%22%2C%20%22CTAO-South%20(Paranal)%22%3A%20%22south%22%7D%2C%0A%20%20%20%20%20%20%20%20value%3D%22CTAO-North%20(La%20Palma)%22%2C%20label%3D%22site%22%2C%20full_width%3DTrue%2C%0A%20%20%20%20)%0A%20%20%20%20night%20%3D%20mo.ui.dropdown(%0A%20%20%20%20%20%20%20%20%5B%222026-03-01%22%2C%20%222026-06-01%22%2C%20%222026-09-01%22%2C%20%222026-12-01%22%5D%2C%0A%20%20%20%20%20%20%20%20value%3D%222026-12-01%22%2C%20label%3D%22night%20of%22%2C%20full_width%3DTrue%2C%0A%20%20%20%20)%0A%20%20%20%20mo.hstack(%5Bsource%2C%20site%2C%20night%5D%2C%20widths%3D%22equal%22)%0A%20%20%20%20return%20SOURCES%2C%20night%2C%20site%2C%20source%0A%0A%0A%40app.cell(hide_code%3DTrue)%0Adef%20_(%0A%20%20%20%20AltAz%2C%0A%20%20%20%20Observation%2C%0A%20%20%20%20SOURCES%2C%0A%20%20%20%20SkyCoord%2C%0A%20%20%20%20get_body%2C%0A%20%20%20%20night%2C%0A%20%20%20%20np%2C%0A%20%20%20%20site%2C%0A%20%20%20%20source%2C%0A%20%20%20%20u%2C%0A)%3A%0A%20%20%20%20TARGET%20%3D%20SkyCoord(%0A%20%20%20%20%20%20%20%20ra%3DSOURCES%5Bsource.value%5D%5B0%5D%20*%20u.deg%2C%20dec%3DSOURCES%5Bsource.value%5D%5B1%5D%20*%20u.deg%0A%20%20%20%20)%0A%0A%20%20%20%20def%20_solar_midnight()%3A%0A%20%20%20%20%20%20%20%20%22%22%22%0A%20%20%20%20%20%20%20%20The%20middle%20of%20the%20night%2C%20not%20midnight%20UTC.%0A%0A%20%20%20%20%20%20%20%20La%20Palma%20sits%201.2%20hours%20west%20of%20Greenwich%2C%20Paranal%204.7%2C%20and%20the%0A%20%20%20%20%20%20%20%20sun's%20lowest%20point%20drifts%20by%20up%20to%20a%20further%20quarter%20hour%20over%20the%0A%20%20%20%20%20%20%20%20year.%20Centring%20the%20plot%20on%2000%3A00%20UTC%20would%20push%20the%20dark%20band%20off%20to%0A%20%20%20%20%20%20%20%20one%20side%2C%20so%20one%20coarse%20ephemeris%20pass%20finds%20where%20the%20sun%20is%0A%20%20%20%20%20%20%20%20actually%20lowest%20instead.%0A%20%20%20%20%20%20%20%20%22%22%22%0A%20%20%20%20%20%20%20%20start%20%3D%20Observation(site%3Dsite.value%2C%20time%3Df%22%7Bnight.value%7DT00%3A00%3A00%22)%0A%20%20%20%20%20%20%20%20coarse%20%3D%20start.time%20%2B%20np.linspace(-12%2C%2012%2C%2049)%20*%20u.hour%0A%20%20%20%20%20%20%20%20sun%20%3D%20get_body(%22sun%22%2C%20coarse%2C%20location%3Dstart.location).transform_to(%0A%20%20%20%20%20%20%20%20%20%20%20%20AltAz(obstime%3Dcoarse%2C%20location%3Dstart.location)%0A%20%20%20%20%20%20%20%20)%0A%20%20%20%20%20%20%20%20return%20Observation(site%3Dsite.value%2C%20time%3Dcoarse%5Bsun.alt.argmin()%5D)%0A%0A%20%20%20%20MIDNIGHT%20%3D%20_solar_midnight()%0A%20%20%20%20return%20MIDNIGHT%2C%20TARGET%0A%0A%0A%40app.cell(hide_code%3DTrue)%0Adef%20_(AltAz%2C%20MIDNIGHT%2C%20TARGET%2C%20get_body%2C%20np%2C%20u)%3A%0A%20%20%20%20HOURS%20%3D%20np.linspace(-12%2C%2012%2C%2097)%20*%20u.hour%0A%20%20%20%20TIMES%20%3D%20MIDNIGHT.time%20%2B%20HOURS%0A%20%20%20%20FRAME%20%3D%20AltAz(obstime%3DTIMES%2C%20location%3DMIDNIGHT.location)%0A%0A%20%20%20%20SUN%20%3D%20get_body(%22sun%22%2C%20TIMES%2C%20location%3DMIDNIGHT.location).transform_to(FRAME)%0A%20%20%20%20MOON%20%3D%20get_body(%22moon%22%2C%20TIMES%2C%20location%3DMIDNIGHT.location).transform_to(FRAME)%0A%20%20%20%20TRACK%20%3D%20TARGET.transform_to(FRAME)%0A%0A%20%20%20%20%23%20Astronomical%20twilight%3A%20the%20sun%2018%20degrees%20below%20the%20horizon%20is%20the%0A%20%20%20%20%23%20conventional%20start%20of%20usable%20dark%20time%20for%20Cherenkov%20astronomy.%0A%20%20%20%20DARK%20%3D%20SUN.alt%20%3C%20-18%20*%20u.deg%0A%20%20%20%20return%20DARK%2C%20HOURS%2C%20MOON%2C%20SUN%2C%20TRACK%0A%0A%0A%40app.cell(hide_code%3DTrue)%0Adef%20_(DARK%2C%20HOURS%2C%20MOON%2C%20SUN%2C%20TRACK%2C%20np%2C%20plt%2C%20source%2C%20u)%3A%0A%20%20%20%20def%20_navigation_plot()%3A%0A%20%20%20%20%20%20%20%20hours%20%3D%20HOURS.to_value(u.hour)%0A%20%20%20%20%20%20%20%20fig%2C%20ax%20%3D%20plt.subplots(figsize%3D(10%2C%205))%0A%0A%20%20%20%20%20%20%20%20%23%20Shade%20twilight%20and%20full%20dark%20behind%20everything%20else.%0A%20%20%20%20%20%20%20%20ax.fill_between(hours%2C%20-90%2C%2090%2C%20where%3DSUN.alt%20%3C%200%20*%20u.deg%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20color%3D%220.75%22%2C%20zorder%3D0)%0A%20%20%20%20%20%20%20%20ax.fill_between(hours%2C%20-90%2C%2090%2C%20where%3DDARK%2C%20color%3D%220.45%22%2C%20zorder%3D0)%0A%0A%20%20%20%20%20%20%20%20ax.plot(hours%2C%20SUN.alt.to_value(u.deg)%2C%20color%3D%22tab%3Aorange%22%2C%20label%3D%22Sun%22)%0A%20%20%20%20%20%20%20%20ax.plot(hours%2C%20MOON.alt.to_value(u.deg)%2C%20color%3D%220.3%22%2C%20ls%3D%22--%22%2C%20label%3D%22Moon%22)%0A%20%20%20%20%20%20%20%20scatter%20%3D%20ax.scatter(hours%2C%20TRACK.alt.to_value(u.deg)%2C%20c%3DTRACK.az.to_value(u.deg)%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20cmap%3D%22twilight%22%2C%20s%3D12%2C%20zorder%3D3%2C%20label%3Dsource.value)%0A%20%20%20%20%20%20%20%20fig.colorbar(scatter%2C%20ax%3Dax%2C%20label%3D%22source%20azimuth%20%5Bdeg%5D%22)%0A%0A%20%20%20%20%20%20%20%20ax.axhline(0%2C%20color%3D%22k%22%2C%20lw%3D0.8)%0A%20%20%20%20%20%20%20%20ax.set_xlim(-12%2C%2012)%0A%20%20%20%20%20%20%20%20ax.set_ylim(-90%2C%2090)%0A%20%20%20%20%20%20%20%20ax.set_xticks(np.arange(-12%2C%2013%2C%203))%0A%20%20%20%20%20%20%20%20ax.set_xlabel(%22hours%20from%20the%20middle%20of%20the%20night%22)%0A%20%20%20%20%20%20%20%20ax.set_ylabel(%22altitude%20%5Bdeg%5D%22)%0A%20%20%20%20%20%20%20%20ax.set_title(%22what%20is%20up%2C%20and%20when%22)%0A%20%20%20%20%20%20%20%20ax.legend(loc%3D%22upper%20left%22%2C%20frameon%3DFalse)%0A%20%20%20%20%20%20%20%20fig.tight_layout()%0A%20%20%20%20%20%20%20%20return%20fig%0A%0A%20%20%20%20_navigation_plot()%0A%20%20%20%20return%0A%0A%0A%40app.cell(hide_code%3DTrue)%0Adef%20_(mo)%3A%0A%20%20%20%20mo.md(r%22%22%22%0A%20%20%20%20The%20plot%20centres%20on%20the%20middle%20of%20the%20night%2C%20not%20midnight%20UTC%20at%0A%20%20%20%20either%20site.%20The%20dark%20band%20is%20astronomical%20night%2C%20sun%20more%20than%2018%0A%20%20%20%20degrees%20down%3B%20the%20light%20band%20is%20twilight.%20The%20source%20track%20is%0A%20%20%20%20coloured%20by%20azimuth%2C%20so%20you%20can%20see%20it%20swing%20across%20the%20sky%20as%20well%0A%20%20%20%20as%20rise%20and%20fall.%0A%0A%20%20%20%20A%20source%20is%20worth%20observing%20when%20it's%20high%20*and*%20the%20sun%20is%20down%2C%0A%20%20%20%20and%20those%20two%20don't%20always%20overlap%3A%20try%20the%20Galactic%20centre%20from%0A%20%20%20%20CTAO-North%2C%20which%20barely%20clears%20the%20horizon%2C%20against%20the%20same%0A%20%20%20%20source%20from%20CTAO-South.%0A%0A%20%20%20%20%23%23%20What%20the%20array%20sees%20while%20it%20tracks%0A%0A%20%20%20%20Now%20point%20the%20array%20at%20it%2C%20with%20%60div%60%20held%20fixed.%20The%20only%20thing%0A%20%20%20%20changing%20is%20where%20the%20source%20has%20got%20to.%0A%20%20%20%20%22%22%22)%0A%20%20%20%20return%0A%0A%0A%40app.cell(hide_code%3DTrue)%0Adef%20_(files%2C%20load_array)%3A%0A%20%20%20%20ARRAY%20%3D%20load_array(%0A%20%20%20%20%20%20%20%20files(%22divtel%22)%20%2F%20%22data%22%20%2F%20%22cta-north-lapalma-alpha-prod6.ecsv%22%0A%20%20%20%20)%0A%20%20%20%20return%20(ARRAY%2C)%0A%0A%0A%40app.cell(hide_code%3DTrue)%0Adef%20_(mo)%3A%0A%20%20%20%20div%20%3D%20mo.ui.slider(%0A%20%20%20%20%20%20%20%200%2C%200.1%2C%20step%3D0.005%2C%20value%3D0.04%2C%20label%3D%22divergence%22%2C%20show_value%3DTrue%2C%0A%20%20%20%20%20%20%20%20full_width%3DTrue%2C%0A%20%20%20%20)%0A%20%20%20%20hour%20%3D%20mo.ui.slider(%0A%20%20%20%20%20%20%20%20-6%2C%206%2C%20step%3D0.5%2C%20value%3D0%2C%20label%3D%22hours%20from%20the%20middle%20of%20the%20night%22%2C%20show_value%3DTrue%2C%0A%20%20%20%20%20%20%20%20full_width%3DTrue%2C%0A%20%20%20%20)%0A%20%20%20%20mo.vstack(%5Bdiv%2C%20hour%5D)%0A%20%20%20%20return%20div%2C%20hour%0A%0A%0A%40app.cell(hide_code%3DTrue)%0Adef%20_(ARRAY%2C%20MIDNIGHT%2C%20TARGET%2C%20display_hyper_fov%2C%20div%2C%20hour%2C%20mo%2C%20plt%2C%20u)%3A%0A%20%20%20%20def%20_snapshot()%3A%0A%20%20%20%20%20%20%20%20moment%20%3D%20MIDNIGHT.after(hour.value%20*%20u.hour)%0A%20%20%20%20%20%20%20%20alt%2C%20az%20%3D%20moment.altaz_of(TARGET)%0A%0A%20%20%20%20%20%20%20%20if%20alt%20%3C%200%20*%20u.deg%3A%0A%20%20%20%20%20%20%20%20%20%20%20%20return%20mo.md(%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20f%22%23%23%23%20The%20source%20is%20below%20the%20horizon%20at%20%22%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20f%22%7Bmoment.time.isot%5B11%3A16%5D%7D%20UTC%20(altitude%20%7Balt%3A.1f%7D).%5Cn%5Cn%22%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%22Move%20the%20time%20slider%2C%20or%20pick%20a%20source%20that%20is%20up%20on%20this%20night.%22%0A%20%20%20%20%20%20%20%20%20%20%20%20)%0A%0A%20%20%20%20%20%20%20%20ARRAY.divergent_pointing(div.value%2C%20alt%2C%20az)%0A%20%20%20%20%20%20%20%20area%20%3D%20ARRAY.hyper_fov()%5B0%5D%0A%20%20%20%20%20%20%20%20mean%2C%20_%20%3D%20ARRAY.multiplicity_moments()%0A%0A%20%20%20%20%20%20%20%20fig%2C%20ax%20%3D%20plt.subplots(figsize%3D(6.5%2C%205.5))%0A%20%20%20%20%20%20%20%20display_hyper_fov(ARRAY%2C%20ax%3Dax)%0A%20%20%20%20%20%20%20%20fig.tight_layout()%0A%0A%20%20%20%20%20%20%20%20return%20mo.hstack(%5B%0A%20%20%20%20%20%20%20%20%20%20%20%20fig%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20mo.md(%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20f%22**%7Bmoment.time.isot%5B11%3A16%5D%7D%20UTC**%5Cn%5Cn%22%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20f%22-%20altitude%20**%7Balt%3A.1f%7D**%5Cn%22%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20f%22-%20azimuth%20**%7Baz%3A.1f%7D**%5Cn%22%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20f%22-%20hyper%20FoV%20**%7Barea%3A.1f%7D**%5Cn%22%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20f%22-%20mean%20multiplicity%20**%7Bmean%3A.2f%7D**%5Cn%22%0A%20%20%20%20%20%20%20%20%20%20%20%20)%2C%0A%20%20%20%20%20%20%20%20%5D%2C%20widths%3D%5B2%2C%201%5D%2C%20align%3D%22center%22)%0A%0A%20%20%20%20_snapshot()%0A%20%20%20%20return%0A%0A%0A%40app.cell(hide_code%3DTrue)%0Adef%20_(ARRAY%2C%20MIDNIGHT%2C%20TARGET%2C%20div%2C%20np%2C%20plt%2C%20source%2C%20u)%3A%0A%20%20%20%20def%20_across_the_night()%3A%0A%20%20%20%20%20%20%20%20hours%20%3D%20np.arange(-6%2C%206.01%2C%200.5)%0A%20%20%20%20%20%20%20%20times%2C%20areas%2C%20means%20%3D%20%5B%5D%2C%20%5B%5D%2C%20%5B%5D%0A%20%20%20%20%20%20%20%20for%20offset%20in%20hours%3A%0A%20%20%20%20%20%20%20%20%20%20%20%20moment%20%3D%20MIDNIGHT.after(offset%20*%20u.hour)%0A%20%20%20%20%20%20%20%20%20%20%20%20alt%2C%20az%20%3D%20moment.altaz_of(TARGET)%0A%20%20%20%20%20%20%20%20%20%20%20%20if%20alt%20%3C%3D%200%20*%20u.deg%3A%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20continue%0A%20%20%20%20%20%20%20%20%20%20%20%20ARRAY.divergent_pointing(div.value%2C%20alt%2C%20az)%0A%20%20%20%20%20%20%20%20%20%20%20%20times.append(offset)%0A%20%20%20%20%20%20%20%20%20%20%20%20areas.append(ARRAY.hyper_fov()%5B0%5D.to_value(u.deg**2))%0A%20%20%20%20%20%20%20%20%20%20%20%20means.append(ARRAY.multiplicity_moments()%5B0%5D)%0A%0A%20%20%20%20%20%20%20%20fig%2C%20ax%20%3D%20plt.subplots(figsize%3D(9%2C%204.5))%0A%20%20%20%20%20%20%20%20if%20not%20times%3A%0A%20%20%20%20%20%20%20%20%20%20%20%20ax.text(0.5%2C%200.5%2C%20f%22%7Bsource.value%7D%20does%20not%20rise%20on%20this%20night%22%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20ha%3D%22center%22%2C%20va%3D%22center%22%2C%20transform%3Dax.transAxes)%0A%20%20%20%20%20%20%20%20%20%20%20%20return%20fig%0A%0A%20%20%20%20%20%20%20%20ax.plot(times%2C%20areas%2C%20marker%3D%22o%22%2C%20color%3D%22tab%3Ablue%22)%0A%20%20%20%20%20%20%20%20ax.set_xlabel(%22hours%20from%20the%20middle%20of%20the%20night%22)%0A%20%20%20%20%20%20%20%20ax.set_ylabel(%22hyper%20FoV%20%5Bdeg%24%5E2%24%5D%22%2C%20color%3D%22tab%3Ablue%22)%0A%20%20%20%20%20%20%20%20ax.tick_params(axis%3D%22y%22%2C%20labelcolor%3D%22tab%3Ablue%22)%0A%0A%20%20%20%20%20%20%20%20twin%20%3D%20ax.twinx()%0A%20%20%20%20%20%20%20%20twin.plot(times%2C%20means%2C%20marker%3D%22s%22%2C%20color%3D%22tab%3Ared%22)%0A%20%20%20%20%20%20%20%20twin.set_ylabel(%22mean%20multiplicity%22%2C%20color%3D%22tab%3Ared%22)%0A%20%20%20%20%20%20%20%20twin.tick_params(axis%3D%22y%22%2C%20labelcolor%3D%22tab%3Ared%22)%0A%0A%20%20%20%20%20%20%20%20ax.set_title(f%22tracking%20%7Bsource.value%7D%20at%20div%20%3D%20%7Bdiv.value%7D%22)%0A%20%20%20%20%20%20%20%20ax.grid(alpha%3D0.3)%0A%20%20%20%20%20%20%20%20fig.tight_layout()%0A%20%20%20%20%20%20%20%20return%20fig%0A%0A%20%20%20%20_across_the_night()%0A%20%20%20%20return%0A%0A%0A%40app.cell(hide_code%3DTrue)%0Adef%20_(mo)%3A%0A%20%20%20%20mo.md(r%22%22%22%0A%20%20%20%20The%20covered%20area%20isn't%20constant%20while%20the%20array%20tracks%2C%20even%20though%0A%20%20%20%20%60div%60%20never%20changes.%0A%0A%20%20%20%20The%20array%20is%20**flat**%2C%20on%20a%20hillside%2C%20pointing%20in%20a%20**tilted**%0A%20%20%20%20direction.%20As%20the%20source%20drops%20toward%20the%20horizon%2C%20the%20array's%0A%20%20%20%20extent%20along%20the%20pointing%20axis%20grows%20relative%20to%20its%20extent%20across%0A%20%20%20%20it%2C%20so%20telescope%20offsets%20from%20the%20barycenter%20increasingly%20lie%0A%20%20%20%20*along*%20the%20line%20of%20sight%20rather%20than%20across%20it.%20Divergence%20is%0A%20%20%20%20driven%20by%20the%20across-axis%20part%2C%20so%20the%20fan%20closes%20up%2C%20the%20discs%20pile%0A%20%20%20%20back%20onto%20each%20other%2C%20the%20area%20shrinks%2C%20and%20multiplicity%20climbs.%0A%0A%20%20%20%20It's%20easy%20to%20assume%20a%20fixed%20%60div%60%20buys%20a%20fixed%20field%20of%20view.%20It%0A%20%20%20%20doesn't.%20Holding%20a%20hyper%20FoV%20through%20an%20observation%20means%20adjusting%0A%20%20%20%20%60div%60%20as%20the%20source%20moves%2C%20the%20question%20the%20next%20notebook%20takes%20up.%0A%20%20%20%20%22%22%22)%0A%20%20%20%20return%0A%0A%0Aif%20__name__%20%3D%3D%20%22__main__%22%3A%0A%20%20%20%20app.run()%0A