PySD - Stella graphical function arrays not supported

[cday41]

Wondering if this is a bug or just unsupported, but when I convert a Stella xmile file that has graphical functions with multiple elements, I get the following error.

IndexError: list index out of range
in xmile_element.py, line 204: subs_list = self.subscripts[subscripts[0]]

Selecting apply-to-all in the Stella UI to collapse the elements to a single function doesn't work either.

Is this something that may get support in the future? Or can you recommend a workaround?

[Okan-26]

Hi cday, I got this same error in another issue I added support for. Could you send me your xmile/stmx file? I'm suspecting it's an error relating to the parsing of subscripts. If so, the edited pysd in my repository should work.

[cday41]

Thanks for the fast response Okan and for taking a look

PotholeCr_ATP_model_15Oct2020_TWOInverts_ForMs_ALT_Cleaned_preproc.zip

!

Here is a file I am working with. If you get past the index error, just a heads up that you will hit a couple other errors (LOOKUP function, and there is also an issue when 'apply-to-all' is selected for the graphical function but this may be resolved with the previous error).

Thanks again, pysd is a great tool.

[cday41]

Hi Okan,
Just checking in to see if you were able to access my file and reproduce the error?
Thanks again

[Okan-26]

Hi Casey, sorry for the late response. My guess is, after a closer look, the problem seems to be the inability of PySD to read the subscripts in ghost converters properly. This is the exact problem I've been working on recently.

These are the lines in your .xmile file that the first bug originates at:

<aux name="Converter 3" access="input">
   <units>days</units>
</aux>

These lines define the ghost converter for 'Converter 3'. There aren't any dimensions stated explicitly, unlike a regular converter, so the parser cannot get the subscripts from the element 'Converter 3'. I didn't see this converter in your model, so it's kinda odd. Anyways, when I delete this converter and run the PySD again, this is the line that produces a bug:

<aux name="average salmon wet mass" access="input">
   <units>grams wet mass/salmon</units>
</aux>

This is the ghost converter for "average salmon wet mass". Again, no dimensions, because I think stella has a way to refer the ghost converters back to the original converter, but PySD doesn't seem to work that way for now. Instead, PySD tries to read it as a regular converter with dimensions, equations, units, etc... and when it cannot retrieve the dimensions, we get the error in line "subs_list = self.subscripts[subscripts[0]]".