def compute_diagnostics(self, num_iter=<NUM_LIT:3>):
for n in range(num_iter):<EOL><INDENT>ignored = self.compute()<EOL><DEDENT>
Compute all tendencies and diagnostics, but don't update model state.
        By default it will call compute() 3 times to make sure all
        subprocess coupling is accounted for. The number of iterations can
        be changed with the input argument.