get rid of E402 and E741

            raise ValueError("Need to load planck dust map first")

        # Convert to galactic coordinates

        c = SkyCoord(ra=ra * u.degree, dec=dec * u.degree, frame="icrs").galactic

        l, b = c.l.radian, c.b.radian

        phi, b = c.l.radian, c.b.radian

        NSIDE = hp.pixelfunc.get_nside(self.ebv_planck)

        pix = hp.pixelfunc.ang2pix(nside=NSIDE, theta=np.pi / 2.0 - b, phi=l, nest=True)

        pix = hp.pixelfunc.ang2pix(nside=NSIDE, theta=np.pi / 2.0 - b, phi=phi, nest=True)

        return self.ebv_planck[pix] * self.Rv

    def apply_extinction(self, spec, Av=1.0):

from observation_sim.ObservationSim import Observation

from observation_sim._util import parse_args, make_run_dirs, generate_pointing_list

from importlib.metadata import version

from pathlib import Path

import os

import sys

import yaml

import shutil

import datetime

import importlib

import gc

gc.enable()

from pathlib import Path

import sys

ROOT = Path(__file__).resolve().parent.parent

sys.path.insert(0, str(ROOT))