Next: Data part, apply nonlinear Up: Input file Previous: Data part, prepare post-processing Contents
control_timestep 10 1. 1.
control_groundflow_nonsaturated_apply 10 -no
control_groundflow_seepage_apply 10 -no
control_print 10 -time_current -post_node_rhside_ratio
This nonsaturated dam calculation is highly nonlinear because of the group_groundflow_nonsaturated_vangenuchten and the groundflow_seepage_geometry records. In order to get a solution a trick is commonly applied to this kind of calculation: first a linear solution is established, and afterwards the full nonlinear solution is searched iteratively with the linear solution as first guess.
We obtain a linear calculation with the extra records control_groundflow_nonsaturated_apply 10 and control_groundflow_seepage_apply 10. These records de-activate the group_groundflow_nonsaturated_vangenuchten and groundflow_seepage_geometry for the time-steps of control_timestep 10. Notice that we print the current time -time_current and the error ratio -post_node_rhside_ratio. The error ratio is the ratio of the maximum out-of-balance flux at nodes without free hydraulic head and the maximum flux at nodes with prescribed hydraulic head.