BISICLES driver executable
BISICLES' driver executable is used to run standalone time-dependent ice sheet simulations. At the start of a simulation, it carries out the following operations- Read a run-time configuration file, specified as the first command line argument
- Compute initial thickness, topography, internal energy, stiffness and basal friction coefficient data over the whole domain on the base level (l = 0) mesh. This data might be loaded from a file through the LevelData interface, computed through the Python interface or through a number of special purpose modules.
- Set L = 0 and solve the stress-balance equations to find the velocity field on the l = 0 mesh
- Generate a new mesh with L + 1 levels on which some regions are more finely resolved according to the refinement criteria
- Compute the thickness, topography, internal energy, damage and basal friction coefficient data on the L + 1 mesh
- Solve the stress-balance equations to find the velocity field on the L + 1 mesh
- Write a plot file, including the thickness, topography, surface elevation, velocity etc
- if L has not reached its maximum value, set L = L + 1 and repeat 3-7
- Write a check-point file, if required
- If time t > simulation time, exit, else continue
- Compute surface fluxes (accumulation and melting) given the geometry and velocity field
- Compute horizontal thickness and internal energy fluxes given the geometry, velocity field and surface fluxes
- Evolve the thickness and internal energy fields for a time determined by the CFL condition
- Solve the stress-balance equations to find a new velocity field
- Write a plot file, if required
- Update the mesh according to the refinement criteria
- If needed, solve the stress-balance equations to find a new velocity field
- repeat from step 1
It is also possible to restart a simulation from a checkpoint file, for example to create a number of simulation branches, having carried out some sort of relaxation