Skip to content

Commit

Permalink
Fix import
Browse files Browse the repository at this point in the history
  • Loading branch information
e-koch committed Sep 6, 2019
1 parent 640ce8c commit 98ab259
Showing 1 changed file with 10 additions and 9 deletions.
19 changes: 10 additions & 9 deletions 14A-235/HI/imaging/check_flux_recovery.py
Original file line number Diff line number Diff line change
Expand Up @@ -5,8 +5,9 @@
'''

import numpy as np
from spectral_cube import (SpectralCube, OneDSpectrum,
VaryingResolutionOneDSpectrum)
from spectral_cube import (SpectralCube, OneDSpectrum,)
from spectral_cube.lower_dimensional_structures import \
VaryingResolutionOneDSpectrum
import os
from astropy.io import fits
import astropy.units as u
Expand Down Expand Up @@ -66,8 +67,8 @@
plt.ylabel("VLA-to-EBHIS Flux Ratio")
plt.xlabel("Velocity (km / s)")
plt.tight_layout()
plt.savefig(allfigs_path("14A-235_imaging/vla_ebhis_14A_flux_recovery_ratio.png"))
plt.savefig(allfigs_path("14A-235_imaging/vla_ebhis_14A_flux_recovery_ratio.pdf"))
plt.savefig(allfigs_path("14A-235_imaging/vla_ebhis_14A_flux_recovery_ratio_v3.png"))
plt.savefig(allfigs_path("14A-235_imaging/vla_ebhis_14A_flux_recovery_ratio_v3.pdf"))
plt.close()

# Plot the total spectra
Expand All @@ -83,8 +84,8 @@
plt.ylabel("Total Flux (Jy)")
plt.xlabel("Velocity (km / s)")
plt.tight_layout()
plt.savefig(allfigs_path("14A-235_imaging/vla_ebhis_14A_flux_recovery.png"))
plt.savefig(allfigs_path("14A-235_imaging/vla_ebhis_14A_flux_recovery.pdf"))
plt.savefig(allfigs_path("14A-235_imaging/vla_ebhis_14A_flux_recovery_v3.png"))
plt.savefig(allfigs_path("14A-235_imaging/vla_ebhis_14A_flux_recovery_v3.pdf"))
plt.close()

plt.plot(vel_axis, total_ebhis_profile,
Expand All @@ -99,8 +100,8 @@
plt.ylabel("Total Flux (Jy)")
plt.xlabel("Velocity (km / s)")
plt.tight_layout()
plt.savefig(allfigs_path("14A-235_imaging/vla_ebhis_14A_flux_recovery_zoom.png"))
plt.savefig(allfigs_path("14A-235_imaging/vla_ebhis_14A_flux_recovery_zoom.pdf"))
plt.savefig(allfigs_path("14A-235_imaging/vla_ebhis_14A_flux_recovery_zoom_v3.png"))
plt.savefig(allfigs_path("14A-235_imaging/vla_ebhis_14A_flux_recovery_zoom_v3.pdf"))
plt.close()

# We've summed up most of the data already. How about a mass estimate?
Expand Down Expand Up @@ -153,6 +154,6 @@
meta=spec.meta,
beam=ebhis_cube.beam)
ebhis_spec.write(ebhis_m31_HI_data_path(
"14A-235_items/CAR_C01_14A235_match_04kms_spectralregrid.total_flux_spec.fits", no_check=True))
"14A-235_items/CAR_C01_14A235_match_04kms_spectralregrid.total_flux_spec.fits", no_check=True))

default_figure()

0 comments on commit 98ab259

Please sign in to comment.