BISMG:SarahS/glimmer-cism-lanl

Basal water in glimmer-cism

 * 0 Local water balance
 * 1 Local water balance + constant flux
 * 2 None

Running the trilinos test

 * There is a test located in the example/TrilinosSolversUnitTest directory
 * Gethin made a Makefile for ifort for this test called Makefile.ifort. This is now in the repository.
 * To run the test case type make -f Makefile.ifort TrilinosSolversUnitTest.exe

Installing glimmer-cism parallel version on dartagnan

 * svn co https://svn.berlios.de/svnroot/repos/glimmer-cism/glimmer-cism-lanl/branches/parallel
 * ./bootstrap
 * ./configure --prefix=/home/paleo/ggsrs/parallel --enable-tests --with-netcdf=/opt/local/CentOS-64/netcdf/4.0/intel_fc_10.1/  --with-lapack=-llapack --with-blas=/usr/lib64/libblas.a --with-trilinos=/opt/local/CentOS-64/trilinos/trilinos-10.0.4-Source/SERIAL/ FC=ifort F77=ifort
 * make

Installing glimmer-cism on dartagnan
Currently using this version of glimmer-cism: https://svn.berlios.de/svnroot/repos/glimmer-cism/glimmer-cism-lanl/trunk


 * Check out code
 * ./bootstrap
 * ./configure --prefix=/home/paleo/ggsrs/lanl_update/trunk/ --with-netcdf=/opt/local/CentOS-64/netcdf/4.0/intel_fc_10.1/ FC=ifort F77=ifort FCFLAGS='-g -traceback -DNO_RESCALE'. The -DNO_RESCALE is necessary if you want to run the Pattyn HO model.
 * make

If you want to setup the model to use different types of parallel solvers, then you need to add a compile flag.
 * For PARDISO use --with-pardiso
 * For UMFPACK use --with-umfpack
 * The default (i.e no flag) uses SLAP library.

How to enable higher order physics
Add [ho_options] to the configuration file

diagnostic_scheme   No time evolution Higher-order velocity computation scheme 0 SIA, no higher order 1 Pattyn unstaggered 2 Pattyn staggered 3 Payne/Price scheme

prognostic_scheme   Time evolution 0 Do not not use higher-order velocities prognostically; compute and output them but use SIA to evolve the ice 1 Pattyn scheme (compute higher-order diffusivities only) 2 Pollard scheme (Not implemented yet) 3 Bueler scheme (Not implemented yet)

basal_stress_input  Flag that indicates how to compute beta, the higher-order basal stress 0 Set to NaN everywhere (no sliding, ice glued to the bed) 1 Set to 1/soft (default) 2 Set to 1/btrc (re-computed at each time step) 3 Set to beta field of input netcdf file basal_stress_type

guess_specified

which_ho_source

include_thin_ice

Price-Payne higher-order (glam) options:

which_ho_babc

which_ho_efvs

which_ho_resid

which_ho_sparse  Flag that indicates method for solving the sparse linear system that arises from the higher-order solver 0 Biconjugate Gradient, Incomplete LU Preconditioner 1 GMRES, Incomplete LU Preconditioner 2 UMFPACK Unsymmetric multifrontal direct solver 3 PARDIOS - Must have pre-compiled libraries installed

which_ho_sparse_fallback

Modifications to Makefile
Remove lines that compile Edinburgh climate driver (i.e eis* )

Remove lines that envoke python:

filename_io.F90 is made on the fly using phthon (generate_ncvars.py) Remove lines in Makefile that remove filename_io.F90 when 'make clean' is used.

Python is still used in the Makefile to check for dependencies.

glimmer-cism-lanl structure
glide.F90 glide_thck.F90 glide_temp.F90 glide_velo.F90 other_glide_files