scalar: write mean concentration to disk // save_fields: do not compute RHS if not saving pressure