mesh {{{
	  x
	  y
	  z
	  elements -> rename index?
	  dim
	  numlayers
	  numberofelements
	  numberofnodes -> rename all nodes vertices? (nodes in the C++ core are not the same as the ones in model)
	  edges

	  uppernodes
	  upperelements
	  lowerelements
	  lowernodes

	  numberofelements2d
	  numberofnodes2d
	  elements2d
	  x2d
	  y2d
	  z2d

	  extrusionexponent -> remove?
	  elementonbed
	  elementonsurface
	  nodeonbed
	  nodeonsurface

	  nodeconnectivity
	  elementconnectivity
	  connectivity

	  segments
	  segmentmarkers
	  nodeonboundary

	  bamg

	  extractednodes
	  extractedelements
}}}
projection{{{
	lat
	  long
	  hemisphere
}}}
geometry {{{
	  surface
	  thickness
	  bed
	  bathymetry

	  bed_slopex
	  bed_slopey
	  surface_slopex
	  surface_slopey

	  thickness_coeff -> delete or rename (See with John)?
}}}
geography{{{ -> rename mask (as in SeaRISE)?
	elementoniceshelf -> rename all as grounded and floating
	  elementonicesheet
	  elementonwater
	  nodeoniceshelf
	  nodeonicesheet
	  nodeonwater
}}}
rifts{{{
	rifts
	numrifts
	riftinfo -> should be removed (constructed by marshall)
	riftproperties
}}}
materials{{{
	rho_ice
	rho_water
	rho_firn -> remove
	heatcapacity
	latentheat
	thermalconductivity
	meltingpoint
	beta
	mixed_layer_capacity
	thermal_exchange_velocity
	rheology_B
	rheology_n
	rheology_law
}}}
friction{{{
	drag_type
	  drag_coefficient
	  drag_p
	  drag_q
}}}
surfaceforcings{{{
	surface_accumulation_rate -> remove all surface_
	surface_ablation_rate
	surface_mass_balance
}}}
basalforcings{{{
	basal_melting_rate -> remove basal_*
	basal_melting_rate_correction
	basal_melting_rate_correction_apply -> remove (correction 0 has the same effect)
	geothermalflux
}}}

flowapproximation{{{
	elements_type -> change all type to approx?
	vertices_type

	ishutter
	ismacayealpattyn
	isstokes

	elementonhutter
	elementonmacayeal
	elementonpattyn
	elementonstokes

	nodeonhutter
	nodeonmacayeal
	nodeonpattyn
	nodeonstokes

	borderstokes

	elements_type2d -> remove
	vertices_type2d -> remove
}}}
diagnostic{{{
	isdiagnostic
	vx
	vy
	vz
	vel
	spcvx
	spcvy
	spcvz
	stokesreconditioning
	shelf_dampening
	viscosity_overshoot
	pressureload -> rename icefront?
	diagnostic_ref -> rename referential?
	penalties
	pressure

	eps_res -> rename following PETSc: restol
	eps_rel -> rename following PETSc: reltol
	eps_abs-> rename following PETSc:  abstol
	max_nonlinear_iterations -> rename maxiter

	penalty_lock -> rename riftpenaltylock or equiv?
	min_mechanical_constraints -> rename (not clear)
}}}
thermal{{{
	isthermal
	temperature
	spctemperature
	penalty_melting
	min_thermal_constraints -> rename (see min_mechanical_constraints)
	waterfraction
	artificial_diffusivity -> rename stabilization
	max_nonlinear_iterations -> rename maxiter
	stabilize_constraints -> rename (not clear)
}}}
prognostic{{{
	isprognostic
	spcthickness
	hydrostatic_adjustment
	artificial_diffusivity -> rename stabilization
	prognostic_DG -> should be merged with stabilization
	penalties
}}}
balancethickness{{{
	spcthickness
	dhdt
	artificial_diffusivity -> rename stabilization
}}}
steadystate{{{
	eps_rel -> rename following PETSc: reltol
	max_steadystate_iterations -> rename maxiter
}}}
goundingline {{{
	groundingline_migration
	groundingline_melting_rate
}}}
hydrology{{{
	watercolumn
	spwatercolumn
	hydro_n
	hydro_CR
	hydro_p
	hydro_q
	hydro_kn
}}}
control{{{
	control_analysis
	num_control_type -> num_controls
	control_type -> rename control_paramters?
	num_cm_responses
	cm_responses -> rename objective_functions
	cm_gradient -> rename gradient_only
	weights

	cm_min -> min_parameters
	cm_max -> max_parameters

	nsteps
	maxiter
	tolx
	optscal
	cm_jump -> jump
	eps_cm  -> remove (never used)
	epsvel  -> hardcoded?
	meanvel -> hardcoded?

	vx_obs
	vy_obs
	vel_obs
	thickness_obs
}}}
qmu{{{
	variables
	responses
	qmu_method -> remove all qmu_
	qmu_params
	dakotaresults
	dakotain
	dakotaout
	dakotadat
	qmu_analysis
	part
	npart
	numberofvariables
	numberofresponses
	variabledescriptors
	responsedescriptors
	qmu_mass_flux_profile_directory
	qmu_mass_flux_profiles
	qmu_mass_flux_segments
	qmu_relax
	qmu_save_femmodel

	adjacency
	vwgt
}}}
flaim{{{
	fm_tracks -> remove all fm_
	  fm_flightreqs
	  fm_criterion
	  fm_gridsatequator
	  fm_usevalueordering
	  fm_split_antimeridian
	  fm_solution
	  fm_quality
}}}

verbose {{{
	-> already exists
}}}
cluster{{{
	waitonlock -> cluster dependent
	-> already exists
}}}
timestepping{{{
	dt
	ndt
	time_adapt -> rename cfl_apply
	adaptative_cfl -> remove
	cfl_coefficient
}}}
results{{{
	requested_outputs
	results_on_vertices
	output_frequency
	outlog
	errlog
	  + all model results ??? (DiagnosticSolution and others?)
}}}
solver{{{ -> or petsc?
petscoptions -> could become the class itself
}}}
debug{{{
	gprof
	mem_debug -> rename valgrind
	memoryleaks -> never used
}}}
private {{{
	counter
	runtimename
	runmefile
	inputfilename -> is it really needed? it is just runtimename.bon
	outputfilename -> same
	solution_type -> rename solution
}}}
constants{{{
	%not material dependent
g
yts
referencetemperature
}}}
computation{{{
io_gather
penalty_offset
sparsity
lowmem
}}}

radarimage{{{
	sarpwr
	  sarxm
	  sarym
}}}
miscelaneous {{{
	notes
	name
	dummy
}}}

mesh_domainoutline -> remove
mesh_riftoutline -> remove
mesh_resolution -> remove
domainoutline -> remove
riftoutline -> remove
iceshelfoutline ->remove
icesheetoutline -> remove
wateroutline -> remove
parameterfile -> remove

vx_obs_raw -> remove (should be in dummy)
vy_obs_raw -> remove (should be in dummy)
vel_obs_raw -> remove (should be in dummy)
vx_bal -> remove (should be in dummy)
vy_bal -> remove (should be in dummy)
vel_bal -> remove (should be in dummy)

forcings -> delete
minh -> remove or really use in the code!
firn_layer -> remove (never used)
deltaH -> remove (unused)
DeltaH -> remove (unused)
deltaT -> remove (unused)
DeltaT -> remove (unused)
timeacc -> remove (unused)
timedec -> remove (unused)
observed_temperature -> remove
repository -> remove (out of core never used)
loadresults -> remove (should be an option of solve: md=solve(md,DiagnosticSolutionEnum,'loadonly',1))

element_debug -> remove
element_debugid -> remove

viscousheating -> remove (that's a model result)
pressure_elem -> remove (that's a model result)
stress-> remove (that's a model result)
stress_surface-> remove (that's a model result)
stress_bed-> remove (that's a model result)
deviatoricstress-> remove (that's a model result)
strainrate-> remove (that's a model result)

alloc_cleanup -> remove
mask -> delete?
