Difference between revisions of "Adding a module to Glimmer I"

From Interactive System for Ice sheet Simulation
Jump to: navigation, search
(Construct and add a new module to Glimmer)
 
(17 intermediate revisions by one user not shown)
Line 1: Line 1:
 
==Construct and add a new module to Glimmer==
 
==Construct and add a new module to Glimmer==
 +
 
*For this exercise, we will create and add a 1st-order, upwind advection scheme to the model
 
*For this exercise, we will create and add a 1st-order, upwind advection scheme to the model
*The new module will be used with velocities from the HO model to evolve a simple ice shelf domain (image below) to a steady-state
 
*The evolution and steady-state using the 1st-order advection scheme will be compared with those from a higher-order advection scheme
 
*For step-by-step notes on adding this module, go to [[adding dH/dt module]].
 
  
 +
*The new module will be used with velocities from the HO model to evolve a simple ice shelf domain (image below).
 +
 +
*For step-by-step notes on adding this module, go to [[adding dH/dt module]].
  
  
 
----
 
----
  
Confined ice shelf solution (diagnostic velocity in m/yr) from Glimmer/CISM HO model (top, color) compared to ice shelf solution from EISMINT-ROSS benchmarks using SSA model (bottom, B&W).
+
Initial thickness, across-flow, and along-flow velocity fields for confined ice shelf (500 m thick, uniform T=253 K). Boundary conditions are 0-flux at top, left, right, and shelf BCs at bottom. Note that this initial condition is essentially the same solution to the unpublished EISMINT ice-shelf experiments [[COMSOL activities#Shallow shelf approximation|that you worked with in COMSOL a few days ago]]. '''NOTE: Due to a plotting error, the velocities shown here will be the ''negative'' of what you will see with NCVIEW for the ice shelf exercise.'''
 +
 
 +
[[Image:Confined-shelf.ini.jpg]]
 +
 
 +
 
 +
 
 +
Thickness and velocity fields after 100 yrs of evolution with no accumulation. Evolution was calculated using the 1st-order advection scheme.
  
[[Image:shelf-benchmark.jpg|350 px]]
+
[[Image:Confined-shelf.100a_no_acab.jpg]]

Latest revision as of 17:42, 13 August 2009

Construct and add a new module to Glimmer

  • For this exercise, we will create and add a 1st-order, upwind advection scheme to the model
  • The new module will be used with velocities from the HO model to evolve a simple ice shelf domain (image below).



Initial thickness, across-flow, and along-flow velocity fields for confined ice shelf (500 m thick, uniform T=253 K). Boundary conditions are 0-flux at top, left, right, and shelf BCs at bottom. Note that this initial condition is essentially the same solution to the unpublished EISMINT ice-shelf experiments that you worked with in COMSOL a few days ago. NOTE: Due to a plotting error, the velocities shown here will be the negative of what you will see with NCVIEW for the ice shelf exercise.

Confined-shelf.ini.jpg


Thickness and velocity fields after 100 yrs of evolution with no accumulation. Evolution was calculated using the 1st-order advection scheme.

Confined-shelf.100a no acab.jpg