Segmentation fault in calc_cp_des

Program terminated with signal SIGSEGV, Segmentation fault.
#0  0x0000557607efcad7 in des_thermo::calc_cp_des (pnp=1)
    at mfix/src/model/des/des_thermo_mod.f:118
118	      lTs = DES_T_s(pNP)
(gdb) p pNP
$1 = 1
(gdb) p DES_T_s
$2 = <not allocated>

Energy equations are not enabled, so this array is not allocated. But I’m not sure why we got to calc_cp_des without energy equations. If I enable energy equations, I get a different error message about no species being defined.

Here’s the full stack trace:

#0  0x0000557607efcad7 in des_thermo::calc_cp_des (pnp=1)
    at mfix/src/model/des/des_thermo_mod.f:118
#1  0x000055760836c64b in calc_thermo_des_mod::calc_thermo_des ()
    at mfix/src/model/des/calc_thermo_des.f:48
#2  0x0000557607f47189 in pic_time_march_mod::pic_time_step (pic_iters=1)
    at mfix/src/model/des/pic/time_march_pic.f90:131
#3  0x0000557607e566b1 in run_mfix::run_pic ()
    at mfix/src/model/mfix.f:239