Coupled.sif

check keywords warn
echo on

$Name = "Coupled"
$gravity = 9.7696e15
$rhoi = 9.1501e-19
$rhow = 1000.0*9.1501e-19/910.0 !1.0055e-18

Header
  Mesh DB "." "cube"
End

Constants
  Water Compressibility = Real 5.04e-4  !MPa-1
End

!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
Simulation
  Coordinate System  = Cartesian 3D 
  Coordinate Mapping(3) = Integer 1 2 3
  Simulation Type = Transient

  Timestepping Method = "bdf"
  BDF Order = 1
  Timestep Intervals = 5
  Timestep Sizes = 0.002

  Steady State Min Iterations = 1
  Steady State Max Iterations = 100

  Output File = "HydroTestCase_$Name".result" !"
  max output level = 4

End

!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
! The ice 
Body 1
  Equation = 1
  Body Force = 1
  Material = 1
  Initial Condition = 1
End

! The Aquifer and the bed Free surface 
Body 2
  Equation = 2
  Body Force = 1
  Material = 1
  Initial Condition = 2
End

!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!

Initial Condition 1
  Pressure = Real 0.0 
End

Initial Condition 2    
  IDSHead = Real 0.0 
  EPLHead = Real 0.0 
End

!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
Body Force 1

  Flow BodyForce 1 = Real 0.0
  Flow BodyForce 2 = Real 0.0
  Flow BodyForce 3 = Real -$gravity
  IDSHead Source Flux = Real 2.0

  EPLToIDS Transfer = Variable Coordinate 1
    Real Procedure "ElmerIceUSF" "EPLToIDS"     

  EPLHead Passive = Equals Open EPL
  Active Element Min Nodes = Integer 1
End
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
Material 1
  Density = Real $rhoi 

! General Hydrology Parameters
  Water Density = Real $rhow
  Leakage Factor = Real 20.0

! IDS Solver
  IDS Transmitivity = Real 5.0e4
  IDS Porosity = Real 0.4
  IDS Thickness = Real 20.0 
  IDS Compressibility = Real 1.0e-2 

  IDSHead Upper Limit = Real 15.0

! EPL Solver
  EPL Transmitivity = Real 2.5e5
  EPL Porosity = Real 0.4
  EPL Thickness = Real 1.0 
  EPL Compressibility = Real 1.0e-2

  EPLHead Upper Limit = Real 15.0

End

!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
Solver 1

  Equation = "Flowdepth"

  Procedure = "ElmerIceSolvers" "FlowDepthSolver"
  Variable = String "Depth"
  Variable DOFs = 1
  Freesurf Name = String "Surface"

  Linear System Solver = "Direct"
  Linear System Direct Method = umfpack
  Gradient = Real -1.0E00

  Calc Free Surface = Logical FALSE
End

!------------------------------------------------------
Solver 2
 Equation = "IDS Equation"

  Procedure = "ElmerIceSolvers" "IDSSolver"
  Variable = IDSHead
  Variable DOFs = 1

  Steady State Convergence Tolerance = Real 1.0E-5

  Linear System Solver = Direct         
  Linear System Direct Method = umfpack
  Linear System Convergence Tolerance = Real 1.0E-7
  Linear System Residual Output = integer 1

  Nonlinear System Max Iterations = Integer 100
  Nonlinear System Convergence Tolerance = Real 1.0E-6
  Nonlinear System Relaxation Factor = Real 1.0

  Apply Dirichlet = Logical TRUE

  Exported Variable 1 = String "IDSHead Residual"
  Exported Variable 2 = String "IDSHead Homologous"
  Exported Variable 3 = String "IDSHead Pressure"
  
End

!------------------------------------------------------
Solver 3
 Equation = "EPL Equation"

  Procedure = "ElmerIceSolvers" "EPLSolver"
  Variable = EPLHead
  Variable DOFs = 1

  Steady State Convergence Tolerance = Real 1.0E-5

  Linear System Solver = Direct         
  Linear System Direct Method = umfpack
  Linear System Convergence Tolerance = Real 1.0E-7
  Linear System Residual Output = integer 1

  Nonlinear System Max Iterations = Integer 100
  Nonlinear System Convergence Tolerance = Real 1.0E-6
  Nonlinear System Relaxation Factor = Real 1.0

  IDS Residual Name = String "IDSHead Residual"
  IDS Load Name = String "IDSHead"

  Apply Dirichlet = Logical TRUE

  Exported Variable 1 = String "EPLHead Homologous"
  Exported Variable 2 = String "Open EPL"
  
End

!------------------------------------------------
Solver 4

  Equation = "Save Line"
   Exec Solver = After TimeStep

  Procedure = File "SaveData" "SaveLine"
  Filename =  "hydro_boundary_$Name".dat" !"
  File Append = Logical False

End

!------------------------------------------------------
Solver 5

  Equation = "Save Scalars"
  Exec Solver = After TimeStep

  Procedure = File "SaveData" "SaveScalars"
  Filename =  results.dat !"hydro_scalars_$Name".dat" !"
  File Append = Logical True     ! For transient simulation

  Variable 1 = String "Time"
  !Operator 2 = String "cpu time"
  !Operator 3 = String "cpu memory"
  
  Variable 2 = String "IDSHead"
  Operator 2 = String "sum"
  Operator 3 = String "max"
  Operator 4 = String "min"

  Variable 5 = String "EPLHead"
  Operator 5 = String "sum"
  Operator 6 = String "Max"
  Operator 7 = String "Min"

End

!------------------------------------------------------
Solver 6
  
  Equation = "result output"
  Exec Solver = After TimeStep

  Procedure = "ResultOutputSolve" "ResultOutputSolver"
  Output File Name = File "HydroTestCase_$Name".vtu" !"
  Output Format = String vtu
  Save Geometry Ids = Logical True
End

!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!

Equation 1
  Active Solvers(4) = 1 4 5 6
End

Equation 2
  Active Solvers(2) = 2 3
End

!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
!Bedrock
!========
Boundary Condition 1
  Target Boundaries  = 5

!!! this BC is equal to body no. 2 !!!
  Body Id = 2

  Save Line = Logical True
End

! Upper Surface
!==============
Boundary Condition 2
  Target Boundaries = 6
!!THE FLOWDEPTH SOLVER
  Depth = Real 0.0

  Save Line = Logical True
End

! Front and sides
!================
Boundary Condition 3
  Target Boundaries (4) = 1 2 3 4

 Save Line = Logical True
End

!Snout
!=====
Boundary Condition 4
  Target Boundaries = 8

!! FOR THE HYDRO                                                                               
   IDSHead = Real 0.0
   EPLHead = Real 0.0
  Save Line = Logical True
End

!Lower frame
!===========

Boundary Condition 5
  Target Boundaries = 7
!! FOR THE HYDRO
   IDSHead Flux BC = Logical True
   IDSHead Water Flux = Real 0.0
   EPLHead Flux BC = Logical True
   EPLHead Water Flux = Real 0.0
  Save Line = Logical True
End

Solver 2 :: Reference Norm = Real 6.4419305
Solver 2 :: Reference Norm Tolerance = Real 1E-06