diff --git a/templates/input.nml b/templates/input.nml
index 25502154b4a4268f65e0f57d7e42136591e4208e..cfca979ace700586cfb0ec4bd02bfacedc00c40b 100644
--- a/templates/input.nml
+++ b/templates/input.nml
@@ -51,7 +51,7 @@
    output_timestamps        = .false.,
    trace_execution          = .false.,
 
-   stages_to_write          = 'preassim', 'postassim'
+   stages_to_write          = 'preassim', 'postassim', 'output'
    output_members           = .true.
    output_mean              = .true.
    output_sd                = .true.