from
equation (37). The model variable WR is a diagnostic variable in the model
and it is not strictly necessary to call WREAL each time step to compute WR.
However, to diagnose the solution and to understand where and when convergence/divergence
occur it is interesting to study the WR field.