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

From greve at lowtem.hokudai.ac.jp  Tue Feb  4 09:16:10 2014
From: greve at lowtem.hokudai.ac.jp (Ralf Greve)
Date: Tue, 04 Feb 2014 17:16:10 +0900
Subject: [SICOPOLIS:29] Enthalpy method now stable
Message-ID: <52F0A1CA.7090202@lowtem.hokudai.ac.jp>

Dear colleagues,

coming back to the enthalpy method in SICOPOLIS (see [SICOPOLIS:27] of 
2013-12-06):

The implementation is now complete and stable. Heinz Blatter and I have 
been doing a lot of testing in order to fine-tune the method, and it 
works really nice now. For the EISMINT Phase 2 Simplied Geometry 
Experiment A, the polythermal run specified by the header file 
sico_specs_v32_emtp2sge25_expA.h (CALCMOD==1, CTS_MELTING_FREEZING==2) 
and a similar run with the enthalpy method (CALCMOD==2, ENTHMOD==1) 
produce temperate ice volumes V_temp that differ by less than 1%.

Please feel free to check out the latest revision of the code and 
experiment with the new method. 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/

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

From greve at lowtem.hokudai.ac.jp  Tue Feb  4 09:18:45 2014
From: greve at lowtem.hokudai.ac.jp (Ralf Greve)
Date: Tue, 04 Feb 2014 17:18:45 +0900
Subject: [SICOPOLIS:30] Enthalpy method now stable
Message-ID: <52F0A265.5060005@lowtem.hokudai.ac.jp>

(Some text was broken in the previous mail, so that I'm resending it. 
Apologies!)

Dear colleagues,

coming back to the enthalpy method in SICOPOLIS (see [SICOPOLIS:27] of 
2013-12-06):

The implementation is now complete and stable. Heinz Blatter and I have 
been doing a lot of testing in order to fine-tune the method, and it 
works really nice now. For the EISMINT Phase 2 Simplified Geometry 
Experiment A, the polythermal run specified by the header file 
sico_specs_v32_emtp2sge25_expA.h (CALCMOD==1, CTS_MELTING_FREEZING==2) 
and a similar run with the enthalpy method (CALCMOD==2, ENTHMOD==1) 
produce temperate ice volumes V_temp that differ by less than 1%.

Please feel free to check out the latest revision of the code and 
experiment with the new method. 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/

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