Thermal Models

ACISpy provides the ability to run Xija thermal models via a special class, ThermalModelRunner, which can input commanded states from a variety of sources. The really nice thing about ThermalModelRunner is that it is actually a Dataset object, so we can look at the different fields, make plots, and create derived fields.

import acispy
import numpy as np

Basic Use of ThermalModelRunner

If a thermal model is being run entirely within the past or another situation where kadi commanded states are fully specified, you can simply run a thermal model by specifying the MSID to be modeled and the start and stop times of the model. The initial temperature will be specified from telemetry unless one sets an initial value using the T_init keyword argument. If we want to take model “bad times” into account, we also need to set mask_bad_times=True:

dea_model = acispy.ThermalModelRunner("1deamzt", "2020:142", "2020:150",
                                      get_msids=True, mask_bad_times=True)

acispy: [INFO     ] 2024-09-11 13:28:55,080 Using model for dea from chandra_models version = 3.54

Fetching msid: aoeclips over 2020:141:23:33:58.816 to 2020:150:00:21:26.816

get_msids=True is set, so the resulting Dataset object will be populated with actual telemetry as well as model data, so that you can examine and plot both:

print(dea_model["model","1deamzt"][:10])
print(dea_model["msids","1deamzt"][:10])

[18.04144287 17.69061317 17.34651811 17.03248643 16.73537596 16.44396592
 16.15814725 15.87781297 15.60285814 15.33317981] deg_C
[17.756287 17.506805 17.328583 17.22165  16.972137 16.687012 16.50879
 16.259277 15.902863 15.653381] deg_C

dp = acispy.DatePlot(dea_model, [("msids", "1deamzt"), ("model", "1deamzt")], plot_bad=True)
dp

/Users/jzuhone/mambaforge/envs/ska-dev/lib/python3.11/site-packages/Ska/Matplotlib/core.py:151: UserWarning: This figure was using a layout engine that is incompatible with subplots_adjust and/or tight_layout; not calling subplots_adjust.
  fig.autofmt_xdate()

Note that for this period of time a segment of it was marked as a “bad time” in the model, and this shows up in the cyan region in the plot.

By default, model specification files will be searched for in the chandra_models package. If your model spec file is not in that package, or you want to specify a different model specification file, you can pass ThermalModelRunner a model_spec optional argument:

dea_model2 = acispy.ThermalModelRunner("1deamzt", "2020:142", "2020:150", model_spec="dea_model_spec.json")

acispy: [INFO     ] 2024-09-11 13:28:58,482 model_spec = /Users/jzuhone/Source/acispy/doc/source/dea_model_spec.json

Fetching msid: aoeclips over 2020:141:23:33:58.816 to 2020:150:00:21:26.816

Running Thermal Models Using States from Various Sources

The ThermalModelRunner class also has various methods which one can use to run thermal models using states from various sources.

Running a Thermal Model from Commands

You can use a set of commands (either a kadi CommandTable or list of command dicts) to run a thermal model using the from_commands method:

from kadi import commands
# commands as a CommandTable
cmds = commands.get_cmds('2018:001:00:00:00', '2018:002:00:00:00')
psmc_model = acispy.ThermalModelRunner.from_commands("1pdeaat", cmds)
# commands as a list of dicts
dict_cmds = cmds.as_list_of_dict()
psmc_model2 = acispy.ThermalModelRunner.from_commands("1pdeaat", dict_cmds)

acispy: [INFO     ] 2024-09-11 13:28:59,044 Using model for psmc from chandra_models version = 3.54

Fetching msid: aoeclips over 2018:001:02:07:50.816 to 2018:002:00:10:46.816

acispy: [INFO     ] 2024-09-11 13:28:59,891 Using model for psmc from chandra_models version = 3.54

Fetching msid: aoeclips over 2018:001:02:07:50.816 to 2018:002:00:10:46.816

Running a Thermal Model from a Backstop File

You can use a backstop file as input to ThermalModelRunner to generate states for a thermal model run. In this case, ThermalModelRunner will go back to the latest telemetry to determine the initial temperature (unless one is specified, see below) and construct states from that point and handle continuity between those and the commands in the backstop file appropriately. To show this, we’ll use the ACA thermal model, which requires that the initial value for the aca0 pseudonode be specified, which we can do using the other_init keyword argument:

tm_aca = acispy.ThermalModelRunner.from_backstop(
    "aacccdpt", "/data/acis/LoadReviews/2020/AUG1720/ofls/CR229_2202.backstop",
    other_init={"aca0": -10})

acispy: [INFO     ] 2024-09-11 13:29:00,654 Using model for aca from chandra_models version = 3.54

Fetching msid: aoeclips over 2020:226:22:27:02.816 to 2020:237:01:59:26.816

dp = acispy.DatePlot(tm_aca, ('model','aacccdpt'), field2="pitch")
dp

/Users/jzuhone/mambaforge/envs/ska-dev/lib/python3.11/site-packages/Ska/Matplotlib/core.py:151: UserWarning: This figure was using a layout engine that is incompatible with subplots_adjust and/or tight_layout; not calling subplots_adjust.
  fig.autofmt_xdate()

Running a Thermal Model Using a states.dat File

One can run a thermal model from a states.dat table file which would be outputted by the thermal model check scripts during a load review, using the from_states_file() method:

dpa_model = acispy.ThermalModelRunner.from_states_file("1dpamzt", "my_states.dat")

acispy: [INFO     ] 2024-09-11 13:29:01,672 Using model for dpa from chandra_models version = 3.54

Fetching msid: aoeclips over 2017:035:18:00:22.816 to 2017:051:06:00:46.816

Running Thermal Models Using States Constructed By Hand

You can create a dictionary of states completely by hand and submit them as a keyword argument to ThermalModelRunner, for running completely hypothetical thermal models. For simplicity, we’ll pick constant states except change the CCD count.

states = {"ccd_count": np.array([5,6,1]),
          "pitch": np.array([150.0]*3),
          "fep_count": np.array([5,6,1]),
          "clocking": np.array([1]*3),
          "vid_board": np.array([1]*3),
          "off_nom_roll": np.array([0.0]*3),
          "simpos": np.array([-99616.0]*3),
          "datestart": np.array(["2015:002:00:00:00","2015:002:12:00:00","2015:003:12:00:00"]),
          "datestop": np.array(["2015:002:12:00:00","2015:003:12:00:00","2015:005:00:00:00"])}

In the previous examples, we never specified an initial temperature, but we can always do so, and must do so if we are past the end value in telemetry. This is not the case here, but we’ll set one anyway:

T_init = 13.0 # in degrees C

Now we can pass in the states dict we created, as well as the T_init:

dea_model = acispy.ThermalModelRunner("1deamzt", "2015:002:00:00:00",
                                      "2015:005:00:00:00", states=states, T_init=T_init)

acispy: [INFO     ] 2024-09-11 13:29:02,454 Using model for dea from chandra_models version = 3.54

Fetching msid: aoeclips over 2015:001:23:34:56.816 to 2015:005:00:17:20.816

dp = acispy.DatePlot(dea_model, ("model","1deamzt"), field2="ccd_count")
dp.set_ylim2(0,7)
dp

/Users/jzuhone/mambaforge/envs/ska-dev/lib/python3.11/site-packages/Ska/Matplotlib/core.py:151: UserWarning: This figure was using a layout engine that is incompatible with subplots_adjust and/or tight_layout; not calling subplots_adjust.
  fig.autofmt_xdate()

We can also dump the results of the model run to disk, both the states and the model components:

dea_model.write_model("model.dat", overwrite=True)
dea_model.write_states("states.dat", overwrite=True)

These files can be loaded in at a later date using ModelDataFromFiles.

Making Dashboard Plots

NOTE: This functionality requires the xijafit package to be installed.

It is possible to use the thermal model objects and the xijafit package to make dashboard plots. For this, use the make_dashboard_plots() method:

dpa_model_long = acispy.ThermalModelRunner("1dpamzt", "2016:200", "2017:200",
                                           mask_bad_times=True, get_msids=True)
dpa_model_long.make_dashboard_plots("1dpamzt", figfile="my_dpa_dash.png")

acispy: [INFO     ] 2024-09-11 13:29:03,334 Using model for dpa from chandra_models version = 3.54

Fetching msid: aoeclips over 2016:199:23:37:59.816 to 2017:200:00:17:34.816

WARNING: imaging routines will not be available,
failed to import sherpa.image.ds9_backend due to
'RuntimeErr: DS9Win unusable: Could not find xpaget on your PATH'

<Figure size 2000x1000 with 4 Axes>

figfile sets the filename to save the dashboard plot to.

One can also use the errorplotlimits and yplotlimits arguments to set the bounds of the temperature and the errors on the plots:

fp_model_long = acispy.ThermalModelRunner("fptemp_11", "2016:200", "2017:200",
                                          mask_bad_times=True, get_msids=True)
fp_model_long.make_dashboard_plots("fptemp_11", yplotlimits=(-120.0, -104.0),
                                   errorplotlimits=(-5.0, 5.0))

acispy: [INFO     ] 2024-09-11 13:29:16,085 Using model for acisfp from chandra_models version = 3.54

Fetching msid: aoeclips over 2016:199:23:37:59.816 to 2017:200:00:17:34.816

<Figure size 2000x1000 with 4 Axes>

Plotting Pitch and State Power Heating Values

All thermal models have a solar heating component. To make a quick solar heating plot, use the make_solarheat_plot method, with the node that the solar heating component acts upon as the first argument.

dpa_model_long.make_solarheat_plot("dpa0", figfile="dpa0_pitches.png")

<Figure size 1500x1000 with 1 Axes>

Similarly, to plot the ACIS state power coefficients for a model (if present), use the make_power_plot method:

# For the 1DEAMZT model set use_ccd_count=True
dpa_model_long.make_power_plot(figfile="acis_state_power.png", use_ccd_count=False)

<Figure size 1000x1000 with 1 Axes>

Simulating Single States

The SimulateSingleState class is a simplified implemenation of ThermalModelRunner which assumes that the spacecraft state is constant over a period of time. This class takes the name of the temperature to simulate, a start and stop date, and a dictionary of states. It is not required to specify all states in this dictionary, but note that states not specified will be assumed to be the following:

• pitch: 90 deg

• off_nom_roll: 0 deg

• ccd_count: 0

• fep_count: 0

• clocking: 0

• vid_board: 0

• simpos: -99616.0

• hetg: "RETR"

• letg: "RETR"

• dh_heater: 0

• q1: 1.0

• q2: 0.0

• q3: 0.0

• q4: 0.0

datestart = "2019:181:01:00:00" # start time of run
datestop = "2019:183:01:00:00" # stop time of run
states = {"pitch": 70.0}
T_init = -11.0 # in degrees F
aca_m = acispy.SimulateSingleState("aacccdpt", datestart, datestop, states, T_init, other_init={"aca0": T_init})

acispy: [INFO     ] 2024-09-11 13:29:23,536 Using model for aca from chandra_models version = 3.54

Fetching msid: aoeclips over 2019:181:00:35:02.816 to 2019:183:01:19:42.816

dp = acispy.DatePlot(aca_m, "aacccdpt")
dp

/Users/jzuhone/mambaforge/envs/ska-dev/lib/python3.11/site-packages/Ska/Matplotlib/core.py:151: UserWarning: This figure was using a layout engine that is incompatible with subplots_adjust and/or tight_layout; not calling subplots_adjust.
  fig.autofmt_xdate()

Simulating ECS Runs

A further special case of running a single state is an ECS run, which may be performed after a safing action. The SimulateECSRun class is a limited version of SimulateSingleState which assumes that the SIM-Z position is HRC-S (-99616) and that ACIS is clocking. The goal is to predict if the temperature will hit the planning limit within the time frame of the ECS run. A start time is specified, as well as the number of hours in the ECS run (the same number of hours as specified on the ECS CAP). Then, an attitude must be specified, which can take one of three forms: * A (pitch, off-nominal roll) combination * An attitude quaternion * The name of a load (for vehicle load ECS runs)

We pick a model to run (most relevant is "1dpamzt") along with a start time a length of the ECS run in hours (24 in this case), and feed the parameters into SimulateECSRun:

datestart = "2015:002:00:00:00" # start time of run
hours = 24 # length of ECS run in hours
attitude = [173.0, 0.0] # pitch, off-nominal roll in degrees
T_init = 7.5 # in degrees C
ccd_count = 6 # number of CCDs
dpa_ecs_run = acispy.SimulateECSRun("1dpamzt", datestart, hours, T_init, attitude,
                                    ccd_count)

acispy: [INFO     ] 2024-09-11 13:29:24,622 Using model for dpa from chandra_models version = 3.54

Fetching msid: aoeclips over 2015:001:23:34:56.816 to 2015:003:12:18:00.816

acispy: [INFO     ] 2024-09-11 13:29:24,859 Run Parameters

acispy: [INFO     ] 2024-09-11 13:29:24,859 --------------

acispy: [INFO     ] 2024-09-11 13:29:24,860 Modeled Temperature: 1dpamzt

acispy: [INFO     ] 2024-09-11 13:29:24,860 Start Datestring: 2015:002:00:00:00.000

acispy: [INFO     ] 2024-09-11 13:29:24,860 Length of state in hours: 24

acispy: [INFO     ] 2024-09-11 13:29:24,861 Stop Datestring: 2015:003:12:00:00.000

acispy: [INFO     ] 2024-09-11 13:29:24,861 Initial Temperature: 7.5 degrees C

acispy: [INFO     ] 2024-09-11 13:29:24,861 CCD/FEP Count: 6

acispy: [INFO     ] 2024-09-11 13:29:24,861 Pitch: 173.0

acispy: [INFO     ] 2024-09-11 13:29:24,862 Off-nominal Roll: 0.0

acispy: [INFO     ] 2024-09-11 13:29:24,862 Detector Housing Heater: OFF

acispy: [INFO     ] 2024-09-11 13:29:24,863 Model Result

acispy: [INFO     ] 2024-09-11 13:29:24,863 ------------

acispy: [INFO     ] 2024-09-11 13:29:24,867 The limit of 38.5 degrees C will be reached at 2015:002:12:47:36.816, after 46.05681599998474 ksec.

acispy: [INFO     ] 2024-09-11 13:29:24,867 The limit is reached before the end of the observation.

acispy: [WARNING  ] 2024-09-11 13:29:24,868 This observation is NOT safe from a thermal perspective.

The run reports back the input parameters and the time when the limit was reached, if it was at all. We can plot the model using the plot_model() method, which shows the limit value as a dashed green line and the time at which the limit was reached as a dashed red line, as well as whether or not this was a safe ECS run:

dpa_ecs_run.plot_model()

/Users/jzuhone/mambaforge/envs/ska-dev/lib/python3.11/site-packages/Ska/Matplotlib/core.py:151: UserWarning: This figure was using a layout engine that is incompatible with subplots_adjust and/or tight_layout; not calling subplots_adjust.
  fig.autofmt_xdate()

On the other hand, if the ECS run had been shorter, the limit would be reached after the ECS run, so this is safe.

hours = 12
dpa_ecs_run = acispy.SimulateECSRun("1dpamzt", datestart, hours, T_init, attitude,
                                    ccd_count)
dpa_ecs_run.plot_model()

acispy: [INFO     ] 2024-09-11 13:29:25,400 Using model for dpa from chandra_models version = 3.54

Fetching msid: aoeclips over 2015:001:23:34:56.816 to 2015:002:18:21:04.816

acispy: [INFO     ] 2024-09-11 13:29:25,606 Run Parameters

acispy: [INFO     ] 2024-09-11 13:29:25,606 --------------

acispy: [INFO     ] 2024-09-11 13:29:25,606 Modeled Temperature: 1dpamzt

acispy: [INFO     ] 2024-09-11 13:29:25,607 Start Datestring: 2015:002:00:00:00.000

acispy: [INFO     ] 2024-09-11 13:29:25,607 Length of state in hours: 12

acispy: [INFO     ] 2024-09-11 13:29:25,607 Stop Datestring: 2015:002:18:00:00.000

acispy: [INFO     ] 2024-09-11 13:29:25,607 Initial Temperature: 7.5 degrees C

acispy: [INFO     ] 2024-09-11 13:29:25,608 CCD/FEP Count: 6

acispy: [INFO     ] 2024-09-11 13:29:25,608 Pitch: 173.0

acispy: [INFO     ] 2024-09-11 13:29:25,608 Off-nominal Roll: 0.0

acispy: [INFO     ] 2024-09-11 13:29:25,608 Detector Housing Heater: OFF

acispy: [INFO     ] 2024-09-11 13:29:25,609 Model Result

acispy: [INFO     ] 2024-09-11 13:29:25,610 ------------

acispy: [INFO     ] 2024-09-11 13:29:25,610 The limit of 38.5 degrees C will be reached at 2015:002:12:47:36.816, after 46.05681599998474 ksec.

acispy: [INFO     ] 2024-09-11 13:29:25,611 The limit is reached after the end of the observation.

acispy: [INFO     ] 2024-09-11 13:29:25,611 This observation is safe from a thermal perspective.

/Users/jzuhone/mambaforge/envs/ska-dev/lib/python3.11/site-packages/Ska/Matplotlib/core.py:151: UserWarning: This figure was using a layout engine that is incompatible with subplots_adjust and/or tight_layout; not calling subplots_adjust.
  fig.autofmt_xdate()

Note that for some combinations of parameters the limit may never be reached. For example, let’s knock the CCD count down to 4:

ccd_count = 4 # only 4 CCDs
dpa_ecs_run = acispy.SimulateECSRun("1dpamzt", datestart, hours, T_init, attitude,
                                    ccd_count)
dpa_ecs_run.plot_model()

acispy: [INFO     ] 2024-09-11 13:29:26,055 Using model for dpa from chandra_models version = 3.54

Fetching msid: aoeclips over 2015:001:23:34:56.816 to 2015:002:18:21:04.816

acispy: [INFO     ] 2024-09-11 13:29:26,233 Run Parameters

acispy: [INFO     ] 2024-09-11 13:29:26,235 --------------

acispy: [INFO     ] 2024-09-11 13:29:26,235 Modeled Temperature: 1dpamzt

acispy: [INFO     ] 2024-09-11 13:29:26,236 Start Datestring: 2015:002:00:00:00.000

acispy: [INFO     ] 2024-09-11 13:29:26,236 Length of state in hours: 12

acispy: [INFO     ] 2024-09-11 13:29:26,236 Stop Datestring: 2015:002:18:00:00.000

acispy: [INFO     ] 2024-09-11 13:29:26,236 Initial Temperature: 7.5 degrees C

acispy: [INFO     ] 2024-09-11 13:29:26,237 CCD/FEP Count: 4

acispy: [INFO     ] 2024-09-11 13:29:26,237 Pitch: 173.0

acispy: [INFO     ] 2024-09-11 13:29:26,237 Off-nominal Roll: 0.0

acispy: [INFO     ] 2024-09-11 13:29:26,238 Detector Housing Heater: OFF

acispy: [INFO     ] 2024-09-11 13:29:26,238 Model Result

acispy: [INFO     ] 2024-09-11 13:29:26,239 ------------

acispy: [INFO     ] 2024-09-11 13:29:26,239 The limit of 38.5 degrees C is never reached.

acispy: [INFO     ] 2024-09-11 13:29:26,239 This observation is safe from a thermal perspective.

/Users/jzuhone/mambaforge/envs/ska-dev/lib/python3.11/site-packages/Ska/Matplotlib/core.py:151: UserWarning: This figure was using a layout engine that is incompatible with subplots_adjust and/or tight_layout; not calling subplots_adjust.
  fig.autofmt_xdate()

Simulating ECS Runs with Vehicle Loads

If the spacecraft executed SCS-107, we may be running an ECS run while the vehicle load is still running, which means that the pitch and off-nominal roll may change during the ECS run. If this is the case, pass the name of the load to the attitude parameter. An example of an ECS run with a vehicle load that is not safe:

datestart = "2017:256:03:20:00"
hours = 24
T_init = 12.0 # in degrees C
attitude = "vehicle" # attitude info will be take from vehicle loads
ccd_count = 6 # number of CCDs
dpa_ecs_run = acispy.SimulateECSRun("1dpamzt", datestart, hours, T_init, attitude, ccd_count)
dpa_ecs_run.plot_model()

acispy: [INFO     ] 2024-09-11 13:29:27,654 Modeling a 6-chip state concurrent with states from the following vehicle loads: {'SEP0917C', 'SEP1317B'}

acispy: [INFO     ] 2024-09-11 13:29:27,919 Using model for dpa from chandra_models version = 3.54

Fetching msid: aoeclips over 2017:256:02:54:54.816 to 2017:257:15:37:58.816

acispy: [INFO     ] 2024-09-11 13:29:28,115 Run Parameters

acispy: [INFO     ] 2024-09-11 13:29:28,115 --------------

acispy: [INFO     ] 2024-09-11 13:29:28,116 Modeled Temperature: 1dpamzt

acispy: [INFO     ] 2024-09-11 13:29:28,116 Start Datestring: 2017:256:03:20:00.000

acispy: [INFO     ] 2024-09-11 13:29:28,116 Length of state in hours: 24

acispy: [INFO     ] 2024-09-11 13:29:28,117 Stop Datestring: 2017:257:15:20:00.000

acispy: [INFO     ] 2024-09-11 13:29:28,117 Initial Temperature: 12.0 degrees C

acispy: [INFO     ] 2024-09-11 13:29:28,117 CCD/FEP Count: 6

acispy: [INFO     ] 2024-09-11 13:29:28,117 Detector Housing Heater: OFF

acispy: [INFO     ] 2024-09-11 13:29:28,118 Model Result

acispy: [INFO     ] 2024-09-11 13:29:28,118 ------------

acispy: [INFO     ] 2024-09-11 13:29:28,119 The limit of 38.5 degrees C will be reached at 2017:256:12:34:22.816, after 33.26281599998474 ksec.

acispy: [INFO     ] 2024-09-11 13:29:28,119 The limit is reached before the end of the observation.

acispy: [WARNING  ] 2024-09-11 13:29:28,120 This observation is NOT safe from a thermal perspective.

/Users/jzuhone/mambaforge/envs/ska-dev/lib/python3.11/site-packages/Ska/Matplotlib/core.py:151: UserWarning: This figure was using a layout engine that is incompatible with subplots_adjust and/or tight_layout; not calling subplots_adjust.
  fig.autofmt_xdate()

But if we drop it down to 5 chips, it is safe:

ccd_count = 5 # number of CCDs
dpa_ecs_run = acispy.SimulateECSRun("1dpamzt", datestart, hours, T_init, attitude, ccd_count)
dpa_ecs_run.plot_model()

acispy: [INFO     ] 2024-09-11 13:29:28,788 Modeling a 5-chip state concurrent with states from the following vehicle loads: {'SEP0917C', 'SEP1317B'}

acispy: [INFO     ] 2024-09-11 13:29:28,984 Using model for dpa from chandra_models version = 3.54

Fetching msid: aoeclips over 2017:256:02:54:54.816 to 2017:257:15:37:58.816

acispy: [INFO     ] 2024-09-11 13:29:29,220 Run Parameters

acispy: [INFO     ] 2024-09-11 13:29:29,221 --------------

acispy: [INFO     ] 2024-09-11 13:29:29,221 Modeled Temperature: 1dpamzt

acispy: [INFO     ] 2024-09-11 13:29:29,221 Start Datestring: 2017:256:03:20:00.000

acispy: [INFO     ] 2024-09-11 13:29:29,221 Length of state in hours: 24

acispy: [INFO     ] 2024-09-11 13:29:29,222 Stop Datestring: 2017:257:15:20:00.000

acispy: [INFO     ] 2024-09-11 13:29:29,222 Initial Temperature: 12.0 degrees C

acispy: [INFO     ] 2024-09-11 13:29:29,223 CCD/FEP Count: 5

acispy: [INFO     ] 2024-09-11 13:29:29,223 Detector Housing Heater: OFF

acispy: [INFO     ] 2024-09-11 13:29:29,224 Model Result

acispy: [INFO     ] 2024-09-11 13:29:29,224 ------------

acispy: [INFO     ] 2024-09-11 13:29:29,231 The limit of 38.5 degrees C is never reached.

acispy: [INFO     ] 2024-09-11 13:29:29,232 This observation is safe from a thermal perspective.

/Users/jzuhone/mambaforge/envs/ska-dev/lib/python3.11/site-packages/Ska/Matplotlib/core.py:151: UserWarning: This figure was using a layout engine that is incompatible with subplots_adjust and/or tight_layout; not calling subplots_adjust.
  fig.autofmt_xdate()

Simulating the Focal Plane Temperature for ECS Runs

The ACIS Focal Plane (FP) temperature can also be modeled for an ECS run. This will track the amount of “cold time” for the ECS measurement. In this case, a quaternion generally should be supplied for the value of attitude, since the FP temperature can depend on the value of the Earth solid angle.

datestart = "2020:148:14:45:00" # start time of run
hours = 24 # length of ECS run in hours
T_init = -115.0 # in degrees C
ccd_count = 4 # number of CCDs
attitude = [-0.04470333, 0.63502552, -0.67575906, 0.37160988] # attitude quaternion required

fp_ecs_run = acispy.SimulateECSRun("fptemp_11", datestart, hours, T_init, attitude, ccd_count)
fp_ecs_run.plot_model()

acispy: [INFO     ] 2024-09-11 13:29:30,031 Using model for acisfp from chandra_models version = 3.54

acispy: [WARNING  ] 2024-09-11 13:29:30,104 Using quaternions to calculate pitch and roll because both were not specified in the states.

Fetching msid: aoeclips over 2020:148:14:22:22.816 to 2020:150:03:05:26.816

acispy: [INFO     ] 2024-09-11 13:29:30,268 Run Parameters

acispy: [INFO     ] 2024-09-11 13:29:30,269 --------------

acispy: [INFO     ] 2024-09-11 13:29:30,269 Modeled Temperature: fptemp_11

acispy: [INFO     ] 2024-09-11 13:29:30,269 Start Datestring: 2020:148:14:45:00.000

acispy: [INFO     ] 2024-09-11 13:29:30,269 Length of state in hours: 24

acispy: [INFO     ] 2024-09-11 13:29:30,270 Stop Datestring: 2020:150:02:45:00.000

acispy: [INFO     ] 2024-09-11 13:29:30,270 Initial Temperature: -115.0 degrees C

acispy: [INFO     ] 2024-09-11 13:29:30,270 CCD/FEP Count: 4

acispy: [INFO     ] 2024-09-11 13:29:30,271 Quaternion: [-0.04470333, 0.63502552, -0.67575906, 0.37160988]

acispy: [INFO     ] 2024-09-11 13:29:30,271 Pitch: 155.13828486037553

acispy: [INFO     ] 2024-09-11 13:29:30,271 Off-nominal Roll: 10.637704660891393

acispy: [INFO     ] 2024-09-11 13:29:30,272 Detector Housing Heater: OFF

acispy: [INFO     ] 2024-09-11 13:29:30,272 Model Result

acispy: [INFO     ] 2024-09-11 13:29:30,273 ------------

acispy: [INFO     ] 2024-09-11 13:29:30,274 The focal plane is never cold for this ECS measurement.

acispy: [INFO     ] 2024-09-11 13:29:30,275 The limit of -86.0 degrees C is never reached.

/Users/jzuhone/mambaforge/envs/ska-dev/lib/python3.11/site-packages/Ska/Matplotlib/core.py:151: UserWarning: This figure was using a layout engine that is incompatible with subplots_adjust and/or tight_layout; not calling subplots_adjust.
  fig.autofmt_xdate()

/Users/jzuhone/mambaforge/envs/ska-dev/lib/python3.11/site-packages/Ska/Matplotlib/core.py:97: UserWarning: constrained_layout not applied because axes sizes collapsed to zero.  Try making figure larger or axes decorations smaller.
  ax.figure.canvas.draw_idle()

This can also be done for a vehicle load as the attitude parameter:

datestart = "2017:256:03:20:00"
hours = 24
T_init = -119.0 # in degrees C
attitude = "vehicle" # attitude info will be take from the vehicle loads
ccd_count = 6 # number of CCDs
fp_ecs_run = acispy.SimulateECSRun("fptemp_11", datestart, hours, T_init, attitude, ccd_count)
fp_ecs_run.plot_model()

acispy: [INFO     ] 2024-09-11 13:29:30,659 Modeling a 6-chip state concurrent with states from the following vehicle loads: {'SEP0917C', 'SEP1317B'}

acispy: [INFO     ] 2024-09-11 13:29:30,851 Using model for acisfp from chandra_models version = 3.54

Fetching msid: aoeclips over 2017:256:02:54:54.816 to 2017:257:15:37:58.816

acispy: [INFO     ] 2024-09-11 13:29:31,029 Run Parameters

acispy: [INFO     ] 2024-09-11 13:29:31,029 --------------

acispy: [INFO     ] 2024-09-11 13:29:31,030 Modeled Temperature: fptemp_11

acispy: [INFO     ] 2024-09-11 13:29:31,030 Start Datestring: 2017:256:03:20:00.000

acispy: [INFO     ] 2024-09-11 13:29:31,030 Length of state in hours: 24

acispy: [INFO     ] 2024-09-11 13:29:31,030 Stop Datestring: 2017:257:15:20:00.000

acispy: [INFO     ] 2024-09-11 13:29:31,030 Initial Temperature: -119.0 degrees C

acispy: [INFO     ] 2024-09-11 13:29:31,031 CCD/FEP Count: 6

acispy: [INFO     ] 2024-09-11 13:29:31,031 Detector Housing Heater: OFF

acispy: [INFO     ] 2024-09-11 13:29:31,032 Model Result

acispy: [INFO     ] 2024-09-11 13:29:31,032 ------------

acispy: [INFO     ] 2024-09-11 13:29:31,032 The focal plane is cold for 48.872 ks.

acispy: [INFO     ] 2024-09-11 13:29:31,033 The limit of -86.0 degrees C is never reached.

/Users/jzuhone/mambaforge/envs/ska-dev/lib/python3.11/site-packages/Ska/Matplotlib/core.py:151: UserWarning: This figure was using a layout engine that is incompatible with subplots_adjust and/or tight_layout; not calling subplots_adjust.
  fig.autofmt_xdate()

/Users/jzuhone/mambaforge/envs/ska-dev/lib/python3.11/site-packages/Ska/Matplotlib/core.py:97: UserWarning: constrained_layout not applied because axes sizes collapsed to zero.  Try making figure larger or axes decorations smaller.
  ax.figure.canvas.draw_idle()