Index: sm/workshop/2012/Talks/10_PIG_UQ/10_PIG_UQ.tex
===================================================================
--- /issm/workshop/2012/Talks/10_PIG_UQ/10_PIG_UQ.tex	(revision 14038)
+++ 	(revision )
@@ -1,325 +1,0 @@
-%Preamble {{{1
-\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 {{{1
-%%Automatically insert outline slide at the begining of each new section (very nice!)
-%\AtBeginSection[]{
-%\begin{frame}<beamer>
-%	\ghostframe                      %do not count slide
-%	\frametitle{Outline}             % title of the slide
-%	\tableofcontents[currentsection]
-%\end{frame}
-%}
-%}}}
-%Presentation properties{{{1
-%between brackets    []: shown on headers and footers of every slide (must be short)
-%between parenthesis (): shown on title page only (long)
-\title[PIG: UQ Application]{Ice Sheet System model}
-\subtitle{Application to Pine Island Glacier: Uncertainty Quantification (ice thickness, ice rigidity and basal
-friction).}
-\author[Larour et al.]{
-\textbf{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},
-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 %{{{1
-\begin{frame}[plain] %No headers or footers
-	%Do not change, this page is automatically generated
-	\ghostframe
-	\titlepage
-\end{frame} %}}}
-\begin{frame}{Overview}%{{{1
-	\tableofcontents
-\end{frame}%}}}
-
-\section{Introduction}
-
-\begin{frame}{Introduction}%{{{1
-
-	Follow on study on Pine Island Glacier. Apply Uncertainty Quantification tools (based on the Dakota 
-	framework) to quantify the propagation of errors in model inputs through a diagnostic solution of the 
-	2D SSA model (MacAyeal), and assess the sensitivity of model outputs (mass flux at different flux gates
-	across Pine Island Glacier) to model inputs.
-	\vspace{1em} 
-
-	Steps:
-	\small{
-	\begin{itemize}
-		\item  Load 09 PIG Application results (starting from the end of the basal friction inversion)
-		\item  Load ice thickness cross-over errors from IceBridge 2009 WAIS campaign
-		\item  Run sampling analysis  using ice thickness cross-over  and mass flux diagnostics.
-		\item  Run sensitivity analysis  using ice thickness, ice rigidity and basal friction as inputs and mass
-		flux diagnostics
-		\item  Plot results
-	\end{itemize}
-	}                        
-	\vspace{2em}                                           
-
-\end{frame}%}}}
-
-\section{Flux Gates}
-
-\begin{frame}{First Step: plot gates}%{{{1
-
-	The flux gates we will rely on are all in \alert{Exp\_Par/MassFluxes}. There are 13 of them. We will compute 
-	mass fluxes (in Gt/yr) for all of these gates (using the depth-averaged ice velocity, ice thickness and ice 
-	density). 
-
-	Run step 1 of the runme.m file to plot the gates, overlayed over the ice velocity (from our model inversion
-	of the basal friction)
-
-	\lstinputlisting[firstnumber=17,firstline=17,lastline=32,basicstyle=\lstbasicfont\scriptsize,numberstyle=\tiny,tabsize=1]{runme.m}
-
-	\vspace{1em} 
-
-	The gates are positioned across PIG at the inset of tributary glaciers.
-
-\end{frame}
-
-\begin{frame}{Load Cross-Over errors for ice thickness}
-
-	Next, adapt your mesh to SeaRISE velocities and thickness.  The data resides in:
-
-	\lstinputlisting[firstnumber=5,firstline=5,lastline=5,basicstyle=\lstbasicfont\scriptsize,numberstyle=\tiny,tabsize=1]{runme.m}
-
-	Adapt:
-	\small {
-	\begin{itemize}
-
-		\item Fill \alert{md.inversion.vx\_obs}, \alert{md.inversion.vy\_obs}, and \alert{md.inversion.vel\_obs} with interpolated velocities 
-		\item Fill \alert{md.geometry.thickness} with interpolated thickness
-		\item Mesh adapt your model (bamg) 
-			\begin{itemize}
-				\item Use the fields \alert{md.inversion.vel\_obs} and \alert{md.geometry.thickness} 
-				\item Set hmax=400000m and hmin=80000m 
-			\end {itemize}
-
-		\item Set model lat/long using SeaRISE projection information (see
-			projection information in \alert{Greenland\_5km\_dev1.2.nc} - Hint: in matlab you can use ndisp)
-		\item Save your model to a file
-	\end {itemize}
-	}
-
-\end{frame}
-
-\begin{frame}{Mesh, 3/3}{Check Your Work}
-	\vspace{1em}
-	Solution:
-	\lstinputlisting[firstnumber=39,firstline=39,lastline=45,basicstyle=\lstbasicfont\scriptsize,numberstyle=\tiny,tabsize=1]{runme.m}
-
-	\vspace{1em}
-	Plot your mesh.  It should look like:
-
-	\center{
-	\includegraphics[width=0.4\textwidth]{Mesh}
-	}
-
-\end{frame}%}}}
-
-\begin{frame}{Mesh}{ Full Solution}%{{{1
-
-	\lstinputlisting[firstnumber=22,firstline=22,lastline=45,basicstyle=\lstbasicfont\scriptsize,numberstyle=\tiny,tabsize=1]{runme.m}
-
-\end{frame}%}}}
-
-\section{Parameterization}
-
-\begin{frame}{Parameterization}%{{{1
-
-	Call the \alert{setmask} function with empty arguments and then parameterize your mesh with file \alert{Exp\_Par/Greenland.par}.
-
-	\vspace{1em}
-	Read though the parameter file \alert{Exp\_Par/Greenland.par}, similar to your Jks par file.
-	\vspace{1em}
-
-	NB:  No function is called to define boundary conditions.  Instead, set \alert{spc} velocities
-	explicitly to zero, ensuring non-singularity.  
-
-	\lstinputlisting[firstnumber=95,firstline=95,lastline=103,basicstyle=\lstbasicfont\scriptsize,numberstyle=\tiny,tabsize=1]{Exp_Par/Greenland.par} 
-
-\end{frame}%}}}
-\begin{frame}{Parameterization}{Solution}%{{{1
-	\lstinputlisting[firstnumber=52,firstline=52,lastline=56,basicstyle=\lstbasicfont\scriptsize,numberstyle=\tiny,tabsize=1]{runme.m}
-
-	\vspace{1em}
-	Plot thickness and velocity.  
-	\lstinputlisting[basicstyle=\lstbasicfont\scriptsize]{Code/plotparam.m}
-	\vspace{1em}
-	They should look like:
-
-	\center{
-	\includegraphics[width=0.4\textwidth]{Init_thickness}\includegraphics[width=0.4\textwidth]{Init_vel}
-	}
-\end{frame}%}}}
-
-\section{Control Method Solution}
-
-\begin{frame}\frametitle{Diagnostic}%{{{1
-	\vspace{1em}
-	Use control methods to inversely solve for Greenland \alert{FrictionCoefficient}
-	\vspace{1em}
-
-	Steps:
-	\begin{itemize}
-
-		\item Set \alert{cost functions} 
-			\begin{itemize} 
-				\item log of surface velocity
-				\item drag coefficient gradient
-			\end{itemize}  
-		\item Set \alert{cost functions coefficient} to 8*10\textsuperscript{-15}
-		\item Specify \alert{max inversion parameter} = 200, \alert{min inversion parameter} = 1
-		\item Solve a 50-step Diagnostic in 2D, Macayeal
-		\item Copy result \alert{Vx}, \alert{Vy}, and \alert{Friction Coefficient} to model values
-		\item Save your model 
-	\end{itemize}
-	\vspace{1em}
-	Startoff Code:
-	\lstinputlisting[firstnumber=62,firstline=62,lastline=69,basicstyle=\lstbasicfont\scriptsize,numberstyle=\tiny,tabsize=1]{runme.m}
-	NB: Remember that \alert{md.inversion} can be called for help!
-
-\end{frame}%}}}
-
-\begin{frame}{Diagnostic}{Solution}%{{{1
-	\lstinputlisting[firstnumber=71,firstline=71,lastline=91,basicstyle=\lstbasicfont\scriptsize,numberstyle=\tiny,tabsize=1]{runme.m}
-
-\end{frame}%}}}
-
-\section{Transient Solution}
-
-\begin{frame}{Transient}{Forcing}%{{{1
-
-	You are ready to run a transient!
-
-	\vspace{1em}
-	Simulate Greenland warming by forcing a temporal decrease in \alert{md.surfaceforcings.mass\_balance}:
-
-	\includegraphics[width=0.9\textwidth]{Forcing_smb}
-
-	\vspace{1em}
-
-	Specify a transient forcing by adding a time value to the end (in the \alert{end+1} position) of the column of forcing variable values.  For example, let \alert{smb} be the initial values of surface mass balance. To impose the forcing above:
-
-	\lstinputlisting[basicstyle=\lstbasicfont\scriptsize]{Code/smb.m}
-
-	\small{
-	NB: Prior to first and after last imposed time, forcing values remain constant.  Between imposed times, forcings are linearly interpolated.
-	}
-\end{frame}%}}}
-
-\begin{frame}{Transient}{Solve}%{{{1
-	Set up your transient
-
-	\vspace{1em}
-	Steps:
-	\begin{itemize}
-		\item Set zero \alert{basal melting rate}
-		\item Interpolate surface mass balance from SeaRISE dataset
-		\item Impose SeaRise surface mass balance for 10 years then linearly decrease to 1 m/yr by year 20
-		\item Set \alert{time step} to 0.2 and \alert{output frequency} to 1 
-		\item Solve a 20 year Transient in 2D, Macayeal
-		\item Save your model 
-	\end{itemize}
-	\vspace{1em}
-	NB: Save the \alert{IceVolume} in your transient results for plotting later:
-	\lstinputlisting[firstnumber=118,firstline=118,lastline=118,basicstyle=\lstbasicfont\scriptsize,numberstyle=\tiny,tabsize=1]{runme.m}
-\end{frame}%}}}
-
-\begin{frame}{Transient}{Solution Setup}%{{{1
-	\lstinputlisting[firstnumber=98,firstline=98,lastline=126,basicstyle=\lstbasicfont\scriptsize,numberstyle=\tiny,tabsize=1]{runme.m}
-
-\end{frame}%}}}
-
-
-\section{Plot Results}
-
-\begin{frame}{Transient Results}{Plot Plan}%{{{1
-
-	\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 in one row of three subplots. 
-
-	\vspace{1em}
-	They should look like:
-
-	\center{
-	\includegraphics[width=1.0\textwidth]{Transient_plan}
-	}
-\end{frame}
-%}}}
-\begin{frame}{Plot}{Solution}%{{{1
-	\lstinputlisting[firstnumber=134,firstline=134,lastline=142,basicstyle=\lstbasicfont\scriptsize,numberstyle=\tiny,tabsize=1]{runme.m}
-
-	\vspace{1em}
-
-	Next, plot a timeseries of mean surface mass balance, mean velocity, and ice volume. Hint to plot mean surface mass balance results:
-
-	\lstinputlisting[firstnumber=145,firstline=145,lastline=150,basicstyle=\lstbasicfont\scriptsize,numberstyle=\tiny,tabsize=1]{runme.m}
-\end{frame}
-\begin{frame}
-
-	\vspace{1em}
-	Your results should look like:
-
-	\center{
-	\includegraphics[width=0.9\textwidth]{Transient_line}
-	}
-\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/10_PIG_UQ/Makefile
===================================================================
--- /issm/workshop/2012/Talks/10_PIG_UQ/Makefile	(revision 14038)
+++ 	(revision )
@@ -1,36 +1,0 @@
-#Makefile for beamer
-TARGET=10_PIG_UQ
-
-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 -rf ${TARGET}.{dvi,ps,pdf,toc,log,aux,out,nav,snm,bbl,blg,vrb} PIG.* PIG-* Temporary* qmu*
-
-allclean: 
-	make clean
-	make unlinks
Index: sm/workshop/2012/Talks/10_PIG_UQ/mcode.sty
===================================================================
--- /issm/workshop/2012/Talks/10_PIG_UQ/mcode.sty	(revision 14038)
+++ 	(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'.
