To overcome these difficulties we follow in this work an approach described recently . The idea is to not select all but a subset of the process variables as optimization variables. Next, the structure of the MESH equations is exploited to calculate values for all process variables and presolve a part of the equations. This way one arrives at a comparably small optimization problem with favorable convergence properties.