GLIDE is a 3D finite-difference ice sheet model. GLIDE is the core system of the GLIMMER ice sheet model. It was originally developed by Tony Payne at the University of Bristol, but has since been modified to form the land-ice component of the GENIE project. This system is a collaborative effort to create a standard ice sheet model which can be used by other projects.
Some features of GLIDE:
These pages document my involvement with GLIMMER and GLIDE and serve as a place where I can put some experiments and musings. You can also get a pdf of the current documentation or browse the HTML Documentation online (complain to me when it is too out of date).
A while back I gave a tutorial of how to use version control systems, CVS and subversion, and GNU automake and autoconf in the context of GLIMMER. The tutorial also included some bits on the GLIMMER configuration files and the code generator. You can get the pdf slides or if you want something more compact the handouts of the tutorial. The presentation was produced with the beamer class for LaTeX. Finally, no more powerpoint!
I've also run some experiments for test D and created some new images. You can find the results here. I consider these tests as successful and will add a module to the GLIMMER cvs repository once I get to go from Ed.
We can now generate HTML Documentation for GLIMMER. The LaTeX sources are translated using TeX4ht.
You can find the initial results of the tests proposed by Ed Bueler and others to verify ice sheet models using exact solutions for the isothermal case here.
Trying to figure out why we get very narrow regions of fast flow when basal sliding is switched on. This effect can be clearly seen in my simulations of the past Fennoscandian ice sheet. I have repeated the EISMINT-2 Exp H with varying sliding parameters to investigate this further.
I've fixed some bugs with the non-linear solver and played with it by varying the timestep size of the EISMINT-2 Experiment H. The result suggest that it doesn't really matter if the non-linear or linear solver are used, at least for this experiment. Obviously, there might also still be a bug lurking.
The additional heat flux due to sliding was indeed missing. I've fixed that now. Both Experiment G and Experiment H look quite good now as well.
I've also written a new PyCF program which will compare EISMINT-2 output with the results reported in the EISMINT-2 paper (e.g. for the EISMINT-2 Experiment H). Again, the results look good.
The EISMINT-2 tests are now also added and look very promising. However, there seems to be a problem when the ice is sliding: basal temperatures are too low and solution is patchy.
I finally managed to reproduce the EISMINT-1 results.