http://websrv.cs.umt.edu/isis/api.php?action=feedcontributions&user=Adamc&feedformat=atomInteractive System for Ice sheet Simulation - User contributions [en]2020-06-04T02:33:43ZUser contributionsMediaWiki 1.21.1http://websrv.cs.umt.edu/isis/index.php/COMSOL_MultiphysicsCOMSOL Multiphysics2011-07-23T00:12:48Z<p>Adamc: added sections for development</p>
<hr />
<div>COMSOL Multiphysics is a commercially available finite element analysis and solver software package.<br />
<br />
==Lessons relating to COMSOL==<br />
Introductory COMSOL lessons developed for the [[Summer Modeling School]]<br />
*[[COMSOL introduction]]<br />
*[[COMSOL activities]]<br />
<br />
==Model Library==<br />
Sample COMSOL models and documentation<br />
<br />
===Grounded Ice Models===<br />
*[[2-D profile with deformation only]]<br />
*[[2-D profile with Budd-type sliding law]]<br />
*[[2-D profile with a stretched vertical coordinate system]]<br />
<br />
===Floating Ice Models===<br />
*[[Simple ice shelf]]<br />
<br />
==External Links==<br />
[http://en.wikipedia.org/wiki/COMSOL_Multiphysics COMSOL Multiphysics] - Wikipedia entry</div>Adamchttp://websrv.cs.umt.edu/isis/index.php/COMSOL_MultiphysicsCOMSOL Multiphysics2011-07-21T16:01:14Z<p>Adamc: </p>
<hr />
<div>COMSOL Multiphysics is a commercially available finite element analysis and solver software package.<br />
<br />
==Lessons relating to COMSOL==<br />
Introductory COMSOL lessons developed for the [[Summer Modeling School]]<br />
*[[COMSOL introduction]]<br />
*[[COMSOL activities]]<br />
<br />
==Model Library==<br />
Sample COMSOL models and documentation<br />
<br />
===Grounded Ice===<br />
*[[2-D deformation only]]<br />
<br />
==External Links==<br />
[http://en.wikipedia.org/wiki/COMSOL_Multiphysics COMSOL Multiphysics] - Wikipedia entry</div>Adamchttp://websrv.cs.umt.edu/isis/index.php/Talk:COMSOL_MultiphysicsTalk:COMSOL Multiphysics2011-07-21T05:42:20Z<p>Adamc: New page: I'm attempting to revamp this page to allow for members of the glaciology community to upload models to demonstrate with simple models the types of physics that have been successfully impl...</p>
<hr />
<div>I'm attempting to revamp this page to allow for members of the glaciology community to upload models to demonstrate with simple models the types of physics that have been successfully implemented in COMSOL.<br />
--[[User:Adamc|Adamc]] 05:42, 21 July 2011 (UTC)</div>Adamchttp://websrv.cs.umt.edu/isis/index.php/COMSOL_MultiphysicsCOMSOL Multiphysics2011-07-21T05:37:39Z<p>Adamc: </p>
<hr />
<div>COMSOL Multiphysics is a commercially available finite element analysis and solver software package.<br />
<br />
==Lessons relating to COMSOL==<br />
Introductory COMSOL lessons developed for the [[Summer Modeling School]]<br />
*[[COMSOL introduction]]<br />
*[[COMSOL activities]]<br />
<br />
==Model Library==<br />
<br />
==External Links==<br />
[http://en.wikipedia.org/wiki/COMSOL_Multiphysics COMSOL Multiphysics] - Wikipedia entry</div>Adamchttp://websrv.cs.umt.edu/isis/index.php/COMSOL_MultiphysicsCOMSOL Multiphysics2011-07-21T05:35:29Z<p>Adamc: edit to wikipedia link</p>
<hr />
<div>COMSOL Multiphysics is a commercially available finite element analysis and solver software package.<br />
<br />
==Lessons relating to COMSOL==<br />
*[[COMSOL introduction]]<br />
*[[COMSOL activities]]<br />
<br />
==External Links==<br />
[http://en.wikipedia.org/wiki/COMSOL_Multiphysics COMSOL Multiphysics] - Wikipedia entry</div>Adamchttp://websrv.cs.umt.edu/isis/index.php/COMSOL_MultiphysicsCOMSOL Multiphysics2011-07-21T05:34:08Z<p>Adamc: add wikipedia entry</p>
<hr />
<div>COMSOL Multiphysics is a commercially available finite element analysis and solver software package.<br />
<br />
==Lessons relating to COMSOL==<br />
*[[COMSOL introduction]]<br />
*[[COMSOL activities]]<br />
<br />
==External Links==<br />
[http://en.wikipedia.org/wiki/COMSOL_Multiphysics link COMSOL on Wikipedia]</div>Adamchttp://websrv.cs.umt.edu/isis/index.php/COMSOL_MultiphysicsCOMSOL Multiphysics2011-07-21T05:32:36Z<p>Adamc: added introductory material</p>
<hr />
<div>COMSOL Multiphysics is a commercially available finite element analysis and solver software package.<br />
<br />
==Lessons relating to COMSOL==<br />
*[[COMSOL introduction]]<br />
*[[COMSOL activities]]</div>Adamchttp://websrv.cs.umt.edu/isis/index.php/SoftwareSoftware2011-07-21T05:27:15Z<p>Adamc: </p>
<hr />
<div>[[Image:mutant_model.jpg|right|thumb|350 px|A community model is a collection of original and borrowed parts. CISM efforts began with the [http://forge.nesc.ac.uk/projects/glimmer/ Glimmer ice sheet model], borrowed components from [http://climate.lanl.gov/Models/CICE/ CICE] and [http://homepages.ulb.ac.be/~fpattyn/3dhomodel.html ICE3D<sup>HO</sup>], and uses infrastructure from [http://www.unidata.ucar.edu/software/netcdf/ NetCDF] and [http://www.cise.ufl.edu/research/sparse/umfpack/ UMPFPACK]. ]]<br />
==Community Ice Sheet Model==<br />
''[http://oceans11.lanl.gov/trac/CISM The Community Ice Sheet Model] is being developed at [http://www.lanl.gov Los Alamos National Laboratories], [http://www.umt.edu University of Montana, Missoula] and elsewhere. The effort began from the [http://forge.nesc.ac.uk/projects/glimmer/ Glimmer ice sheet model] developed in Bristol by Tony Payne, Ian Rutt, and Magnus Hagdorn. ''<br />
<br />
The core improvements that CISM provides are:<br />
===[[Higher Order Physics]]===<br />
===[[Parallel Solvers]]===<br />
===[[Validation and Verification]]===<br />
===[[Climate Drivers]]===<br />
===[[Basal Water Modeling]]===<br />
===[[Evaluating Model Output]]===<br />
===[[Code Refactoring]]===<br />
===[[Grounding Line Migration]]===<br />
<br />
==ISIS==<br />
===ISIS===<br />
[[Image:ISISViz.jpg|thumb|right|400 px|The ISIS Visualization provides quick viewing of model output, even while the simulation is running.]]<br />
[[Image:ISISConfig.png|thumb|right|400 px| ISIS provides an intuitive interface for setting the values of important simulation parameters and model settings.]]<br />
<br />
'''ISIS''' is an acronym for ''Interactive System for Ice sheet Simulation''. [[Wikipedia:Isis| Isis]] also happens to be the Egyptian goddess of fertility. This is an appropriate name in that our objective is to produce more modelers. ISIS attempts to spawn more modelers by lowering the barrier to trying ice sheet modeling by <br />
* Providing a compiled version of an ice sheet model on Mac, Linux, and Windows platforms.<br />
* Provide interesting data sets and scenarios that can be run with the model.<br />
* Provide some simple tools for visualizing the model results.<br />
<br />
====Download====<br />
Versions of ISIS have been prepared for Linux, OS X, and Windows operating systems. Download them here.<br />
<br />
{|<br />
|-valign="top"<br />
|[[Image:Mac.jpg|thumb|100 px|[[Media:ISISMac.zip]]<br>Mac OS X version of ISIS, Intel Mac only.]]<br />
|[[Image:Windows.jpg|thumb|150 px|[[Media:ISISWindows.zip]]<br>Windows version of ISIS]]<br />
|[[Image:Linux.jpg|thumb|150 px|[[Media:ISISLinux.tar.gz]]<br>Linux version of ISIS]]<br />
|}<br />
<br><br><br />
<br />
===ETSU ISIS===<br />
The [http://www.cs.etsu.edu/ Computer Science Department at East Tennessee State University] has made several changes to ISIS. ETSU ISIS provides a<br />
* command line option for enabling educational mode, <br />
* a generic model class (PISM and Glimmer included), and<br />
* various bug fixes<br />
<br />
For more information on ETSU ISIS see [[Software/ETSUISIS/Model]]<br />
<br />
====Download====<br />
ETSU ISIS has been tested on Windows and OS X<br><br />
Downloads will be available here soon<br />
<br />
===eduISIS===<br />
There is potential for Earth systems models to provide instruction for students in a range of subjects and at a range of levels. Seeking to target high school students, eduISIS was developed as an extension of ISIS with a simplified user interface. <br />
<br />
<br />
{|<br />
|-valign="top"<br />
|[[Image:Mac.jpg|thumb|100 px|[[Media:eduISISMac.dmg]]<br>Mac OS X version of eduISIS, Intel Mac only.]]<br />
|[[Image:Windows.jpg|thumb|100 px|[[Media:eduISISWindows.zip]]<br>Windows version of eduISIS]]<br />
<!--|[[Image:Linux.jpg|thumb|100 px|[[Media:ISISLinux.tar.gz]]<br>Linux version of eduISIS]] --><br />
|}<br />
<br />
==Ice Sheet Models==<br />
===PISM===<br />
The [https://gna.org/projects/pism/ Parallel Ice Sheet Model] being developed at the [http://www.uaf.edu/ University of Alaska, Fairbanks] has provided CISM with many fresh ideas for model verification and data interchange. Ultimately, it will be desirable to have a common interface to both models, allowing for rapid cross model comparisons. For the time being, the first step is a common data format. This is largely completed, but remains to be tested in more depth.<br />
<br />
===Glimmer===<br />
CISM efforts began from the [http://forge.nesc.ac.uk/projects/glimmer/ Glimmer ice sheet model] developed in Bristol by Tony Payne, Ian Rutt, and Magnus Hagdorn. This is a thermo-mechanically coupled shallow ice model that uses an implicit finite difference scheme to solve the field equations. Beyond this level of simulation, Glimmer has a number of attractive software features including;<br />
<br />
*NetCDF support,<br />
*Automake build system,<br />
*modular architecture created in FORTRAN 90,<br />
*Open Source, GPL license,<br />
*designed with interface to climate models in mind, and<br />
*set of test suites<br />
<br />
===COMSOL Multiphysics===<br />
[[COMSOL Multiphysics]] is a commercially available finite element solver. <br />
<br />
==Glimmer Netcdf Data Conversion Program==<br />
<br />
===Data Conversion===<br />
Data used by Glimmer is initially stored in separate data files with a ".dat" file extension. The advantage to the data conversion is to eliminate having multiple files by merging them into one Netcdf file. The data converting program is a Java application that allows users to open data files and convert them to Netcdf format. The conversion program examines each data file and reads the format instruction at the beginning of each file. This information is read based on a parse tree that is created using a data file with the rules for the instruction's grammar. The file's data is then read and stored in a multi-dimensional ArrayList. This multi-dimensional ArrayList is then written to a Netcdf file.<br />
<br />
==Glimmer Netcdf Data Portability Analysis==<br />
<br />
===ESML===<br />
ESML stands for Earth Science Markup Language. "ESML is an interchange technology that enables data (both structural and semantic) interoperability with applications without enforcing a standard format within the Earth Science community".[1] ESML is basically a metadata file written in the XML grammar syntax which describes the contents of another data file. ESML metadata file can be used by different data readers with the help of ESML library. This way data readers can read any data format that has an ESML file associated with it. ESML library is currently supported in C++.<br />
An ESML file for the glimmer netcdf data can be created making the glimmer data more portable.<br />
<br />
===ncmL, ncML-GML===<br />
ncML is an XML based language for netCDF data files. It is also a metadata language which explains the contents of an associated netCDF file. <br />
<br />
ncML-GML is a hybrid technology of ncML and GML. ncML-GML uses both ncML and GML grammer. Study of this technology is in progress and more information will be updated very soon.<br />
<br />
== Other (Post Processing / Visualization) ==<br />
<br />
See the [[Data#NetCDF_Viewers]] section of the [[Data]] page.<br />
<br />
==References==<br />
[1]http://esml.itsc.uah.edu/</div>Adamchttp://websrv.cs.umt.edu/isis/index.php/SoftwareSoftware2011-07-21T05:23:32Z<p>Adamc: added COMSOL section</p>
<hr />
<div>[[Image:mutant_model.jpg|right|thumb|350 px|A community model is a collection of original and borrowed parts. CISM efforts began with the [http://forge.nesc.ac.uk/projects/glimmer/ Glimmer ice sheet model], borrowed components from [http://climate.lanl.gov/Models/CICE/ CICE] and [http://homepages.ulb.ac.be/~fpattyn/3dhomodel.html ICE3D<sup>HO</sup>], and uses infrastructure from [http://www.unidata.ucar.edu/software/netcdf/ NetCDF] and [http://www.cise.ufl.edu/research/sparse/umfpack/ UMPFPACK]. ]]<br />
==Community Ice Sheet Model==<br />
''[http://oceans11.lanl.gov/trac/CISM The Community Ice Sheet Model] is being developed at [http://www.lanl.gov Los Alamos National Laboratories], [http://www.umt.edu University of Montana, Missoula] and elsewhere. The effort began from the [http://forge.nesc.ac.uk/projects/glimmer/ Glimmer ice sheet model] developed in Bristol by Tony Payne, Ian Rutt, and Magnus Hagdorn. ''<br />
<br />
The core improvements that CISM provides are:<br />
===[[Higher Order Physics]]===<br />
===[[Parallel Solvers]]===<br />
===[[Validation and Verification]]===<br />
===[[Climate Drivers]]===<br />
===[[Basal Water Modeling]]===<br />
===[[Evaluating Model Output]]===<br />
===[[Code Refactoring]]===<br />
===[[Grounding Line Migration]]===<br />
<br />
==ISIS==<br />
===ISIS===<br />
[[Image:ISISViz.jpg|thumb|right|400 px|The ISIS Visualization provides quick viewing of model output, even while the simulation is running.]]<br />
[[Image:ISISConfig.png|thumb|right|400 px| ISIS provides an intuitive interface for setting the values of important simulation parameters and model settings.]]<br />
<br />
'''ISIS''' is an acronym for ''Interactive System for Ice sheet Simulation''. [[Wikipedia:Isis| Isis]] also happens to be the Egyptian goddess of fertility. This is an appropriate name in that our objective is to produce more modelers. ISIS attempts to spawn more modelers by lowering the barrier to trying ice sheet modeling by <br />
* Providing a compiled version of an ice sheet model on Mac, Linux, and Windows platforms.<br />
* Provide interesting data sets and scenarios that can be run with the model.<br />
* Provide some simple tools for visualizing the model results.<br />
<br />
====Download====<br />
Versions of ISIS have been prepared for Linux, OS X, and Windows operating systems. Download them here.<br />
<br />
{|<br />
|-valign="top"<br />
|[[Image:Mac.jpg|thumb|100 px|[[Media:ISISMac.zip]]<br>Mac OS X version of ISIS, Intel Mac only.]]<br />
|[[Image:Windows.jpg|thumb|150 px|[[Media:ISISWindows.zip]]<br>Windows version of ISIS]]<br />
|[[Image:Linux.jpg|thumb|150 px|[[Media:ISISLinux.tar.gz]]<br>Linux version of ISIS]]<br />
|}<br />
<br><br><br />
<br />
===ETSU ISIS===<br />
The [http://www.cs.etsu.edu/ Computer Science Department at East Tennessee State University] has made several changes to ISIS. ETSU ISIS provides a<br />
* command line option for enabling educational mode, <br />
* a generic model class (PISM and Glimmer included), and<br />
* various bug fixes<br />
<br />
For more information on ETSU ISIS see [[Software/ETSUISIS/Model]]<br />
<br />
====Download====<br />
ETSU ISIS has been tested on Windows and OS X<br><br />
Downloads will be available here soon<br />
<br />
===eduISIS===<br />
There is potential for Earth systems models to provide instruction for students in a range of subjects and at a range of levels. Seeking to target high school students, eduISIS was developed as an extension of ISIS with a simplified user interface. <br />
<br />
<br />
{|<br />
|-valign="top"<br />
|[[Image:Mac.jpg|thumb|100 px|[[Media:eduISISMac.dmg]]<br>Mac OS X version of eduISIS, Intel Mac only.]]<br />
|[[Image:Windows.jpg|thumb|100 px|[[Media:eduISISWindows.zip]]<br>Windows version of eduISIS]]<br />
<!--|[[Image:Linux.jpg|thumb|100 px|[[Media:ISISLinux.tar.gz]]<br>Linux version of eduISIS]] --><br />
|}<br />
<br />
==Ice Sheet Models==<br />
===PISM===<br />
The [https://gna.org/projects/pism/ Parallel Ice Sheet Model] being developed at the [http://www.uaf.edu/ University of Alaska, Fairbanks] has provided CISM with many fresh ideas for model verification and data interchange. Ultimately, it will be desirable to have a common interface to both models, allowing for rapid cross model comparisons. For the time being, the first step is a common data format. This is largely completed, but remains to be tested in more depth.<br />
<br />
===Glimmer===<br />
CISM efforts began from the [http://forge.nesc.ac.uk/projects/glimmer/ Glimmer ice sheet model] developed in Bristol by Tony Payne, Ian Rutt, and Magnus Hagdorn. This is a thermo-mechanically coupled shallow ice model that uses an implicit finite difference scheme to solve the field equations. Beyond this level of simulation, Glimmer has a number of attractive software features including;<br />
<br />
*NetCDF support,<br />
*Automake build system,<br />
*modular architecture created in FORTRAN 90,<br />
*Open Source, GPL license,<br />
*designed with interface to climate models in mind, and<br />
*set of test suites<br />
<br />
===COMSOL Multiphysics===<br />
<br />
==Glimmer Netcdf Data Conversion Program==<br />
<br />
===Data Conversion===<br />
Data used by Glimmer is initially stored in separate data files with a ".dat" file extension. The advantage to the data conversion is to eliminate having multiple files by merging them into one Netcdf file. The data converting program is a Java application that allows users to open data files and convert them to Netcdf format. The conversion program examines each data file and reads the format instruction at the beginning of each file. This information is read based on a parse tree that is created using a data file with the rules for the instruction's grammar. The file's data is then read and stored in a multi-dimensional ArrayList. This multi-dimensional ArrayList is then written to a Netcdf file.<br />
<br />
==Glimmer Netcdf Data Portability Analysis==<br />
<br />
===ESML===<br />
ESML stands for Earth Science Markup Language. "ESML is an interchange technology that enables data (both structural and semantic) interoperability with applications without enforcing a standard format within the Earth Science community".[1] ESML is basically a metadata file written in the XML grammar syntax which describes the contents of another data file. ESML metadata file can be used by different data readers with the help of ESML library. This way data readers can read any data format that has an ESML file associated with it. ESML library is currently supported in C++.<br />
An ESML file for the glimmer netcdf data can be created making the glimmer data more portable.<br />
<br />
===ncmL, ncML-GML===<br />
ncML is an XML based language for netCDF data files. It is also a metadata language which explains the contents of an associated netCDF file. <br />
<br />
ncML-GML is a hybrid technology of ncML and GML. ncML-GML uses both ncML and GML grammer. Study of this technology is in progress and more information will be updated very soon.<br />
<br />
== Other (Post Processing / Visualization) ==<br />
<br />
See the [[Data#NetCDF_Viewers]] section of the [[Data]] page.<br />
<br />
==References==<br />
[1]http://esml.itsc.uah.edu/</div>Adamchttp://websrv.cs.umt.edu/isis/index.php/Adam_and_Kristin%27s_UW_BrownbagAdam and Kristin's UW Brownbag2009-10-28T02:54:37Z<p>Adamc: </p>
<hr />
<div>[[Image:GLunch.pdf]]</div>Adamchttp://websrv.cs.umt.edu/isis/index.php/Adam_and_Kristin%27s_UW_BrownbagAdam and Kristin's UW Brownbag2009-10-28T02:54:06Z<p>Adamc: New page: Image:Glunch.pdf</p>
<hr />
<div>[[Image:Glunch.pdf]]</div>Adamchttp://websrv.cs.umt.edu/isis/index.php/File:GLunch.pdfFile:GLunch.pdf2009-10-28T02:52:44Z<p>Adamc: Presentation Adam Campbell and Kristin Poinar gave at the University of Washington Glaciology Group Lunch, 13 Oct 2009.</p>
<hr />
<div>Presentation Adam Campbell and Kristin Poinar gave at the University of Washington Glaciology Group Lunch, 13 Oct 2009.</div>Adamchttp://websrv.cs.umt.edu/isis/index.php/Summer_Modeling_SchoolSummer Modeling School2009-10-28T02:50:41Z<p>Adamc: /* Student Participants */ added section heading</p>
<hr />
<div>[[Image:Portland.jpg|thumb|right|400 px|The summer ice sheet modeling school will be held in Portland Oregon, August 3-14, 2009]]<br />
<br />
==Overview==<br />
Ice Sheet Models for the 21st Century Summer School will be an intensive short course that brings current and future ice-sheet scientists together to develop better models for the projection of future sea-level rise (slr). The IPCC Fourth Assessment Report [http://www.ipcc.ch/ipccreports/ar4-syr.htm] acknowledged that current models do not adequately treat the dynamic response of ice sheets to climate change, and that this is the largest uncertainty in assessing potential rapid sea-level rise. Recognizing this, an ice-sheet modelling Workshop was held during the July 2008 SCAR/IASC [https://www.comnap.aq/content/events/osc2008] meeting, in St. Petersburg, Russia. This meeting developed a community strategy on how best to (i) improve the physical understanding of ice-sheet processes responsible for rapid change; (ii) incorporate improved physical understanding into numerical models; (iii) assimilate appropriate data into the models for calibration and validation; and (iv) develop prognostic whole ice-sheet models that better incorporate non-linear ice-sheet response to environmental forcing (such as change in surface mass balance, loss of buttressing from floating ice shelves and ice tongues, and rising sea level). <br />
<br />
The two-week Summer School is a first step towards implementing this strategy. It will bring scientists from differing backgrounds together and allow more extensive and in-depth interactions between the relevant scientific research communities. A series of general background lectures as well as discussions of more specialized and advanced topics during this Summer School will provide the foundation for cross-disciplinary research, particularly for early career scientists. We anticipate publication of lecture notes both in hard copy and on a dedicated home page, to provide the glaciological community with an up-to-date overview of the science and observational techniques that will serve to guide further research efforts. Direct beneficiaries will be young researchers; indirect beneficiaries will be coastal zone communities who will gain improved sea level change forecasts to underpin their plans for sustainable development.<br />
<br />
===Venue===<br />
The modeling school was held on the campus of [[Wikipedia:Portland State University|Portland State University]] in [[Wikipedia:Portland, Oregon|Portland, Oregon]] August 3-14, 2009.<br />
<br />
* [http://maps.google.com/maps?f=d&source=s_d&saddr=Portland+Airport&daddr=310+SW+Lincoln+St,+Portland,+OR+97201-5007+(University+Place-Portland)&geocode=&hl=en&mra=ls&dirflg=r&date=07%2F28%2F09&time=8:59am&ttype=dep&noexp=0&noal=0&sort=&tline=&sll=45.54878,-122.629155&sspn=0.092445,0.144367&ie=UTF8&ll=45.548679,-122.619438&spn=0.092445,0.144367&z=13&start=0 Map] from airport to [http://cegs.pdx.edu/stay/upl/ University Place Hotel] using public transport (note that the directions in your travel letter are better than the Google generated instructions here).<br />
<br />
* [http://maps.google.com/maps?f=d&source=s_d&saddr=310+SW+Lincoln+St,+Portland,+OR+97201-5007+(University+Place-Portland)&daddr=1721+SW+Broadway,+Portland,+OR+97201+(Cramer+Hall)&hl=en&geocode=FdVhtgIdZwqw-CHO0mMQPCwi0Q%3BFRN3tgIdvP2v-CHxCBg32xEzXA&mra=ls&dirflg=w&sll=45.51029,-122.681675&sspn=0.005782,0.009023&ie=UTF8&ll=45.510091,-122.68232&spn=0.005782,0.009023&z=17 Map] from [http://cegs.pdx.edu/stay/upl/ University Place Hotel] to [http://www.pdx.edu/campus-map Cramer Hall].<br />
<br />
=== Student Participants ===<br />
<br />
*[[Nominations for Student/Instructor Awards]]<br />
*[[Student Bios]]<br />
*[[Student Presentation Development]]<br />
*[[Groups]] example of [[connections in groups]]<br />
*[[Terminology]]<br />
*[[Questions]]<br />
*[[Adam and Kristin's UW Brownbag]]<br />
<br />
===Lectures and Planned Activities===<br />
<br />
For information about editing this page, see [[Wikipedia:How to edit]].<br />
<br />
You can view a summary of the subjects covered at the Summer Modeling School by clicking on the following link. [[Subject Overview]]<br />
<br />
{| border="1" cellpadding="5" cellspacing="0"<br />
|-valign="top" style="background:RoyalBlue"<br />
!width="20%"|Dates<br />
!width="25%"|Lecture Topics<br />
!width="15%"|Lecturers<br />
!width="25%"|Laboratory Topics<br />
!width="15%"|Laboratory Instructors <br />
|-valign="top" style="background:AliceBlue"<br />
| [[4-5 August]]<br />
| Introduction to and theoretical basis for ice sheet modeling. <br />
| Kees van der Veen, [[Nina Kirchner]] <br />
| [[Finite differencing|Finite differencing]] and [[Pragmatic Programming|pragmatic programming]] using Fortran[http://en.wikipedia.org/wiki/Fortran] 95...<br />
computing divergence and gradient...<br />
from conservation equation to matrix algebra...<br />
rheology and that which makes ice ice...<br />
simple, ideal models...<br />
that which makes ice-sheet modeling hard...<br />
| Gethin Williams, [[Ian Rutt]], [[Jesse Johnson]]<br />
|-valign="top" style="background:PowderBlue"<br />
| 6 August <br />
| [[Basal Conditions]], [[Data sets for ice sheet modeling]]<br />
| Alan Rempel, Slawek Tulaczyk and Ken Jezek<br />
| [[COMSOL Multiphysics]]<br />
| Olga Sergienko and Jesse Johnson<br />
|-valign="top" style="background:AliceBlue"<br />
| 7 August<br />
| The world of [[ice shelves]] and 'distributed stress-field solutions'. [[Modelling mountain glaciers]].<br />
| Todd Dupont, Olga Sergienko, and Brian Anderson<br />
| Linear Algebra of ice-sheet modeling, relaxation methods, finite-element methodology, solution of Laplace equation in arbitrary domain, creation of an ice-shelf flow-field model (snap shot of flow field), Models of the Ross Ice Shelf<br />
| Olga Sergienko and Todd Dupont<br />
<br />
|-valign="top" style="background:PowderBlue"<br />
| 8 August<br />
| [[Student Presentation]]<br />
| [[Modeling School Students]]<br />
| open work day with breakfast at 8 am & student presentation at 9 am<br />
| go to the farmer's market<br />
|-valign="top" style="background:AliceBlue"<br />
| 9 August<br />
| Free day; possible PDX tour<br />
|<br />
|<br />
|<br />
|-valign="top" style="background:PowderBlue"<br />
| 10 August<br />
| Excursion to Mt. Hood and [[Eliot Glacier field trip]] <br> Meet at 7:00 am for breakfast in Cramer Hall. Bring a jacket, water, and sensible (supportive) shoes.<br> Box lunches are provided and we will stop for dinner on the return trip (bring cash). <br />
| Guided by [http://web.pdx.edu/~basagic/ Hassan Basagic]<br />
|<br />
|<br />
|-valign="top" style="background:AliceBlue"<br />
| [[Quantifying_model_uncertainty |11 August]]<br />
| [[Quantifying model uncertainty]]<br />
| Charles Jackson and [http://www.mit.edu/~heimbach Patrick Heimbach]<br />
| Uncertain lab, [[Dynamic response to the enhanced basal flow in the Greenland ice sheet]] Weili Wang<br />
| Charles Jackson, Patrick Heimbach, and Weili Wang<br />
|-valign="top" style="background:PowderBlue"<br />
| [[12-13 August]]<br />
| Introduction to Glimmer-CISM ([[Introduction to Glimmer I|Part I]], [[Introduction to Glimmer II|Part II]]); [[Higher order velocity schemes|Higher-order models]]<br />
| [[Ian Rutt]], [[Magnus Hagdorn]], [[Stephen Price]], Bill Lipscomb, [[Jesse Johnson]]<br />
| Software development and [[Adding a module to Glimmer I|creating a module for Glimmer]], [[representing and manipulating data]]. [[Grounding line treatments]], presented by Sophie Nowicki. [[Verifying ice sheet models]], presented by Aitbala Sargent<br />
| [[Ian Rutt]], [[Magnus Hagdorn]], Gethin Williams, Stephen Price, Bill Lipscomb, [[Jesse Johnson]]<br />
|-valign="top" style="background:AliceBlue"<br />
| 14 August<br />
| Coupling the Cryosphere to other Earth systems: [[Coupling the Cryosphere to other Earth systems, part I|part I]], [[Coupling the Cryosphere to other Earth systems, part II|part II]] ; [[Course Evaluation]]<br />
| Bill Lipscomb and [[Ian Rutt]]<br />
| Community Climate System Model (CCSM) Lab<br />
| Bill Lipscomb, [[Jesse Johnson]], Stephen Price and [[Ian Rutt]]<br />
|}<br />
<br />
====[[Typical Daily Schedule]]====<br />
<br />
===Resources===<br />
<br />
Additional student/instructor resources for the Summer School:<br />
* List of [[Computing Resources and Room Description]]<br />
* [[Cryptocards]]<br />
* Details of [[Eliot Glacier field trip]]<br />
* An outline [[Reading List]]<br />
* [[Notes]] from daily lectures<br />
* Portland [[dining and brewpub suggestions]]<br />
* [[age verification|PDX afterhours]]<br />
* [[ideas for Portland extracurricular activities]]<br />
* [[Newsletter Article]]<br />
<br />
===Application and Registration===<br />
''The window for receipt of student applications has closed. Thank you for your interest in the program. ''<br />
<br />
The registration fee for the course is US $350.<br />
<br />
===Funding Agencies===<br />
<br />
<br />
{|<br />
|-valign="top"<br />
|[[Image:iscu.jpg|300 px]]<br />
|[[Image:scar.jpg|150 px]]<br />
|-valign="top"<br />
|[[Image:wcrp.jpg|200 px]]<br />
|[[Image:nsf_logo.gif|300px]]<br />
|-valign="top"<br />
|[[Image:cresis.jpg|100 px]]<br />
|[[Image:cires.jpg|350 px]]<br />
|-valign="top"<br />
|[[Image:IASC_logo_07_RGB.jpg|100 px]]<br />
|}<br />
<br />
===Organizing Committee===<br />
Christina Hulbe, Jesse Johnson, Cornelis van der Veen</div>Adamchttp://websrv.cs.umt.edu/isis/index.php/Talk:Newsletter_ArticleTalk:Newsletter Article2009-09-10T00:17:54Z<p>Adamc: added disusion about modeling software</p>
<hr />
<div>"like COMSOL, GLIMMER, and even a clunky finite differences script we wrote ourselves - to solve problems."<br />
<br />
COMSOL is not a model in the same that GLIMMER is a model. The way it reads now it suggests that COMSOL is kind of glacier model.</div>Adamchttp://websrv.cs.umt.edu/isis/index.php/Newsletter_ArticleNewsletter Article2009-09-10T00:13:53Z<p>Adamc: /* Afterword */ corrected my own mistake</p>
<hr />
<div>Hey Kristin,<br />
<br />
Why don't you add something here?<br />
<br />
-love Adam<br />
<br />
<br />
<br />
==Afterword==<br />
<br />
Why is a full-stress model important? Which area should I focus on in the next step of my career? What's with having to declare everything in Fortran? <br />
I learned the answers to all of these questions this summer at Portland State University at the excellent summer school that Christina Hulbe, Jesse Johnson, and Kees van der Veen convened. Summer schools are important to the glaciology community because - well, let's face it, there just aren't enough of us to offer a satisfyingly diverse array of ice classes at each institution. This summer school taught us those basics, from force balance and finite differences on up, but it quickly <br />
accelerated into specialized topics like basal processes, uncertainty in models, and using models - like COMSOL, GLIMMER, and even a clunky finite differences script we wrote ourselves - to solve problems. <br />
<br />
Actually, not nearly every group's model was as clunky as mine, and some were downright elegant. This was another strength of the modeling school - the nineteen students and nineteen instructor-students came from a broad splay of math, engineering, geology, physics, and computer science backgrounds. Some people could code in eight languages but encountered the term "grounding line" for the first time at the school, while others were just the opposite. Since everyone had a different area of expertise, we were able to teach each other. I learned at least as much from my five-person group as I did from the lectures and derivations. <br />
<br />
The makeup of the groups was the modeling school's cleverest element. The school grouped three students and two instructor-students - early-career researchers or professors, often who had entered glaciology recently from another field - together to work on the daily assignments. I tend to think of "networking" as a stressful, prim-and-proper activity, but it turns out that the whole time I was installing GLIMMER, turning model knobs, brow-scratching, and laughing with my group partners, I was building professional relationships. We've got some really sharp, friendly people in glaciology, and I see a future of happy collaboration with the many of them that I met at the summer school.<br />
<br />
368 words!!!!!</div>Adamchttp://websrv.cs.umt.edu/isis/index.php/Newsletter_ArticleNewsletter Article2009-09-10T00:11:36Z<p>Adamc: /* Afterword */ changed 'modelling' to 'modeling'</p>
<hr />
<div>Hey Kristin,<br />
<br />
Why don't you add something here?<br />
<br />
-love Adam<br />
<br />
<br />
<br />
==Afterword==<br />
<br />
Why is a full-stress model important? Which area should I focus on in the next step of my career? What's with having to declare everything in Fortran? <br />
I learned the answers to all of these questions this summer at Portland State University at the excellent summer school that Christina Hulbe, Jesse Johnson and, Kees van der Veen convened. Summer schools are important to the glaciology community because - well, let's face it, there just aren't enough of us to offer a satisfyingly diverse array of ice classes at each institution. This summer school taught us those basics, from force balance and finite differences on up, but it quickly <br />
accelerated into specialized topics like basal processes, uncertainty in models, and using models - like COMSOL, GLIMMER, and even a clunky finite differences script we wrote ourselves - to solve problems. <br />
<br />
Actually, not nearly every group's model was as clunky as mine, and some were downright elegant. This was another strength of the modeling school - the nineteen students and nineteen instructor-students came from a broad splay of math, engineering, geology, physics, and computer science backgrounds. Some people could code in eight languages but encountered the term "grounding line" for the first time at the school, while others were just the opposite. Since everyone had a different area of expertise, we were able to teach each other. I learned at least as much from my five-person group as I did from the lectures and derivations. <br />
<br />
The makeup of the groups was the modeling school's cleverest element. The school grouped three students and two instructor-students - early-career researchers or professors, often who had entered glaciology recently from another field - together to work on the daily assignments. I tend to think of "networking" as a stressful, prim-and-proper activity, but it turns out that the whole time I was installing GLIMMER, turning model knobs, brow-scratching, and laughing with my group partners, I was building professional relationships. We've got some really sharp, friendly people in glaciology, and I see a future of happy collaboration with the many of them that I met at the summer school.<br />
<br />
368 words!!!!!</div>Adamchttp://websrv.cs.umt.edu/isis/index.php/Newsletter_ArticleNewsletter Article2009-09-10T00:11:00Z<p>Adamc: /* Afterword */ changed conjection</p>
<hr />
<div>Hey Kristin,<br />
<br />
Why don't you add something here?<br />
<br />
-love Adam<br />
<br />
<br />
<br />
==Afterword==<br />
<br />
Why is a full-stress model important? Which area should I focus on in the next step of my career? What's with having to declare everything in Fortran? <br />
I learned the answers to all of these questions this summer at Portland State University at the excellent summer school that Christina Hulbe, Jesse Johnson and, Kees van der Veen convened. Summer schools are important to the glaciology community because - well, let's face it, there just aren't enough of us to offer a satisfyingly diverse array of ice classes at each institution. This summer school taught us those basics, from force balance and finite differences on up, but it quickly <br />
accelerated into specialized topics like basal processes, uncertainty in models, and using models - like COMSOL, GLIMMER, and even a clunky finite differences script we wrote ourselves - to solve problems. <br />
<br />
Actually, not nearly every group's model was as clunky as mine, and some were downright elegant. This was another strength of the modelling school - the nineteen students and nineteen instructor-students came from a broad splay of math, engineering, geology, physics, and computer science backgrounds. Some people could code in eight languages but encountered the term "grounding line" for the first time at the school, while others were just the opposite. Since everyone had a different area of expertise, we were able to teach each other. I learned at least as much from my five-person group as I did from the lectures and derivations. <br />
<br />
The makeup of the groups was the modelling school's cleverest element. The school grouped three students and two instructor-students - early-career researchers or professors, often who had entered glaciology recently from another field - together to work on the daily assignments. I tend to think of "networking" as a stressful, prim-and-proper activity, but it turns out that the whole time I was installing GLIMMER, turning model knobs, brow-scratching, and laughing with my group partners, I was building professional relationships. We've got some really sharp, friendly people in glaciology, and I see a future of happy collaboration with the many of them that I met at the summer school.<br />
<br />
368 words!!!!!</div>Adamchttp://websrv.cs.umt.edu/isis/index.php/Newsletter_ArticleNewsletter Article2009-09-10T00:07:57Z<p>Adamc: /* Afterword */ changed location to Portland State University</p>
<hr />
<div>Hey Kristin,<br />
<br />
Why don't you add something here?<br />
<br />
-love Adam<br />
<br />
<br />
<br />
==Afterword==<br />
<br />
Why is a full-stress model important? Which area should I focus on in the next step of my career? What's with having to declare everything in Fortran? <br />
I learned the answers to all of these questions this summer at Portland State University at the excellent summer school that Christina Hulbe, Jesse Johnson and, Kees van der Veen convened. Summer schools are important to the glaciology community because - well, let's face it, there just aren't enough of us to offer a satisfyingly diverse array of ice classes at each institution. This summer school taught us those basics, from force balance and finite differences on up, but it quickly <br />
accelerated into specialized topics like basal processes, uncertainty in models, and using models - like COMSOL, GLIMMER, and even a clunky finite differences script we wrote ourselves - to solve problems. <br />
<br />
Actually, not nearly every group's model was as clunky as mine, and some were downright elegant. This was another strength of the modelling school - the nineteen students and nineteen instructor-students came from a broad splay of math, engineering, geology, physics, and computer science backgrounds. Some people could code in eight languages but encountered the term "grounding line" for the first time at the school, and others were just the opposite. Since everyone had a different area of expertise, we were able to teach each other. I learned at least as much from my five-person group as I did from the lectures and derivations. <br />
<br />
The makeup of the groups was the modelling school's cleverest element. The school grouped three students and two instructor-students - early-career researchers or professors, often who had entered glaciology recently from another field - together to work on the daily assignments. I tend to think of "networking" as a stressful, prim-and-proper activity, but it turns out that the whole time I was installing GLIMMER, turning model knobs, brow-scratching, and laughing with my group partners, I was building professional relationships. We've got some really sharp, friendly people in glaciology, and I see a future of happy collaboration with the many of them that I met at the summer school.<br />
<br />
368 words!!!!!</div>Adamchttp://websrv.cs.umt.edu/isis/index.php/Newsletter_ArticleNewsletter Article2009-09-10T00:07:11Z<p>Adamc: /* Afterword */ alphabetized name order</p>
<hr />
<div>Hey Kristin,<br />
<br />
Why don't you add something here?<br />
<br />
-love Adam<br />
<br />
<br />
<br />
==Afterword==<br />
<br />
Why is a full-stress model important? Which area should I focus on in the next step of my career? What's with having to declare everything in Fortran? <br />
I learned the answers to all of these questions this summer in Portland, Oregon, at the excellent summer school that Christina Hulbe, Jesse Johnson and, Kees van der Veen convened. Summer schools are important to the glaciology community because - well, let's face it, there just aren't enough of us to offer a satisfyingly diverse array of ice classes at each institution. This summer school taught us those basics, from force balance and finite differences on up, but it quickly <br />
accelerated into specialized topics like basal processes, uncertainty in models, and using models - like COMSOL, GLIMMER, and even a clunky finite differences script we wrote ourselves - to solve problems. <br />
<br />
Actually, not nearly every group's model was as clunky as mine, and some were downright elegant. This was another strength of the modelling school - the nineteen students and nineteen instructor-students came from a broad splay of math, engineering, geology, physics, and computer science backgrounds. Some people could code in eight languages but encountered the term "grounding line" for the first time at the school, and others were just the opposite. Since everyone had a different area of expertise, we were able to teach each other. I learned at least as much from my five-person group as I did from the lectures and derivations. <br />
<br />
The makeup of the groups was the modelling school's cleverest element. The school grouped three students and two instructor-students - early-career researchers or professors, often who had entered glaciology recently from another field - together to work on the daily assignments. I tend to think of "networking" as a stressful, prim-and-proper activity, but it turns out that the whole time I was installing GLIMMER, turning model knobs, brow-scratching, and laughing with my group partners, I was building professional relationships. We've got some really sharp, friendly people in glaciology, and I see a future of happy collaboration with the many of them that I met at the summer school.<br />
<br />
368 words!!!!!</div>Adamchttp://websrv.cs.umt.edu/isis/index.php/Newsletter_ArticleNewsletter Article2009-08-27T18:59:56Z<p>Adamc: Getting 'er started</p>
<hr />
<div>Hey Kristin,<br />
<br />
Why don't you add something here?<br />
<br />
-love Adam</div>Adamchttp://websrv.cs.umt.edu/isis/index.php/Summer_Modeling_SchoolSummer Modeling School2009-08-27T18:58:34Z<p>Adamc: Added Section for working on</p>
<hr />
<div>[[Image:Portland.jpg|thumb|right|400 px|The summer ice sheet modeling school will be held in Portland Oregon, August 3-14, 2009]]<br />
<br />
==Overview==<br />
The Summer Modeling School will be an intensive Summer School that will bring current and future ice-sheet scientists together to develop better models for the projection of future sea-level rise (slr). The IPCC Fourth Assessment Report [http://www.ipcc.ch/ipccreports/ar4-syr.htm] acknowledged that current models do not adequately treat the dynamic response of ice sheets to climate change, and that this is the largest uncertainty in assessing potential rapid sea-level rise. Recognizing this, an ice-sheet modelling Workshop was held during the July 2008 SCAR/IASC [https://www.comnap.aq/content/events/osc2008] meeting, in St. Petersburg, Russia. This meeting developed a community strategy on how best to (i) improve the physical understanding of ice-sheet processes responsible for rapid change; (ii) incorporate improved physical understanding into numerical models; (iii) assimilate appropriate data into the models for calibration and validation; and (iv) develop prognostic whole ice-sheet models that better incorporate non-linear ice-sheet response to environmental forcing (such as change in surface mass balance, loss of buttressing from floating ice shelves and ice tongues, and rising sea level). <br />
<br />
The two-week Summer School is a first step towards implementing this strategy. It will bring scientists from differing backgrounds together and allow more extensive and in-depth interactions between the relevant scientific research communities. A series of general background lectures as well as discussions of more specialized and advanced topics during this Summer School will provide the foundation for cross-disciplinary research, particularly for early career scientists. We anticipate publication of lecture notes both in hard copy and on a dedicated home page, to provide the glaciological community with an up-to-date overview of the science and observational techniques that will serve to guide further research efforts. Direct beneficiaries will be young researchers; indirect beneficiaries will be coastal zone communities who will gain improved sea level change forecasts to underpin their plans for sustainable development.<br />
<br />
===Venue===<br />
The modeling school was held on the campus of [[Wikipedia:Portland State University|Portland State University]] in [[Wikipedia:Portland, Oregon|Portland, Oregon]] August 3-14, 2009.<br />
<br />
* [http://maps.google.com/maps?f=d&source=s_d&saddr=Portland+Airport&daddr=310+SW+Lincoln+St,+Portland,+OR+97201-5007+(University+Place-Portland)&geocode=&hl=en&mra=ls&dirflg=r&date=07%2F28%2F09&time=8:59am&ttype=dep&noexp=0&noal=0&sort=&tline=&sll=45.54878,-122.629155&sspn=0.092445,0.144367&ie=UTF8&ll=45.548679,-122.619438&spn=0.092445,0.144367&z=13&start=0 Map] from airport to [http://cegs.pdx.edu/stay/upl/ University Place Hotel] using public transport (note that the directions in your travel letter are better than the Google generated instructions here).<br />
<br />
* [http://maps.google.com/maps?f=d&source=s_d&saddr=310+SW+Lincoln+St,+Portland,+OR+97201-5007+(University+Place-Portland)&daddr=1721+SW+Broadway,+Portland,+OR+97201+(Cramer+Hall)&hl=en&geocode=FdVhtgIdZwqw-CHO0mMQPCwi0Q%3BFRN3tgIdvP2v-CHxCBg32xEzXA&mra=ls&dirflg=w&sll=45.51029,-122.681675&sspn=0.005782,0.009023&ie=UTF8&ll=45.510091,-122.68232&spn=0.005782,0.009023&z=17 Map] from [http://cegs.pdx.edu/stay/upl/ University Place Hotel] to [http://www.pdx.edu/campus-map Cramer Hall].<br />
<br />
=== Student Participants ===<br />
<br />
*[[Nominations for Student/Instructor Awards]]<br />
*[[Student Bios]]<br />
*[[Student Presentation Development]]<br />
*[[Groups]] example of [[connections in groups]]<br />
*[[Terminology]]<br />
*[[Questions]]<br />
<br />
===Lectures and Planned Activities===<br />
<br />
For information about editing this page, see [[Wikipedia:How to edit]].<br />
<br />
{| border="1" cellpadding="5" cellspacing="0"<br />
|-valign="top" style="background:RoyalBlue"<br />
!width="20%"|Dates<br />
!width="25%"|Lecture Topics<br />
!width="15%"|Lecturers<br />
!width="25%"|Laboratory Topics<br />
!width="15%"|Laboratory Instructors <br />
|-valign="top" style="background:AliceBlue"<br />
| [[4-5 August]]<br />
| Introduction to and theoretical basis for ice sheet modeling. <br />
| Kees van der Veen, [[Nina Kirchner]] <br />
| [[Finite differencing|Finite differencing]] and [[Pragmatic Programming|pragmatic programming]] using Fortran[http://en.wikipedia.org/wiki/Fortran] 95...<br />
computing divergence and gradient...<br />
from conservation equation to matrix algebra...<br />
rheology and that which makes ice ice...<br />
simple, ideal models...<br />
that which makes ice-sheet modeling hard...<br />
| Gethin Williams, [[Ian Rutt]], [[Jesse Johnson]]<br />
|-valign="top" style="background:PowderBlue"<br />
| 6 August <br />
| [[Basal Conditions]], [[Data sets for ice sheet modeling]]<br />
| Alan Rempel, Slawek Tulaczyk and Ken Jezek<br />
| [[COMSOL Multiphysics]]<br />
| Olga Sergienko and Jesse Johnson<br />
|-valign="top" style="background:AliceBlue"<br />
| 7 August<br />
| The world of [[ice shelves]] and 'distributed stress-field solutions'. [[Modelling mountain glaciers]].<br />
| Todd Dupont, Olga Sergienko, and Brian Anderson<br />
| Linear Algebra of ice-sheet modeling, relaxation methods, finite-element methodology, solution of Laplace equation in arbitrary domain, creation of an ice-shelf flow-field model (snap shot of flow field), Models of the Ross Ice Shelf<br />
| Olga Sergienko and Todd Dupont<br />
<br />
|-valign="top" style="background:PowderBlue"<br />
| 8 August<br />
| [[Student Presentation]]<br />
| [[Modeling School Students]]<br />
| open work day with breakfast at 8 am & student presentation at 9 am<br />
| go to the farmer's market<br />
|-valign="top" style="background:AliceBlue"<br />
| 9 August<br />
| Free day; possible PDX tour<br />
|<br />
|<br />
|<br />
|-valign="top" style="background:PowderBlue"<br />
| 10 August<br />
| Excursion to Mt. Hood and [[Eliot Glacier field trip]] <br> Meet at 7:00 am for breakfast in Cramer Hall. Bring a jacket, water, and sensible (supportive) shoes.<br> Box lunches are provided and we will stop for dinner on the return trip (bring cash). <br />
| Guided by [http://web.pdx.edu/~basagic/ Hassan Basagic]<br />
|<br />
|<br />
|-valign="top" style="background:AliceBlue"<br />
| [[Quantifying_model_uncertainty |11 August]]<br />
| [[Quantifying model uncertainty]]<br />
| Charles Jackson and [http://www.mit.edu/~heimbach Patrick Heimbach]<br />
| Uncertain lab, [[Dynamic response to the enhanced basal flow in the Greenland ice sheet]] Weili Wang<br />
| Charles Jackson, Patrick Heimbach, and Weili Wang<br />
|-valign="top" style="background:PowderBlue"<br />
| [[12-13 August]]<br />
| Introduction to Glimmer-CISM ([[Introduction to Glimmer I|Part I]], [[Introduction to Glimmer II|Part II]]); [[Higher order velocity schemes|Higher-order models]]<br />
| [[Ian Rutt]], [[Magnus Hagdorn]], [[Stephen Price]], Bill Lipscomb, [[Jesse Johnson]]<br />
| Software development and [[Adding a module to Glimmer I|creating a module for Glimmer]], [[representing and manipulating data]]. [[Grounding line treatments]], presented by Sophie Nowicki. [[Verifying ice sheet models]], presented by Aitbala Sargent<br />
| [[Ian Rutt]], [[Magnus Hagdorn]], Gethin Williams, Stephen Price, Bill Lipscomb, [[Jesse Johnson]]<br />
|-valign="top" style="background:AliceBlue"<br />
| 14 August<br />
| Coupling the Cryosphere to other Earth systems: [[Coupling the Cryosphere to other Earth systems, part I|part I]], [[Coupling the Cryosphere to other Earth systems, part II|part II]] ; [[Course Evaluation]]<br />
| Bill Lipscomb and [[Ian Rutt]]<br />
| Community Climate System Model (CCSM) Lab<br />
| Bill Lipscomb, [[Jesse Johnson]], Stephen Price and [[Ian Rutt]]<br />
|}<br />
<br />
====[[Typical Daily Schedule]]====<br />
<br />
===Resources===<br />
<br />
Additional student/instructor resources for the Summer School:<br />
* List of [[Computing Resources and Room Description]]<br />
* [[Cryptocards]]<br />
* Details of [[Eliot Glacier field trip]]<br />
* An outline [[Reading List]]<br />
* [[Notes]] from daily lectures<br />
* Portland [[dining and brewpub suggestions]]<br />
* [[age verification|PDX afterhours]]<br />
* [[ideas for Portland extracurricular activities]]<br />
* [[Newsletter Article]]<br />
<br />
===Application and Registration===<br />
''The window for receipt of student applications has closed. Thank you for your interest in the program. ''<br />
<br />
The registration fee for the course is US $350.<br />
<br />
===Funding Agencies===<br />
<br />
<br />
{|<br />
|-valign="top"<br />
|[[Image:iscu.jpg|300 px]]<br />
|[[Image:scar.jpg|150 px]]<br />
|-valign="top"<br />
|[[Image:wcrp.jpg|200 px]]<br />
|[[Image:nsf_logo.gif|300px]]<br />
|-valign="top"<br />
|[[Image:cresis.jpg|100 px]]<br />
|[[Image:cires.jpg|350 px]]<br />
|-valign="top"<br />
|[[Image:IASC_logo_07_RGB.jpg|100 px]]<br />
|}<br />
<br />
===Organizing Committee===<br />
Christina Hulbe, Jesse Johnson, Cornelis van der Veen</div>Adamchttp://websrv.cs.umt.edu/isis/index.php/PDX_afterhoursPDX afterhours2009-08-15T18:55:07Z<p>Adamc: /* Thursday Aug 13, 2009 */ added the thing</p>
<hr />
<div>[[Image:Erin_adam.JPG|thumb|right|300px|The PDX Afterhours king & queen in Tube Bar, home of Wednesday night $1 Miller High Life.]]<br />
<br />
==tuesday Aug 4, 2009==<br />
* meet at [http://paccinirestaurant.com/ Paccini's] pub at 7:30<br />
<br />
==wednesday Aug 5, 2009==<br />
* after FORTRAN session, leave PSU to go to [http://maps.google.com/maps?hl=en&client=firefox-a&rls=com.ubuntu:en-US:unofficial&hs=zuU&um=1&ie=UTF-8&q=sushi+ichiban+portland&fb=1&split=1&gl=us&view=text&latlng=899833397715679289 Sushi Ichiban]. Adam Campbell will guide you.<br />
* go to [http://www.groundkontrol.com Ground Kontrol] a retro arcade with beer<br />
* go to [http://www.voodoodoughnut.com Voodoo Doughnut], please someone buy Ian Rutt the $5 doughnut.<br />
<br />
<br />
==thursday Aug 6, 2009==<br />
[http://amontobin.com/field/ Amon Tobin] and [http://www.pitchblack.co.nz/?s1=index Pitch Black], along with two opening bands, are playing at the Roseland Theatre (8 NW 6th Ave) Thursday night (starting at 9:00 pm or thereabouts). Tickets are $26 (available online at [http://ticketswest.rdln.com/Venue.aspx?ven=ROS TicketsWest]). The music is best described as sampled electronica (Amon Tobin) and Kiwi-style dub (Pitch Black). I'd expect a late night of electronic music: an evening nap may be in order! Jeremy's already got his ticket and can fill you in with more, including a sample of the music.<br />
<br />
*other local music recommendations from Adam<br />
<br />
'''Boy Eats Drum Machine, French Miami, Southern Belle and Electric Opera Company''' - Indie Rock -<br />
Thu., Aug. 6, 9 p.m.<br />
$6-8<br />
Berbati's Pan<br />
10 SW 3rd Ave.<br />
Downtown<br />
<br />
'''Nurses, Inside Voices and Slaves''' - Indie Rock -<br />
Thu., Aug. 6, 8:30 p.m.<br />
$7<br />
Holocene <br />
1001 SE Morrison<br />
Southeast<br />
<br />
==friday Aug 7, 2009==<br />
[http://www.biteoforegon.com/ The Bite of Oregon] is a food festival that takes place at Tom McCall Waterfront Park. Featuring food, wine, beer and entertainment from Oregon. Entry is $8, food and beverages are extra.<br />
<br />
[http://www.pioneercourthousesquare.org/calendar_august.htm Flicks on the Bricks] will be showing Jurassic Park at dusk outside at Pioneer Courthouse Square, FREE (including popcorn). (10 minute walk)<br />
<br />
[http://www.portlandtwilight.com/ Portland Downtown Twilight Bike Criterium]: Professional (by U.S. standards) bike race through downtown Portland. Complete with beer garden, food, and an expected 15,000 spectators. Pro race starts at 7:30.<br />
<br />
==saturday Aug 8, 2009==<br />
[http://www.biteoforegon.com/ The Bite of Oregon] is a food festival that takes place at Tom McCall Waterfront Park. Featuring food, wine, beer and entertainment from Oregon. Entry is $8, food and beverages are extra.<br />
<br />
==sunday Aug 9, 2009==<br />
The [http://providence.org/bridgepedal/ Portland Bridge Pedal] is a fun event where you can go on 14, 24, or 37 mi. bike ride over Portland's Bridges. Adam is trying to assemble a group to go on Sunday morning. Please speak with him if you are interested in going to this. I tentatively have 4 bikes I can get ahold of.<br />
<br />
RIDERS: Register [http://providence.org/bridgepedal/ here]. Make sure to click the team rider option and the 11 bridge option! Ok so here is the deal I tried to register a team:<br />
Team Name The Icepocalypse<br />
Team password BPT943<br />
<br />
HEY It works now!<br />
<br />
[http://www.biteoforegon.com/ The Bite of Oregon] is a food festival that takes place at Tom McCall Waterfront Park. Featuring food, wine, beer and entertainment from Oregon. Entry is $8, food and beverages are extra.<br />
<br />
==tuesday Aug 11, 2009==<br />
[http://www.rootsorganicbrewing.com/ Roots Brewery] has a $2.50 imperial pint night. The beer is great! It's on the eastside but it's not very far. There are carts and more beer nearby.<br />
<br />
==wednesday Aug 12, 2009==<br />
[[Image:420.png|thumb|right|300px|The Glimmer model output gives suggestions for illicit afterhours activities in Portland. I suggest looking in the North park blocks.]]<br />
<br />
Quiet night. Kristin ate a Whole Foods Meal Pod in her room and Toby went for pho by himself. Good work!<br />
<br />
==Thursday Aug 13, 2009==<br />
We attempted to make an expedition to the movie theater last Tuesday but I forgot my ID (:D). Tonight, we will make another try... So let's meet at the hotel at 6:30 pm. Don't forget to bring your ID!!!!!<br />
<br />
[http://www.freedomridersthemovie.com/ Freedom Riders] is a mountain bike documentary about the evolution of freeriding near Jackson Hole, WY (I think). The show is at the [http://www.clintonsttheater.com/ Clinton St. Theater]. It starts at 7 pm.<br />
<br />
Summary of the evening:<br />
* Star Trek discussion / planning of future paper about the Red Black Hole<br />
* Three people from Seattle get touch-and-go jello shots with one very cool Portlander<br />
* Kristin tries to order Munster cheese at the Portland Timbers bar HQ. The retort: "This is a sportsbar, hon, we have cheddar or swiss."<br />
* Gethin is a pool shark, but it's Doug who holds his own against the expert-shooting hipster who usurped our table<br />
* None of us is a shuffleboard shark<br />
* Todd has the best time of all - measured with <math>N_{drinks}</math> data and a degree-of-swaying interpolation scheme - and has not showed up at the workshop this morning<br />
<br />
==Friday and Saturday==<br />
The Thing at [http://5thavenuecinema.org/ Fifth Avenue Cinema]. Showtimes: FRI. / SAT. - 7:00 P.M. & 9:30 P.M SUN. - 3:00 P.M.<br />
<br />
[[How_to_make_a_model|50c beer all night long]]<br />
<br />
<br />
==For the end...==<br />
<br />
CEILING POPCORN!!!!! Did you check your bedroom ceiling?<br />
Do you want to take a bubble bath?<br />
<br />
+++ We are doing a survey about the ceiling popcorn and the bubble bath.+++ <br />
<br />
Ask Kristin for more details.<br />
<br />
<br />
If you want the recipe of the Mushroon Linguini ask Flo for more details...<br />
<br />
Ciao! Ciao!<br />
<br />
SOOOOO SSAAAAADDDDD!!!! I invite you all in italy in one year to make an after school!!!!!<br />
<3 <3 <3 <3 <3</div>Adamchttp://websrv.cs.umt.edu/isis/index.php/Age_verificationAge verification2009-08-15T18:44:52Z<p>Adamc: added verification</p>
<hr />
<div>WARNING: You must be 21 years of age or greater to enter [[PDX afterhours]]</div>Adamchttp://websrv.cs.umt.edu/isis/index.php/Summer_Modeling_SchoolSummer Modeling School2009-08-15T18:43:52Z<p>Adamc: /* Resources */ added link to age verification</p>
<hr />
<div>[[Image:Portland.jpg|thumb|right|400 px|The summer ice sheet modeling school will be held in Portland Oregon, August 3-14, 2009]]<br />
<br />
==Overview==<br />
The Summer Modeling School will be an intensive Summer School that will bring current and future ice-sheet scientists together to develop better models for the projection of future sea-level rise (slr). The IPCC Fourth Assessment Report [http://www.ipcc.ch/ipccreports/ar4-syr.htm] acknowledged that current models do not adequately treat the dynamic response of ice sheets to climate change, and that this is the largest uncertainty in assessing potential rapid sea-level rise. Recognizing this, an ice-sheet modelling Workshop was held during the July 2008 SCAR/IASC [https://www.comnap.aq/content/events/osc2008] meeting, in St. Petersburg, Russia. This meeting developed a community strategy on how best to (i) improve the physical understanding of ice-sheet processes responsible for rapid change; (ii) incorporate improved physical understanding into numerical models; (iii) assimilate appropriate data into the models for calibration and validation; and (iv) develop prognostic whole ice-sheet models that better incorporate non-linear ice-sheet response to environmental forcing (such as change in surface mass balance, loss of buttressing from floating ice shelves and ice tongues, and rising sea level). <br />
<br />
The two-week Summer School is a first step towards implementing this strategy. It will bring scientists from differing backgrounds together and allow more extensive and in-depth interactions between the relevant scientific research communities. A series of general background lectures as well as discussions of more specialized and advanced topics during this Summer School will provide the foundation for cross-disciplinary research, particularly for early career scientists. We anticipate publication of lecture notes both in hard copy and on a dedicated home page, to provide the glaciological community with an up-to-date overview of the science and observational techniques that will serve to guide further research efforts. Direct beneficiaries will be young researchers; indirect beneficiaries will be coastal zone communities who will gain improved sea level change forecasts to underpin their plans for sustainable development.<br />
<br />
===Venue===<br />
The modeling school was held on the campus of [[Wikipedia:Portland State University|Portland State University]] in [[Wikipedia:Portland, Oregon|Portland, Oregon]] August 3-14, 2009.<br />
<br />
* [http://maps.google.com/maps?f=d&source=s_d&saddr=Portland+Airport&daddr=310+SW+Lincoln+St,+Portland,+OR+97201-5007+(University+Place-Portland)&geocode=&hl=en&mra=ls&dirflg=r&date=07%2F28%2F09&time=8:59am&ttype=dep&noexp=0&noal=0&sort=&tline=&sll=45.54878,-122.629155&sspn=0.092445,0.144367&ie=UTF8&ll=45.548679,-122.619438&spn=0.092445,0.144367&z=13&start=0 Map] from airport to [http://cegs.pdx.edu/stay/upl/ University Place Hotel] using public transport (note that the directions in your travel letter are better than the Google generated instructions here).<br />
<br />
* [http://maps.google.com/maps?f=d&source=s_d&saddr=310+SW+Lincoln+St,+Portland,+OR+97201-5007+(University+Place-Portland)&daddr=1721+SW+Broadway,+Portland,+OR+97201+(Cramer+Hall)&hl=en&geocode=FdVhtgIdZwqw-CHO0mMQPCwi0Q%3BFRN3tgIdvP2v-CHxCBg32xEzXA&mra=ls&dirflg=w&sll=45.51029,-122.681675&sspn=0.005782,0.009023&ie=UTF8&ll=45.510091,-122.68232&spn=0.005782,0.009023&z=17 Map] from [http://cegs.pdx.edu/stay/upl/ University Place Hotel] to [http://www.pdx.edu/campus-map Cramer Hall].<br />
<br />
=== Student Participants ===<br />
<br />
*[[Nominations for Student/Instructor Awards]]<br />
*[[Student Bios]]<br />
*[[Student Presentation Development]]<br />
*[[Groups]] example of [[connections in groups]]<br />
*[[Terminology]]<br />
*[[Questions]]<br />
<br />
===Lectures and Planned Activities===<br />
<br />
For information about editing this page, see [[Wikipedia:How to edit]].<br />
<br />
{| border="1" cellpadding="5" cellspacing="0"<br />
|-valign="top" style="background:RoyalBlue"<br />
!width="20%"|Dates<br />
!width="25%"|Lecture Topics<br />
!width="15%"|Lecturers<br />
!width="25%"|Laboratory Topics<br />
!width="15%"|Laboratory Instructors <br />
|-valign="top" style="background:AliceBlue"<br />
| [[4-5 August]]<br />
| Introduction to and theoretical basis for ice sheet modeling. <br />
| Kees van der Veen, [[Nina Kirchner]] <br />
| [[Finite differencing|Finite differencing]] and [[Pragmatic Programming|pragmatic programming]] using Fortran[http://en.wikipedia.org/wiki/Fortran] 95...<br />
computing divergence and gradient...<br />
from conservation equation to matrix algebra...<br />
rheology and that which makes ice ice...<br />
simple, ideal models...<br />
that which makes ice-sheet modeling hard...<br />
| Gethin Williams, [[Ian Rutt]], [[Jesse Johnson]]<br />
|-valign="top" style="background:PowderBlue"<br />
| 6 August <br />
| [[Basal Conditions]], [[Data sets for ice sheet modeling]]<br />
| Alan Rempel, Slawek Tulaczyk and Ken Jezek<br />
| [[COMSOL Multiphysics]]<br />
| Olga Sergienko and Jesse Johnson<br />
|-valign="top" style="background:AliceBlue"<br />
| 7 August<br />
| The world of [[ice shelves]] and 'distributed stress-field solutions'. [[Modelling mountain glaciers]].<br />
| Todd Dupont, Olga Sergienko, and Brian Anderson<br />
| Linear Algebra of ice-sheet modeling, relaxation methods, finite-element methodology, solution of Laplace equation in arbitrary domain, creation of an ice-shelf flow-field model (snap shot of flow field), Models of the Ross Ice Shelf<br />
| Olga Sergienko and Todd Dupont<br />
<br />
|-valign="top" style="background:PowderBlue"<br />
| 8 August<br />
| [[Student Presentation]]<br />
| [[Modeling School Students]]<br />
| open work day with breakfast at 8 am & student presentation at 9 am<br />
| go to the farmer's market<br />
|-valign="top" style="background:AliceBlue"<br />
| 9 August<br />
| Free day; possible PDX tour<br />
|<br />
|<br />
|<br />
|-valign="top" style="background:PowderBlue"<br />
| 10 August<br />
| Excursion to Mt. Hood and [[Eliot Glacier field trip]] <br> Meet at 7:00 am for breakfast in Cramer Hall. Bring a jacket, water, and sensible (supportive) shoes.<br> Box lunches are provided and we will stop for dinner on the return trip (bring cash). <br />
| Guided by [http://web.pdx.edu/~basagic/ Hassan Basagic]<br />
|<br />
|<br />
|-valign="top" style="background:AliceBlue"<br />
| [[Quantifying_model_uncertainty |11 August]]<br />
| [[Quantifying model uncertainty]]<br />
| Charles Jackson and [http://www.mit.edu/~heimbach Patrick Heimbach]<br />
| Uncertain lab, [[Dynamic response to the enhanced basal flow in the Greenland ice sheet]] Weili Wang<br />
| Charles Jackson, Patrick Heimbach, and Weili Wang<br />
|-valign="top" style="background:PowderBlue"<br />
| [[12-13 August]]<br />
| Introduction to Glimmer-CISM ([[Introduction to Glimmer I|Part I]], [[Introduction to Glimmer II|Part II]]); [[Higher order velocity schemes|Higher-order models]]<br />
| [[Ian Rutt]], [[Magnus Hagdorn]], [[Stephen Price]], Bill Lipscomb, [[Jesse Johnson]]<br />
| Software development and [[Adding a module to Glimmer I|creating a module for Glimmer]], [[representing and manipulating data]]. [[Grounding line treatments]], presented by Sophie Nowicki. [[Verifying ice sheet models]], presented by Aitbala Sargent<br />
| [[Ian Rutt]], [[Magnus Hagdorn]], Gethin Williams, Stephen Price, Bill Lipscomb, [[Jesse Johnson]]<br />
|-valign="top" style="background:AliceBlue"<br />
| 14 August<br />
| Coupling the Cryosphere to other Earth systems: [[Coupling the Cryosphere to other Earth systems, part I|part I]], [[Coupling the Cryosphere to other Earth systems, part II|part II]] ; [[Course Evaluation]]<br />
| Bill Lipscomb and [[Ian Rutt]]<br />
| Community Climate System Model (CCSM) Lab<br />
| Bill Lipscomb, [[Jesse Johnson]], Stephen Price and [[Ian Rutt]]<br />
|}<br />
<br />
====[[Typical Daily Schedule]]====<br />
<br />
===Resources===<br />
<br />
Additional student/instructor resources for the Summer School:<br />
* List of [[Computing Resources and Room Description]]<br />
* [[Cryptocards]]<br />
* Details of [[Eliot Glacier field trip]]<br />
* An outline [[Reading List]]<br />
* [[Notes]] from daily lectures<br />
* Portland [[dining and brewpub suggestions]]<br />
* [[age verification|PDX afterhours]]<br />
* [[ideas for Portland extracurricular activities]]<br />
<br />
===Application and Registration===<br />
''The window for receipt of student applications has closed. Thank you for your interest in the program. ''<br />
<br />
The registration fee for the course is US $350.<br />
<br />
===Funding Agencies===<br />
<br />
<br />
{|<br />
|-valign="top"<br />
|[[Image:iscu.jpg|300 px]]<br />
|[[Image:scar.jpg|150 px]]<br />
|-valign="top"<br />
|[[Image:wcrp.jpg|200 px]]<br />
|[[Image:nsf_logo.gif|300px]]<br />
|-valign="top"<br />
|[[Image:cresis.jpg|100 px]]<br />
|[[Image:cires.jpg|350 px]]<br />
|-valign="top"<br />
|[[Image:IASC_logo_07_RGB.jpg|100 px]]<br />
|}<br />
<br />
===Organizing Committee===<br />
Christina Hulbe, Jesse Johnson, Cornelis van der Veen</div>Adamchttp://websrv.cs.umt.edu/isis/index.php/Summer_Modeling_SchoolSummer Modeling School2009-08-15T18:39:56Z<p>Adamc: /* Venue */ corrected verb tense</p>
<hr />
<div>[[Image:Portland.jpg|thumb|right|400 px|The summer ice sheet modeling school will be held in Portland Oregon, August 3-14, 2009]]<br />
<br />
==Overview==<br />
The Summer Modeling School will be an intensive Summer School that will bring current and future ice-sheet scientists together to develop better models for the projection of future sea-level rise (slr). The IPCC Fourth Assessment Report [http://www.ipcc.ch/ipccreports/ar4-syr.htm] acknowledged that current models do not adequately treat the dynamic response of ice sheets to climate change, and that this is the largest uncertainty in assessing potential rapid sea-level rise. Recognizing this, an ice-sheet modelling Workshop was held during the July 2008 SCAR/IASC [https://www.comnap.aq/content/events/osc2008] meeting, in St. Petersburg, Russia. This meeting developed a community strategy on how best to (i) improve the physical understanding of ice-sheet processes responsible for rapid change; (ii) incorporate improved physical understanding into numerical models; (iii) assimilate appropriate data into the models for calibration and validation; and (iv) develop prognostic whole ice-sheet models that better incorporate non-linear ice-sheet response to environmental forcing (such as change in surface mass balance, loss of buttressing from floating ice shelves and ice tongues, and rising sea level). <br />
<br />
The two-week Summer School is a first step towards implementing this strategy. It will bring scientists from differing backgrounds together and allow more extensive and in-depth interactions between the relevant scientific research communities. A series of general background lectures as well as discussions of more specialized and advanced topics during this Summer School will provide the foundation for cross-disciplinary research, particularly for early career scientists. We anticipate publication of lecture notes both in hard copy and on a dedicated home page, to provide the glaciological community with an up-to-date overview of the science and observational techniques that will serve to guide further research efforts. Direct beneficiaries will be young researchers; indirect beneficiaries will be coastal zone communities who will gain improved sea level change forecasts to underpin their plans for sustainable development.<br />
<br />
===Venue===<br />
The modeling school was held on the campus of [[Wikipedia:Portland State University|Portland State University]] in [[Wikipedia:Portland, Oregon|Portland, Oregon]] August 3-14, 2009.<br />
<br />
* [http://maps.google.com/maps?f=d&source=s_d&saddr=Portland+Airport&daddr=310+SW+Lincoln+St,+Portland,+OR+97201-5007+(University+Place-Portland)&geocode=&hl=en&mra=ls&dirflg=r&date=07%2F28%2F09&time=8:59am&ttype=dep&noexp=0&noal=0&sort=&tline=&sll=45.54878,-122.629155&sspn=0.092445,0.144367&ie=UTF8&ll=45.548679,-122.619438&spn=0.092445,0.144367&z=13&start=0 Map] from airport to [http://cegs.pdx.edu/stay/upl/ University Place Hotel] using public transport (note that the directions in your travel letter are better than the Google generated instructions here).<br />
<br />
* [http://maps.google.com/maps?f=d&source=s_d&saddr=310+SW+Lincoln+St,+Portland,+OR+97201-5007+(University+Place-Portland)&daddr=1721+SW+Broadway,+Portland,+OR+97201+(Cramer+Hall)&hl=en&geocode=FdVhtgIdZwqw-CHO0mMQPCwi0Q%3BFRN3tgIdvP2v-CHxCBg32xEzXA&mra=ls&dirflg=w&sll=45.51029,-122.681675&sspn=0.005782,0.009023&ie=UTF8&ll=45.510091,-122.68232&spn=0.005782,0.009023&z=17 Map] from [http://cegs.pdx.edu/stay/upl/ University Place Hotel] to [http://www.pdx.edu/campus-map Cramer Hall].<br />
<br />
=== Student Participants ===<br />
<br />
*[[Nominations for Student/Instructor Awards]]<br />
*[[Student Bios]]<br />
*[[Student Presentation Development]]<br />
*[[Groups]] example of [[connections in groups]]<br />
*[[Terminology]]<br />
*[[Questions]]<br />
<br />
===Lectures and Planned Activities===<br />
<br />
For information about editing this page, see [[Wikipedia:How to edit]].<br />
<br />
{| border="1" cellpadding="5" cellspacing="0"<br />
|-valign="top" style="background:RoyalBlue"<br />
!width="20%"|Dates<br />
!width="25%"|Lecture Topics<br />
!width="15%"|Lecturers<br />
!width="25%"|Laboratory Topics<br />
!width="15%"|Laboratory Instructors <br />
|-valign="top" style="background:AliceBlue"<br />
| [[4-5 August]]<br />
| Introduction to and theoretical basis for ice sheet modeling. <br />
| Kees van der Veen, [[Nina Kirchner]] <br />
| [[Finite differencing|Finite differencing]] and [[Pragmatic Programming|pragmatic programming]] using Fortran[http://en.wikipedia.org/wiki/Fortran] 95...<br />
computing divergence and gradient...<br />
from conservation equation to matrix algebra...<br />
rheology and that which makes ice ice...<br />
simple, ideal models...<br />
that which makes ice-sheet modeling hard...<br />
| Gethin Williams, [[Ian Rutt]], [[Jesse Johnson]]<br />
|-valign="top" style="background:PowderBlue"<br />
| 6 August <br />
| [[Basal Conditions]], [[Data sets for ice sheet modeling]]<br />
| Alan Rempel, Slawek Tulaczyk and Ken Jezek<br />
| [[COMSOL Multiphysics]]<br />
| Olga Sergienko and Jesse Johnson<br />
|-valign="top" style="background:AliceBlue"<br />
| 7 August<br />
| The world of [[ice shelves]] and 'distributed stress-field solutions'. [[Modelling mountain glaciers]].<br />
| Todd Dupont, Olga Sergienko, and Brian Anderson<br />
| Linear Algebra of ice-sheet modeling, relaxation methods, finite-element methodology, solution of Laplace equation in arbitrary domain, creation of an ice-shelf flow-field model (snap shot of flow field), Models of the Ross Ice Shelf<br />
| Olga Sergienko and Todd Dupont<br />
<br />
|-valign="top" style="background:PowderBlue"<br />
| 8 August<br />
| [[Student Presentation]]<br />
| [[Modeling School Students]]<br />
| open work day with breakfast at 8 am & student presentation at 9 am<br />
| go to the farmer's market<br />
|-valign="top" style="background:AliceBlue"<br />
| 9 August<br />
| Free day; possible PDX tour<br />
|<br />
|<br />
|<br />
|-valign="top" style="background:PowderBlue"<br />
| 10 August<br />
| Excursion to Mt. Hood and [[Eliot Glacier field trip]] <br> Meet at 7:00 am for breakfast in Cramer Hall. Bring a jacket, water, and sensible (supportive) shoes.<br> Box lunches are provided and we will stop for dinner on the return trip (bring cash). <br />
| Guided by [http://web.pdx.edu/~basagic/ Hassan Basagic]<br />
|<br />
|<br />
|-valign="top" style="background:AliceBlue"<br />
| [[Quantifying_model_uncertainty |11 August]]<br />
| [[Quantifying model uncertainty]]<br />
| Charles Jackson and [http://www.mit.edu/~heimbach Patrick Heimbach]<br />
| Uncertain lab, [[Dynamic response to the enhanced basal flow in the Greenland ice sheet]] Weili Wang<br />
| Charles Jackson, Patrick Heimbach, and Weili Wang<br />
|-valign="top" style="background:PowderBlue"<br />
| [[12-13 August]]<br />
| Introduction to Glimmer-CISM ([[Introduction to Glimmer I|Part I]], [[Introduction to Glimmer II|Part II]]); [[Higher order velocity schemes|Higher-order models]]<br />
| [[Ian Rutt]], [[Magnus Hagdorn]], [[Stephen Price]], Bill Lipscomb, [[Jesse Johnson]]<br />
| Software development and [[Adding a module to Glimmer I|creating a module for Glimmer]], [[representing and manipulating data]]. [[Grounding line treatments]], presented by Sophie Nowicki. [[Verifying ice sheet models]], presented by Aitbala Sargent<br />
| [[Ian Rutt]], [[Magnus Hagdorn]], Gethin Williams, Stephen Price, Bill Lipscomb, [[Jesse Johnson]]<br />
|-valign="top" style="background:AliceBlue"<br />
| 14 August<br />
| Coupling the Cryosphere to other Earth systems: [[Coupling the Cryosphere to other Earth systems, part I|part I]], [[Coupling the Cryosphere to other Earth systems, part II|part II]] ; [[Course Evaluation]]<br />
| Bill Lipscomb and [[Ian Rutt]]<br />
| Community Climate System Model (CCSM) Lab<br />
| Bill Lipscomb, [[Jesse Johnson]], Stephen Price and [[Ian Rutt]]<br />
|}<br />
<br />
====[[Typical Daily Schedule]]====<br />
<br />
===Resources===<br />
<br />
Additional student/instructor resources for the Summer School:<br />
* List of [[Computing Resources and Room Description]]<br />
* [[Cryptocards]]<br />
* Details of [[Eliot Glacier field trip]]<br />
* An outline [[Reading List]]<br />
* [[Notes]] from daily lectures<br />
* Portland [[dining and brewpub suggestions]]<br />
* [[PDX afterhours]]<br />
* [[ideas for Portland extracurricular activities]]<br />
<br />
===Application and Registration===<br />
''The window for receipt of student applications has closed. Thank you for your interest in the program. ''<br />
<br />
The registration fee for the course is US $350.<br />
<br />
===Funding Agencies===<br />
<br />
<br />
{|<br />
|-valign="top"<br />
|[[Image:iscu.jpg|300 px]]<br />
|[[Image:scar.jpg|150 px]]<br />
|-valign="top"<br />
|[[Image:wcrp.jpg|200 px]]<br />
|[[Image:nsf_logo.gif|300px]]<br />
|-valign="top"<br />
|[[Image:cresis.jpg|100 px]]<br />
|[[Image:cires.jpg|350 px]]<br />
|-valign="top"<br />
|[[Image:IASC_logo_07_RGB.jpg|100 px]]<br />
|}<br />
<br />
===Organizing Committee===<br />
Christina Hulbe, Jesse Johnson, Cornelis van der Veen</div>Adamchttp://websrv.cs.umt.edu/isis/index.php/Cryptocard_holdersCryptocard holders2009-08-14T22:48:05Z<p>Adamc: Add my contact info</p>
<hr />
<div>Jesse Johnson <br />
<br />
*417 Social Science Building, 32 Ave Missoula MT, 59812, jesse.v.johnson@cs.umt.edu<br />
<br />
<br />
Florence Colleoni<br />
<br />
*Centro Euro-Mediterraneo per i Cambiamenti Climatici, Via Aldo Moro 44, 40127 Bologna, Italy, flocolleoni@gmail.com<br />
<br />
<br />
Kristin Poinar<br />
* University of Washington Dept. of Earth and Space Sciences, Box 351310, 4000 15th Ave. NE, Seattle, WA 98125. kpoinar@u.washington.edu<br />
<br />
<br />
Gethin Williams<br />
<br />
* School of Geographical Sciences, University of Bristol, University Road, Bristol BS8 1SS. United Kingdom. gethin.williams@bristol.ac.uk<br />
<br />
<br />
Saffia Hossainzadeh<br />
* University of California, Santa Cruz, Dept. of Earth and Planetary Sciences, 1156 High St., Santa Cruz, CA 95064, hoss@uchicago.edu<br />
<br />
<br />
Adam Campbell<br />
* University of Washington Dept. of Earth and Space Sciences, Box 351310, 4000 15th Ave. NE, Seattle, WA 98125. campbead@u.washington.edu</div>Adamchttp://websrv.cs.umt.edu/isis/index.php/PDX_afterhoursPDX afterhours2009-08-12T18:35:34Z<p>Adamc: added glimmer output</p>
<hr />
<div>[[Image:Erin_adam.JPG|thumb|right|300px|The PDX Afterhours king & queen in Tube Bar, home of Wednesday night $1 Miller High Life.]]<br />
<br />
==tuesday Aug 4, 2009==<br />
* meet at [http://paccinirestaurant.com/ Paccini's] pub at 7:30<br />
<br />
==wednesday Aug 5, 2009==<br />
* after FORTRAN session, leave PSU to go to [http://maps.google.com/maps?hl=en&client=firefox-a&rls=com.ubuntu:en-US:unofficial&hs=zuU&um=1&ie=UTF-8&q=sushi+ichiban+portland&fb=1&split=1&gl=us&view=text&latlng=899833397715679289 Sushi Ichiban]. Adam Campbell will guide you.<br />
* go to [http://www.groundkontrol.com Ground Kontrol] a retro arcade with beer<br />
* go to [http://www.voodoodoughnut.com Voodoo Doughnut], please someone buy Ian Rutt the $5 doughnut.<br />
<br />
<br />
==thursday Aug 6, 2009==<br />
[http://amontobin.com/field/ Amon Tobin] and [http://www.pitchblack.co.nz/?s1=index Pitch Black], along with two opening bands, are playing at the Roseland Theatre (8 NW 6th Ave) Thursday night (starting at 9:00 pm or thereabouts). Tickets are $26 (available online at [http://ticketswest.rdln.com/Venue.aspx?ven=ROS TicketsWest]). The music is best described as sampled electronica (Amon Tobin) and Kiwi-style dub (Pitch Black). I'd expect a late night of electronic music: an evening nap may be in order! Jeremy's already got his ticket and can fill you in with more, including a sample of the music.<br />
<br />
*other local music recommendations from Adam<br />
<br />
'''Boy Eats Drum Machine, French Miami, Southern Belle and Electric Opera Company''' - Indie Rock -<br />
Thu., Aug. 6, 9 p.m.<br />
$6-8<br />
Berbati's Pan<br />
10 SW 3rd Ave.<br />
Downtown<br />
<br />
'''Nurses, Inside Voices and Slaves''' - Indie Rock -<br />
Thu., Aug. 6, 8:30 p.m.<br />
$7<br />
Holocene <br />
1001 SE Morrison<br />
Southeast<br />
<br />
==friday Aug 7, 2009==<br />
[http://www.biteoforegon.com/ The Bite of Oregon] is a food festival that takes place at Tom McCall Waterfront Park. Featuring food, wine, beer and entertainment from Oregon. Entry is $8, food and beverages are extra.<br />
<br />
[http://www.pioneercourthousesquare.org/calendar_august.htm Flicks on the Bricks] will be showing Jurassic Park at dusk outside at Pioneer Courthouse Square, FREE (including popcorn). (10 minute walk)<br />
<br />
[http://www.portlandtwilight.com/ Portland Downtown Twilight Bike Criterium]: Professional (by U.S. standards) bike race through downtown Portland. Complete with beer garden, food, and an expected 15,000 spectators. Pro race starts at 7:30.<br />
<br />
==saturday Aug 8, 2009==<br />
[http://www.biteoforegon.com/ The Bite of Oregon] is a food festival that takes place at Tom McCall Waterfront Park. Featuring food, wine, beer and entertainment from Oregon. Entry is $8, food and beverages are extra.<br />
<br />
==sunday Aug 9, 2009==<br />
The [http://providence.org/bridgepedal/ Portland Bridge Pedal] is a fun event where you can go on 14, 24, or 37 mi. bike ride over Portland's Bridges. Adam is trying to assemble a group to go on Sunday morning. Please speak with him if you are interested in going to this. I tentatively have 4 bikes I can get ahold of.<br />
<br />
RIDERS: Register [http://providence.org/bridgepedal/ here]. Make sure to click the team rider option and the 11 bridge option! Ok so here is the deal I tried to register a team:<br />
Team Name The Icepocalypse<br />
Team password BPT943<br />
<br />
HEY It works now!<br />
<br />
[http://www.biteoforegon.com/ The Bite of Oregon] is a food festival that takes place at Tom McCall Waterfront Park. Featuring food, wine, beer and entertainment from Oregon. Entry is $8, food and beverages are extra.<br />
<br />
==tuesday Aug 11, 2009==<br />
[http://www.rootsorganicbrewing.com/ Roots Brewery] has a $2.50 imperial pint night. The beer is great! It's on the eastside but it's not very far. There are carts and more beer nearby.<br />
<br />
==wednesday Aug 12, 2009==<br />
[[Image:420.png|thumb|right|300px|The Glimmer model output gives suggestions for illicit afterhours activities in Portland. I suggest looking in the North park blocks.]]</div>Adamchttp://websrv.cs.umt.edu/isis/index.php/File:420.pngFile:420.png2009-08-12T18:33:02Z<p>Adamc: interesting model output of horizontal velocity profile on Glimmer EISMINT-2 experiment F</p>
<hr />
<div>interesting model output of horizontal velocity profile on Glimmer EISMINT-2 experiment F</div>Adamchttp://websrv.cs.umt.edu/isis/index.php/PDX_afterhoursPDX afterhours2009-08-11T16:14:35Z<p>Adamc: added roots info</p>
<hr />
<div>[[Image:Erin_adam.JPG|thumb|right|300px|The PDX Afterhours king & queen in Tube Bar, home of Wednesday night $1 Miller High Life.]]<br />
<br />
==tuesday Aug 4, 2009==<br />
* meet at [http://paccinirestaurant.com/ Paccini's] pub at 7:30<br />
<br />
==wednesday Aug 5, 2009==<br />
* after FORTRAN session, leave PSU to go to [http://maps.google.com/maps?hl=en&client=firefox-a&rls=com.ubuntu:en-US:unofficial&hs=zuU&um=1&ie=UTF-8&q=sushi+ichiban+portland&fb=1&split=1&gl=us&view=text&latlng=899833397715679289 Sushi Ichiban]. Adam Campbell will guide you.<br />
* go to [http://www.groundkontrol.com Ground Kontrol] a retro arcade with beer<br />
* go to [http://www.voodoodoughnut.com Voodoo Doughnut], please someone buy Ian Rutt the $5 doughnut.<br />
<br />
<br />
==thursday Aug 6, 2009==<br />
[http://amontobin.com/field/ Amon Tobin] and [http://www.pitchblack.co.nz/?s1=index Pitch Black], along with two opening bands, are playing at the Roseland Theatre (8 NW 6th Ave) Thursday night (starting at 9:00 pm or thereabouts). Tickets are $26 (available online at [http://ticketswest.rdln.com/Venue.aspx?ven=ROS TicketsWest]). The music is best described as sampled electronica (Amon Tobin) and Kiwi-style dub (Pitch Black). I'd expect a late night of electronic music: an evening nap may be in order! Jeremy's already got his ticket and can fill you in with more, including a sample of the music.<br />
<br />
*other local music recommendations from Adam<br />
<br />
'''Boy Eats Drum Machine, French Miami, Southern Belle and Electric Opera Company''' - Indie Rock -<br />
Thu., Aug. 6, 9 p.m.<br />
$6-8<br />
Berbati's Pan<br />
10 SW 3rd Ave.<br />
Downtown<br />
<br />
'''Nurses, Inside Voices and Slaves''' - Indie Rock -<br />
Thu., Aug. 6, 8:30 p.m.<br />
$7<br />
Holocene <br />
1001 SE Morrison<br />
Southeast<br />
<br />
==friday Aug 7, 2009==<br />
[http://www.biteoforegon.com/ The Bite of Oregon] is a food festival that takes place at Tom McCall Waterfront Park. Featuring food, wine, beer and entertainment from Oregon. Entry is $8, food and beverages are extra.<br />
<br />
[http://www.pioneercourthousesquare.org/calendar_august.htm Flicks on the Bricks] will be showing Jurassic Park at dusk outside at Pioneer Courthouse Square, FREE (including popcorn). (10 minute walk)<br />
<br />
[http://www.portlandtwilight.com/ Portland Downtown Twilight Bike Criterium]: Professional (by U.S. standards) bike race through downtown Portland. Complete with beer garden, food, and an expected 15,000 spectators. Pro race starts at 7:30.<br />
<br />
==saturday Aug 8, 2009==<br />
[http://www.biteoforegon.com/ The Bite of Oregon] is a food festival that takes place at Tom McCall Waterfront Park. Featuring food, wine, beer and entertainment from Oregon. Entry is $8, food and beverages are extra.<br />
<br />
==sunday Aug 9, 2009==<br />
The [http://providence.org/bridgepedal/ Portland Bridge Pedal] is a fun event where you can go on 14, 24, or 37 mi. bike ride over Portland's Bridges. Adam is trying to assemble a group to go on Sunday morning. Please speak with him if you are interested in going to this. I tentatively have 4 bikes I can get ahold of.<br />
<br />
RIDERS: Register [http://providence.org/bridgepedal/ here]. Make sure to click the team rider option and the 11 bridge option! Ok so here is the deal I tried to register a team:<br />
Team Name The Icepocalypse<br />
Team password BPT943<br />
<br />
HEY It works now!<br />
<br />
[http://www.biteoforegon.com/ The Bite of Oregon] is a food festival that takes place at Tom McCall Waterfront Park. Featuring food, wine, beer and entertainment from Oregon. Entry is $8, food and beverages are extra.<br />
<br />
==tuesday Aug 11, 2009==<br />
[http://www.rootsorganicbrewing.com/ Roots Brewery] has a $2.50 imperial pint night. The beer is great! It's on the eastside but it's not very far. There are carts and more beer nearby.</div>Adamchttp://websrv.cs.umt.edu/isis/index.php/PDX_afterhoursPDX afterhours2009-08-08T20:59:22Z<p>Adamc: /* sunday Aug 9, 2009 */ added bridgepedal team info</p>
<hr />
<div>[[Image:Erin_adam.JPG|thumb|right|300px|The PDX Afterhours king & queen in Tube Bar, home of Wednesday night $1 Miller High Life.]]<br />
<br />
==tuesday Aug 4, 2009==<br />
* meet at [http://paccinirestaurant.com/ Paccini's] pub at 7:30<br />
<br />
==wednesday Aug 5, 2009==<br />
* after FORTRAN session, leave PSU to go to [http://maps.google.com/maps?hl=en&client=firefox-a&rls=com.ubuntu:en-US:unofficial&hs=zuU&um=1&ie=UTF-8&q=sushi+ichiban+portland&fb=1&split=1&gl=us&view=text&latlng=899833397715679289 Sushi Ichiban]. Adam Campbell will guide you.<br />
* go to [http://www.groundkontrol.com Ground Kontrol] a retro arcade with beer<br />
* go to [http://www.voodoodoughnut.com Voodoo Doughnut], please someone buy Ian Rutt the $5 doughnut.<br />
<br />
<br />
==thursday Aug 6, 2009==<br />
[http://amontobin.com/field/ Amon Tobin] and [http://www.pitchblack.co.nz/?s1=index Pitch Black], along with two opening bands, are playing at the Roseland Theatre (8 NW 6th Ave) Thursday night (starting at 9:00 pm or thereabouts). Tickets are $26 (available online at [http://ticketswest.rdln.com/Venue.aspx?ven=ROS TicketsWest]). The music is best described as sampled electronica (Amon Tobin) and Kiwi-style dub (Pitch Black). I'd expect a late night of electronic music: an evening nap may be in order! Jeremy's already got his ticket and can fill you in with more, including a sample of the music.<br />
<br />
*other local music recommendations from Adam<br />
<br />
'''Boy Eats Drum Machine, French Miami, Southern Belle and Electric Opera Company''' - Indie Rock -<br />
Thu., Aug. 6, 9 p.m.<br />
$6-8<br />
Berbati's Pan<br />
10 SW 3rd Ave.<br />
Downtown<br />
<br />
'''Nurses, Inside Voices and Slaves''' - Indie Rock -<br />
Thu., Aug. 6, 8:30 p.m.<br />
$7<br />
Holocene <br />
1001 SE Morrison<br />
Southeast<br />
<br />
==friday Aug 7, 2009==<br />
[http://www.biteoforegon.com/ The Bite of Oregon] is a food festival that takes place at Tom McCall Waterfront Park. Featuring food, wine, beer and entertainment from Oregon. Entry is $8, food and beverages are extra.<br />
<br />
[http://www.pioneercourthousesquare.org/calendar_august.htm Flicks on the Bricks] will be showing Jurassic Park at dusk outside at Pioneer Courthouse Square, FREE (including popcorn). (10 minute walk)<br />
<br />
[http://www.portlandtwilight.com/ Portland Downtown Twilight Bike Criterium]: Professional (by U.S. standards) bike race through downtown Portland. Complete with beer garden, food, and an expected 15,000 spectators. Pro race starts at 7:30.<br />
<br />
==saturday Aug 8, 2009==<br />
[http://www.biteoforegon.com/ The Bite of Oregon] is a food festival that takes place at Tom McCall Waterfront Park. Featuring food, wine, beer and entertainment from Oregon. Entry is $8, food and beverages are extra.<br />
<br />
==sunday Aug 9, 2009==<br />
The [http://providence.org/bridgepedal/ Portland Bridge Pedal] is a fun event where you can go on 14, 24, or 37 mi. bike ride over Portland's Bridges. Adam is trying to assemble a group to go on Sunday morning. Please speak with him if you are interested in going to this. I tentatively have 4 bikes I can get ahold of.<br />
<br />
RIDERS: Register [http://providence.org/bridgepedal/ here]. Make sure to click the team rider option and the 11 bridge option! Ok so here is the deal I tried to register a team:<br />
Team Name The Icepocalypse<br />
Team password BPT943<br />
<br />
But the form will not accept this password. It might be that the website takes a little while to process this but I'm not sure. I'd say if you can't use the password successfully by 4pm then just register as an individual rider... it cost $5 extra but whatever. I sent an email about this trying to get it figured out.<br />
<br />
[http://www.biteoforegon.com/ The Bite of Oregon] is a food festival that takes place at Tom McCall Waterfront Park. Featuring food, wine, beer and entertainment from Oregon. Entry is $8, food and beverages are extra.</div>Adamchttp://websrv.cs.umt.edu/isis/index.php/Erin_BurkeErin Burke2009-08-08T18:25:12Z<p>Adamc: </p>
<hr />
<div>i like creating awkward situations.</div>Adamchttp://websrv.cs.umt.edu/isis/index.php/Student_Presentation_DevelopmentStudent Presentation Development2009-08-07T21:40:31Z<p>Adamc: </p>
<hr />
<div>Students: Please email [[User:Mankoff|Ken Mankoff]] the Lat/Long of your institution and study area by the end of the day. (You can click 'email this user' in the lower left toolbox from his user page.<br />
<br />
==Roundtable 1: Study regions==<br />
* Antarctica. 6 students: [[User:Hoffman|Matt]], [[User:Adamc,|Adam]], [[User:mankoff|Ken Mankoff]]...<br />
* Greenland. 6 students: [[User:Kpoinar|Kristin]], [[User:meierbtw|Toby]]...<br />
* Mountain glaciers. 4 students: <br />
* Other/Global. 2 students:<br />
Possible discussion questions:<br />
* What questions are the climate change community pressuring us to answer?<br />
* What do we know now that would have been a big surprise 10 years ago?<br />
* How important is field data to your research? If you could collect any field data/observations to progress your work, what would it/they be? <br />
* <br />
*<br />
...<br />
<br />
==Roundtable 2: Background==<br />
* Geology. 5 students: [[User:meierbtw,|Toby]],...<br />
* Physics. 5 students: [[User:Adamc,|Adam]], [[User:Kpoinar|Kristin]]...<br />
* Math(s). 4 students:<br />
* Engineering/CS/Other. 4 students: [[User:Hoffman|Matt]], [[User:mankoff|Ken Mankoff]], ...<br />
Possible discussion questions:<br />
*<br />
* Given your background in XX, what have you done/would you do outside of your area of expertise to make yourself a better cryosphere scientist?<br />
* Why have we come to glaciology, over the other careers / topics available to people with our background?<br />
* Glaciology is 'hot' now... how hot will it stay? What are possible exit ramps once it cools down? Will data gaps between the current and next generation of satellite sensors affect what work can be done?<br />
* What should we do immediately post-PhD? travel, stay at home institution, ship off to postdoc abroad, ...<br />
* How can I continue to deal with people who say, "Oh, you study glaciers, huh? Better hurry up, hohnk hohnk"<br />
...<br />
<br />
==Other Ways to Sort Ourselves==<br />
Seneca - notes from yesterday?<br />
<br />
<br />
==Lists of valuable skills/knowledge learned, connections made, etc.?==</div>Adamchttp://websrv.cs.umt.edu/isis/index.php/Student_Presentation_DevelopmentStudent Presentation Development2009-08-07T21:31:50Z<p>Adamc: </p>
<hr />
<div>Students: Please email [[User:Mankoff|Ken Mankoff]] the Lat/Long of your institution and study area by the end of the day. (You can click 'email this user' in the lower left toolbox from his user page.<br />
<br />
==Roundtable 1: Study regions==<br />
* Antarctica. 6 students: [[User:Hoffman|Matt]], [[User:adamc,|Adam]], [[User:mankoff|Ken Mankoff]]...<br />
* Greenland. 6 students: [[User:Kpoinar|Kristin]], [[User:meierbtw|Toby]]...<br />
* Mountain glaciers. 4 students: <br />
* Other/Global. 2 students:<br />
Possible discussion questions:<br />
* What questions are the climate change community pressuring us to answer?<br />
* What do we know now that would have been a big surprise 10 years ago?<br />
* How important is field data to your research? If you could collect any field data/observations to progress your work, what would it/they be? <br />
* <br />
*<br />
...<br />
<br />
==Roundtable 2: Background==<br />
* Geology. 5 students: [[User:meierbtw,|Toby]],...<br />
* Physics. 5 students: [[User:Adamc,|Adam]], [[User:Kpoinar|Kristin]]...<br />
* Math(s). 4 students:<br />
* Engineering/CS/Other. 4 students: [[User:Hoffman|Matt]], [[User:mankoff|Ken Mankoff]], ...<br />
Possible discussion questions:<br />
*<br />
* Given your background in XX, what have you done/would you do outside of your area of expertise to make yourself a better cryosphere scientist?<br />
* Why have we come to glaciology, over the other careers / topics available to people with our background?<br />
* Glaciology is 'hot' now... how hot will it stay? What are possible exit ramps once it cools down?<br />
* What should we do immediately post-PhD? travel, stay at home institution, ship off to postdoc abroad, ...<br />
* How can I continue to deal with people who say, "Oh, you study glaciers, huh? Better hurry up, hohnk hohnk"<br />
...<br />
<br />
==Other Ways to Sort Ourselves==<br />
Seneca - notes from yesterday?<br />
<br />
<br />
==Lists of valuable skills/knowledge learned, connections made, etc.?==</div>Adamchttp://websrv.cs.umt.edu/isis/index.php/Student_Presentation_DevelopmentStudent Presentation Development2009-08-07T21:21:55Z<p>Adamc: added my info</p>
<hr />
<div>Students: Please email [[User:Mankoff|Ken Mankoff]] the Lat/Long of your institution and study area by the end of the day. (You can click 'email this user' in the lower left toolbox from his user page.<br />
<br />
==Roundtable 1: Study regions==<br />
* Antarctica. 6 students: [[User:Hoffman|Matt]], , [[User:Adamc,|Adam]],...<br />
* Greenland. 6 students: [[User:Kpoinar|Kristin]], ...<br />
* Mountain glaciers. 4 students: <br />
* Other/Global. 2 students:<br />
Possible discussion questions:<br />
* What questions are the climate change community pressuring us to answer?<br />
* What do we know now that would have been a big surprise 10 years ago?<br />
* <br />
* <br />
*<br />
...<br />
<br />
==Roundtable 2: Background==<br />
* Geology. 5 students:<br />
* Physics. 5 students: [[User:Adamc,|Adam]],...<br />
* Math(s). 4 students:<br />
* Engineering/CS/Other. 4 students: [[User:Hoffman|Matt]], ...<br />
Possible discussion questions:<br />
*<br />
*<br />
...<br />
<br />
==Other Ways to Sort Ourselves==<br />
Seneca - notes from yesterday?<br />
<br />
<br />
==Lists of valuable skills/knowledge learned, connections made, etc.?==</div>Adamchttp://websrv.cs.umt.edu/isis/index.php/COMSOL_activitiesCOMSOL activities2009-08-07T20:57:00Z<p>Adamc: </p>
<hr />
<div>==Overview==<br />
Now, let's see if COMSOL can be used to solve problems of glaciological relevance. We'll look at shallow ice approximation flow, and shallow shelf approximation flows.<br />
<br />
==Isothermal Shallow Ice Approximation==<br />
<br />
Begin with the often used shallow ice form for ice thickness evolution, which casts evolution as a non-linear diffusion problem<br />
<br />
:<math>\frac{\partial H}{\partial t} = - \nabla D \nabla H + M</math><br />
where<br />
:<math>D = \frac{2A(\rho g)^n}{n+2} H^{n+2} \left[\nabla H \cdot \nabla H \right]^{(n-1)/2}</math><br />
<br />
with boundary condition <math>H=0</math> on the edge of the computational domain.<br />
<br />
===Comsol Modeling===<br />
We will use the '''PDE, General Form''' transient mode to solve this equation. For convenience, make the dependent variable <math>H</math>.<br />
<br />
====Geometry====<br />
You should not find it difficult to create a unit square. Once it's made, you can double click it to change it's size and do other transformations. Read below to find the appropriate domain.<br />
<br />
====Field equations====<br />
This equation mode solves equations of the form<br />
:<math>e_a\frac{\partial^2 H}{\partial t^2} + d_a \frac{\partial H}{\partial t} + \nabla \cdot \Gamma = F</math><br />
<br />
Which is just what we want if we recognize that in our system <math>e_a</math>=0, <math>d_a</math>=1, <math>F</math>=<math>M</math>, and<br />
:<math>\Gamma_x = -\frac{2A(\rho g)^n}{n+2} H^{n+2} \left[\nabla H \cdot \nabla H \right]^{(n-1)/2} \frac{\partial H}{\partial x} </math><br />
<br />
:<math>\Gamma_y = -\frac{2A(\rho g)^n}{n+2} H^{n+2} \left[\nabla H \cdot \nabla H \right]^{(n-1)/2} \frac{\partial H}{\partial y} </math><br />
<br />
Now the problem has been reduced to one of typing. It will make the COMSOL model easier to read if you create a scalar expression for <math>D</math> . Then your <math>\Gamma_x = -D \frac{\partial H}{\partial x}</math> and <math>\Gamma_y =- D\frac{\partial H}{\partial y}</math> are very clear.<br />
<br />
====Boundary conditions====<br />
This type of problem requires a Dirchlet boundary condition. Set <math>H</math> = 0 on all four sides.<br />
<br />
====Other====<br />
You'll also be needing to know how to tell COMSOL to use a derivative. That is '''Hx''', '''Hy''', and '''Ht''' for <math>\frac{\partial H}{\partial x}</math>, <math>\frac{\partial H}{\partial y}</math>, and <math>\frac{\partial H}{\partial t}</math> respectively.<br />
<br />
===Exercises===<br />
#Complete the model, and do the isothermal ''fixed margin'' experiment Huybrechts (1996)<ref name="Huybrechts">Huybrechts et al. The EISMINT Benchmarks for Testing Ice--Sheet Models. Ann. Glaciol. (1996) vol. 23 pp. 1-12 [http://homepages.vub.ac.be/~phuybrec/pdf/Huyb.Ann.Glac.23.pdf pdf]</ref>. You'll find all values of constants there as well. Verify that your model is providing results consistent with those reported in the paper.<br />
#Now alter your model (the accumulation field) to do the isothermal ''moving margin'', again verify that it's at least just as wrong as the other models. You're going to have to come up with something to deal with the negative values of thickness that you'll get...<br />
<br />
==Shallow shelf approximation==<br />
===Field equations===<br />
Now, consider the equations describing a flow that is vertically integrated. The equations are<br />
:<math>\frac{\partial}{\partial x}\left ( 2 \eta H <br />
\left(2\frac{\partial u}{\partial x}+\frac{\partial v}{\partial y}\right)\right)<br />
+\frac{\partial}{\partial y}\left(\eta H\left(<br />
\frac{\partial u}{\partial y}+\frac{\partial v}{\partial x}\right)\right)<br />
=\rho gH \frac{\partial s}{\partial x}<br />
</math><br />
<br />
:<math><br />
\frac{\partial}{\partial y}\left ( 2 \eta H <br />
\left(2\frac{\partial v}{\partial y}+\frac{\partial u}{\partial x}\right)\right)<br />
+\frac{\partial}{\partial x}\left(\eta H\left(<br />
\frac{\partial u}{\partial y}+\frac{\partial v}{\partial x}\right)\right)<br />
=\rho gH \frac{\partial s}{\partial y}<br />
</math><br />
<br />
<math> \eta</math> is the non-linear, vertically averaged viscosity. <br />
It will need to be entered as a '''scalar expression''', and is written<br />
<br />
:<math>\eta = \frac{B}{2}\left[ \left(\frac{\partial u}{\partial x}\right)^2 + \left(\frac{\partial v}{\partial y}\right)^2 + \frac{1}{4} \left(\frac{\partial u}{\partial x} + \frac{\partial v}{\partial y}\right)^2 + \frac{\partial u}{\partial x}\frac{\partial v}{\partial y}\right]^{-1/n}</math><br />
===Boundary conditions===<br />
There are two flavors of boundary conditions applied in the typical ice sheet model. <br />
====Kinematic====<br />
First, the Dirichlet, or ''kinematic'' boundary condition specify the velocity<br />
:<math> \mathbf{u} = \mathbf{u_b},~\forall \partial \Omega_k \in \partial \Omega.</math><br />
This is the boundary condition where the ice moving across the grounding line. In order to model an ice shelf, one determines (or estimates), and specifies that velocity.<br />
====Dynamic====<br />
This is the Nuemann, or ''dynamic'' boundary condition that is applied along the ice front,<br />
<br />
:<math>- 2 \eta H \left(2\frac{\partial u}{\partial x}+\frac{\partial v}{\partial y}\right)\mathbf{n_x}<br />
<br />
-\eta H\left(<br />
\frac{\partial u}{\partial y}+\frac{\partial v}{\partial x}\right)\mathbf{n_y}<br />
= <br />
-\rho g H^2\left(1-\frac{\rho_i}{\rho_w}\right) \mathbf{n_x} </math><br />
<br />
and<br />
<br />
<math>- \eta H \left(\frac{\partial u}{\partial y}+\frac{\partial v}{\partial x}\right)\mathbf{n_x}<br />
<br />
-2\eta H\left(<br />
\frac{\partial u}{\partial x}+2\frac{\partial v}{\partial y}\right)\mathbf{n_y}<br />
= <br />
-\rho g H^2\left(1-\frac{\rho_i}{\rho_w}\right) \mathbf{n_y} </math><br />
<br />
<br />
:<math> ~\forall \partial \Omega_d \in \partial \Omega.</math><br />
<br />
Note that in case you will need to know how to tell Comsol the normal vectors. These are '''nx''', and '''ny''''.<br />
<br />
===Exercises===<br />
#As a first exercise in solving these equations, try the experiments described in the EISMINT ice shelf models, but never published [http://homepages.vub.ac.be/~phuybrec/eismint/iceshelf.html]. Get the '''self-descr.pdf''', or the first hyperlink on the page. Let's do experiments 3-4 on page 6 of the document (Note that we will see and work with the solution to these experiments again when we [[Adding a module to Glimmer I|do some exercises with the higher-order dynamics routines in Glimmer/CISM]]).<br />
# The above is neat, but ultimately not that useful because the relation realistic geometry is weak. Try out [[Media:Ross.tar.gz|this]] model, that has the geometry and boundary conditions for the Ross Ice shelf, but otherwise the same equations solved in the previous exercise. See the utility of this? Is the solution dependent on the mesh? Can you do anything with the solver to improve the time required for a solution.<br />
<br />
==References==<br />
<references/></div>Adamchttp://websrv.cs.umt.edu/isis/index.php/PDX_afterhoursPDX afterhours2009-08-06T00:57:04Z<p>Adamc: /* sunday Aug 9, 2009 */ updated bicycle count</p>
<hr />
<div>==tuesday Aug 4, 2009==<br />
* meet at [http://paccinirestaurant.com/ Paccini's] pub at 7:30<br />
<br />
==wednesday Aug 5, 2009==<br />
* after FORTRAN session, leave PSU to go to [http://maps.google.com/maps?hl=en&client=firefox-a&rls=com.ubuntu:en-US:unofficial&hs=zuU&um=1&ie=UTF-8&q=sushi+ichiban+portland&fb=1&split=1&gl=us&view=text&latlng=899833397715679289 Sushi Ichiban]. Adam Campbell will guide you.<br />
* go to [http://www.groundkontrol.com Ground Kontrol] a retro arcade with beer<br />
* go to [http://www.voodoodoughnut.com Voodoo Doughnut], please someone buy Ian Rutt the $5 doughnut.<br />
<br />
==thursday Aug 6, 2009==<br />
[http://amontobin.com/field/ Amon Tobin] and [http://www.pitchblack.co.nz/?s1=index Pitch Black], along with two opening bands, are playing at the Roseland Theatre (8 NW 6th Ave) Thursday night (starting at 9:00 pm or thereabouts). Tickets are $26 (available online at [http://ticketswest.rdln.com/Venue.aspx?ven=ROS TicketsWest]). The music is best described as sampled electronica (Amon Tobin) and Kiwi-style dub (Pitch Black). I'd expect a late night of electronic music: an evening nap may be in order! Jeremy's already got his ticket and can fill you in with more, including a sample of the music.<br />
<br />
*other local music recommendations from Adam<br />
<br />
'''Boy Eats Drum Machine, French Miami, Southern Belle and Electric Opera Company''' - Indie Rock -<br />
Thu., Aug. 6, 9 p.m.<br />
$6-8<br />
Berbati's Pan<br />
10 SW 3rd Ave.<br />
Downtown<br />
<br />
'''Nurses, Inside Voices and Slaves''' - Indie Rock -<br />
Thu., Aug. 6, 8:30 p.m.<br />
$7<br />
Holocene <br />
1001 SE Morrison<br />
Southeast<br />
<br />
==friday Aug 7, 2009==<br />
[http://www.biteoforegon.com/ The Bite of Oregon] is a food festival that takes place at Tom McCall Waterfront Park. Featuring food, wine, beer and entertainment from Oregon. Entry is $8, food and beverages are extra.<br />
<br />
==saturday Aug 8, 2009==<br />
[http://www.biteoforegon.com/ The Bite of Oregon] is a food festival that takes place at Tom McCall Waterfront Park. Featuring food, wine, beer and entertainment from Oregon. Entry is $8, food and beverages are extra.<br />
<br />
==sunday Aug 9, 2009==<br />
The [http://providence.org/bridgepedal/ Portland Bridge Pedal] is a fun event where you can go on 14, 24, or 37 mi. bike ride over Portland's Bridges. Adam is trying to assemble a group to go on Sunday morning. Please speak with him if you are interested in going to this. I tentatively have 4 bikes I can get ahold of.<br />
<br />
[http://www.biteoforegon.com/ The Bite of Oregon] is a food festival that takes place at Tom McCall Waterfront Park. Featuring food, wine, beer and entertainment from Oregon. Entry is $8, food and beverages are extra.</div>Adamchttp://websrv.cs.umt.edu/isis/index.php/Typical_Daily_ScheduleTypical Daily Schedule2009-08-06T00:16:22Z<p>Adamc: /* Special Events */ clarified location of Nancy's home</p>
<hr />
<div>==A typical day==<br />
<br />
{| border="1" cellpadding="5" cellspacing="0"<br />
|-valign="top" style="background:RoyalBlue"<br />
!width="15%"|Time<br />
!width="15%"|What's happening?<br />
!width="15%"|Location <br />
<br />
|-valign="top" style="background:AliceBlue"<br />
| 7:30-8:00<br />
| Breakfast<br />
| Cramer Hall 17<br />
|-valign="top" style="background:PowderBlue"<br />
| 8:00-10:00<br />
| Morning Lectures I<br />
| Cramer Hall, Room 1<br />
|-valign="top" style="background:AliceBlue"<br />
| 10:00-10:30<br />
| Coffee Break<br />
| Cramer Hall outside Room 1<br />
|-valign="top" style="background:PowderBlue"<br />
| 10:30-12:00<br />
| Morning Lectures II<br />
| Cramer Hall 1<br />
|-valign="top" style="background:AliceBlue"<br />
| 12:00-13:30<br />
| Lunch<br />
| <br />
|-valign="top" style="background:PowderBlue"<br />
| 13:30-15:30<br />
| Afternoon Activities I<br />
| Cramer Hall 1<br />
|-valign="top" style="background:AliceBlue"<br />
| 15:30-16:00<br />
| Coffee Break<br />
| Cramer Hall outside Room 1<br />
|-valign="top" style="background:PowderBlue"<br />
| 16:00-17:30<br />
| Afternoon Activities II<br />
| Cramer Hall 1<br />
|}<br />
<br />
==Special Events==<br />
<br />
* Please join us at 6 pm August 3rd for a welcome dinner at Hot Lips Pizza. The address is 1909 SW 6th Avenue, between College and Hall Streets, just a few blocks from University Place. <br />
* Breakfast will be provided in the Geology Department lounge every class meeting day.<br />
* Morning and afternoon coffee breaks will be provided each class day. If you want coffee before this, we recommend Cafe Ono on 5th between Hall and Harrison (where you can find Stumptown Coffee).<br />
* Catered lunches will be provided on Tuesday the 4th and Tuesday the 11th.<br />
* Saturday will be a working day, with the computer room open.<br />
* Sunday will be a day off with a PDX tour if there is interest.<br />
* Monday will be an excursion to Mt. Hood, Eliot Glacier.<br />
* There will be a final group dinner on Friday the 14th at Nancy's home (along the blue MAX line; must backtrack slightly to get red MAX to airport).</div>Adamchttp://websrv.cs.umt.edu/isis/index.php/PDX_afterhoursPDX afterhours2009-08-06T00:02:57Z<p>Adamc: /* sunday Aug 9, 2009 */ added the Bite of Oregon</p>
<hr />
<div>==tuesday Aug 4, 2009==<br />
* meet at [http://paccinirestaurant.com/ Paccini's] pub at 7:30<br />
<br />
==wednesday Aug 5, 2009==<br />
* after FORTRAN session, leave PSU to go to [http://maps.google.com/maps?hl=en&client=firefox-a&rls=com.ubuntu:en-US:unofficial&hs=zuU&um=1&ie=UTF-8&q=sushi+ichiban+portland&fb=1&split=1&gl=us&view=text&latlng=899833397715679289 Sushi Ichiban]. Adam Campbell will guide you.<br />
* go to [http://www.groundkontrol.com Ground Kontrol] a retro arcade with beer<br />
* go to [http://www.voodoodoughnut.com Voodoo Doughnut], please someone buy Ian Rutt the $5 doughnut.<br />
<br />
==thursday Aug 6, 2009==<br />
[http://amontobin.com/field/ Amon Tobin] and [http://www.pitchblack.co.nz/?s1=index Pitch Black], along with two opening bands, are playing at the Roseland Theatre (8 NW 6th Ave) Thursday night (starting at 9:00 pm or thereabouts). Tickets are $26 (available online at [http://ticketswest.rdln.com/Venue.aspx?ven=ROS TicketsWest]). The music is best described as sampled electronica (Amon Tobin) and Kiwi-style dub (Pitch Black). I'd expect a late night of electronic music: an evening nap may be in order! Jeremy's already got his ticket and can fill you in with more, including a sample of the music.<br />
<br />
*other local music recommendations from Adam<br />
<br />
'''Boy Eats Drum Machine, French Miami, Southern Belle and Electric Opera Company''' - Indie Rock -<br />
Thu., Aug. 6, 9 p.m.<br />
$6-8<br />
Berbati's Pan<br />
10 SW 3rd Ave.<br />
Downtown<br />
<br />
'''Nurses, Inside Voices and Slaves''' - Indie Rock -<br />
Thu., Aug. 6, 8:30 p.m.<br />
$7<br />
Holocene <br />
1001 SE Morrison<br />
Southeast<br />
<br />
==friday Aug 7, 2009==<br />
[http://www.biteoforegon.com/ The Bite of Oregon] is a food festival that takes place at Tom McCall Waterfront Park. Featuring food, wine, beer and entertainment from Oregon. Entry is $8, food and beverages are extra.<br />
<br />
==saturday Aug 8, 2009==<br />
[http://www.biteoforegon.com/ The Bite of Oregon] is a food festival that takes place at Tom McCall Waterfront Park. Featuring food, wine, beer and entertainment from Oregon. Entry is $8, food and beverages are extra.<br />
<br />
==sunday Aug 9, 2009==<br />
The [http://providence.org/bridgepedal/ Portland Bridge Pedal] is a fun event where you can go on 14, 24, or 37 mi. bike ride over Portland's Bridges. Adam is trying to assemble a group to go on Sunday morning. Please speak with him if you are interested in going to this.<br />
<br />
[http://www.biteoforegon.com/ The Bite of Oregon] is a food festival that takes place at Tom McCall Waterfront Park. Featuring food, wine, beer and entertainment from Oregon. Entry is $8, food and beverages are extra.</div>Adamchttp://websrv.cs.umt.edu/isis/index.php/PDX_afterhoursPDX afterhours2009-08-06T00:02:38Z<p>Adamc: /* saturday Aug 8, 2009 */ add the Bite of Oregon</p>
<hr />
<div>==tuesday Aug 4, 2009==<br />
* meet at [http://paccinirestaurant.com/ Paccini's] pub at 7:30<br />
<br />
==wednesday Aug 5, 2009==<br />
* after FORTRAN session, leave PSU to go to [http://maps.google.com/maps?hl=en&client=firefox-a&rls=com.ubuntu:en-US:unofficial&hs=zuU&um=1&ie=UTF-8&q=sushi+ichiban+portland&fb=1&split=1&gl=us&view=text&latlng=899833397715679289 Sushi Ichiban]. Adam Campbell will guide you.<br />
* go to [http://www.groundkontrol.com Ground Kontrol] a retro arcade with beer<br />
* go to [http://www.voodoodoughnut.com Voodoo Doughnut], please someone buy Ian Rutt the $5 doughnut.<br />
<br />
==thursday Aug 6, 2009==<br />
[http://amontobin.com/field/ Amon Tobin] and [http://www.pitchblack.co.nz/?s1=index Pitch Black], along with two opening bands, are playing at the Roseland Theatre (8 NW 6th Ave) Thursday night (starting at 9:00 pm or thereabouts). Tickets are $26 (available online at [http://ticketswest.rdln.com/Venue.aspx?ven=ROS TicketsWest]). The music is best described as sampled electronica (Amon Tobin) and Kiwi-style dub (Pitch Black). I'd expect a late night of electronic music: an evening nap may be in order! Jeremy's already got his ticket and can fill you in with more, including a sample of the music.<br />
<br />
*other local music recommendations from Adam<br />
<br />
'''Boy Eats Drum Machine, French Miami, Southern Belle and Electric Opera Company''' - Indie Rock -<br />
Thu., Aug. 6, 9 p.m.<br />
$6-8<br />
Berbati's Pan<br />
10 SW 3rd Ave.<br />
Downtown<br />
<br />
'''Nurses, Inside Voices and Slaves''' - Indie Rock -<br />
Thu., Aug. 6, 8:30 p.m.<br />
$7<br />
Holocene <br />
1001 SE Morrison<br />
Southeast<br />
<br />
==friday Aug 7, 2009==<br />
[http://www.biteoforegon.com/ The Bite of Oregon] is a food festival that takes place at Tom McCall Waterfront Park. Featuring food, wine, beer and entertainment from Oregon. Entry is $8, food and beverages are extra.<br />
<br />
==saturday Aug 8, 2009==<br />
[http://www.biteoforegon.com/ The Bite of Oregon] is a food festival that takes place at Tom McCall Waterfront Park. Featuring food, wine, beer and entertainment from Oregon. Entry is $8, food and beverages are extra.<br />
<br />
==sunday Aug 9, 2009==<br />
The [http://providence.org/bridgepedal/ Portland Bridge Pedal] is a fun event where you can go on 14, 24, or 37 mi. bike ride over Portland's Bridges. Adam is trying to assemble a group to go on Sunday morning. Please speak with him if you are interested in going to this.</div>Adamchttp://websrv.cs.umt.edu/isis/index.php/PDX_afterhoursPDX afterhours2009-08-06T00:02:26Z<p>Adamc: /* friday Aug 7, 2009 */ added the Bite of Oregon</p>
<hr />
<div>==tuesday Aug 4, 2009==<br />
* meet at [http://paccinirestaurant.com/ Paccini's] pub at 7:30<br />
<br />
==wednesday Aug 5, 2009==<br />
* after FORTRAN session, leave PSU to go to [http://maps.google.com/maps?hl=en&client=firefox-a&rls=com.ubuntu:en-US:unofficial&hs=zuU&um=1&ie=UTF-8&q=sushi+ichiban+portland&fb=1&split=1&gl=us&view=text&latlng=899833397715679289 Sushi Ichiban]. Adam Campbell will guide you.<br />
* go to [http://www.groundkontrol.com Ground Kontrol] a retro arcade with beer<br />
* go to [http://www.voodoodoughnut.com Voodoo Doughnut], please someone buy Ian Rutt the $5 doughnut.<br />
<br />
==thursday Aug 6, 2009==<br />
[http://amontobin.com/field/ Amon Tobin] and [http://www.pitchblack.co.nz/?s1=index Pitch Black], along with two opening bands, are playing at the Roseland Theatre (8 NW 6th Ave) Thursday night (starting at 9:00 pm or thereabouts). Tickets are $26 (available online at [http://ticketswest.rdln.com/Venue.aspx?ven=ROS TicketsWest]). The music is best described as sampled electronica (Amon Tobin) and Kiwi-style dub (Pitch Black). I'd expect a late night of electronic music: an evening nap may be in order! Jeremy's already got his ticket and can fill you in with more, including a sample of the music.<br />
<br />
*other local music recommendations from Adam<br />
<br />
'''Boy Eats Drum Machine, French Miami, Southern Belle and Electric Opera Company''' - Indie Rock -<br />
Thu., Aug. 6, 9 p.m.<br />
$6-8<br />
Berbati's Pan<br />
10 SW 3rd Ave.<br />
Downtown<br />
<br />
'''Nurses, Inside Voices and Slaves''' - Indie Rock -<br />
Thu., Aug. 6, 8:30 p.m.<br />
$7<br />
Holocene <br />
1001 SE Morrison<br />
Southeast<br />
<br />
==friday Aug 7, 2009==<br />
[http://www.biteoforegon.com/ The Bite of Oregon] is a food festival that takes place at Tom McCall Waterfront Park. Featuring food, wine, beer and entertainment from Oregon. Entry is $8, food and beverages are extra.<br />
<br />
==saturday Aug 8, 2009==<br />
<br />
==sunday Aug 9, 2009==<br />
The [http://providence.org/bridgepedal/ Portland Bridge Pedal] is a fun event where you can go on 14, 24, or 37 mi. bike ride over Portland's Bridges. Adam is trying to assemble a group to go on Sunday morning. Please speak with him if you are interested in going to this.</div>Adamchttp://websrv.cs.umt.edu/isis/index.php/PDX_afterhoursPDX afterhours2009-08-05T22:59:00Z<p>Adamc: /* thursday Aug 6, 2009 */ added possible shows</p>
<hr />
<div>==tuesday Aug 4, 2009==<br />
* meet at [http://paccinirestaurant.com/ Paccini's] pub at 7:30<br />
<br />
==wednesday Aug 5, 2009==<br />
* after FORTRAN session, leave PSU to go to [http://maps.google.com/maps?hl=en&client=firefox-a&rls=com.ubuntu:en-US:unofficial&hs=zuU&um=1&ie=UTF-8&q=sushi+ichiban+portland&fb=1&split=1&gl=us&view=text&latlng=899833397715679289 Sushi Ichiban]. Adam Campbell will guide you.<br />
* go to [http://www.groundkontrol.com Ground Kontrol] a retro arcade with beer<br />
* go to [http://www.voodoodoughnut.com Voodoo Doughnut], please someone buy Ian Rutt the $5 doughnut.<br />
<br />
==thursday Aug 6, 2009==<br />
[http://amontobin.com/field/ Amon Tobin] and [http://www.pitchblack.co.nz/?s1=index Pitch Black], along with two opening bands, are playing at the Roseland Theatre (8 NW 6th Ave) Thursday night (starting at 9:00 pm or thereabouts). Tickets are $26 (available online at [http://ticketswest.rdln.com/Venue.aspx?ven=ROS TicketsWest]). The music is best described as sampled electronica (Amon Tobin) and Kiwi-style dub (Pitch Black). I'd expect a late night of electronic music: an evening nap may be in order! Jeremy's already got his ticket and can fill you in with more, including a sample of the music.<br />
<br />
*other local music recommendations from Adam<br />
<br />
'''Boy Eats Drum Machine, French Miami, Southern Belle and Electric Opera Company''' - Indie Rock -<br />
Thu., Aug. 6, 9 p.m.<br />
$6-8<br />
Berbati's Pan<br />
10 SW 3rd Ave.<br />
Downtown<br />
<br />
'''Nurses, Inside Voices and Slaves''' - Indie Rock -<br />
Thu., Aug. 6, 8:30 p.m.<br />
$7<br />
Holocene <br />
1001 SE Morrison<br />
Southeast<br />
<br />
==friday Aug 7, 2009==<br />
<br />
==saturday Aug 8, 2009==<br />
<br />
==sunday Aug 9, 2009==<br />
The [http://providence.org/bridgepedal/ Portland Bridge Pedal] is a fun event where you can go on 14, 24, or 37 mi. bike ride over Portland's Bridges. Adam is trying to assemble a group to go on Sunday morning. Please speak with him if you are interested in going to this.</div>Adamchttp://websrv.cs.umt.edu/isis/index.php/PDX_afterhoursPDX afterhours2009-08-05T22:51:38Z<p>Adamc: /* wednesday Aug 5, 2009 */ updated agenda</p>
<hr />
<div>==tuesday Aug 4, 2009==<br />
* meet at [http://paccinirestaurant.com/ Paccini's] pub at 7:30<br />
<br />
==wednesday Aug 5, 2009==<br />
* after FORTRAN session, leave PSU to go to [http://maps.google.com/maps?hl=en&client=firefox-a&rls=com.ubuntu:en-US:unofficial&hs=zuU&um=1&ie=UTF-8&q=sushi+ichiban+portland&fb=1&split=1&gl=us&view=text&latlng=899833397715679289 Sushi Ichiban]. Adam Campbell will guide you.<br />
* go to [http://www.groundkontrol.com Ground Kontrol] a retro arcade with beer<br />
* go to [http://www.voodoodoughnut.com Voodoo Doughnut], please someone buy Ian Rutt the $5 doughnut.<br />
<br />
==thursday Aug 6, 2009==<br />
[http://amontobin.com/field/ Amon Tobin] and [http://www.pitchblack.co.nz/?s1=index Pitch Black], along with two opening bands, are playing at the Roseland Theatre (8 NW 6th Ave) Thursday night (starting at 9:00 pm or thereabouts). Tickets are $26 (available online at [http://ticketswest.rdln.com/Venue.aspx?ven=ROS TicketsWest]). The music is best described as sampled electronica (Amon Tobin) and Kiwi-style dub (Pitch Black). I'd expect a late night of electronic music: an evening nap may be in order! Jeremy's already got his ticket and can fill you in with more, including a sample of the music.<br />
<br />
==friday Aug 7, 2009==<br />
<br />
==saturday Aug 8, 2009==<br />
<br />
==sunday Aug 9, 2009==<br />
The [http://providence.org/bridgepedal/ Portland Bridge Pedal] is a fun event where you can go on 14, 24, or 37 mi. bike ride over Portland's Bridges. Adam is trying to assemble a group to go on Sunday morning. Please speak with him if you are interested in going to this.</div>Adamchttp://websrv.cs.umt.edu/isis/index.php/Team_2_SolutionTeam 2 Solution2009-08-05T22:21:35Z<p>Adamc: updated code and added python script</p>
<hr />
<div>* FORTRAN mountain glacier<br />
<br />
<source lang="fortran"><br />
program Kees1<br />
<br />
implicit none<br />
real, parameter :: dt = 1/12.0<br />
real, parameter :: dx = 1<br />
real, parameter :: tend = 5000<br />
real, parameter :: xend = 30<br />
real, parameter :: n = 3<br />
real, parameter :: rho = 920<br />
real, parameter :: g = 9.2<br />
real, parameter :: A = 1.0e-16<br />
real, parameter :: Hunnecessary = 6378.1 * sqrt(9.8/9.2-1) !height where g =9.2m/s^2<br />
real, parameter :: C = 2*A*(rho*g)**n/(n+2)*(1.e3)**n<br />
integer :: nt = int(tend/dt)<br />
integer :: nx = int(xend/dx)<br />
real, dimension(int(tend/dt)+1,int(xend/dx)+1) :: S<br />
integer :: t,i<br />
real, dimension(int(xend/dx)+1)::xarr=(/((i-1)*dx,i=1,&<br />
int(xend/dx)+1)/), MB, H, B<br />
real, dimension(int(xend/dx)) :: Diff<br />
<br />
B = 1-0.0001*xarr <br />
MB = 0.004 - xarr * 0.0002<br />
<br />
S(1,:) = B<br />
<br />
do t=2, nt+1<br />
S(t,1)=0<br />
H = S(t-1,:) - B<br />
Diff=C*((H(1:nx)+H(2:nx+1))/2)**(n+2)&<br />
*abs((S(t-1,2:nx+1)-S(t-1,1:nx))/dx)**(n-1)<br />
S(t,2:nx)=S(t-1,2:nx)+(dt/dx**2)&<br />
*(Diff(2:nx)*(S(t-1,3:nx+1)-S(t-1,2:nx))&<br />
- Diff(1:nx-1)*(S(t-1,2:nx)-S(t-1,1:nx-1)))&<br />
+MB(2:nx)*dt<br />
where(S(t,:)<B)<br />
S(t,:)=B<br />
end where<br />
if (mod(t,100) .eq.0) then<br />
write (*,*) S(t,:)+Hunnecessary <br />
end if<br />
end do<br />
end program<br />
<br />
</source><br />
<br />
*python plotting script<br />
<br />
<source lang="python"><br />
#!/usr/bin/env python<br />
<br />
# Import only what is needed<br />
from numpy import loadtxt,shape,linspace<br />
from pylab import plot,show,clf,show,ion,ylim<br />
<br />
# Import data file, called 'data'<br />
d=loadtxt('data')<br />
<br />
#Determine how much data came in<br />
dims = shape(d)<br />
<br />
clf() # Clears the screen<br />
ion() # Interactive plot mode, critical for animation<br />
<br />
# x data, note that this must correspond to program's domain<br />
x = linspace(0,30,dims[1]) <br />
<br />
# Initial plot, very Matlab(ish), note return of plot handle that allows plot to<br />
# be altered elsewhere in code.<br />
ph,=plot(x,d[0,:],'k') <br />
ph.figure.show() # matplot lib requires show to be called<br />
ylim(1629.7,1630.5)<br />
<br />
# Loop to plot each time step<br />
for i in range(1,dims[0]):<br />
ph.set_ydata(d[i,:]) # Only update y data (faster than replot)<br />
ph.figure.show()<br />
<br />
</source></div>Adamchttp://websrv.cs.umt.edu/isis/index.php/Team_2_SolutionTeam 2 Solution2009-08-05T22:08:52Z<p>Adamc: deleted unused line</p>
<hr />
<div><source lang="fortran"><br />
program Kees1<br />
<br />
implicit none<br />
real, parameter :: dt = 1/12.0<br />
real, parameter :: dx = 1<br />
real, parameter :: tend = 5000<br />
real, parameter :: xend = 30<br />
real, parameter :: n = 3<br />
real, parameter :: rho = 920<br />
real, parameter :: g = 9.2<br />
real, parameter :: A = 1.0e-16<br />
real, parameter :: C = 2*A*(rho*g)**n/(n+2)*(1.e3)**n<br />
integer :: nt = int(tend/dt)<br />
integer :: nx = int(xend/dx)<br />
real, dimension(int(tend/dt)+1,int(xend/dx)+1) :: S<br />
integer :: t,i<br />
real, dimension(int(xend/dx)+1)::xarr=(/((i-1)*dx,i=1,&<br />
int(xend/dx)+1)/), MB, H, B<br />
real, dimension(int(xend/dx)) :: Diff<br />
<br />
B = 1-0.0001*xarr <br />
MB = 0.004 - xarr * 0.0002<br />
<br />
S(1,:) = B<br />
<br />
do t=2, nt+1<br />
S(t,1)=0<br />
H = S(t-1,:) - B<br />
Diff=C*((H(1:nx)+H(2:nx+1))/2)**(n+2)&<br />
*abs((S(t-1,2:nx+1)-S(t-1,1:nx))/dx)**(n-1)<br />
S(t,2:nx)=S(t-1,2:nx)+(dt/dx**2)&<br />
*(Diff(2:nx)*(S(t-1,3:nx+1)-S(t-1,2:nx))&<br />
- Diff(1:nx-1)*(S(t-1,2:nx)-S(t-1,1:nx-1)))&<br />
+MB(2:nx)*dt<br />
where(S(t,:)<B)<br />
S(t,:)=B<br />
end where<br />
if (mod(t,100) .eq.0) then<br />
write (*,*) S(t,:)+1633.42 !height at which g=9.2<br />
end if<br />
end do<br />
end program<br />
</source></div>Adamchttp://websrv.cs.umt.edu/isis/index.php/Team_2_SolutionTeam 2 Solution2009-08-05T22:07:06Z<p>Adamc: added source code</p>
<hr />
<div><source lang="fortran"><br />
program Kees1<br />
<br />
implicit none<br />
real, parameter :: dt = 1/12.0<br />
real, parameter :: dx = 1<br />
real, parameter :: tend = 5000<br />
real, parameter :: xend = 30<br />
real, parameter :: n = 3<br />
real, parameter :: rho = 920<br />
real, parameter :: g = 9.2<br />
real, parameter :: A = 1.0e-16<br />
real, parameter :: C = 2*A*(rho*g)**n/(n+2)*(1.e3)**n<br />
integer :: nt = int(tend/dt)<br />
integer :: nx = int(xend/dx)<br />
real, dimension(int(tend/dt)+1,int(xend/dx)+1) :: S<br />
integer :: t,i<br />
real, dimension(int(xend/dx)+1)::xarr=(/((i-1)*dx,i=1,&<br />
int(xend/dx)+1)/), MB, H, B<br />
real, dimension(int(xend/dx)) :: Diff<br />
<br />
B = 1-0.0001*xarr <br />
MB = 0.004 - xarr * 0.0002<br />
<br />
S(1,:) = B<br />
<br />
do t=2, nt+1<br />
S(t,1)=S(t-1,2)<br />
S(t,1)=0<br />
H = S(t-1,:) - B<br />
Diff=C*((H(1:nx)+H(2:nx+1))/2)**(n+2)&<br />
*abs((S(t-1,2:nx+1)-S(t-1,1:nx))/dx)**(n-1)<br />
S(t,2:nx)=S(t-1,2:nx)+(dt/dx**2)&<br />
*(Diff(2:nx)*(S(t-1,3:nx+1)-S(t-1,2:nx))&<br />
- Diff(1:nx-1)*(S(t-1,2:nx)-S(t-1,1:nx-1)))&<br />
+MB(2:nx)*dt<br />
where(S(t,:)<B)<br />
S(t,:)=B<br />
end where<br />
if (mod(t,100) .eq.0) then<br />
write (*,*) S(t,:)+1633.42 !height at which g=9.2<br />
end if<br />
end do<br />
end program<br />
</source></div>Adamchttp://websrv.cs.umt.edu/isis/index.php/Kees%27_assignmentKees' assignment2009-08-05T22:04:52Z<p>Adamc: added solution section</p>
<hr />
<div>==Model equation==<br />
:<math>\frac{\partial H}{\partial t} = - \frac{\partial}{\partial x}\left(-D(x) \frac{\partial h}{\partial x}\right) + M</math><br />
<br />
where<br />
<br />
<math>D(x) = C H^{n+2}\left|\frac{\partial h} {\partial x}\right| ^{n-1},</math><br />
<br />
and <br />
<br />
:<math>C = \frac{2 A}{n+2} \left(\rho g\right)^n</math><br />
<br />
==Model parameters==<br />
*<math>\frac{\partial b}{\partial x} = -0.1</math><br />
* <math>M(x) = M_0 - x M_1 = 4.0\times 10^{-3} - 0.2 x</math> km/yr<br />
* <math>\rho</math> = 920 <math>kg/m^3</math><br />
*g=9.8 <math>m/s^2</math><br />
*A = 1e-16 <math> Pa^{-3} a^{-1}</math><br />
* n=3<br />
* dx=1.0 km<br />
<br />
==Boundary conditions==<br />
<br />
* <math>H_l = 0 </math> (left boundary)<br />
<br />
* <math> H_r>0</math> (right boundary)<br />
<br />
==Numerical tips==<br />
<br />
Use a staggered grid such that the <math>D(x_{j+1/2})</math> are computed at the '''centers''' of the grid (as opposed to the vertices, as we have been doing), so<br />
<br />
:<math>D(x_{j+1/2}) = C \left(\frac{H_j + H_{j+1}}{2}\right)^{n+2} \left(\frac{h_{j+1} - h_j}{\Delta x}\right)^{n-1}.</math><br />
<br />
From the diffusivity, the flux is computed<br />
<br />
:<math>\phi_{i+1/2} = D(x_{i+1/2}) \frac{\partial h}{\partial x}</math>,<br />
<br />
where <br />
<br />
:<math>\frac{\partial h}{\partial x} = \frac{h_{i+1}-h_{i}}{\Delta x}</math><br />
<br />
and then the flux (<math>\phi_i</math>) can be used to compute the rate of change of the surface from<br />
<br />
:<math>\frac{\partial H}{\partial t} = -\frac{\partial }{\partial x} H\bar u + M = - \frac{\phi_{i+1/2} - \phi_{i-1/2} }{\Delta x} + M</math><br />
<br />
==Solutions==<br />
* [[Team 1 Solution]]<br />
* [[Team 2 Solution]]<br />
* [[Team 3 Solution]]<br />
* [[Team 4 Solution]]<br />
* [[Team 5 Solution]]<br />
* [[Team 6 Solution]]</div>Adamchttp://websrv.cs.umt.edu/isis/index.php/Kees%27_assignmentKees' assignment2009-08-05T20:13:25Z<p>Adamc: /* Numerical tips */ generalized to arbitrary n</p>
<hr />
<div>==Model equation==<br />
:<math>\frac{\partial H}{\partial t} = - \frac{\partial}{\partial x}\left(-D(x) \frac{\partial h}{\partial x}\right) + M</math><br />
<br />
where<br />
<br />
<math>D(x) = C H^{n+2}\left|\frac{\partial h} {\partial x}\right| ^{n-1},</math><br />
<br />
and <br />
<br />
:<math>C = \frac{2 A}{n+2} \left(\rho g\right)^n</math><br />
<br />
==Model parameters==<br />
*<math>\frac{\partial b}{\partial x} = -0.1</math><br />
* <math>M(x) = M_0 - x M_1 = 4.0\times 10^{-3} - 0.2 x</math> km/yr<br />
* <math>\rho</math> = 920 <math>kg/m^3</math><br />
*g=9.8 <math>m/s^2</math><br />
*A = 1e-16 <math> Pa^{-3} a^{-1}</math><br />
* n=3<br />
<br />
==Boundary conditions==<br />
<br />
* <math>H_l = 0 </math> (left boundary)<br />
<br />
* <math> H_r>0</math> (right boundary)<br />
<br />
==Numerical tips==<br />
<br />
Use a staggered grid such that the <math>D(x_{j+1/2})</math> are computed at the '''centers''' of the grid (as opposed to the vertices, as we have been doing), so<br />
<br />
:<math>D(x_{j+1/2}) = C \left(\frac{H_j + H_{j+1}}{2}\right)^{n+2} \left(\frac{h_{j+1} - h_j}{\Delta x}\right)^{n-1}.</math><br />
<br />
From the diffusivity, the flux is computed<br />
<br />
:<math>\phi_{i+1/2} = D(x_{i+1/2}) \frac{\partial h}{\partial x}</math>,<br />
<br />
where <br />
<br />
:<math>\frac{\partial h}{\partial x} = \frac{h_{i+1}-h_{i}}{\Delta x}</math><br />
<br />
and then the flux (<math>\phi_i</math>) can be used to compute the rate of change of the surface from<br />
<br />
:<math>\frac{\partial H}{\partial t} = -\frac{\partial }{\partial x} H\bar u + M = - \frac{\phi_{i+1/2} - \phi_{i-1/2} }{\Delta x} + M</math></div>Adamchttp://websrv.cs.umt.edu/isis/index.php/Kees%27_assignmentKees' assignment2009-08-05T20:09:48Z<p>Adamc: /* Model equation */ fixed absolute value issue</p>
<hr />
<div>==Model equation==<br />
:<math>\frac{\partial H}{\partial t} = - \frac{\partial}{\partial x}\left(-D(x) \frac{\partial h}{\partial x}\right) + M</math><br />
<br />
where<br />
<br />
<math>D(x) = C H^{n+2}\left|\frac{\partial h} {\partial x}\right| ^{n-1},</math><br />
<br />
and <br />
<br />
:<math>C = \frac{2 A}{n+2} \left(\rho g\right)^n</math><br />
<br />
==Model parameters==<br />
*<math>\frac{\partial b}{\partial x} = -0.1</math><br />
* <math>M(x) = M_0 - x M_1 = 4.0\times 10^{-3} - 0.2 x</math> km/yr<br />
* <math>\rho</math> = 920 <math>kg/m^3</math><br />
*g=9.8 <math>m/s^2</math><br />
*A = 1e-16 <math> Pa^{-3} a^{-1}</math><br />
* n=3<br />
<br />
==Boundary conditions==<br />
<br />
* <math>H_l = 0 </math> (left boundary)<br />
<br />
* <math> H_r>0</math> (right boundary)<br />
<br />
==Numerical tips==<br />
<br />
Use a staggered grid such that the <math>D(x_{j+1/2})</math> are computed at the '''centers''' of the grid (as opposed to the vertices, as we have been doing), so<br />
<br />
:<math>D(x_{j+1/2}) = C \left(\frac{H_j + H_{j+1}^5}{2}\right) \left(\frac{h_{j+1} - h_j}{\Delta x}\right)^2.</math><br />
<br />
From the diffusivity, the flux is computed<br />
<br />
:<math>\phi_{i+1/2} = D(x_{i+1/2}) \frac{\partial h}{\partial x}</math>,<br />
<br />
where <br />
<br />
:<math>\frac{\partial h}{\partial x} = \frac{h_{i+1}-h_{i}}{\Delta x}</math><br />
<br />
and then the flux (<math>\phi_i</math>) can be used to compute the rate of change of the surface from<br />
<br />
:<math>\frac{\partial H}{\partial t} = -\frac{\partial }{\partial x} H\bar u + M = - \frac{\phi_{i+1/2} - \phi_{i-1/2} }{\Delta x} + M</math></div>Adamchttp://websrv.cs.umt.edu/isis/index.php/Kees%27_assignmentKees' assignment2009-08-05T20:08:09Z<p>Adamc: </p>
<hr />
<div>==Model equation==<br />
:<math>\frac{\partial H}{\partial t} = - \frac{\partial}{\partial x}\left(-D(x) \frac{\partial h}{\partial x}\right) + M</math><br />
<br />
where<br />
<br />
<math>D(x) = C H^{n+2}\frac{\partial h} {\partial x} ^{n-1},</math><br />
<br />
and <br />
<br />
:<math>C = \frac{2 A}{n+2} \left(\rho g\right)^n</math><br />
<br />
==Model parameters==<br />
*<math>\frac{\partial b}{\partial x} = -0.1</math><br />
* <math>M(x) = M_0 - x M_1 = 4.0\times 10^{-3} - 0.2 x</math> km/yr<br />
* <math>\rho</math> = 920 <math>kg/m^3</math><br />
*g=9.8 <math>m/s^2</math><br />
*A = 1e-16 <math> Pa^{-3} a^{-1}</math><br />
* n=3<br />
<br />
==Boundary conditions==<br />
<br />
* <math>H_l = 0 </math> (left boundary)<br />
<br />
* <math> H_r>0</math> (right boundary)<br />
<br />
==Numerical tips==<br />
<br />
Use a staggered grid such that the <math>D(x_{j+1/2})</math> are computed at the '''centers''' of the grid (as opposed to the vertices, as we have been doing), so<br />
<br />
:<math>D(x_{j+1/2}) = C \left(\frac{H_j + H_{j+1}^5}{2}\right) \left(\frac{h_{j+1} - h_j}{\Delta x}\right)^2.</math><br />
<br />
From the diffusivity, the flux is computed<br />
<br />
:<math>\phi_{i+1/2} = D(x_{i+1/2}) \frac{\partial h}{\partial x}</math>,<br />
<br />
where <br />
<br />
:<math>\frac{\partial h}{\partial x} = \frac{h_{i+1}-h_{i}}{\Delta x}</math><br />
<br />
and then the flux (<math>\phi_i</math>) can be used to compute the rate of change of the surface from<br />
<br />
:<math>\frac{\partial H}{\partial t} = -\frac{\partial }{\partial x} H\bar u + M = - \frac{\phi_{i+1/2} - \phi_{i-1/2} }{\Delta x} + M</math></div>Adamchttp://websrv.cs.umt.edu/isis/index.php/PDX_afterhoursPDX afterhours2009-08-05T15:10:40Z<p>Adamc: Added Bridge Pedal Event on Sunday</p>
<hr />
<div>==tuesday Aug 4, 2009==<br />
* meet at [http://paccinirestaurant.com/ Paccini's] pub at 7:30<br />
<br />
==wednesday Aug 5, 2009==<br />
* after class leave PSU to go to [http://maps.google.com/maps?hl=en&client=firefox-a&rls=com.ubuntu:en-US:unofficial&hs=zuU&um=1&ie=UTF-8&q=sushi+ichiban+portland&fb=1&split=1&gl=us&view=text&latlng=899833397715679289 Sushi Ichiban]<br />
* possibly head to a pub??<br />
* go to [http://www.groundkontrol.com Ground Kontrol] a retro arcade with beer<br />
* go to [http://www.voodoodoughnut.com Voodoo Doughnut], please someone buy Ian Rutt the $5 doughnut.<br />
* ???<br />
<br />
==thursday Aug 6, 2009==<br />
<br />
==friday Aug 7, 2009==<br />
<br />
==saturday Aug 8, 2009==<br />
<br />
==sunday Aug 9, 2009==<br />
The [http://providence.org/bridgepedal/ Portland Bridge Pedal] is a fun event where you can go on 14, 24, or 37 mi. bike ride over Portland's Bridges. Adam is trying to assemble a group to go on Sunday morning. Please speak with him if you are interested in going to this.</div>Adamchttp://websrv.cs.umt.edu/isis/index.php/Finite_differencing_IFinite differencing I2009-08-05T01:33:36Z<p>Adamc: /* Final program and exercise */</p>
<hr />
<div>==Overview==<br />
<br />
In one dimension, the general form of the [http://en.wikipedia.org/wiki/Convection_diffusion_equation convection diffusion equation] is<br />
<br />
:<math>\frac{\partial u(x,t) }{\partial t} - \frac{\partial}{\partial x} D(x) \frac{\partial}{\partial x} u(x,t) - C(x)\frac{\partial}{\partial x} u(x,t) = S(x,t), </math><br />
<br />
<math>u</math> is a general variable, <math>D</math> is a spatially-varying diffusivity, <math>C</math> is a spatially-varying convection rate, and <math>S</math> is a source term. The second term on the left represents diffusion of a solute or other material property, the third term represent convection.<br />
<br />
This equation can be used to model a wide range of phenomena, including the distribution of temperatures (or energy conservation) in an ice sheet. It also bears similarity to the equations expressing conservation of momentum, and analysis of the numerical solutions to this equation are representative of the analysis of numerous numerical treatments in computational fluid dynamics. For these reasons, we will make convection diffusion our first "model problem", or problem to solve in order to strengthen intuition.<br />
<br />
We will take a step-wise approach to solving this equation, first solving for the diffusion, or parabolic equation, then solving for the convection portion. Finally, we will solve the complete equation. Through this process, we will be looking at the stability of the numerical methods used to solve the equation. <br />
<br />
You should think of this as a starting point for both your learning to program, as well as your learning to solve PDEs with programs.<br />
<br />
==Diffusion and explicit solution==<br />
First, we will solve a simplified version of the equation ''explicitly''. Explicit here refers to the way (or what) the differentiation operators are applied to. In this situation they are directly applied to the solution at the present time step in order to determine the next time. <br />
<br />
[[Image:Explicit_method-stencil.png|left|thumb|500 px|The Stencil for the most common explicit method for the parabolic equation.]]<br />
<br />
To better understand, apply the idea to what is called the parabolic, diffusion, or sometimes heat equation. In terms of convection diffusion this is <math>D(x,t) = 1</math>, <math>C(x,t) = 0</math> and <math>S(x,t) = 0</math>,<br />
<br />
:<math> \frac{\partial u(x,t) }{\partial t} = \frac{\partial ^2 u(x,t)}{\partial x^2},</math><br />
<br />
The finite difference approximation of the equation is<br />
<br />
:<math> \frac{u(x,t+\Delta t) - u(x,t)}{\Delta t} = \frac{u(x+\Delta x,t) - 2u(x,t) + u(x-\Delta x,t)}{\Delta x^2}.</math><br />
<br />
Where both derivative approximations are known from the previous lesson. One is called the 'forward Euler' approximation of the time derivative, and the other is the second order accurate, centered second derivative. <br />
<br />
The equation is then algebraically solved for <math> u(x,t)</math> <br />
<br />
:<math>u(x,t + \Delta t) = u(x,t) + \Delta t \frac{u(x+\Delta x,t) - 2u(x,t) +<br />
u(x-\Delta x,t)}{\Delta x^2}</math><br />
<br />
There you have it, a way to compute the future, using the present. Your first task will be to change this into an algorithm.<br />
<br />
One final note, the stencil at the left makes great sense, and understanding it will help make other algorithms clear. However, to understand it we must modify the notation a little. Make super-scripts <math>n</math> refer to time and subscripts <math>j</math> refer to space. The previous equation becomes<br />
<br />
: <math>u_j^{n+1} = u_j^n + \frac{\Delta t}{\Delta x ^2} \left( u_{j-1}^{n} - 2 u_j^{n} + u_{j+1}^{n} \right).</math><br />
<br />
Make sure you recognize how this corresponds with the diagram of the explicit stencil.<br />
<br />
==Numerical solution==<br />
One can see from the above equation that one way to numerically solve the parabolic equation is to use stencils or operators for computing second derivatives. Once the derivative is computed, finding the solution corresponding to the next time step, <math>u(x,t+\Delta t)</math> is just a matter of multiplying the derivative by <math>\Delta t</math> and adding <math> u(x,t)</math>. <br />
<br />
In psuedocode, the solution looks something like<br />
<br />
<source lang=text><br />
<br />
Initialize variables<br />
<br />
loop t over time:<br />
loop i over space:<br />
u(t,i) = u(t-1,i) + delta_t * (u(t-1,i-1) - 2*u(t-1,i) + u(t-1,i+1) ) / delta_x**2<br />
store solution as needed<br />
end loop over space<br />
end loop over time<br />
<br />
</source><br />
<br />
Note that in psuedo-code the final result is an T by N array (T is time steps, N space points). This is typical, the data structures used to store output is often as complex as the algorithm. Sometimes more. <br />
<br />
You'll need to do this in fortran 90, and you'll find Gethin's [[Pragmatic Programming]] very helpful. Plotting is also a key to understanding simulation output. Again, Gethin's got what you need, although you'll also find the [http://matplotlib.sourceforge.net/ matplotlib] documentation helpful (should be familiar if you've used Matlab).<br />
<br />
===Plotting results===<br />
The above line <br />
<source lang=text><br />
store solution as needed<br />
</source><br />
<br />
is perhaps, frustratingly vague. Gethin's [http://source.ggy.bris.ac.uk/wiki/Fortran1 Fortran1: Fortran for beginners] discusses input/output, IO, so the writing should not be a problem. Even<br />
<br />
<source lang=fortran><br />
write(*,*) u<br />
</source><br />
<br />
will suffice, if you can re-direct program output to a file. Assuming your program is called '''cd_prg''', this is done with<br />
<br />
<source lang=bash><br />
./cd_prg > data<br />
</source><br />
<br />
where the '''>''' operator will replace whatever was previously in the file '''data'''.<br />
<br />
As for reading it into Python and plotting it, the main challenges are reading in the data and animating the time series data. Here is some code to do that<br />
<br />
<source lang=python><br />
#!/usr/bin/env python<br />
<br />
# Import only what is needed<br />
from numpy import loadtxt,shape,linspace<br />
from pylab import plot,show,clf,show,ion<br />
<br />
# Import data file, called 'data'<br />
d=loadtxt('data')<br />
<br />
#Determine how much data came in<br />
dims = shape(d)<br />
<br />
clf() # Clears the screen<br />
ion() # Interactive plot mode, critical for animation<br />
<br />
# x data, note that this must correspond to program's domain<br />
x = linspace(0,1,dims[1]) <br />
<br />
# Initial plot, very Matlab(ish), note return of plot handle that allows plot to<br />
# be altered elsewhere in code.<br />
ph,=plot(x,d[0,:],'k') <br />
ph.figure.show() # matplot lib requires show to be called<br />
<br />
# Loop to plot each time step<br />
for i in range(1,dims[0]):<br />
ph.set_ydata(d[i,:]) # Only update y data (faster than replot)<br />
ph.figure.show()<br />
</source><br />
<br />
==Exercise==<br />
#Using the algorithm for the 'explicit' method, find a numerical solution to this heat conduction problem:<br />
<br />
:<math> \frac{\partial u(x,t) }{\partial t} = \frac{\partial ^2 u(x,t)}{\partial x^2},</math><br />
:<math> u(x,0) = sin(\pi x)</math><br />
:<math> u(0,t) = u(1,t) = 0</math><br />
<br />
Use <math>\Delta x</math> = 0.1, <math>\Delta t</math> = 0.005125, and <math>t_{end}</math> = 1.025. Compare the computed solution to the exact solution <math> u(x,t) = \exp(-\pi^2 t) \sin(\pi x)</math>. Repeat the experiment with <math>\Delta t</math> = 0.006 and <math>t_{end}</math> = 1.026.<br />
* [[Group one, parabolic, explicit]]<br />
* [[Group two, parabolic, explicit]]<br />
* [[Group three, parabolic, explicit]]<br />
* [[Group four, parabolic, explicit]]<br />
* [[Group five, parabolic, explicit]]<br />
* [[Group six, parabolic, explicit]]<br />
<br />
==Convection and numerical stability==<br />
For this unit consider the first-order hyperbolic PDE<br />
:<math> \frac{\partial u}{\partial t} = - v \frac{\partial u}{\partial x}.</math><br />
<br />
Mathematically, this statement is saying that a quantity <math>u(x,t)</math> exists on some grid, and is being carried along by a wind with velocity <math>v</math>. Before applying finite difference operators, clean up the notation so that super-scripts (<math>u^n</math>) refer to time and subscripts (<math>u_j</math>) refer to space.<br />
<br />
Continuing to work with explicit schemes, the machinery of descritization allows us to quickly move to the form<br />
:<math> \frac{u_j^{n+1} - u_j^n}{\Delta t} = -v \left( \frac{u^n_{j+1} - u_{j-1}^n}{2\Delta x} \right ) </math><br />
<br />
and solve to give a recurrence relation<br />
<br />
:<math> u_j^{n+1} =u_j^n - \frac{v \Delta t}{2\Delta x} \left(u^n_{j+1} - u_{j-1}^n \right ) </math><br />
<br />
===von Neumann Stability Analysis===<br />
Before implementing this, consider the stability of the solutions by assuming a very generic form of solution<br />
:<math> u_j^n = A(k)^n e^{ijk}</math>.<br />
<br />
Maybe you can recall a course in differential equations where you spent the better part of a semester making similar substitutions into equations to find solutions? This [[Wikipedia:Euler's formula | complex exponential]] is the Swiss Army knife of functions, and satisfies many equations.<br />
<br />
In our assumed solutions the amplitude is <math>A^n(k)</math> (exponentiated to higher powers with time) and the ''wave number'' is <math>k = \frac{2 \pi}{\lambda}</math>. Said in words, we assume that the solution will be oscillatory (recall <math> e^{ikx} = cos(kx) + i sin(kx)</math>) and that the solution's amplitude will depend on the frequency, or <math>k</math>. In our discrete case <math>j</math> serves as a proxy for space, <math>x</math>.<br />
<br />
:<math> A^{n+1} e^{ijk} = A^{n} e^{ijk} - \frac{v \Delta t}{2\Delta x} \left( A^{n} e^{i(j+1)k} - A^{n} e^{i(j-1)k} \right ) </math><br />
<br />
divide through by <math> A^n e^{ijk}</math><br />
<br />
:<math> A = 1 - \frac{v \Delta t}{2\Delta x} (e^{ik} - e^{-ik}) = 1 - i\frac{v \Delta t}{\Delta x} sin k</math><br />
<br />
''What does that mean?'' If <math>|A|^2>1</math>, the solution grows without bound in time, because each time step applies an higher exponent to <math> A(k)^n</math>. So, this solution scheme is unstable for all time steps, and space steps. Bummer.<br />
<br />
===Getting stability===<br />
Now, let's try that again, and when discretizing do what will become a favorite trick, average or smear the values of the function. A new discretization will be <br />
<br />
:<math> u_j^{n+1} = \frac{1}{2} (u_{j+1}^n + u_{j-1}^n) - \frac{v \Delta t}{2\Delta x} \left(u^n_{j+1} - u_{j-1}^n \right ), </math><br />
<br />
this is called the ''Lax method''. Now consider stability in the same way. Omitting some algebra<br />
:<math> A = cos k - i \frac{v \Delta t}{2\Delta x} sin k, </math><br />
is the amplitude. Requiring <br />
:<math> |A|^2 = cos^2 k + \left(\frac{v \Delta t}{2\Delta x}\right)^2 sin^2 k \leq 1, </math><br />
to avoid unbound growth, yields<br />
:<math> \frac{|v|\Delta t}{\Delta x} \leq 1.</math><br />
<br />
This is called the ''[[Wikipedia:Courant–Friedrichs–Lewy condition |Courant-Friedrichs-Levy]]'' stability criterion. It states that the information on a grid has a velocity of <math> \frac{\Delta x}{\Delta t} </math> and that the velocity in the system can not be exceeded by it (causing the ratio to exceed one). For such a thing to happen would be completely unphysical. Consider what happens when an object exceeds the velocity of waves in the media that carries it, a sonic boom. This is a "numerical boom".<br />
<br />
==Exercises==<br />
#Implement the Lax method for a linear system. Is this method explicit or implicit? Use a 10 unit domain and begin with a height 1.0 square wave between 4.5 <math>\leq x \leq</math> 5.5. Fix the ends at 0. let v be 1.0. Also track the sum of the solution before and after the simulation. End the simulation after 4 seconds. Report the behavior with and without the CFL being satisfied. If the CFL is very small, do things improve. What about when it's just under 1.0. What's going on here. Try subtracting <math>u^j_n</math> from both sides of the discretization and inspect for differences between the original discretization and Lax. See an extra term?<br />
#Try the ''leapfrog method'' for descretization<br />
:<math>u^{n+1}_j = u^{n-1}_j - \frac{v\Delta t}{\Delta x} (u^n_{j+1} - u^n_{j-1}).</math><br />
:Does this improve the numerical diffusion?<br />
<br />
* [[Group2 Convection Exercise]]<br />
<br />
==Final program and exercise==<br />
Bring the descretization schemes for diffusion and convection to solve the convection-diffusion equation explicitly, with finite differences. Consider a non-dimensional form of the equation.<br />
<br />
:<math>\frac{\partial \phi}{\partial t} + u \frac{\partial \phi}{\partial x} - \frac{\partial }{\partial x}\mathrm{Pe}^{-1} \frac{\partial }{\partial x} \phi = q </math><br />
<br />
The <math>\mathrm{Pe}^{-1} </math>, the inverse of the [[Wikipedia:Peclet number|Peclet number]], the ratio of the velocity scale <math>U</math> times the length scale <math>L</math> to the diffusivity <math>D</math>, <br />
:<math> \mathrm{Pe} = \frac{UL}{D}.</math><br />
*On a unit domain, specify <math>\phi_(t,0)</math>= 0, <math>\phi_(t,1)</math>=1, and Pe = 10.<br />
*Compare the solution to <math>c_a(x)</math>, the analytic solution of the equation, is<br />
:<math> c_a(x) = a + (b-a)*\frac{\exp((x-1)\mathrm{Pe}) - \exp(\mathrm{-Pe})}{1-\exp(\mathrm{-Pe})} </math><br />
* Experiment with the Peclet number and mesh resolution to determine how stable your numerical scheme is.<br />
<br />
* [[Group2 Convection + Diffusion Exercise]]</div>Adamc