Initialization of optional values depending on other system models

As in #539 and other PRs it is a problem to initialise other systems based on the model selection:

I suggest adding in semidiscretization:

for_each(semi) do system_A 
   for_each(semi) do system_B 
      # pseudocode
      if system_A == system_B continue end
      initialize!(system_A, system_B) 
   end
end

And handling all initialization that is based on other systems in the dispatch of initialize!(system_A, system_B).

Examples of this are:

• Initialization concerning surface normals and surface tension
• Corrections

Start at the semidiscretization logic and the initialize!(system_A, system_B) dispatch described in the issue. Determine how dependent initialization should cover surface normals, surface tension, and corrections, then verify that initialization works for every pair of distinct systems.