nsb 0.1.0

Night-sky background evaluation and observation planning for ground-based astronomy.
Documentation 
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40

// SPDX-License-Identifier: AGPL-3.0-only
// Copyright (C) 2026 Vallés Puig, Ramon
use chrono::{DateTime, Utc};
use nsb::{ComponentMask, NsbEvaluator, PointQuery, Target, DEG};
use siderust::catalogs::observatories;
use tempoch::{Time, UTC};

fn main() -> nsb::Result<()> {
    let time = DateTime::parse_from_rfc3339("2023-09-04T01:48:00Z")
        .expect("valid example timestamp")
        .with_timezone(&Utc);

    let query = PointQuery::new(
        observatories::EL_PARANAL.geodetic(),
        Time::<UTC>::from_chrono(time),
        Target::new(266.41683 * DEG, -29.00781 * DEG),
    )
    .with_components(ComponentMask::ZODIACAL | ComponentMask::AIRGLOW);

    let result = NsbEvaluator::new()?.evaluate(&query)?;

    for component in &result.components {
        println!(
            "{:>10}: integrated = {:.6e}  ph/(cm² ns sr)  | B = {:.3} S10 | V = {:.3} S10",
            component.name,
            component.integrated.value(),
            component.b_flux_s10.value(),
            component.v_flux_s10.value()
        );
    }
    println!("--------------------");
    println!(
        "    total: {:.6e} ph/(cm² ns sr)",
        result.integrated.value()
    );
    println!("       B = {:.3} mag/arcsec²", result.b_mag.value());
    println!("       V = {:.3} mag/arcsec²", result.v_mag.value());

    Ok(())
}