Input gacode

[jmcclena]

This pull request adds reading and writing of input.gacode files.

Example reading and writing an input.gacode file:

using GACODE
input_gacode = GACODE.expro_read("input.gacode")
GACODE.expro_write(input_gacode,"input.gacode_new")

[TimSlendebroek]

@jmcclena

Looks like it works I did some minor bug fixes and made the naming more uniform with the rest of the code base

This is now the API:

using GACODE

input_gacode = GACODE.InputGACODE(dd)
GACODE.save(ip,"<your filename>")
input_gacode = GACODE.load("<your filename>")

[TimSlendebroek]

I think we should test once reading an external input.gacode and export one to OMFIT to see if this works @jmcclena and then merge @orso82

[TimSlendebroek]

@orso82 or @bclyons12 it looks like I don't have merge rights on this repo can either of you merge this?

[orso82]

I still don't understand why one would need the special handling of DT in input.gacode

If the dd has a single hydrogenic ion with mass of 2.5 then that's what should also go in the input.gacode that is derived from dd

[jmcclena]

Are there IMAS functions to bundle and unbundle DT?

[orso82]

@bclyons12 I wish we did not yet merge this, since I don't believe this is the right way of handling DT.

The dd --> input.gacode should not be in the business of modifying the ion species. The translation from one plasma statefile to the other should be done verbatim.

Instead, I see the ability to modify the ion species (like split DT into separate D and T, or combine separate D and T into DT) as functionalities that should be added to IMAS.jl. Consider for example what we do with quasineutrality. We have IMAS.jl functions for those. We don't add a quasineutrality function for different file formats that we could export dd to.

To answer @jmcclena question, IMAS.jl has a lump_ions_as_bulk_and_impurity function, and we could add a bundle_DT and unbundle_DT functions there as well.

[bclyons12]

@orso82 There was no release made, so it should be simple to revert this if desired.