BISMG:SarahS/weekly reports

Aims
Testing the ability of glimmer-cism to predict changes in ice thickness and velocity for Antarctica over the period 1980 to 2004.

Simulating changes in the Antarctic ice sheet for the next 200 years using data from regional climate models for the A1B and E1 SRES scenarios.

A sensitivity study of key model parameters in glimmer-cism will be carried out and I will also help support an ice model inter-comparison project with Tamsin Edwards.

This work forms part of the ice2sea project.

Week ending 13th November 2009
Ran glimmer-cism using the glint driver. Forcing with one year NCEP data repeated for 10 years. Initial conditions 1) Ice surface elevation 5km provided by Anne Le Brocq (5km resolution) and 2)  Bed Topography from BEDMAP1_plus  (5km resolution).

I had to run on the 20km resolution because I encountered a run time error about insufficient virtual memory when I ran on a 5km resolution. Perhaps 5km is too fine a resolution to run the model or I need to look for a bug. Prepared ERA-40 monthly 2m temp and total precipitation data for input into glimmer-cism for the years 1958-2002. Changed the driver to read in ERA-40 data.

Week ending 20th November 2009
I made changes to the floation equation in glimmer-cism to account for firn snow. These are the changes

1.	In glide_io.F90 : read in firn snow from netCDF file

2.	In glide_vars.def : add extra variable data%geometry%firn to the list

3.	Change density of ice to value used in Anne's data(original value 910kgm-3, new rhoi=918kgm-3)

4.	In glide_setup.F90 insert new subroutine called glide_calclsrf_firn(thck,topg,eus,firn,lsrf)

!*FD Calculates the elevation of the lower surface of the ice, !*FD by considering whether it is floating or not, modify to   !*FD include firn

use glimmer_global, only : dp   use glimmer_physcon, only : rhoi, rhoo implicit none

real(dp), intent(in), dimension :: thck !*FD Ice thickness real(dp), intent(in), dimension :: topg !*FD Bedrock topography elevation real, intent(in)                     :: eus  !*FD global sea level real(dp), intent(in), dimension :: firn !*FD Firn depth real(dp), intent(out), dimension :: lsrf !*FD Lower ice surface elevation

real(dp), parameter :: con = - rhoi / rhoo

where (topg-eus < con * thck + firn) lsrf = con * thck + firn elsewhere lsrf = topg end where

end subroutine glide_calclsrf_firn

5.	In glide_mask.F90 modify subroutine glide_set_mask

!Identify points where the ice is floating or where there is open ocean including firn

where (topg - eus < con * thck + firn) mask = ior(mask, GLIDE_MASK_OCEAN) elsewhere mask = ior(mask, GLIDE_MASK_LAND) endwhere

Test pixels to check if calculated usrf and lsrf equals data provided by Anne.

Looks ok!

Week ending 27th Nov 2009
Installed an updated version of glimmer-cism.

Modified the code to read in fixed temperature and accumulation data from Anne's dataset.

I’ve been comparing the re-factored temperature solver to the original code to get a handle on what changes have been implemented.

Week ending 4th Dec 2009
Started a new wiki section for glimmer-cism notes.

Made modifications to glimmer-cism Makefile to remove python.

Some of the week has been spent writing a paper on dust modeling and tuning.