Jumping in here (as the advisor!), our 1st step will be to run the pure trajectory code (no dispersion etc.) So what we are trying to do is figure out what to edit and where to edit it! For example: (1) in searching the code (in /library) for the planet's radius, we have found that to be set to 6371 km in several locations. Does that seem right (as opposed to being defined in just one location)? Either way, we can edit that for Mars. (2) If we keep sedimentation, we'll have to do the save for gravity. (3) The next thing is terrain: is there ONE location/routine in which terrain is read in (so that we can override and input Mars' terrain data)? And (4): is there ONE location in the code at which meteorological data is read in? Ideally, we can then slot in our own read statements to ingest wind output from the Mars GCM (as opposed to convincing the Mars GCM to output data in a form hysplit already knows).
Any help with these source code questions would be great! Final example: is there a subset of the code which is for trajectory calculation ONLY (so we could ignore the rest for now)??