Since the vast majority of the CPU time in a typical BISICLES run is spent in the velocity solver, it can be useful to run the solver as a standalone code to examine performance and convergence issues. The code in BISICLES/code/benchmark is designed to read in a BISICLES plotfile written out as a part of a standard BISICLES run and perform the momentum balance velocity solve based on the configuration in the plotfile. In order for this to work, the option
amr.write_solver_rhs = 1
must be set in the input file used by the driver run. Then, the option
main.filename =hdf5 filename
must be set in the inputs file used to run the benchmark code, which tells the code where to look for the problem setup. To see how this all works (assuming that Chombo, etc are all set up already), try the following test: