ADD: md.groundingline in class, no implementation for now

Author: enigne

src/core/elements.jl

@@ -579,20 +579,34 @@ function InputUpdateFromVector(element::Tria, vector::Vector{Float64}, enum::Iss
 	return nothing
 end #}}}
 function IsAllFloating(element::Tria) #{{{
-	input = GetInput(element, MaskOceanLevelsetEnum)
 
-	# by default use none migration
-	if GetInputMin(input) > 0. 
-		return false;
+	input=GetInput(element, MaskIceLevelsetEnum)
+	migration_style = FindParam(String, element, GroundinglineMigrationEnum)
+	# get the Enum
+	migration_style_enum = StringToEnum(migration_style)
+
+	if (migration_style_enum==SubelementMigrationEnum)
+		if GetInputMax(input) <= 0.
+			return true;
+		else
+			return false
+		end
+	elseif (migration_style_enum==ContactEnum)
+		if GetInputMin(input) < 0.
+			return true;
+		else
+			return false
+		end
+	elseif (migration_style_enum==NoneEnum || migration_style_enum==AggressiveMigrationEnum || migration_style_enum==SoftMigrationEnum || migration_style_enum==GroundingOnlyEnum)
+		# by default use none migration
+		if GetInputMin(input) > 0. 
+			return false;
+		else
+			return true
+		end
 	else
-		return true
+		error("migration_style not implemented yet")
 	end
-	# TODO: add subelement migration
-	#if GetInputMax(input) <= 0. 
-	#	return true;
-	#else
-	#	return false
-	#end
 end#}}}
 function IsIcefront(element::Tria) #{{{
 

src/core/modules.jl

@@ -136,6 +136,9 @@ function CreateParameters(parameters::Parameters,md::model,solutionstring::Symbo
 
 	#Is moving front 
 	AddParam(parameters,md.transient.ismovingfront,TransientIsmovingfrontEnum)
+	#Is grounding line
+	AddParam(parameters,md.transient.isgroundingline,TransientIsgroundinglineEnum)
+	AddParam(parameters,md.groundingline.migration,GroundinglineMigrationEnum)
 
 	return nothing
 end# }}}

src/usr/classes.jl

@@ -321,6 +321,17 @@ function Base.show(io::IO, this::Frontalforcings)# {{{
 	IssmStructDisp(io, this)
 end# }}}
 # }}}
+#Groundingline{{{
+mutable struct Groundingline
+	migration::String
+end
+function Groundingline() #{{{
+	return Groundingline("None")
+end# }}}
+function Base.show(io::IO, this::Groundingline)# {{{
+	IssmStructDisp(io, this)
+end# }}}
+# }}}
 
 #Model structure
 mutable struct model{Mesh<:AbstractMesh, Friction<:AbstractFriction, Basalforcings<:AbstractBasalforcings, Calving<:AbstractCalving}
@@ -342,20 +353,21 @@ mutable struct model{Mesh<:AbstractMesh, Friction<:AbstractFriction, Basalforcin
 	calving::Calving
 	levelset::Levelset
 	frontalforcings::Frontalforcings
+	groundingline::Groundingline
 end
 function model() #{{{
       return model( Mesh2dTriangle(), Geometry(), Mask(), Materials(),
                                        Initialization(),Stressbalance(), Constants(), Dict(),
                                        BuddFriction(), DefaultBasalforcings(), SMBforcings(), DefaultTimestepping(),
 													Masstransport(), Transient(), Inversion(), DefaultCalving(), 
-													Levelset(), Frontalforcings())
+													Levelset(), Frontalforcings(), Groundingline())
 end#}}}
 function model(md::model; mesh::AbstractMesh=md.mesh, friction::AbstractFriction=md.friction, calving::AbstractCalving=md.calving, basalforcings::AbstractBasalforcings=md.basalforcings) #{{{
 	return model(mesh, md.geometry, md.mask, md.materials, 
 					 md.initialization, md.stressbalance, md.constants, md.results, 
 					 friction, basalforcings, md.smb, md.timestepping, 
 					 md.masstransport, md.transient, md.inversion, md.calving, 
-					 md.levelset, md.frontalforcings)
+					 md.levelset, md.frontalforcings, md.groundingline)
 end#}}}
 function model(matmd::Dict; verbose::Bool=true, friction::AbstractFriction=BuddFriction(), basalforcings::AbstractBasalforcings=DefaultBasalforcings()) #{{{
 
@@ -436,6 +448,7 @@ function Base.show(io::IO, md::model)# {{{
 	@printf "%19s: %-26s -- %s\n" "calving" typeof(md.calving) "parameters for calving"
 	@printf "%19s: %-26s -- %s\n" "levelset" typeof(md.levelset) "parameters for moving boundaries (level-set method)"
 	@printf "%19s: %-26s -- %s\n" "frontalforcings" typeof(md.frontalforcings) "parameters for frontalforcings"
+	@printf "%19s: %-26s -- %s\n" "groundingline" typeof(md.groundingline) "parameters for groundingline"
 	@printf "%19s: %-26s -- %s\n" "results" typeof(md.results) "model results"
 
 end# }}}