WIP: Debug windows CI for YAC coupling #691

BenjaminRodenberg · 2025-08-28T06:26:19Z

I open this PR to trigger the CI pipeline since I do not have access to a Windows machine for local testing.

I try to find the reason for the test 6 - Contact failing in #682 @ ac2d8cb

…mpi-handshake as fallback if working without YAC.

This reverts commit 994152c.

…solution)

…n devel)

BenjaminRodenberg · 2025-09-02T08:41:13Z

I reduced the diff to a minimum compared to https://github.com/ElmerCSC/elmerfem/tree/devel @ f6798c3 and the test 6 - Contact still fails.

tzwinger · 2025-09-03T05:45:48Z

@raback we would need output from the Windows runner to identify the issue. Don't merge, if fixed, since this is just a dummy PR to get rid of failing windows testing

mmuetzel · 2025-10-07T08:50:45Z

Looking at recent CI runs for the build-windows-mingw workflow, the test "6 - Contact" is failing intermittently with or without the changes related to YAC. It is passing (most of the time or always?) when the test is re-run.
In case that should matter: These re-runs are serial (compared to the initial tests which are executed in parallel).
The fact that the test sometimes fails and sometimes passes hints that there might be an issue with threading or memory layout. But it might also be something else.

@tzwinger: Which kind of output are you looking for. If tests fail, the step "re-run tests" lists the directory content of the affected tests and the content of the test-stderr*.log and test-stdout*.log files before re-running the previously failed tests.
Does that include the output you are looking for?

tzwinger · 2025-10-13T10:44:49Z

@mmuetzel that is what @raback reported on that issue #682 (comment)

mmuetzel · 2025-10-13T10:48:10Z

Thanks for clarifying.

So, is there anything missing in this PR? Or can it be closed?

tzwinger · 2025-10-13T19:03:14Z

Well, we still have no real clue why this one test appears to have latency in writing the TEST_PASSED file (as suggested by @raback as the cause of it to fail in a random way)

mmuetzel · 2025-10-14T09:40:55Z

Maybe, @raback was looking at a different CI run than the one here.
At least in the last run here, the test did not fail because the TEST_PASSED file was written too late. But the test actually segmentation faulted:
https://github.com/ElmerCSC/elmerfem/actions/runs/17396535189/job/49379617033#step:10:71

Copying the output to stderr and stdout here in case it is difficult to find in the CI logs:

  ---- Content of elmerice/Tests/Contact/test-stderr.log ----
  
  Program received signal SIGSEGV: Segmentation fault - invalid memory reference.
  
  Backtrace for this error:
  #0  0x7aa18b24
  #1  0x7a981e33
  #2  0x31f52201
  #3  0xa198ad47
  #4  0xa1c448ee
  #5  0xa1bd192d
  #6  0xa1c438ed
  #7  0x723aaf7f
  #8  0x723c7fbc
  #9  0x7180324c
  #10  0x723d92e8
  #11  0x72628b76
  #12  0x71a88f9e
  #13  0x723ecee2
  #14  0x72406351
  #15  0x72408502
  #16  0x724098df
  #17  0x727b9238
  #18  0x727c266b
  #19  0x31f515f8
  #20  0x31f52b3f
  #21  0x31f510c8
  #22  0x31f51415
  #23  0xa02a4caf
  #24  0xa1c1edca
  #25  0xffffffff
  ---- Content of elmerice/Tests/Contact/test-stdout.log ----
  ELMER SOLVER (v 9.0) STARTED AT: 2025/09/02 07:47:03
  ParCommInit: Initialize #PEs: 1
  MAIN: OMP_NUM_THREADS not set. Using only 1 thread per task.
  MAIN: 
  MAIN: =============================================================
  MAIN: ElmerSolver finite element software, Welcome!
  MAIN: This program is free software licensed under (L)GPL
  MAIN: Copyright 1st April 1995 - , CSC - IT Center for Science Ltd.
  MAIN: Webpage http://www.csc.fi/elmer, Email elmeradm@csc.fi
  MAIN: Version: 9.0 (Rev: babb0ed, Compiled: 2025-09-02)
  MAIN:  Running one task without MPI parallelization.
  MAIN:  Running with just one thread per task.
  MAIN:  MUMPS library linked in.
  MAIN:  CHOLMOD library linked in.
  MAIN:  Lua interpreter linked in.
  MAIN:  Zoltan library linked in.
  MAIN: =============================================================
  MAIN: 
  MAIN: 
  MAIN: -------------------------------------
  MAIN: Reading Model: cavity.sif
  LoadInputFile: Scanning input file: cavity.sif
  LoadInputFile: Scanning only size info
  LoadInputFile: First time visiting
  LoadInputFile: Reading base load of sif file
  Header
  Mesh DB "." "ice"
  End
  Constants
  Water Density = Real 1.00413490931e-18
   Real 1.00413490931e-18
  End
  Simulation
  Coordinate System  = Cartesian 2D
   Cartesian 2D
  Simulation Type = transient
   transient
  Timestepping Method = "bdf"
   "bdf"
  BDF Order = 1
   1
  Timestep Intervals = 2
   2
  Output Intervals = 10
   10
  Timestep Sizes = 0.005
   0.005
  Steady State Max Iterations = 1 
   1 
  Steady State Min Iterations = 1
   1
  Initialize Dirichlet Conditions = Logical False
   Logical False
  Output File = "cavity.result
   "cavity.result
  Post File = "cavity.vtu
   "cavity.vtu
  max output level = 5
   5
  End
  Body 1
  Name = "ice"
   "ice"
  Equation = 1
   1
  Body Force = 1
   1
  Material = 1
   1
  Initial Condition = 1
   1
  End
  Body 2
  Name= "cavity free surface"
   "cavity free surface"
  Equation = 2
   2
  Material = 1
   1
  Body Force = 2
   2
  Initial Condition = 2
   2
  End
  Initial Condition 1
  Pressure = Real 0.0
   Real 0.0
  Velocity 1 = Real 0.0
   Real 0.0
  Velocity 2 = Real 0.0
   Real 0.0
  Mesh Update 1 = Real 0.0
   Real 0.0
  Mesh Update 2 = real 0.0
   real 0.0
  End
  Initial Condition 2
  Zs Bottom = Variable Coordinate 1
   Variable Coordinate 1
  Real Procedure "ElmerIceUSF" "ZsBottomIni"
  Bedrock = Variable Coordinate 1
   Variable Coordinate 1
  Real
  0.0 0.0
  2.0 0.1
  4.0 -0.3
  10.0 0.0
  End
  End
  Body Force 1
  Flow BodyForce 1 = Real 0.0
   Real 0.0
  Flow BodyForce 2 = Real 0.0
   Real 0.0
  Flow Solver Name = String Flow Solution
   String Flow Solution
  Mesh Update 1 = real 0.0
   real 0.0
  End
  Body Force 2
  Zs Bottom Accumulation = Real 0.0
   Real 0.0
  End
  Material 1
  Density = Real 9.03721418379e-19
   Real 9.03721418379e-19
  Viscosity Model = String "glen"
   String "glen"
  Viscosity = 1.0 
   1.0 
  Glen Exponent = Real 3.0
   Real 3.0
  Limit Temperature = Real -10.0
   Real -10.0
  Rate Factor 1 = Real 9.125353404e+12
   Real 9.125353404e+12
  Rate Factor 2 = Real 7.6601655936e+23
   Real 7.6601655936e+23
  Activation Energy 1 = Real 60e3
   Real 60e3
  Activation Energy 2 = Real 115e3
   Real 115e3
  Glen Enhancement Factor = Real 1.0
   Real 1.0
  Critical Shear Rate = Real 1.0e-15
   Real 1.0e-15
  Constant Temperature = Real -1.0
   Real -1.0
  Min Zs Bottom = Equals Bedrock
   Equals Bedrock
  Max Zs Bottom = Real 1.0e6
   Real 1.0e6
  Mesh Youngs Modulus = Real 1.0
   Real 1.0
  Mesh Poisson Ratio = real 0.3
   real 0.3
  Cauchy = Logical True
   Logical True
  End
  Solver 1
  Exec Solver = Before All
   Before All
  Procedure = "ElmerIceSolvers" "GroundedSolver"
   "ElmerIceSolvers" "GroundedSolver"
  Equation = "GroundedMaskInit"
   "GroundedMaskInit"
  Variable = "GroundedMask"
   "GroundedMask"
  Variable DOFs = 1
   1
  Toler = Real 0.0001
   Real 0.0001
  Bedrock Variable = String "Bedrock"
   String "Bedrock"
  End
  Solver2
  Equation = "Mesh Update"
   "Mesh Update"
  Variable = "Mesh Update"
   "Mesh Update"
  Variable DOFs = 2
   2
  Linear System Solver = "Direct"
   "Direct"
  Linear System Direct Method = umfpack
   umfpack
  Steady State Convergence Tolerance = 1.0e-04
   1.0e-04
  End
  Solver 3
  Equation = "GroundedMask"
   "GroundedMask"
  Procedure = "ElmerIceSolvers" "GroundedSolver"
   "ElmerIceSolvers" "GroundedSolver"
  Variable = "GroundedMask"
   "GroundedMask"
  Variable DOFs = 1
   1
  Toler = Real 0.0001
   Real 0.0001
  Bedrock Variable = String "Bedrock"
   String "Bedrock"
  End
  Solver 4
  Equation = "NormalVector"
   "NormalVector"
  Procedure = "ElmerIceSolvers" "ComputeNormalSolver"
   "ElmerIceSolvers" "ComputeNormalSolver"
  Variable = String "Normal Vector"
   String "Normal Vector"
  Variable DOFs = 2
   2
  Exported Variable 1 = BodyMask
   BodyMask
  Exported Variable 1 DOFs = 1
   1
  ComputeAll = Logical False
   Logical False
  Optimize Bandwidth = Logical False
   Logical False
  End
  Solver 5
  Equation = Fw
   Fw
  Procedure = "ElmerIceSolvers" "GetHydrostaticLoads"
   "ElmerIceSolvers" "GetHydrostaticLoads"
  Variable = Fw[Fwater:2]
   Fw[Fwater:2]
  Variable DOFs = 2
   2
  End
  Solver 6
  Equation = "Navier-Stokes"
   "Navier-Stokes"
  Linear System Solver = Direct
   Direct
  Linear System Direct Method = umfpack
   umfpack
  Nonlinear System Max Iterations = 50
   50
  Nonlinear System Convergence Tolerance  = 1.0e-6
   1.0e-6
  Nonlinear System Newton After Iterations = 51
   51
  Nonlinear System Newton After Tolerance = 1.0e-07
   1.0e-07
  Nonlinear System Relaxation Factor = 1.00
   1.00
  Steady State Convergence Tolerance = Real 1.0e-4
   Real 1.0e-4
  Stabilization Method = String Bubbles
   String Bubbles
  Exported Variable 1 = Flow Solution Loads[Stress Vector:2 CEQ Residual:1]
   Flow Solution Loads[Stress Vector:2 CEQ Residual:1]
  Calculate Loads = Logical True
   Logical True
  Flow Model = String 'Stokes'
   String 'Stokes'
  End
  Solver 7
  Equation = Sij
   Sij
  Variable = -nooutput "Sij"
   -nooutput "Sij"
  Variable DOFs = 1
   1
  Exported Variable 1 = ViscousStress
   ViscousStress
  Exported Variable 1 DOFs = 4
   4
  Procedure = "ElmerIceSolvers" "ComputeDevStress"
   "ElmerIceSolvers" "ComputeDevStress"
  Flow Solver Name = String "Flow Solution"
   String "Flow Solution"
  Stress Variable Name = String "ViscousStress"
   String "ViscousStress"
  Linear System Solver = Direct
   Direct
  Linear System Direct Method = umfpack
   umfpack
  End
  Solver 8
  Equation = "Free Surface cavity"
   "Free Surface cavity"
  Variable = "Zs Bottom"
   "Zs Bottom"
  Variable DOFS =  1
    1
  Exported Variable 1 = "Zs Bottom Residual"
   "Zs Bottom Residual"
  Exported Variable 1 DOFs = 1
   1
  Exported Variable 2 = "Bedrock"
   "Bedrock"
  CheckKeyword: Found keyword type assuming suffix 1 for: exported variable 2
  CheckKeyword: Found keyword type assuming suffix 1 for: exported variable 2
  Exported Variable 2 DOFs = 1
   1
  Procedure = "FreeSurfaceSolver" "FreeSurfaceSolver"
   "FreeSurfaceSolver" "FreeSurfaceSolver"
  Before Linsolve = "EliminateDirichlet" "EliminateDirichlet"
   "EliminateDirichlet" "EliminateDirichlet"
  Linear System Solver = Iterative
   Iterative
  Linear System Max Iterations = 100
   100
  Linear System Iterative Method = BiCGStab
   BiCGStab
  Linear System Preconditioning = ILU0
   ILU0
  Linear System Convergence Tolerance = Real 1.0e-9
   Real 1.0e-9
  Linear System Abort Not Converged = False
   False
  Linear System Residual Output = 1
   1
  Nonlinear System Max Iterations = 100
   100
  Nonlinear System Convergence Tolerance  = 1.0e-6
   1.0e-6
  Nonlinear System Relaxation Factor = 1.00
   1.00
  Steady State Convergence Tolerance = 1.0e-03
   1.0e-03
  Stabilization Method = Stabilized
   Stabilized
  Apply Dirichlet = Logical True
   Logical True
  Relaxation Factor = Real 1.0
   Real 1.0
  Use Accumulation = Logical True
   Logical True
  Normal Flux = Logical False
   Logical False
  End
  Solver 9
  Exec Solver = After TimeStep
   After TimeStep
  Equation = SaveScalars
   SaveScalars
  Procedure = File "SaveData" "SaveScalars"
   File "SaveData" "SaveScalars"
  Filename = "results.dat"
   "results.dat"
  Show Norm = True
   True
  Show Norm Index = 3
   3
  Variable 1 = "Stress Vector 2"
   "Stress Vector 2"
  Operator 1 = "boundary min"
   "boundary min"
  Operator 2 = "boundary max"
   "boundary max"
  CheckKeyword: Found keyword type assuming suffix 1 for: operator 2
  CheckKeyword: Found keyword type assuming suffix 1 for: operator 2
  Operator 3 = "boundary mean"
   "boundary mean"
  CheckKeyword: Found keyword type assuming suffix 1 for: operator 3
  CheckKeyword: Found keyword type assuming suffix 1 for: operator 3
  End
  Equation 1
  Active Solvers (5) = 2 4 6 7 9
  Size 5  = 2 4 6 7 9
   2 4 6 7 9
  End
  Equation 2
  Active Solvers(4) = 1 3 5 8
  Size 4  = 1 3 5 8
   1 3 5 8
  Flow Solution Name = String "Flow Solution"
   String "Flow Solution"
  Convection = String Computed
   String Computed
  End
  Boundary Condition 1
  Name = "cavity"
   "cavity"
  Target Boundaries = 1
   1
  Body Id = 2
   2
  Normal-Tangential Velocity = Logical True
   Logical True
  Flow Force BC = Logical True
   Logical True
  Test Contact Tolerance = Real 1.0e-3
   Real 1.0e-3
  Grounding Line Definition = String "Last Grounded"
   String "Last Grounded"
  Zs Bottom = Equals Bedrock
   Equals Bedrock
  Zs Bottom Condition = Variable Coordinate 1
   Variable Coordinate 1
  Real MATC "2.0 - tx"
  Slip Coefficient 2 = Variable Coordinate 1
   Variable Coordinate 1
  Real Procedure "ElmerIceUSF" "SlidCoef_Contact"
  Sliding Law = String weertman
   String weertman
  Weertman Friction Coefficient = Real 0.0 
   Real 0.0 
  Weertman Exponent = Real 1.0
   Real 1.0
  Weertman Linear Velocity = Real 1.0
   Real 1.0
  Velocity 1 = Real 0.0
   Real 0.0
  Velocity 1 Condition = Variable GroundedMask
   Variable GroundedMask
  Real MATC "tx + 0.5"
  External Pressure = Real -1
   Real -1
  ComputeNormal = Logical True
   Logical True
  ComputeNormal Condition = Variable GroundedMask
   Variable GroundedMask
  Real MATC "tx + 0.5"
  Mesh Update 1 = Real 0.0e0
   Real 0.0e0
  Mesh Update 2 = Variable Zs Bottom
   Variable Zs Bottom
  Real Procedure "ElmerIceUSF" "ZsBottomMzsIni"
  Save Scalars = Logical True
   Logical True
  End
  Boundary Condition 2
  Name = "ice right"
   "ice right"
  Target Boundaries = 2
   2
  Periodic BC = 4
   4
  Periodic BC Translate(2) = Real 10.0 0.0
  Size 2  = Real 10.0 0.0
   Real 10.0 0.0
  Periodic BC Velocity 1  = Logical True
   Logical True
  Periodic BC Velocity 2  = Logical True
   Logical True
  Periodic BC Pressure  = Logical True
   Logical True
  Periodic BC Zs Bottom = Logical True
   Logical True
  Periodic BC Mesh Update 2 = Logical True
   Logical True
  End
  Boundary Condition 3
  Name = "ice top"
   "ice top"
  Target Boundaries = 3
   3
  Velocity 1 = Real 10
   Real 10
  Flow Force BC = Logical True
   Logical True
  External Pressure = Real -1.2
   Real -1.2
  Mesh Update 1 = real 0.0
   real 0.0
  Mesh Update 2 = real 0.0
   real 0.0
  End
  Boundary Condition 4
  Name = "ice left"
   "ice left"
  Target Boundaries = 4
   4
  End
  Solver 6 :: Reference Norm = 7.0439455
   Reference Norm = 7.0439455
   7.0439455
  Solver 6 :: Reference Norm Tolerance = 1.0E-05
   Reference Norm Tolerance = 1.0E-05
   1.0E-05
  Solver 9 :: Reference Norm = 1.19528960E-01
   Reference Norm = 1.19528960E-01
   1.19528960E-01
  Solver 9 :: Reference Norm Tolerance = 1.0E-03
   Reference Norm Tolerance = 1.0E-03
   1.0E-03
  LoadInputFile: Loading input file: cavity.sif
  LoadInputFile: Reading base load of sif file
  check keywords warn
  echo on
  Header
  Mesh DB "." "ice"
  End
  Constants
  Water Density = Real 1.00413490931e-18
   Real 1.00413490931e-18
  CheckKeyword:  Unlisted keyword: [water density] in section: [constants]
  End
  Simulation
  Coordinate System  = Cartesian 2D
   Cartesian 2D
  Simulation Type = transient
   transient
  Timestepping Method = "bdf"
   "bdf"
  BDF Order = 1
   1
  Timestep Intervals = 2
   2
  Output Intervals = 10
   10
  Timestep Sizes = 0.005
   0.005
  Steady State Max Iterations = 1 
   1 
  Steady State Min Iterations = 1
   1
  Initialize Dirichlet Conditions = Logical False
   Logical False
  Output File = "cavity.result
   "cavity.result
  Post File = "cavity.vtu
   "cavity.vtu
  max output level = 5
   5
  End
  Body 1
  Name = "ice"
   "ice"
  Equation = 1
   1
  Body Force = 1
   1
  Material = 1
   1
  Initial Condition = 1
   1
  End
  Body 2
  Name= "cavity free surface"
   "cavity free surface"
  Equation = 2
   2
  Material = 1
   1
  Body Force = 2
   2
  Initial Condition = 2
   2
  End
  Initial Condition 1
  Pressure = Real 0.0
   Real 0.0
  Velocity 1 = Real 0.0
   Real 0.0
  Velocity 2 = Real 0.0
   Real 0.0
  Mesh Update 1 = Real 0.0
   Real 0.0
  Mesh Update 2 = real 0.0
   real 0.0
  End
  Initial Condition 2
  Zs Bottom = Variable Coordinate 1
   Variable Coordinate 1
  Real Procedure "ElmerIceUSF" "ZsBottomIni"
  CheckKeyword:  Unlisted keyword: [zs bottom] in section: [initial condition 2]
  Loading user function library: [ElmerIceUSF]...[ZsBottomIni]
  Bedrock = Variable Coordinate 1
   Variable Coordinate 1
  Real
  CheckKeyword:  Unlisted keyword: [bedrock] in section: [initial condition 2]
  0.0 0.0
  2.0 0.1
  4.0 -0.3
  10.0 0.0
  End
  End
  Body Force 1
  Flow BodyForce 1 = Real 0.0
   Real 0.0
  Flow BodyForce 2 = Real 0.0
   Real 0.0
  Flow Solver Name = String Flow Solution
   String Flow Solution
  CheckKeyword:  Unlisted keyword: [flow solver name] in section: [body force 1]
  Mesh Update 1 = real 0.0
   real 0.0
  CheckKeyword:  Unlisted keyword: [mesh update 1] in section: [body force 1]
  End
  Body Force 2
  Zs Bottom Accumulation = Real 0.0
   Real 0.0
  CheckKeyword:  Unlisted keyword: [zs bottom accumulation] in section: [body force 2]
  End
  Material 1
  Density = Real 9.03721418379e-19
   Real 9.03721418379e-19
  Viscosity Model = String "glen"
   String "glen"
  Viscosity = 1.0 
   1.0 
  Glen Exponent = Real 3.0
   Real 3.0
  Limit Temperature = Real -10.0
   Real -10.0
  CheckKeyword:  Unlisted keyword: [limit temperature] in section: [material 1]
  Rate Factor 1 = Real 9.125353404e+12
   Real 9.125353404e+12
  Rate Factor 2 = Real 7.6601655936e+23
   Real 7.6601655936e+23
  Activation Energy 1 = Real 60e3
   Real 60e3
  Activation Energy 2 = Real 115e3
   Real 115e3
  Glen Enhancement Factor = Real 1.0
   Real 1.0
  Critical Shear Rate = Real 1.0e-15
   Real 1.0e-15
  Constant Temperature = Real -1.0
   Real -1.0
  CheckKeyword:  Unlisted keyword: [constant temperature] in section: [material 1]
  Min Zs Bottom = Equals Bedrock
   Equals Bedrock
  Max Zs Bottom = Real 1.0e6
   Real 1.0e6
  CheckKeyword:  Unlisted keyword: [max zs bottom] in section: [material 1]
  Mesh Youngs Modulus = Real 1.0
   Real 1.0
  Mesh Poisson Ratio = real 0.3
   real 0.3
  Cauchy = Logical True
   Logical True
  CheckKeyword:  Unlisted keyword: [cauchy] in section: [material 1]
  End
  Solver 1
  Exec Solver = Before All
   Before All
  Procedure = "ElmerIceSolvers" "GroundedSolver"
   "ElmerIceSolvers" "GroundedSolver"
  Equation = "GroundedMaskInit"
   "GroundedMaskInit"
  Variable = "GroundedMask"
   "GroundedMask"
  Variable DOFs = 1
   1
  Toler = Real 0.0001
   Real 0.0001
  CheckKeyword:  Unlisted keyword: [toler] in section: [solver 1]
  Bedrock Variable = String "Bedrock"
   String "Bedrock"
  CheckKeyword:  Unlisted keyword: [bedrock variable] in section: [solver 1]
  End
  Solver2
  Equation = "Mesh Update"
   "Mesh Update"
  Variable = "Mesh Update"
   "Mesh Update"
  Variable DOFs = 2
   2
  Linear System Solver = "Direct"
   "Direct"
  Linear System Direct Method = umfpack
   umfpack
  Steady State Convergence Tolerance = 1.0e-04
   1.0e-04
  End
  Solver 3
  Equation = "GroundedMask"
   "GroundedMask"
  Procedure = "ElmerIceSolvers" "GroundedSolver"
   "ElmerIceSolvers" "GroundedSolver"
  Variable = "GroundedMask"
   "GroundedMask"
  Variable DOFs = 1
   1
  Toler = Real 0.0001
   Real 0.0001
  CheckKeyword:  Unlisted keyword: [toler] in section: [solver 3]
  Bedrock Variable = String "Bedrock"
   String "Bedrock"
  CheckKeyword:  Unlisted keyword: [bedrock variable] in section: [solver 3]
  End
  Solver 4
  Equation = "NormalVector"
   "NormalVector"
  Procedure = "ElmerIceSolvers" "ComputeNormalSolver"
   "ElmerIceSolvers" "ComputeNormalSolver"
  Variable = String "Normal Vector"
   String "Normal Vector"
  Variable DOFs = 2
   2
  Exported Variable 1 = BodyMask
   BodyMask
  Exported Variable 1 DOFs = 1
   1
  ComputeAll = Logical False
   Logical False
  CheckKeyword:  Unlisted keyword: [computeall] in section: [solver 4]
  Optimize Bandwidth = Logical False
   Logical False
  End
  Solver 5
  Equation = Fw
   Fw
  Procedure = "ElmerIceSolvers" "GetHydrostaticLoads"
   "ElmerIceSolvers" "GetHydrostaticLoads"
  Variable = Fw[Fwater:2]
   Fw[Fwater:2]
  Variable DOFs = 2
   2
  End
  Solver 6
  Equation = "Navier-Stokes"
   "Navier-Stokes"
  Linear System Solver = Direct
   Direct
  Linear System Direct Method = umfpack
   umfpack
  Nonlinear System Max Iterations = 50
   50
  Nonlinear System Convergence Tolerance  = 1.0e-6
   1.0e-6
  Nonlinear System Newton After Iterations = 51
   51
  Nonlinear System Newton After Tolerance = 1.0e-07
   1.0e-07
  Nonlinear System Relaxation Factor = 1.00
   1.00
  Steady State Convergence Tolerance = Real 1.0e-4
   Real 1.0e-4
  Stabilization Method = String Bubbles
   String Bubbles
  Exported Variable 1 = Flow Solution Loads[Stress Vector:2 CEQ Residual:1]
   Flow Solution Loads[Stress Vector:2 CEQ Residual:1]
  Calculate Loads = Logical True
   Logical True
  Flow Model = String 'Stokes'
   String 'Stokes'
  End
  Solver 7
  Equation = Sij
   Sij
  Variable = -nooutput "Sij"
   -nooutput "Sij"
  Variable DOFs = 1
   1
  Exported Variable 1 = ViscousStress
   ViscousStress
  Exported Variable 1 DOFs = 4
   4
  Procedure = "ElmerIceSolvers" "ComputeDevStress"
   "ElmerIceSolvers" "ComputeDevStress"
  Flow Solver Name = String "Flow Solution"
   String "Flow Solution"
  CheckKeyword:  Unlisted keyword: [flow solver name] in section: [solver 7]
  Stress Variable Name = String "ViscousStress"
   String "ViscousStress"
  CheckKeyword:  Unlisted keyword: [stress variable name] in section: [solver 7]
  Linear System Solver = Direct
   Direct
  Linear System Direct Method = umfpack
   umfpack
  End
  Solver 8
  Equation = "Free Surface cavity"
   "Free Surface cavity"
  Variable = "Zs Bottom"
   "Zs Bottom"
  Variable DOFS =  1
    1
  Exported Variable 1 = "Zs Bottom Residual"
   "Zs Bottom Residual"
  Exported Variable 1 DOFs = 1
   1
  Exported Variable 2 = "Bedrock"
   "Bedrock"
  CheckKeyword: Found keyword type assuming suffix 1 for: exported variable 2
  CheckKeyword: Found keyword type assuming suffix 1 for: exported variable 2
  Exported Variable 2 DOFs = 1
   1
  Procedure = "FreeSurfaceSolver" "FreeSurfaceSolver"
   "FreeSurfaceSolver" "FreeSurfaceSolver"
  Before Linsolve = "EliminateDirichlet" "EliminateDirichlet"
   "EliminateDirichlet" "EliminateDirichlet"
  Linear System Solver = Iterative
   Iterative
  Linear System Max Iterations = 100
   100
  Linear System Iterative Method = BiCGStab
   BiCGStab
  Linear System Preconditioning = ILU0
   ILU0
  Linear System Convergence Tolerance = Real 1.0e-9
   Real 1.0e-9
  Linear System Abort Not Converged = False
   False
  Linear System Residual Output = 1
   1
  Nonlinear System Max Iterations = 100
   100
  Nonlinear System Convergence Tolerance  = 1.0e-6
   1.0e-6
  Nonlinear System Relaxation Factor = 1.00
   1.00
  Steady State Convergence Tolerance = 1.0e-03
   1.0e-03
  Stabilization Method = Stabilized
   Stabilized
  Apply Dirichlet = Logical True
   Logical True
  Relaxation Factor = Real 1.0
   Real 1.0
  CheckKeyword:  Unlisted keyword: [relaxation factor] in section: [solver 8]
  Use Accumulation = Logical True
   Logical True
  Normal Flux = Logical False
   Logical False
  End
  Solver 9
  Exec Solver = After TimeStep
   After TimeStep
  Equation = SaveScalars
   SaveScalars
  Procedure = File "SaveData" "SaveScalars"
   File "SaveData" "SaveScalars"
  Filename = "results.dat"
   "results.dat"
  Show Norm = True
   True
  Show Norm Index = 3
   3
  Variable 1 = "Stress Vector 2"
   "Stress Vector 2"
  Operator 1 = "boundary min"
   "boundary min"
  Operator 2 = "boundary max"
   "boundary max"
  CheckKeyword: Found keyword type assuming suffix 1 for: operator 2
  CheckKeyword: Found keyword type assuming suffix 1 for: operator 2
  Operator 3 = "boundary mean"
   "boundary mean"
  CheckKeyword: Found keyword type assuming suffix 1 for: operator 3
  CheckKeyword: Found keyword type assuming suffix 1 for: operator 3
  End
  Equation 1
  Active Solvers (5) = 2 4 6 7 9
  Size 5  = 2 4 6 7 9
   2 4 6 7 9
  End
  Equation 2
  Active Solvers(4) = 1 3 5 8
  Size 4  = 1 3 5 8
   1 3 5 8
  Flow Solution Name = String "Flow Solution"
   String "Flow Solution"
  Convection = String Computed
   String Computed
  End
  Boundary Condition 1
  Name = "cavity"
   "cavity"
  Target Boundaries = 1
   1
  Body Id = 2
   2
  Normal-Tangential Velocity = Logical True
   Logical True
  Flow Force BC = Logical True
   Logical True
  Test Contact Tolerance = Real 1.0e-3
   Real 1.0e-3
  CheckKeyword:  Unlisted keyword: [test contact tolerance] in section: [boundary condition 1]
  Grounding Line Definition = String "Last Grounded"
   String "Last Grounded"
  CheckKeyword:  Unlisted keyword: [grounding line definition] in section: [boundary condition 1]
  Zs Bottom = Equals Bedrock
   Equals Bedrock
  Zs Bottom Condition = Variable Coordinate 1
   Variable Coordinate 1
  Real MATC "2.0 - tx"
  CheckKeyword:  Unlisted keyword: [zs bottom condition] in section: [boundary condition 1]
  Loading user function library: [ElmerIceUSF]...[SlidCoef_Contact]
  Slip Coefficient 2 = Variable Coordinate 1
   Variable Coordinate 1
  Real Procedure "ElmerIceUSF" "SlidCoef_Contact"
  Sliding Law = String weertman
   String weertman
  CheckKeyword:  Unlisted keyword: [sliding law] in section: [boundary condition 1]
  Weertman Friction Coefficient = Real 0.0 
   Real 0.0 
  CheckKeyword:  Unlisted keyword: [weertman friction coefficient] in section: [boundary condition 1]
  Weertman Exponent = Real 1.0
   Real 1.0
  CheckKeyword:  Unlisted keyword: [weertman exponent] in section: [boundary condition 1]
  Weertman Linear Velocity = Real 1.0
   Real 1.0
  CheckKeyword:  Unlisted keyword: [weertman linear velocity] in section: [boundary condition 1]
  Velocity 1 = Real 0.0
   Real 0.0
  Velocity 1 Condition = Variable GroundedMask
   Variable GroundedMask
  Real MATC "tx + 0.5"
  External Pressure = Real -1
   Real -1
  ComputeNormal = Logical True
   Logical True
  CheckKeyword:  Unlisted keyword: [computenormal] in section: [boundary condition 1]
  ComputeNormal Condition = Variable GroundedMask
   Variable GroundedMask
  Real MATC "tx + 0.5"
  CheckKeyword:  Unlisted keyword: [computenormal condition] in section: [boundary condition 1]
  Loading user function library: [ElmerIceUSF]...[ZsBottomMzsIni]
  Mesh Update 1 = Real 0.0e0
   Real 0.0e0
  Mesh Update 2 = Variable Zs Bottom
   Variable Zs Bottom
  Real Procedure "ElmerIceUSF" "ZsBottomMzsIni"
  Save Scalars = Logical True
   Logical True
  End
  Boundary Condition 2
  Name = "ice right"
   "ice right"
  Target Boundaries = 2
   2
  Periodic BC = 4
   4
  Periodic BC Translate(2) = Real 10.0 0.0
  Size 2  = Real 10.0 0.0
   Real 10.0 0.0
  Periodic BC Velocity 1  = Logical True
   Logical True
  Periodic BC Velocity 2  = Logical True
   Logical True
  Periodic BC Pressure  = Logical True
   Logical True
  Periodic BC Zs Bottom = Logical True
   Logical True
  CheckKeyword:  Unlisted keyword: [periodic bc zs bottom] in section: [boundary condition 2]
  Periodic BC Mesh Update 2 = Logical True
   Logical True
  End
  Boundary Condition 3
  Name = "ice top"
   "ice top"
  Target Boundaries = 3
   3
  Velocity 1 = Real 10
   Real 10
  Flow Force BC = Logical True
   Logical True
  External Pressure = Real -1.2
   Real -1.2
  Mesh Update 1 = real 0.0
   real 0.0
  Mesh Update 2 = real 0.0
   real 0.0
  End
  Boundary Condition 4
  Name = "ice left"
   "ice left"
  Target Boundaries = 4
   4
  End
  Solver 6 :: Reference Norm = 7.0439455
   Reference Norm = 7.0439455
   7.0439455
  Solver 6 :: Reference Norm Tolerance = 1.0E-05
   Reference Norm Tolerance = 1.0E-05
   1.0E-05
  Solver 9 :: Reference Norm = 1.19528960E-01
   Reference Norm = 1.19528960E-01
   1.19528960E-01
  Solver 9 :: Reference Norm Tolerance = 1.0E-03
   Reference Norm Tolerance = 1.0E-03
   1.0E-03
  LoadInputFile: Number of BCs: 4
  LoadInputFile: Number of Body Forces: 2
  LoadInputFile: Number of Initial Conditions: 2
  LoadInputFile: Number of Materials: 1
  LoadInputFile: Number of Equations: 2
  LoadInputFile: Number of Solvers: 9
  LoadInputFile: Number of Bodies: 2
  ListTagKeywords: Setting weight for keywords!
  ListTagKeywords: No parameters width suffix: normalize by area
  ListTagKeywords: Setting weight for keywords!
  ListTagKeywords: No parameters width suffix: normalize by volume
  ElmerAsciiMesh: Base mesh name: ./ice
  OverlayInterfaceMeshes: Scaling with:    1.000000E+00   1.000000E+00   1.000000E+00
  OverlayInterfaceMeshes: Translation:    1.000000E+01   0.000000E+00   0.000000E+00
  BuildQuandrantTree: Start
  BuildQuandrantTree: Ready
  PeriodicProjector: Elapsed REAL time:     0.0000 (s)
  LoadMesh: Elapsed REAL time:     0.0160 (s)
  MAIN: -------------------------------------
  AddVtuOutputSolverHack: Adding ResultOutputSolver to write VTU output in file: cavity
  OptimizeBandwidth: Initial bandwidth for groundedmaskinit: 100
  OptimizeBandwidth: Optimized bandwidth for groundedmaskinit: 4
  OptimizeBandwidth: Initial bandwidth for mesh update: 303
  OptimizeBandwidth: Optimized bandwidth for mesh update: 66
  OptimizeBandwidth: Initial bandwidth for groundedmask: 100
  OptimizeBandwidth: Optimized bandwidth for groundedmask: 4
  OptimizeBandwidth: Initial bandwidth for normalvector: 303
  OptimizeBandwidth: Initial bandwidth for fw: 100
  OptimizeBandwidth: Optimized bandwidth for fw: 4
  OptimizeBandwidth: Initial bandwidth for navier-stokes: 303
  OptimizeBandwidth: Optimized bandwidth for navier-stokes: 66
  OptimizeBandwidth: Initial bandwidth for sij: 303
  OptimizeBandwidth: Optimized bandwidth for sij: 66
  OptimizeBandwidth: Initial bandwidth for free surface cavity: 100
  OptimizeBandwidth: Optimized bandwidth for free surface cavity: 4
  MAIN: Number of timesteps to be saved: 2
  GroundedSolver: Computing grounded mask from geometry
  GroundedSolver: Memory allocation done.
  GroundedSolver: No front variable defined. Some basal frontal nodes may be left with GroundedMask=1
  GroundedSolver:  Grounding Line, x   2.0000000000000000
  GroundedSolver:  Grounding Line, x   8.0000000000000000
  GroundedSolver: Done
  MAIN: 
  MAIN: -------------------------------------
  MAIN: Time: 1/2:   5.000E-03
  MAIN: -------------------------------------
  MAIN: 
  MeshSolve: 
  MeshSolve: -------------------------------------
  MeshSolve: MESH UPDATE SOLVER:
  MeshSolve: -------------------------------------
  MeshSolve: 
  MeshSolve: Starting assembly...
  MeshSolve: Assembly done
  MeshSolve: Set boundaries done
  SolveLinearSystem: Solution trivially zero!
  GroundedSolver: Computing grounded mask from geometry
  GroundedSolver: No front variable defined. Some basal frontal nodes may be left with GroundedMask=1
  GroundedSolver:  Grounding Line, x   2.0000000000000000
  GroundedSolver:  Grounding Line, x   8.0000000000000000
  GroundedSolver: Done
  ComputeNormalSolver: Computing Normal Vector for Nodes
  ComputeNormalSolver: End
  GetHydrostaticLoads: Memory allocation done.
  Update Periodic Nodes for: fwater 1: Start
  UpdatePeriodicNodes: All Done
  Update Periodic Nodes for: fwater 2: Start
  UpdatePeriodicNodes: All Done
  GetHydrostaticLoads: End
  FlowSolve: 
  FlowSolve: 
  FlowSolve: -------------------------------------
  FlowSolve:  NAVIER-STOKES ITERATION           1
  FlowSolve: -------------------------------------
  FlowSolve: 
  FlowSolve: Starting Assembly...
  FlowSolve: Assembly:
  FlowSolve: Assembly done
  SlidCoef_Contact: Using default name >Flow Solution<
  SlidCoef_Contact: Grounding Line Defined as last grounded
  SlidCoef_Contact: GL may move by default
  SlidCoef_Contact: If you want to fix the Grounding Line, put the keyword "Grounding line moves" to False
  SlidCoef_Contact: Contact tested for a tolerance of 0.10000000E-02
  SlidCoef_Contact: far inland nodes have the possibility to detach by default
  SlidCoef_Contact: to avoid detachment (when bedrock is well below sea level),
  SlidCoef_Contact: use the keyword "non detachment inland distance" to the distance you wish
  SlidCoef_Contact: This works with the DistanceSolver
  Sliding_Weertman:  Flow Solver Name:Flow Solution SSA: F
  FlowSolve: Dirichlet conditions done
  ComputeChange: NS (ITER=1) (NRM,RELC): (  7.0514976      2.0000000     ) :: navier-stokes
  FlowSolve: iter:    1 Assembly: (s)    0.33    0.33
  FlowSolve: iter:    1 Solve:    (s)    0.06    0.06
  FlowSolve:  Result Norm     :    7.0514975700813389
  FlowSolve:  Relative Change :    2.0000000000000000
  FlowSolve: 
  FlowSolve: 
  FlowSolve: -------------------------------------
  FlowSolve:  NAVIER-STOKES ITERATION           2
  FlowSolve: -------------------------------------
  FlowSolve: 
  FlowSolve: Starting Assembly...
  FlowSolve: Assembly:
  FlowSolve: Assembly done
  FlowSolve: Dirichlet conditions done
  ComputeChange: NS (ITER=2) (NRM,RELC): (  7.0594487     0.11269428E-02 ) :: navier-stokes
  FlowSolve: iter:    2 Assembly: (s)    0.33    0.66
  FlowSolve: iter:    2 Solve:    (s)    0.06    0.13
  FlowSolve:  Result Norm     :    7.0594486848851901
  FlowSolve:  Relative Change :    1.1269428229949692E-003
  FlowSolve: 
  FlowSolve: 
  FlowSolve: -------------------------------------
  FlowSolve:  NAVIER-STOKES ITERATION           3
  FlowSolve: -------------------------------------
  FlowSolve: 
  FlowSolve: Starting Assembly...
  FlowSolve: Assembly:
  FlowSolve: Assembly done
  FlowSolve: Dirichlet conditions done
  ComputeChange: NS (ITER=3) (NRM,RELC): (  7.0540673     0.76258202E-03 ) :: navier-stokes
  FlowSolve: iter:    3 Assembly: (s)    0.32    0.98
  FlowSolve: iter:    3 Solve:    (s)    0.05    0.18
  FlowSolve:  Result Norm     :    7.0540673280817314
  FlowSolve:  Relative Change :    7.6258202399947194E-004
  FlowSolve: 
  FlowSolve: 
  FlowSolve: -------------------------------------
  FlowSolve:  NAVIER-STOKES ITERATION           4
  FlowSolve: -------------------------------------
  FlowSolve: 
  FlowSolve: Starting Assembly...
  FlowSolve: Assembly:
  FlowSolve: Assembly done
  SlidCoef_Contact: FLOW SOLVER HAS SLIGHTLY CONVERGED: look for new basal conditions
  SlidCoef_Contact:  Retreat of the Grounding Line :    2.0000000000000000       0.10000000000000001
  SlidCoef_Contact:  Retreat of the Grounding Line :    8.0000000000000000      -0.10000000000000001
  FlowSolve: Dirichlet conditions done
  ComputeChange: NS (ITER=4) (NRM,RELC): (  7.0509845     0.43712550E-03 ) :: navier-stokes
  FlowSolve: iter:    4 Assembly: (s)    0.33    1.31
  FlowSolve: iter:    4 Solve:    (s)    0.05    0.23
  FlowSolve:  Result Norm     :    7.0509844891890392
  FlowSolve:  Relative Change :    4.3712549696804139E-004
  FlowSolve: 
  FlowSolve: 
  FlowSolve: -------------------------------------
  FlowSolve:  NAVIER-STOKES ITERATION           5
  FlowSolve: -------------------------------------
  FlowSolve: 
  FlowSolve: Starting Assembly...
  FlowSolve: Assembly:
  FlowSolve: Assembly done
  FlowSolve: Dirichlet conditions done
  ComputeChange: NS (ITER=5) (NRM,RELC): (  7.0487638     0.31500385E-03 ) :: navier-stokes
  FlowSolve: iter:    5 Assembly: (s)    0.34    1.65
  FlowSolve: iter:    5 Solve:    (s)    0.05    0.27
  FlowSolve:  Result Norm     :    7.0487637517102186
  FlowSolve:  Relative Change :    3.1500384842033632E-004
  FlowSolve: 
  FlowSolve: 
  FlowSolve: -------------------------------------
  FlowSolve:  NAVIER-STOKES ITERATION           6
  FlowSolve: -------------------------------------
  FlowSolve: 
  FlowSolve: Starting Assembly...
  FlowSolve: Assembly:
  FlowSolve: Assembly done
  FlowSolve: Dirichlet conditions done
  ComputeChange: NS (ITER=6) (NRM,RELC): (  7.0472813     0.21033071E-03 ) :: navier-stokes
  FlowSolve: iter:    6 Assembly: (s)    0.33    1.98
  FlowSolve: iter:    6 Solve:    (s)    0.05    0.32
  FlowSolve:  Result Norm     :    7.0472813361197080
  FlowSolve:  Relative Change :    2.1033071067437563E-004
  FlowSolve: 
  FlowSolve: 
  FlowSolve: -------------------------------------
  FlowSolve:  NAVIER-STOKES ITERATION           7
  FlowSolve: -------------------------------------
  FlowSolve: 
  FlowSolve: Starting Assembly...
  FlowSolve: Assembly:
  FlowSolve: Assembly done
  FlowSolve: Dirichlet conditions done
  ComputeChange: NS (ITER=7) (NRM,RELC): (  7.0462950     0.13996970E-03 ) :: navier-stokes
  FlowSolve: iter:    7 Assembly: (s)    0.31    2.29
  FlowSolve: iter:    7 Solve:    (s)    0.06    0.38
  FlowSolve:  Result Norm     :    7.0462949992761272
  FlowSolve:  Relative Change :    1.3996970252377533E-004
  FlowSolve: 
  FlowSolve: 
  FlowSolve: -------------------------------------
  FlowSolve:  NAVIER-STOKES ITERATION           8
  FlowSolve: -------------------------------------
  FlowSolve: 
  FlowSolve: Starting Assembly...
  FlowSolve: Assembly:
  FlowSolve: Assembly done
  FlowSolve: Dirichlet conditions done
  ComputeChange: NS (ITER=8) (NRM,RELC): (  7.0456376     0.93305139E-04 ) :: navier-stokes
  FlowSolve: iter:    8 Assembly: (s)    0.31    2.60
  FlowSolve: iter:    8 Solve:    (s)    0.06    0.45
  FlowSolve:  Result Norm     :    7.0456375744140169
  FlowSolve:  Relative Change :    9.3305138762544466E-005
  FlowSolve: 
  FlowSolve: 
  FlowSolve: -------------------------------------
  FlowSolve:  NAVIER-STOKES ITERATION           9
  FlowSolve: -------------------------------------
  FlowSolve: 
  FlowSolve: Starting Assembly...
  FlowSolve: Assembly:
  FlowSolve: Assembly done
  FlowSolve: Dirichlet conditions done
  ComputeChange: NS (ITER=9) (NRM,RELC): (  7.0451985     0.62313923E-04 ) :: navier-stokes
  FlowSolve: iter:    9 Assembly: (s)    0.31    2.92
  FlowSolve: iter:    9 Solve:    (s)    0.06    0.51
  FlowSolve:  Result Norm     :    7.0451985467781615
  FlowSolve:  Relative Change :    6.2313922620273215E-005
  FlowSolve: 
  FlowSolve: 
  FlowSolve: -------------------------------------
  FlowSolve:  NAVIER-STOKES ITERATION          10
  FlowSolve: -------------------------------------
  FlowSolve: 
  FlowSolve: Starting Assembly...
  FlowSolve: Assembly:
  FlowSolve: Assembly done
  FlowSolve: Dirichlet conditions done
  ComputeChange: NS (ITER=10) (NRM,RELC): (  7.0449052     0.41641666E-04 ) :: navier-stokes
  FlowSolve: iter:   10 Assembly: (s)    0.33    3.24
  FlowSolve: iter:   10 Solve:    (s)    0.06    0.57
  FlowSolve:  Result Norm     :    7.0449051790821748
  FlowSolve:  Relative Change :    4.1641665908848106E-005
  FlowSolve: 
  FlowSolve: 
  FlowSolve: -------------------------------------
  FlowSolve:  NAVIER-STOKES ITERATION          11
  FlowSolve: -------------------------------------
  FlowSolve: 
  FlowSolve: Starting Assembly...
  FlowSolve: Assembly:
  FlowSolve: Assembly done
  FlowSolve: Dirichlet conditions done
  ComputeChange: NS (ITER=11) (NRM,RELC): (  7.0447092     0.27817515E-04 ) :: navier-stokes
  FlowSolve: iter:   11 Assembly: (s)    0.32    3.57
  FlowSolve: iter:   11 Solve:    (s)    0.05    0.62
  FlowSolve:  Result Norm     :    7.0447092100531812
  FlowSolve:  Relative Change :    2.7817514884548234E-005
  FlowSolve: 
  FlowSolve: 
  FlowSolve: -------------------------------------
  FlowSolve:  NAVIER-STOKES ITERATION          12
  FlowSolve: -------------------------------------
  FlowSolve: 
  FlowSolve: Starting Assembly...
  FlowSolve: Assembly:
  FlowSolve: Assembly done
  FlowSolve: Dirichlet conditions done
  ComputeChange: NS (ITER=12) (NRM,RELC): (  7.0445784     0.18568469E-04 ) :: navier-stokes
  FlowSolve: iter:   12 Assembly: (s)    0.33    3.90
  FlowSolve: iter:   12 Solve:    (s)    0.05    0.67
  FlowSolve:  Result Norm     :    7.0445784018006385
  FlowSolve:  Relative Change :    1.8568469342989646E-005
  FlowSolve: 
  FlowSolve: 
  FlowSolve: -------------------------------------
  FlowSolve:  NAVIER-STOKES ITERATION          13
  FlowSolve: -------------------------------------
  FlowSolve: 
  FlowSolve: Starting Assembly...
  FlowSolve: Assembly:
  FlowSolve: Assembly done
  FlowSolve: Dirichlet conditions done
  ComputeChange: NS (ITER=13) (NRM,RELC): (  7.0444912     0.12384643E-04 ) :: navier-stokes
  FlowSolve: iter:   13 Assembly: (s)    0.33    4.22
  FlowSolve: iter:   13 Solve:    (s)    0.05    0.72
  FlowSolve:  Result Norm     :    7.0444911577546474
  FlowSolve:  Relative Change :    1.2384642665345102E-005
  FlowSolve: 
  FlowSolve: 
  FlowSolve: -------------------------------------
  FlowSolve:  NAVIER-STOKES ITERATION          14
  FlowSolve: -------------------------------------
  FlowSolve: 
  FlowSolve: Starting Assembly...
  FlowSolve: Assembly:
  FlowSolve: Assembly done
  FlowSolve: Dirichlet conditions done
  ComputeChange: NS (ITER=14) (NRM,RELC): (  7.0444330     0.82546795E-05 ) :: navier-stokes
  FlowSolve: iter:   14 Assembly: (s)    0.34    4.57
  FlowSolve: iter:   14 Solve:    (s)    0.05    0.76
  FlowSolve:  Result Norm     :    7.0444330079778537
  FlowSolve:  Relative Change :    8.2546795070671733E-006
  FlowSolve: 
  FlowSolve: 
  FlowSolve: -------------------------------------
  FlowSolve:  NAVIER-STOKES ITERATION          15
  FlowSolve: -------------------------------------
  FlowSolve: 
  FlowSolve: Starting Assembly...
  FlowSolve: Assembly:
  FlowSolve: Assembly done
  FlowSolve: Dirichlet conditions done
  ComputeChange: NS (ITER=15) (NRM,RELC): (  7.0443943     0.54993260E-05 ) :: navier-stokes
  FlowSolve: iter:   15 Assembly: (s)    0.33    4.90
  FlowSolve: iter:   15 Solve:    (s)    0.06    0.83
  FlowSolve:  Result Norm     :    7.0443942684508842
  FlowSolve:  Relative Change :    5.4993259849576591E-006
  FlowSolve: 
  FlowSolve: 
  FlowSolve: -------------------------------------
  FlowSolve:  NAVIER-STOKES ITERATION          16
  FlowSolve: -------------------------------------
  FlowSolve: 
  FlowSolve: Starting Assembly...
  FlowSolve: Assembly:
  FlowSolve: Assembly done
  FlowSolve: Dirichlet conditions done
  ComputeChange: NS (ITER=16) (NRM,RELC): (  7.0443685     0.36626062E-05 ) :: navier-stokes
  FlowSolve: iter:   16 Assembly: (s)    0.33    5.23
  FlowSolve: iter:   16 Solve:    (s)    0.05    0.87
  FlowSolve:  Result Norm     :    7.0443684676559917
  FlowSolve:  Relative Change :    3.6626062026636725E-006
  FlowSolve: 
  FlowSolve: 
  FlowSolve: -------------------------------------
  FlowSolve:  NAVIER-STOKES ITERATION          17
  FlowSolve: -------------------------------------
  FlowSolve: 
  FlowSolve: Starting Assembly...
  FlowSolve: Assembly:
  FlowSolve: Assembly done
  FlowSolve: Dirichlet conditions done
  ComputeChange: NS (ITER=17) (NRM,RELC): (  7.0443513     0.24389716E-05 ) :: navier-stokes
  FlowSolve: iter:   17 Assembly: (s)    0.33    5.56
  FlowSolve: iter:   17 Solve:    (s)    0.05    0.92
  FlowSolve:  Result Norm     :    7.0443512866625202
  FlowSolve:  Relative Change :    2.4389715703175822E-006
  FlowSolve: 
  FlowSolve: 
  FlowSolve: -------------------------------------
  FlowSolve:  NAVIER-STOKES ITERATION          18
  FlowSolve: -------------------------------------
  FlowSolve: 
  FlowSolve: Starting Assembly...
  FlowSolve: Assembly:
  FlowSolve: Assembly done
  FlowSolve: Dirichlet conditions done
  ComputeChange: NS (ITER=18) (NRM,RELC): (  7.0443398     0.16240771E-05 ) :: navier-stokes
  FlowSolve: iter:   18 Assembly: (s)    0.33    5.88
  FlowSolve: iter:   18 Solve:    (s)    0.06    0.98
  FlowSolve:  Result Norm     :    7.0443398461019884
  FlowSolve:  Relative Change :    1.6240771302326036E-006
  FlowSolve: 
  FlowSolve: 
  FlowSolve: -------------------------------------
  FlowSolve:  NAVIER-STOKES ITERATION          19
  FlowSolve: -------------------------------------
  FlowSolve: 
  FlowSolve: Starting Assembly...
  FlowSolve: Assembly:
  FlowSolve: Assembly done
  FlowSolve: Dirichlet conditions done
  ComputeChange: NS (ITER=19) (NRM,RELC): (  7.0443322     0.10814906E-05 ) :: navier-stokes
  FlowSolve: iter:   19 Assembly: (s)    0.33    6.21
  FlowSolve: iter:   19 Solve:    (s)    0.05    1.03
  FlowSolve:  Result Norm     :    7.0443322277185283
  FlowSolve:  Relative Change :    1.0814906359046379E-006
  FlowSolve: 
  FlowSolve: 
  FlowSolve: -------------------------------------
  FlowSolve:  NAVIER-STOKES ITERATION          20
  FlowSolve: -------------------------------------
  FlowSolve: 
  FlowSolve: Starting Assembly...
  FlowSolve: Assembly:
  FlowSolve: Assembly done
  FlowSolve: Dirichlet conditions done
  ComputeChange: NS (ITER=20) (NRM,RELC): (  7.0443272     0.72023826E-06 ) :: navier-stokes
  FlowSolve: iter:   20 Assembly: (s)    0.34    6.56
  FlowSolve: iter:   20 Solve:    (s)    0.05    1.08
  FlowSolve:  Result Norm     :    7.0443271541227590
  FlowSolve:  Relative Change :    7.2023826139606169E-007
  ComputeDevStress: 
  ComputeDevStress: 
  ComputeDevStress: 
  ComputeDevStress: 
  ComputeDevStress: Starting assembly...
  ComputeDevStress:  Component :   1
  ComputeDevStress: Assembly done
  ComputeDevStress: Set boundaries done
  ComputeChange: NS (ITER=1) (NRM,RELC): (  1.2144601      2.0000000     ) :: sij
  ComputeDevStress:  Component :   2
  ComputeDevStress: Assembly done
  ComputeDevStress: Set boundaries done
  ComputeChange: NS (ITER=2) (NRM,RELC): (  1.2124247     0.16773815E-02 ) :: sij
  ComputeDevStress:  Component :   3
  ComputeDevStress: Assembly done
  ComputeDevStress: Set boundaries done
  ComputeChange: NS (ITER=3) (NRM,RELC): (  1.2117236     0.57838410E-03 ) :: sij
  ComputeDevStress:  Component :   4
  ComputeDevStress: Assembly done
  ComputeDevStress: Set boundaries done
  ComputeChange: NS (ITER=4) (NRM,RELC): ( 0.44539613E-01  1.8581838     ) :: sij
  ComputeDevStress:  Result Norm   :    1.0506122799431621        1.2117236081674525
  ComputeDevStress:  Relative Change :   0.14242918487125467
  FreeSurfaceSolver (zs bottom): Solving for free surface
  FreeSurfaceSolver (zs bottom): Memory allocations done
  FreeSurfaceSolver (zs bottom): Non-linear Iteration 1 out of max. 100
  Dirichlet:  Eliminated 21 unknowns out of 101
         1 0.2749E-07
         2 0.1510E-14
         2 0.1510E-14

raback · 2025-10-17T09:27:49Z

Thanx for the attention to the issue! Nothing to merge anymore.

tzwinger and others added 25 commits June 2, 2025 11:09

added YAC initialization and finalization

c9a94bf

added YAC to setup process

8115eba

added new search functionality to cmake-part

1df8b90

added coupling code to Elmer (contrib altered from M. Hanke, DKRZ)

17d0c36

first compiling version on levante

2526b5a

fixed path issue for Elmer grid file

cd891c5

corrected wrong partition directory for loading Elmer grid

7e50bcf

ranming variables

cd261ff

changed order of module declaration in elmer_coupling.F90

c09780f

Use MPI_Handshake for Comm split

12913a7

Check YAC version and add missing function.

f96b267

Add MPI handshake implementation from https://gitlab.dkrz.de/dkrz-sw/…

8195b2f

…mpi-handshake as fallback if working without YAC.

Minor consistency fixes and checks.

17620f8

debugging yac-xios integration

2ded22f

Add new files to CMakeLists.txt.

cca5384

Add include guard.

4825f03

Revert "Add include guard."

13e5c48

This reverts commit 994152c.

Move up.

3eed4f7

Try with guards for mpi_handshake_stubs.

856de11

Put guard around mpi_handshake.c.

1dab640

Try with a different guard and provide stub implementation.

97a9548

Rename file.

ec50b7b

Merge branch 'devel' into yac-coupling

9dc173c

Improve CMake module for XIOS.

4677b0e

Make imports more specific.

a2298a4

BenjaminRodenberg mentioned this pull request Aug 28, 2025

WIP: YAC-coupling integration in Elmer (colouring and mpi_handshake) #682

Draft

BenjaminRodenberg added 2 commits August 28, 2025 16:06

Avoid destructing YAC if it has not been initialized. (WIP and hacky …

958f5c6

…solution)

Try completely without the dirty fix.

a2b3c1b

BenjaminRodenberg force-pushed the yac-coupling-debug-windows-ci branch from ffea4c9 to a2b3c1b Compare August 28, 2025 14:11

BenjaminRodenberg changed the base branch from yac-coupling to devel September 1, 2025 15:29

BenjaminRodenberg added 3 commits September 1, 2025 17:47

Revert most changes.

a8486d6

Remove grid reader elmer_grid.{c,h} and corresponding calls.

4b2d69b

Revert CMakeLists.txt. (this should effectively lead to same src as o…

996331f

…n devel)

BenjaminRodenberg force-pushed the yac-coupling-debug-windows-ci branch from 218b773 to 996331f Compare September 1, 2025 17:32

Reduce diff to just a file containing a comment.

f1d4c36

tzwinger requested a review from raback September 3, 2025 05:44

tzwinger assigned raback Sep 3, 2025

raback closed this Oct 17, 2025

BenjaminRodenberg deleted the yac-coupling-debug-windows-ci branch October 17, 2025 13:36

Provide feedback

Saved searches

Use saved searches to filter your results more quickly

Uh oh!

WIP: Debug windows CI for YAC coupling #691

WIP: Debug windows CI for YAC coupling #691

Uh oh!

BenjaminRodenberg commented Aug 28, 2025

Uh oh!

BenjaminRodenberg commented Sep 2, 2025

Uh oh!

tzwinger commented Sep 3, 2025

Uh oh!

mmuetzel commented Oct 7, 2025

Uh oh!

tzwinger commented Oct 13, 2025

Uh oh!

mmuetzel commented Oct 13, 2025

Uh oh!

tzwinger commented Oct 13, 2025 •

edited

Loading

Uh oh!

mmuetzel commented Oct 14, 2025

Uh oh!

raback commented Oct 17, 2025

Uh oh!

Reviewers

Assignees

Labels

Projects

Milestone

Development

Uh oh!

5 participants

WIP: Debug windows CI for YAC coupling #691

WIP: Debug windows CI for YAC coupling #691

Uh oh!

Conversation

BenjaminRodenberg commented Aug 28, 2025

Uh oh!

BenjaminRodenberg commented Sep 2, 2025

Uh oh!

tzwinger commented Sep 3, 2025

Uh oh!

mmuetzel commented Oct 7, 2025

Uh oh!

tzwinger commented Oct 13, 2025

Uh oh!

mmuetzel commented Oct 13, 2025

Uh oh!

tzwinger commented Oct 13, 2025 • edited Loading Uh oh! There was an error while loading. Please reload this page.

Uh oh!

Uh oh!

mmuetzel commented Oct 14, 2025

Uh oh!

raback commented Oct 17, 2025

Uh oh!

Reviewers

Assignees

Labels

Projects

Milestone

Development

Uh oh!

5 participants

tzwinger commented Oct 13, 2025 •

edited

Loading