Index: sm/workshop/2012/Talks/13_IceBridge/13_IceBridge.tex =================================================================== --- /issm/workshop/2012/Talks/13_IceBridge/13_IceBridge.tex (revision 14042) +++ (revision ) @@ -1,307 +1,0 @@ -%Preamble {{{ -\documentclass[9pt,notes=hide]{beamer} % Document class: beamer %8 pts is also a good option -\usetheme[jpl,headstyle=institute,fnolabel=,sidebar=.15\paperwidth,footstyle=low,footsep={ $\cdot$ },TPomitframeno]{Pasadena} -\useinnertheme{MLHacks} % ISSM bonuses (ghostframes for example) -\usepackage[T1]{fontenc} % T1 encoding is better for foreign langages -\usepackage{aeguill,aecompl} % Virtual fonts -\usepackage{amsmath,amssymb,amsfonts}% Mathematical symbols (useful for equations) -\usepackage{bm} % maths letters, real = IR -\usepackage{xcolor} % We can define new colors -\usepackage{graphicx} % Include pictures -\usepackage{times} % Use Times Font instead of default -\usepackage{appendixnumberbeamer} -\usepackage[framed,numbered,autolinebreaks,useliterate]{mcode} -\usepackage{booktabs} %nice tables (toprule, bottomrule,...) -\graphicspath{ % directories where the pictures are located -{Images/} -} -\hypersetup{pdfstartview={Fit}} % Acrobat reader: open fit to screen - -\setbeamercovered{transparent} %invisible, transparent, dynamic, highly dynamic -%}}} -%Newcommands {{{ -%%Automatically insert outline slide at the begining of each new section (very nice!) -%\AtBeginSection[]{ -%\begin{frame} -% \ghostframe %do not count slide -% \frametitle{Outline} % title of the slide -% \tableofcontents[currentsection] -%\end{frame} -%} -%}}} -%Presentation properties{{{ -%between brackets []: shown on headers and footers of every slide (must be short) -%between parenthesis (): shown on title page only (long) -\title[Application IceBridge]{Ice Sheet System model} -\subtitle{Application to IceBridge dataset, Greenland Ice Sheet} -\author[Larour et al.]{ -Eric \textsc{Larour}\inst{1}, -Eric \textsc{Rignot}\inst{1,2}, -Mathieu \textsc{Morlighem}\inst{2}, -H\'el\`ene \textsc{Seroussi}\inst{1} -Nicole \textsc{Schlegel}\inst{1}, -Chris \textsc{Borstad}\inst{1}, -Feras \textsc{Habbal}\inst{1,2}, -Daria \textsc{Halkides}\inst{1,3}, -Gilberto \textsc{Perez}\inst{1}, -John \textsc{Schiermeier}\inst{1}, -\textbf{Michiel \textsc{Helsen}}\inst{2} -\vspace{1em} -} -\institute[Jet Propulsion Laboratory]{ -\inst{1}Jet Propulsion Laboratory - California Institute of Technology\\ -\inst{2}University of California, Irvine\\ -\inst{3}Joint Institute for Regional Earth System Science \& Engineering, UCLA -} -\conference[ISSM Workshop 2012]{ISSM Workshop 2012} -\date[]{\hspace{-15em}December 2012\hspace{3em}\copyright Copyright 2012. All rights reserved} -\logo{\includegraphics[width=8em]{ISSMlogo}} -%}}} - -\begin{document} - -%Title slide %{{{ -\begin{frame}[plain] %No headers or footers - %Do not change, this page is automatically generated - \ghostframe - \titlepage -\end{frame} %}}} -\begin{frame}{Overview}%{{{ - \tableofcontents -\end{frame}%}}} - -\section{Introduction} - -%Everything that follows is an example talk... Nicole's from last year - -\begin{frame}{Goals}%{{{ - - Use your new ISSM skills to adjust a coarse Greenland model by adding IceBridge data. - - We refine the mesh in the Jakobshavn basin and add higher resolution bedrock and surface - elevation data from IceBridge in this area. - \vspace{1em} - - Steps: - \small{ - \begin{itemize} - \item Refine Greenland mesh using given Jakobshavn outline - \item Parameterize, include the IceBridge/CReSIS bedrock and surface data - \item Plot bedrock and surface data - \item Diagnostic: run 2 inverse method runs to solve for control drag (20 steps recommended) - \item Transient: run 20 year runs, with coarse and refined bedrock and surface elevation data - \item Plot transient results - \end{itemize} - } - \vspace{2em} - -\end{frame}%}}} - -\section{Mesh} - -\begin{frame}{First Run Step: Mesh}%{{{ - We modify the experiment from the Greenland SeaRISE talk, and improve from there. - First, run the first step in \alert{runme.m} in directory \alert{13\_IceBridge} to mesh the Greenland - domain as done in the previous talk, and plot the result. - \vspace{1em} - - Step 1 is interrupted after making the default mesh. - \vspace{1em} - - We want to refine the mesh in the Jakobshavn area. An outline of this area - \alert{Jak\_outline.exp} can be found in directory \alert{Exp\_Par}. - \vspace{1em} - - Try using \alert{exptool} to view this outline: - - \lstinputlisting[basicstyle=\lstbasicfont\scriptsize]{Code/exptool.m} -\end{frame}%}}} - -\begin{frame}{Mesh, 2/6}{Adapt}%{{{ - It should look like this: - \center{ - \includegraphics[width=0.4\textwidth]{13_Mesh1.pdf} - } -\end{frame}%}}} - -\begin{frame}{Mesh, 3/6}{Adapt}%{{{ - Next, modify the bamg command by imposing a 1 km resolution within the Jakobshavn area: use \alert{hmaxVertices}. - \vspace{1em} - - Note that you need to deactivate the previous \alert{bamg} command. - \vspace{1em} - - Check your results using \alert{plotmodel}. -\end{frame}%}}} - -\begin{frame}{Mesh, 4/6}{Check your work}%{{{ - Solution: - - \lstinputlisting[firstnumber=27,firstline=27,lastline=36,basicstyle=\lstbasicfont\scriptsize,numberstyle=\tiny,tabsize=1]{runme.m} -\end{frame}%}}} - -\begin{frame}{Mesh, 5/6}{Adapt}%{{{ - Solution: - \center{ - \includegraphics[width=0.4\textwidth]{13_Mesh2.pdf} - } - \vspace{1em} - - Try \alert{zoom} to make a close-up of the Jakobshavn domain. -\end{frame}%}}} - -\begin{frame}{Mesh, 6/6}{ Full Solution}%{{{ - - \lstinputlisting[firstnumber=22,firstline=22,lastline=42,basicstyle=\lstbasicfont\scriptsize,numberstyle=\tiny,tabsize=1]{runme.m} - -\end{frame}%}}} - -\section{Parameterization} - -\begin{frame}{Parameterization}%{{{ - - We want to include high-resolution bedrock and surface elevation data acquired in the IceBridge - mission. Accessible on: - \vspace{1em} - - https://data.cresis.ku.edu/data/grids/ - \alert{Jakobshavn\_2008\_2011\_Composite\_XYZGrid.txt} - \vspace{1em} - - \includegraphics[width=0.6\textwidth]{Jakobshavn_2008_2011_Composite_Bottom_Preview.png} -\end{frame}%}}} - -\begin{frame}{Parameterization, 2/6}%{{{ - Bedrock data is read, transformed into a usable grid, and interpolated to the mesh in the parameter - file \alert{Exp\_Par/Greenland.par}: - \lstinputlisting[firstnumber=30,firstline=30,lastline=47,basicstyle=\lstbasicfont\scriptsize,numberstyle=\tiny,tabsize=1]{Exp_Par/Greenland.par} - - Modify \alert{Greenland.par} such that the surface elevation data is also included for the - Jakobshavn area. - -\end{frame}%}}} - -\begin{frame}{Parameterization, 3/6}{Solution}%{{{ - \lstinputlisting[firstnumber=30,firstline=30,lastline=52,basicstyle=\lstbasicfont\scriptsize,numberstyle=\tiny,tabsize=1]{Exp_Par/Greenland_solution.par} - Plot surface elevation, thickness and bedrock. -\end{frame}%}}} - -\begin{frame}{Parameterization, 4/6}{Solution}%{{{ - \lstinputlisting[basicstyle=\lstbasicfont\scriptsize]{Code/plotparam.m} - They should look like: - - \center{ - \includegraphics[width=0.3\textwidth]{13_init_surface.pdf}\includegraphics[width=0.3\textwidth]{13_init_thickness.pdf}\includegraphics[width=0.3\textwidth]{13_init_bed.pdf} - } -\end{frame}%}}} - -\begin{frame}{Parameterization, 5/6}{Solution}%{{{ - To plot the difference in bed topography between the SeaRISE and IceBridge datasets: - modify the parameterization step in your runme and save the model under a different name. - - a difference in fields can be plotted using: - \lstinputlisting[basicstyle=\lstbasicfont\scriptsize]{Code/plotdifbed.m} - - -\end{frame}%}}} - -\begin{frame}{Parameterization, 6/6}{Solution}%{{{ - - a difference in fields can be plotted using: - \lstinputlisting[basicstyle=\lstbasicfont\scriptsize]{Code/plotdifbed.m} - - \includegraphics[width=0.4\textwidth]{13_Jakdifbed.pdf} - -\end{frame}%}}} - -\section{Control Method Solution} - -\begin{frame}\frametitle{Diagnostic}%{{{ - \vspace{1em} - Use control methods to inversely solve for Greenland friction coefficient - \vspace{1em} - Comparable to 12\_SeaRISE experiment - - \includegraphics[width=0.4\textwidth]{13_vel_obs.pdf} - - The observed velocity map contains some gaps; exclude these from the inversion by creating - a new exp file that outlines the gaps in the velocity data, using \alert{exptool}: - \lstinputlisting[basicstyle=\lstbasicfont\scriptsize]{Code/exptool2.m} -\end{frame}%}}} - -\begin{frame}\frametitle{Diagnostic, 2/3}%{{{ - \lstinputlisting[firstnumber=78,firstline=78,lastline=80,basicstyle=\lstbasicfont\scriptsize,numberstyle=\tiny,tabsize=1]{runme.m} - - Next, run the model and plot the resulting friction coefficient and velocity pattern - - \lstinputlisting[basicstyle=\lstbasicfont\scriptsize]{Code/plotinversion.m} - -\end{frame}%}}} - -\begin{frame}\frametitle{Diagnostic, 3/3}%{{{ - \lstinputlisting[basicstyle=\lstbasicfont\scriptsize]{Code/plotinversion.m} - \includegraphics[width=0.4\textwidth]{13_FrictCoef.pdf}\includegraphics[width=0.4\textwidth]{13_Vel.pdf} - -\end{frame}%}}} - -\section{Transient Solution} - -\begin{frame}{Transient}{Forcing}%{{{ - - Do a transient run.. - - \vspace{1em} - -% \includegraphics[width=0.9\textwidth]{Forcing_smb} - -\end{frame}%}}} - -\section{Plot Results} - -\begin{frame}{Transient Results}{Plot Plan}%{{{ - - \vspace{1em} - - Your results are located in \alert{md.results.TransientSolution}. Plot your results. - - \vspace{1em} - First, plot the initial plan view of velocity, surface mass balance, thickness, and surface. - - They should look like: - - \includegraphics[width=0.6\textwidth]{13_Transient_plan.pdf} -\end{frame} -%}}} - -\begin{frame}{Transient Results}{Plot time series}%{{{ - They should look like: - \includegraphics[width=0.8\textwidth]{13_time_series.pdf} -\end{frame} -%}}} - -\begin{frame}{Transient Results}{}%{{{ - Some suggestions what to explore further: - \begin{itemize} - \item How would you make a plot of time series of results from the SeaRISE and - IceBridge experiments? - \item How would you make a plot of the difference between final and initial ice thickness? - \item ... - \end{itemize} - We can help you to implement your own ideas in the code. -\end{frame} -%}}} - -%Thanks slide{{{ -\usebackgroundtemplate{\includegraphics[width=\paperwidth,height=\paperheight]{Thanksbg}} -\begin{frame}[plain] - \ghostframe - \begin{center} - \vspace*{-0.1\paperheight} - \hspace*{-0.15\paperwidth}\Huge{Thanks!} - \end{center} - \note{Thank you very much for your attention} -\end{frame} -\usebackgroundtemplate{} -%}}} - -\end{document} Index: sm/workshop/2012/Talks/13_IceBridge/Makefile =================================================================== --- /issm/workshop/2012/Talks/13_IceBridge/Makefile (revision 14042) +++ (revision ) @@ -1,36 +1,0 @@ -#Makefile for beamer -TARGET=13_IceBridge - -export TEXINPUTS=$(JPL_SVN)/publications/templates/beamer/Theme/Pasadena//: -all: one - -one: - pdflatex -halt-on-error -file-line-error $(TARGET).tex - open $(TARGET).pdf - -complete: - pdflatex -halt-on-error -file-line-error -draftmode $(TARGET).tex - #bibtex $(TARGET) - pdflatex -halt-on-error -file-line-error -draftmode $(TARGET).tex - pdflatex -halt-on-error -file-line-error $(TARGET).tex - open $(TARGET).pdf - -bib: - bibtex $(TARGET) - -links: - ln -s $(JPL_SVN)/publications/templates/beamer/Theme/Pasadena/*.sty . - ln -s $(JPL_SVN)/publications/templates/beamer/Theme/Pasadena/general/*.sty . - ln -s $(JPL_SVN)/publications/templates/beamer/Theme/Pasadena/subtheme/*.sty . - ln -s $(JPL_SVN)/publications/templates/beamer/Theme/Pasadena/ML/*.sty . - ln -s $(JPL_SVN)/publications/templates/beamer/Theme/Pasadena/graphics/JPL/*.pdf . - -unlinks: - find . -maxdepth 1 -type l -print | xargs rm - -clean: - rm -f ${TARGET}.{dvi,ps,pdf,toc,log,aux,out,nav,snm,bbl,blg,vrb} - -allclean: - make clean - make unlinks Index: sm/workshop/2012/Talks/13_IceBridge/mcode.sty =================================================================== --- /issm/workshop/2012/Talks/13_IceBridge/mcode.sty (revision 14042) +++ (revision ) @@ -1,280 +1,0 @@ -%% -%% This is file `mcode.sty' -%% -%% It is supposed to help you easily include MATLAB source code -%% into LaTeX document, but have it nicely highlighted, using -%% the great listings package. -%% -%% PLEASE NOTE that this package does nothing but save you from -%% figuring out some configurations in setting up the LISTINGS -%% package. ALL the work is done by that package! Thus, please -%% refer your questions to the listings package documentation. -%% -%% Usage: You have three ways of including your MATLAB code. As -%% environment, as inline object and directly from an external -%% file. -%% -%% 1) Environment: -%% -%% \begin{lstlisting} -%% YOUR CODE HERE -%% \end{lstlisting} -%% -%% -%% 2) Inline object: -%% -%% Bla bla \mcode{CODEFRAGMENT} bla bla. -%% -%% -%% 3) Include external file (in environment form) -%% -%% \lstinputlisting{YOUR-FILE.m} -%% -%% -%% For your convenience this package has the following options: -%% -%% - bw if you intend to print the document (highlighting done -%% via text formatting (bold, italic) and shades of gray) -%% -%% - numbered if you want line numbers -%% -%% - autolinebreaks if you want the package to automatically -%% wrap your code. This is buggy as it may well break -%% break syntax and it doesn't work well with comments. -%% You REALLY should wrap your code manually. -%% -%% - useliterate if you want some characters / relations in -%% your code to be replace with something more readable. -%% Example: ~= becomes $\neq$, >= becomes $\geq$, delta -%% becomes $\delta$ and so on. -%% -%% - framed if you want a frame around the source code blocks -%% -%% - final if you have ``gloablly'' set the draft option, the -%% listings package will not output the code at all. to -%% force it to do so anyway, load this package with the -%% final option (passes the ``final'' on to listings). -%% -%% For example, you may use \usepackage[numbered,framed]{mcode} -%% in your document preamble. -%% -%% Note: Inside code blocks you can escape to LaTeX text mode -%% using §...§. For ex. §text and some math: $x^2$§, which is -%% especially useful in comments for putting nicely typeset -%% equations etc. To get the same colour/style as in the rest -%% of the comment use \mcommentfont, i.e. §\mcommentfont $x^2$§ -%% -%% To change the font used, edit the first line in the "custo- -%% mise below" section. And feel free to edit other things as -%% well. Refer to the documentation of the listings package to -%% see what else you could do. If an extra small font is re- -%% quired, use {\fontfamily{pcr}\fontsize{3}{4.6}\selectfont} -%% in the definition of \lstbasicfont. -%% -%% Author: -%% Florian Knorn | florian@knorn.org | www.florian-knorn.com -%% -%% Version history: -%% 2.2 -- Bugfix (thanks Willi Gerbig!) -%% 2.1 -- Finally automatic detection between end and end -%% 2.0 -- New options for line breaking and literate prog. -%% 1.8 -- Fixed typo in documentation regarding §...§ -%% 1.7 -- Added MATLAB block comment syntax %{ ...... %} -%% 1.6 -- Added some infos, dealing with keyword ``end'' -%% 1.5 -- Tweaked check to see wether textcomp is loaded -%% 1.4 -- Fixed misconfig (mathescape now set to false) -%% 1.3 -- Purely cosmetic (tabs replaced by spaces) -%% 1.2 -- Added \lstset{showstringspaces=false} -%% 1.1 -- Added \mcode command and [final] option -%% 1.0 -- Release - - -%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% -% D O N ' T T O U C H T H I S % -%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% -\def\fileversion{2.2} -\def\filedate{2011/08/26} - -\typeout{-- Package: `mcode' \fileversion\space <\filedate> --} -\NeedsTeXFormat{LaTeX2e} -\ProvidesPackage{mcode}[\filedate\space\fileversion] - -% for bw-option -\newif\ifbw -\DeclareOption{bw}{\bwtrue} - -% numbered option -\newif\ifnumbered -\DeclareOption{numbered}{\numberedtrue} - -% final option -\newif\iffinal -\DeclareOption{final}{\finaltrue} - -% autolinebreaks option -\newif\ifautolinebreaks -\DeclareOption{autolinebreaks}{\autolinebreakstrue} - -% literate programming (replace certain characters/relations -\newif\ifuseliterate -\DeclareOption{useliterate}{\useliteratetrue} - -% framed option -\newif\ifframed -\DeclareOption{framed}{\framedtrue} - -\DeclareOption*{% default - \PackageWarning{mcode}{Unknown option `\CurrentOption' !}% -} -\ProcessOptions - -\ifbw\typeout{ - settings optimized for printing (bw formating)} -\else\typeout{ - settings optimized for display (colour formating)}\fi -\ifnumbered\typeout{ - line numbering enabled}\else\fi -\ifuseliterate\typeout{ - literate programming (character replacements) enabled}\else\fi -\ifautolinebreaks\typeout{ - automatic line breaking enabled (careful, buggy!)}\else\fi -\ifframed\typeout{ - framed listings}\else\fi - -% This command allows you to typeset syntax highlighted Matlab -% code ``inline''. The font size \small seems to look best... -\newcommand{\mcode}[1]{\lstinline[basicstyle=\lstbasicfont\small]|#1|} - -% check if color command exists -\ifx\color\undefined% - \RequirePackage{xcolor}% -\fi - -% check if listings has been loaded -\ifx\lstset\undefined% - \iffinal - \RequirePackage[final]{listings} - \else - \RequirePackage{listings} - \fi -\fi - -% Check if textcomp has been loaded (this package is needed for -% upright quotes '' (instead of typographic ones `´)... -\ifx\textquotesingle\undefined% - \RequirePackage{textcomp}% -\fi - -%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% -% C U S T O M I S E B E L O W % -%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% - -% --------------------------------------------------------------------------------- -% default font -\def\lstbasicfont{\fontfamily{pcr}\selectfont\footnotesize} - -% --------------------------------------------------------------------------------- -% matlat languate definition -\lstdefinelanguage{matlabfloz}{% - alsoletter={...},% - morekeywords={% % keywords - break,case,catch,continue,elseif,else,end,for,function,global,% - if,otherwise,persistent,return,switch,try,while,...},% - comment=[l]\%, % comments - morecomment=[l]..., % comments - morecomment=[s]{\%\{}{\%\}}, % block comments - morestring=[m]' % strings -}[keywords,comments,strings]% - -% --------------------------------------------------------------------------------- -% general definitions -\lstset{% - basicstyle={\lstbasicfont}, % set font - showstringspaces=false, % do not emphasize spaces in strings - tabsize=3, % number of spaces of a TAB - mathescape=false,escapechar=§, % escape to latex with §...§ - upquote=true, % upright quotes - aboveskip={1.5\baselineskip}, % a bit of space above listings - columns=fixed % nice spacing -} - -% --------------------------------------------------------------------------------- -% define colours and styles -\ifbw % use font formating and gray 'colors' - \def\mcommentfont{\color[gray]{.75}\itshape} %comments light gray and italic - \lstset{language=matlabfloz, % use our version of highlighting - keywordstyle=\bfseries, % keywords in bold - commentstyle=\mcommentfont, % comments - stringstyle=\color[gray]{0.5} % strings darker gray - } -\else% notbw => use colors : ) - \def\mcommentfont{\color[rgb]{.133,.545,.133}} %comments in green - \lstset{language=matlabfloz, % use our version of highlighting - keywordstyle=\color[rgb]{0,0,1}, % keywords in blue - commentstyle=\mcommentfont, % comments - stringstyle=\color[rgb]{.627,.126,.941} % strings in purple - } -\fi%bw - -% --------------------------------------------------------------------------------- -% automatic line breaking --- warning, this is buggy and -% doesn't break comments correctly! -\ifautolinebreaks - \newsavebox{\lbreakdots}\sbox{\lbreakdots}{\lstbasicfont\mcommentfont...} - \lstset{breaklines=true,breakatwhitespace=true,prebreak=\usebox{\lbreakdots}} -\fi - -% --------------------------------------------------------------------------------- -% literate replacements -% the following is for replacing some matlab relations like >= or ~= -% by the corresponding LaTeX symbols, which are much easier to read ... -\ifuseliterate - \lstset{% - literate=% - {~}{{$\neg$}}1 % \neg - {<=}{{\tiny$\leq$}}1 % \leq - {>=}{{\tiny$\geq$}}1 % \geq - {~=}{{\tiny$\neq$}}1 % \neq - {delta}{{\tiny$\Delta$}}1 % \Delta - {(end)}{\lstbasicfont (end)}{5} % black ``end'' when indexing last vector element - {({ }end)}{\lstbasicfont ({ }end)}{6} - {(end{ })}{\lstbasicfont (end{ })}{6} - {({ }end{ })}{\lstbasicfont ({ }end{ })}{7} - {:end}{\lstbasicfont :end}{4} - {:{ }end}{\lstbasicfont :{ }end}{5} - {end:}{\lstbasicfont end:}{4} - {end{ }:}{\lstbasicfont end{ }:}{5} - {,end}{\lstbasicfont ,end}{4} - {,{ }end}{\lstbasicfont ,{ }end}{5} - } -\else - \lstset{% - literate=% - {(end)}{\lstbasicfont (end)}{5} % black ``end'' when indexing last vector element - {({ }end)}{\lstbasicfont ({ }end)}{6} - {(end{ })}{\lstbasicfont (end{ })}{6} - {({ }end{ })}{\lstbasicfont ({ }end{ })}{7} - {:end}{\lstbasicfont :end}{4} - {:{ }end}{\lstbasicfont :{ }end}{5} - {end:}{\lstbasicfont end:}{4} - {end{ }:}{\lstbasicfont end{ }:}{5} - {,end}{\lstbasicfont ,end}{4} - {,{ }end}{\lstbasicfont ,{ }end}{5} - } -\fi%literates - -% --------------------------------------------------------------------------------- -% line numbering -\ifnumbered% numbered option - \lstset{% - numbersep=3mm, numbers=left, numberstyle=\tiny, % number style - } -\fi - -\ifframed% framed option - \lstset{% - frame=single, % frame - } - \ifnumbered% - \lstset{% - framexleftmargin=6mm, xleftmargin=6mm % tweak margins - } - \fi -\fi - -\endinput -%% End of file `mcode.sty'.