Index: ../trunk-jpl/examples/Pig/Pig.par
===================================================================
--- ../trunk-jpl/examples/Pig/Pig.par	(revision 19042)
+++ ../trunk-jpl/examples/Pig/Pig.par	(revision 19043)
@@ -6,7 +6,7 @@
 md.miscellaneous.name='PIG';
 md.mesh.hemisphere='s';
 
-% {{{ NetCdf Loading
+%NetCdf Loading
 disp('   Loading SeaRISE data from NetCDF');
 ncdata='../Data/Antarctica_5km_withshelves_v0.75.nc';
 x1    = ncread(ncdata,'x1');
@@ -46,8 +46,8 @@
 
 y2=(ymax-ny*spacing)+(0:1:ny)'*spacing; 
 y2=double(y2);
-% }}}
-% {{{ Geometry
+
+%Geometry
 disp('   Interpolating surface and ice base');
 md.geometry.base    = InterpFromGridToMesh(x1,y1,topg,md.mesh.x,md.mesh.y,0);
 md.geometry.surface = InterpFromGridToMesh(x1,y1,usrf,md.mesh.x,md.mesh.y,0);
@@ -72,8 +72,8 @@
 pos0=find(md.geometry.thickness<=0);
 md.geometry.thickness(pos0)=1;
 md.geometry.surface=md.geometry.thickness+md.geometry.base;
-% }}} 
-% {{{ Initialization parameters
+
+%Initialization parameters
 disp('   Interpolating temperatures');
 md.initialization.temperature=InterpFromGridToMesh(x1,y1,temp,md.mesh.x,md.mesh.y,0)+273.15+Temp_change;
 clear temp;
@@ -92,13 +92,11 @@
 disp('   Set Pressure');
 md.initialization.pressure=md.materials.rho_ice*md.constants.g*md.geometry.thickness;
 
-% }}}
-
 disp('   Construct ice rheological properties');
 md.materials.rheology_n=3*ones(md.mesh.numberofelements,1);
 md.materials.rheology_B=paterson(md.initialization.temperature);
 
-% {{{ Forcings
+%Forcings
 disp('   Interpolating surface mass balance');
 md.surfaceforcings.mass_balance=InterpFromGridToMesh(x1,y1,smb,md.mesh.x,md.mesh.y,0);
 md.surfaceforcings.mass_balance=md.surfaceforcings.mass_balance*md.materials.rho_water/md.materials.rho_ice;
@@ -107,9 +105,8 @@
 disp('   Set geothermal heat flux');
 md.basalforcings.geothermalflux=InterpFromGridToMesh(x1,y1,gflux,md.mesh.x,md.mesh.y,0);
 clear gflux;
-% }}}
 
-% {{{ Friction and inversion set up
+%Friction and inversion set up
 disp('   Construct basal friction parameters');
 md.friction.coefficient=friction_coefficient*ones(md.mesh.numberofvertices,1);
 md.friction.p=ones(md.mesh.numberofelements,1);
@@ -123,13 +120,10 @@
 md.inversion.vx_obs=vx_obs;
 md.inversion.vy_obs=vy_obs;
 md.inversion.vel_obs=vel_obs;
-% }}}
 
-
 disp('   Set boundary conditions');
 md=SetMarineIceSheetBC(md);
 md.basalforcings.floatingice_melting_rate = zeros(md.mesh.numberofvertices,1);
 md.basalforcings.groundedice_melting_rate = zeros(md.mesh.numberofvertices,1);
 md.thermal.spctemperature     = [md.initialization.temperature;1]; %impose observed temperature on surface
 md.masstransport.spcthickness    = NaN*ones(md.mesh.numberofvertices,1);
-
Index: ../trunk-jpl/examples/Pig/runme.m
===================================================================
--- ../trunk-jpl/examples/Pig/runme.m	(revision 19042)
+++ ../trunk-jpl/examples/Pig/runme.m	(revision 19043)
@@ -3,7 +3,7 @@
 
 %Run Steps
 
-% {{{ Mesh Generation #1
+%Mesh Generation #1
 if any(steps==1)
 
 	md.miscellaneous.name='PIG.Mesh_generation';
@@ -72,9 +72,8 @@
 	% Save model
 	save ./Models/PIG.Mesh_generation md;
 end
-% }}}
 
-% {{{ Masks #2
+%Masks #2
 if any(steps==2) 
 
 	md = loadmodel('./Models/PIG.Mesh_generation');	
@@ -103,9 +102,8 @@
 	% Save model
 	save ./Models/PIG.SetMask md;
 end
-% }}}
 
-% {{{ Parameterization #3
+%Parameterization #3
 if any(steps==3) 
 
 	md = loadmodel('./Models/PIG.SetMask');
@@ -117,9 +115,8 @@
 	% Save model
 	save ./Models/PIG.Parameterization md;
 end
-% }}}
 
-% {{{ Control Method #4
+%Control Method #4
 if any(steps==4)
 
 	md = loadmodel('./Models/PIG.Parameterization');
@@ -162,9 +159,8 @@
 	% Save model
 	save ./Models/PIG.Control_drag md;
 end
-% }}}
 
-% {{{ Plot #5
+%Plot #5
 if any(steps==5)
 
 	md = loadmodel('./Models/PIG.Control_drag');
@@ -181,9 +177,8 @@
 		'caxis#1-2',([1.5,4000]),...
 		'colorbartitle#3','[m]', 'log#1-2',10);
 end
-% }}}
 
-% {{{ HO #6
+%HO #6
 if any(steps==6)
 
 	% Load Model
@@ -199,9 +194,8 @@
 	% Save Model
 
 end
-% }}}
 
-% {{{ Plot #7
+%Plot #7
 if any(steps==7)
 
 	mdHO = loadmodel('./Models/PIG.ModelHO');
@@ -225,4 +219,3 @@
 						'colorbartitle#all','[m/yr]',...
 						'layer#5',1, 'log#1', 10,'log#3', 10,'log#5', 10);
 end
-% }}}
Index: ../trunk-jpl/examples/Pig/PigRegion.m
===================================================================
--- ../trunk-jpl/examples/Pig/PigRegion.m	(revision 19042)
+++ ../trunk-jpl/examples/Pig/PigRegion.m	(revision 19043)
@@ -1,4 +1,4 @@
-% {{{ Getting the velocity in PIG vicinity for the ExpDraw
+%Getting the velocity in PIG vicinity for the ExpDraw
 
 % Load Velocities
 % http://nsidc.org/data/nsidc-0484.html
@@ -50,5 +50,3 @@
 yred=y(id1y:id2y);
 vel_obs=sqrt(vx_obs.^2.+vy_obs.^2.);
 imagesc(xred,yred,vel_obs)
-
-%}}}
Index: ../trunk-jpl/examples/ISMIP/CheatyRunme.m
===================================================================
--- ../trunk-jpl/examples/ISMIP/CheatyRunme.m	(revision 19042)
+++ ../trunk-jpl/examples/ISMIP/CheatyRunme.m	(revision 19043)
@@ -9,7 +9,7 @@
 
 %Run Steps
 
-% {{{ Mesh Generation #1
+%Mesh Generation #1
 if any(steps==1)
 
 	%initialize md as a new model #help model
@@ -30,9 +30,8 @@
 	%->
 	save ./Models/ISMIP.Mesh_generation md;
 end
-% }}}
 
-% {{{ Masks #2
+%Masks #2
 if any(steps==2)
 
 	% load the preceding step #help loadmodel
@@ -50,9 +49,8 @@
 	%->
 	save ./Models/ISMIP.SetMask md;
 end
-% }}}
 
-% {{{ Parameterization #3
+%Parameterization #3
 if any(steps==3)
 
 	% load the preceding step #help loadmodel
@@ -68,9 +66,8 @@
 	%->
 	save ./Models/ISMIP.Parameterization md;
 end
-% }}}
 
-% {{{ Extrusion #4
+%Extrusion #4
 if any(steps==4)
 	
 	% load the preceding step #help loadmodel
@@ -88,9 +85,8 @@
 	%->
 	save ./Models/ISMIP.Extrusion md;
 end
-% }}}
 
-% {{{ Set the flow computing method #5
+%Set the flow computing method #5
 if any(steps==5)
 
 	% load the preceding step #help loadmodel
@@ -105,9 +101,8 @@
 	%->
 	save ./Models/ISMIP.SetFlow md;
 end
-% }}}
 
-% {{{ Set Boundary Conditions #6
+%Set Boundary Conditions #6
 if any(steps==6)
 
 	% load the preceding step #help loadmodel
@@ -160,9 +155,8 @@
 	%->
 	save ./Models/ISMIP.BoundaryCondition md;
 end
-% }}}
 
-% {{{ Solving #7
+%Solving #7
 if any(steps==7)
 	% load the preceding step #help loadmodel
 	% path is given by the organizer with the name of the given step
@@ -187,9 +181,8 @@
 	%->
 	plotmodel(md,'data',md.results.StressbalanceSolution.Vel)
 end
-% }}}
 
-% {{{ Solving #8
+%Solving #8
 if any(steps==8)
 	% load the preceding step #help loadmodel
 	% path is given by the organizer with the name of the given step
@@ -226,4 +219,3 @@
 	%->
 	plotmodel(md,'data',md.results.TransientSolution(20).Vel)
 end
-% }}}
Index: ../trunk-jpl/examples/ISMIP/runme.m
===================================================================
--- ../trunk-jpl/examples/ISMIP/runme.m	(revision 19042)
+++ ../trunk-jpl/examples/ISMIP/runme.m	(revision 19043)
@@ -9,7 +9,7 @@
 
 %Run Steps
 
-% {{{ Mesh Generation #1
+%Mesh Generation #1
 if any(steps==1) 
 
 	%initialize md as a new model #help model
@@ -26,9 +26,8 @@
 	%->
 
 end
-% }}}
 
-% {{{ Masks #2
+%Masks #2
 if any(steps==2) 
 
 	% load the preceding step #help loadmodel
@@ -46,9 +45,8 @@
 	%->
 
 end
-% }}}
 
-% {{{ Parameterization #3
+%Parameterization #3
 if any(steps==3) 
 
 	% load the preceding step #help loadmodel
@@ -64,9 +62,8 @@
 	%->
 
 end
-% }}}
 
-% {{{ Extrusion #4
+%Extrusion #4
 if any(steps==4)
 	
 	% load the preceding step #help loadmodel
@@ -84,10 +81,8 @@
 	%->
 
 end
-% }}}
 
-% {{{ Set the flow computing method #5
-
+%Set the flow computing method #5
 if any(steps==5)
 
 	% load the preceding step #help loadmodel
@@ -102,9 +97,8 @@
 	%->
 
 end
-% }}}
 
-% {{{ Set Boundary Conditions #6
+%Set Boundary Conditions #6
 if any(steps==6)
 
 	% load the preceding step #help loadmodel
@@ -156,9 +150,8 @@
 	%->
 
 end
-% }}}
 
-% {{{ Solving #7
+%Solving #7
 if any(steps==7)
 	% load the preceding step #help loadmodel
 	% path is given by the organizer with the name of the given step
@@ -182,9 +175,8 @@
 	% plot the surface velocities #plotdoc
 	%->
 end
-% }}}
 
-% {{{ Solving #8
+%Solving #8
 if any(steps==8)
 	% load the preceding step #help loadmodel
 	% path is given by the organizer with the name of the given step
@@ -221,4 +213,3 @@
 	%->
 
 end
-% }}}
