Code we could consider deleting

• rifts
• tetras
• flowband models
• ice model coupling

Things to improve/address

• too many memory allocations
• get rid of Esa solver, adapt newly developed solid Earth capabilities instead
• make model a dynamic class (do not set all fields by default). This would help integrate new solutions that may not be ice oriented into ISSM (e.g. GIA, Kriging, etc)

New capabilities

• Julia interface?