From 60cba749d78e432f774963565f9da3eb92b2a380 Mon Sep 17 00:00:00 2001
From: "Karl D. Gordon" <kgordon@stsci.edu>
Date: Tue, 28 Nov 2023 09:18:35 -0500
Subject: [PATCH] adding to docs and tests

---
 docs/dust_extinction/choose_model.rst  |  2 ++
 docs/dust_extinction/model_flavors.rst | 16 ++++++++++------
 dust_extinction/grain_models.py        |  2 +-
 3 files changed, 13 insertions(+), 7 deletions(-)

diff --git a/docs/dust_extinction/choose_model.rst b/docs/dust_extinction/choose_model.rst
index d388bd5..64cf017 100644
--- a/docs/dust_extinction/choose_model.rst
+++ b/docs/dust_extinction/choose_model.rst
@@ -172,6 +172,8 @@ observationally (e.g., in the extreme UV below 912 A).
 +--------------+----------------+------------------+--------------+
 |   J13 MWRV31 |   0.00001 - 25 |    0.04 - 100000 |  MW R(V)=3.1 |
 +--------------+----------------+------------------+--------------+
+|  HD23 MWRV31 |  0.000033 - 10 |      0.1 - 30000 |  MW R(V)=3.1 |
++--------------+----------------+------------------+--------------+
 
 Shape Models
 ============
diff --git a/docs/dust_extinction/model_flavors.rst b/docs/dust_extinction/model_flavors.rst
index 0a6b129..251180d 100644
--- a/docs/dust_extinction/model_flavors.rst
+++ b/docs/dust_extinction/model_flavors.rst
@@ -351,7 +351,7 @@ Grain models
    import matplotlib.pyplot as plt
    import astropy.units as u
 
-   from dust_extinction.grain_models import DBP90, WD01, D03, ZDA04, C11, J13
+   from dust_extinction.grain_models import DBP90, WD01, D03, ZDA04, C11, J13, HD23
 
    fig, ax = plt.subplots()
 
@@ -363,12 +363,14 @@ Grain models
              WD01, WD01, WD01,
              D03, D03, D03,
              ZDA04,
-             C11, J13]
+             C11, J13,
+             HD23]
    modelnames = ["MWRV31",
                  "MWRV31", "MWRV40", "MWRV55",
                  "MWRV31", "MWRV40", "MWRV55",
                  "BARE-GR-S",
-                 "MWRV31", "MWRV31"]
+                 "MWRV31", "MWRV31",
+                 "MWRV31"]
 
    for cmodel, cname in zip(models, modelnames):
       ext_model = cmodel(cname)
@@ -400,7 +402,7 @@ Grain models
    from matplotlib.ticker import ScalarFormatter
    import astropy.units as u
 
-   from dust_extinction.grain_models import DBP90, WD01, D03, ZDA04, C11, J13
+   from dust_extinction.grain_models import DBP90, WD01, D03, ZDA04, C11, J13, HD23
 
    fig, ax = plt.subplots()
 
@@ -412,12 +414,14 @@ Grain models
              WD01, WD01, WD01,
              D03, D03, D03,
              ZDA04,
-             C11, J13]
+             C11, J13,
+             HD23]
    modelnames = ["MWRV31",
                  "MWRV31", "MWRV40", "MWRV55",
                  "MWRV31", "MWRV40", "MWRV55",
                  "BARE-GR-S",
-                 "MWRV31", "MWRV31"]
+                 "MWRV31", "MWRV31",
+                 "MWRV31"]
 
    for cmodel, cname in zip(models, modelnames):
       ext_model = cmodel(cname)
diff --git a/dust_extinction/grain_models.py b/dust_extinction/grain_models.py
index bda8073..f9e04fb 100644
--- a/dust_extinction/grain_models.py
+++ b/dust_extinction/grain_models.py
@@ -628,7 +628,7 @@ class HD23(GMBase):
         plt.show()
     """
 
-    x_range = [1.0 / 1e5, 1.0 / 4e-2]
+    x_range = [1.0 / 3e4, 1.0 / 0.1]
 
     possnames = {"MWRV31": ("astrodust+PAH_MW_RV3.1.fits", 3.1)}