Draft: C13 c12 script
A script that uses prodimopy routines to read and plot the output from the new flag estimate_13CO_12CO_ratio.
@c.rab Hi I would like to create a directory where we can put scripts related to specific ProDiMo output or scripts used for papers that uses ProDiMo and prodimopy. Users can create their own sub-directories for their own projects.
Thank you, Wing-Fai, for the contribution.
I have to think about your comment concerning a script directory. prodimopy is designed as a package, e.g. it can be installed via a package manager. That means (user)scripts cannot be easily edited or used. There are some standard ones, but one has to include them in the package manually, and they cannot be easily edited.
If you have a new output file and want a read routine, it is best to include it in read.py. There are already several read routines that can be used as examples of how to code a read routine. Even if the output is not standard (e.g. not included in every model), it should be there, as it does not harm. This ensures that there is always only one read routine for an output file in ProDiMo. If those read routines are distributed in various scripts, maintaining those scripts will be difficult (and some of them will break over time).
A similar approach is used for plotting (e.g. see plot.py).
However, I understand that there might be project-specific plots or tools that would be useful to share. I am unsure if prodimopy is the right place to store them as it is designed as a general package.
If I understand your case, the C13_C12 is a new feature in ProDiMo, and I don't see why we should not have general read and plotting routines for that in prodimopy.