The function model_output() allows you to generate outputs specific to your model. This could be used, for example, to look for topological defects, to output different combinations of the fields, or for whatever else you wanted. The function takes two inputs, the integer flush and the string ext_. The former tells the function whether or not infrequent calculations such as checkpointing and spectra are being done ( and for yes and no respectively). This way you can choose for your output tasks to be done each time a save is performed or only when infrequent calculations are done. In the latter case you can simply have model_output() return without doing anything except when . The second variable is the extension being used for output files, so that if you wish you can create other output files with the same extension.
This function can also be used to set the global variable rescaling, which will be multiplied by all field values output by meansvars(), histograms(), histograms2d(), and slices(). If you just want to set this value once for the entire run you should probably do that in model_initialize(), but if you do it in model_output() then you can rescale dynamically as the run continues.