From 3ae76e4c921f72b71ff7522462740e904136f428 Mon Sep 17 00:00:00 2001
From: Stephen Bailey <StephenBailey@lbl.gov>
Date: Tue, 21 Apr 2020 14:54:00 -0700
Subject: [PATCH] more freeze_iers

---
 doc/changes.rst      | 2 +-
 py/desisim/simexp.py | 5 +++--
 2 files changed, 4 insertions(+), 3 deletions(-)

diff --git a/doc/changes.rst b/doc/changes.rst
index 76e5e9d1c..e2e1735ca 100644
--- a/doc/changes.rst
+++ b/doc/changes.rst
@@ -5,7 +5,7 @@ desisim change log
 0.35.2 (unreleased)
 -------------------
 
-* No changes yet.
+* Add freeze_iers to more functions in simexp (direct to master).
 
 0.35.1 (2020-04-15)
 -------------------
diff --git a/py/desisim/simexp.py b/py/desisim/simexp.py
index af240f011..418afb1d7 100644
--- a/py/desisim/simexp.py
+++ b/py/desisim/simexp.py
@@ -21,6 +21,7 @@
 import desimodel.io
 from desimodel.focalplane import fiber_area_arcsec2
 import desiutil.depend
+from desiutil.iers import freeze_iers
 import desispec.interpolation
 import desisim.io
 import desisim.specsim
@@ -131,6 +132,7 @@ def simflat(flatfile, nspec=5000, nonuniform=False, exptime=10, testslit=False,
     import specsim.simulator
     from desiutil.log import get_logger
     log = get_logger()
+    freeze_iers()
 
     log.info('Reading flat lamp spectrum from {}'.format(flatfile))
     sbflux, hdr = fits.getdata(flatfile, header=True)
@@ -243,6 +245,7 @@ def simscience(targets, fiberassign, obsconditions='DARK', expid=None,
     '''
     from desiutil.log import get_logger
     log = get_logger()
+    freeze_iers()
 
     flux, wave, meta = targets
 
@@ -393,8 +396,6 @@ def simulate_spectra(wave, flux, fibermap=None, obsconditions=None, redshift=Non
 
     from desiutil.log import get_logger
     log = get_logger('DEBUG')
-
-    from desiutil.iers import freeze_iers
     freeze_iers()
 
     # Input cosmology to calculate the angular diameter distance of the galaxy's redshift