++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++

From greve at lowtem.hokudai.ac.jp  Fri Dec  6 04:42:56 2013
From: greve at lowtem.hokudai.ac.jp (Ralf Greve)
Date: Fri, 06 Dec 2013 12:42:56 +0900
Subject: [SICOPOLIS:27] Implementation of the enthalpy method in SICOPOLIS
Message-ID: <52A147C0.50109@lowtem.hokudai.ac.jp>

Dear colleagues,

in the recent weeks, Heinz Blatter and I have been working quite 
intensively on the enthalpy method [Aschwanden et al., J. Glaciol. 
58(209), 2012 doi: 10.3189/2012JoG11J088] and its implementation in 
SICOPOLIS. The latter has not been finished yet, there are still some 
things to be done (proper tracking of the CTS = cold-temperate 
transition surface, different diffusivities in cold and temperate ice, 
some I/O and documentation), but in principle it already works. This 
means that SICOPOLIS now features three different thermodynamics options:

   CALCMOD==0: Cold-ice method
   CALCMOD==1: Polythermal method with separate numerical domains for
               cold ice (kc=0...KCMAX) and temperate ice (kt=0...KTMAX)
   CALCMOD==2: Enthalpy method; polythermal, but with only one domain
               (kc=0...KCMAX) for both cold and temperate ice

Since I'll be at the AGU Fall Meeting next week, the remaining work on 
the implementation of the enthalpy method must wait until after that. 
Meanwhile, please feel free to check out the "bleeding edge" of the code 
and play a bit with the new method (e.g., by setting CALCMOD to 2 in the 
header files in the repository).

If you come across any issues, please let me know.

Thanks & best regards

Ralf

-- 
Dr. Ralf Greve
Professor, Glacier and Ice Sheet Research
Institute of Low Temperature Science
Hokkaido University
Kita-19, Nishi-8, Kita-ku
Sapporo 060-0819, Japan
Phone: (+81)-(0)11-706-6891
greve at lowtem.hokudai.ac.jp
http://wwwice.lowtem.hokudai.ac.jp/~greve/

++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++