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)
(Construct and add a new module to Glimmer)
 
(9 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).
 
*The new module will be used with velocities from the HO model to evolve a simple ice shelf domain (image below).
*Evolution using the 1st-order advection scheme will be compared with that from a higher-order advection scheme.
+
 
 
*For step-by-step notes on adding this module, go to [[adding dH/dt module]].
 
*For step-by-step notes on adding this module, go to [[adding dH/dt module]].
 
  
  
 
----
 
----
  
 
+
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.'''
 
+
 
+
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.  
+
  
 
[[Image:Confined-shelf.ini.jpg]]
 
[[Image:Confined-shelf.ini.jpg]]
 
 
  
  
Line 23: Line 19:
  
 
[[Image:Confined-shelf.100a_no_acab.jpg]]
 
[[Image:Confined-shelf.100a_no_acab.jpg]]
 
 
 
 
 
Thickness and velocity fields after 100 yrs of evolution with uniform accumulation of 0.5 m/a. Evolution was calculated using the 1st-order advection scheme.
 
 
[[Image:Confined-shelf.100a_5_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