Index: /issm/trunk/test/Verification/test01_IceShelfIceFrontM2d/DomainOutline.exp
===================================================================
--- /issm/trunk/test/Verification/test01_IceShelfIceFrontM2d/DomainOutline.exp	(revision 2237)
+++ /issm/trunk/test/Verification/test01_IceShelfIceFrontM2d/DomainOutline.exp	(revision 2237)
@@ -0,0 +1,10 @@
+## Name:domainoutline
+## Icon:0
+# Points Count  Value
+5 1.
+# X pos Y pos
+0 0
+1000000 0
+1000000 1000000
+0 1000000
+0 0
Index: /issm/trunk/test/Verification/test01_IceShelfIceFrontM2d/Front.exp
===================================================================
--- /issm/trunk/test/Verification/test01_IceShelfIceFrontM2d/Front.exp	(revision 2237)
+++ /issm/trunk/test/Verification/test01_IceShelfIceFrontM2d/Front.exp	(revision 2237)
@@ -0,0 +1,10 @@
+## Name:icefront
+## Icon:0
+# Points Count  Value
+5 1.
+# X pos Y pos
+-1000 900000
+-1000 1100000
+1100000 1100000
+1100000 900000
+-1000 900000
Index: /issm/trunk/test/Verification/test01_IceShelfIceFrontM2d/Square.par
===================================================================
--- /issm/trunk/test/Verification/test01_IceShelfIceFrontM2d/Square.par	(revision 2237)
+++ /issm/trunk/test/Verification/test01_IceShelfIceFrontM2d/Square.par	(revision 2237)
@@ -0,0 +1,37 @@
+%Start defining model parameters here
+
+%dynamics
+md.debug=0;
+md.dt=1; %1 year
+md.ndt=md.dt*10; 
+md.artificial_diffusivity=1;
+
+hmin=300;
+hmax=1000;
+ymin=min(md.y);
+ymax=max(md.y);
+md.thickness=hmax+(hmin-hmax)*(md.y-ymin)/(ymax-ymin);
+md.bed=-md.rho_ice/md.rho_water*md.thickness;
+md.surface=md.bed+md.thickness;
+
+md.drag_type=2; %0 none 1 plastic 2 viscous
+md.drag=200*ones(md.numberofgrids,1); %q=1.
+%Take care of iceshelves: no basal drag
+pos=find(md.elementoniceshelf);
+md.drag(md.elements(pos,:))=0;
+md.p=ones(md.numberofelements,1);
+md.q=ones(md.numberofelements,1);
+md.viscosity_overshoot=0.3;
+
+md.observed_temperature=(273-20)*ones(md.numberofgrids,1);
+
+md.B=paterson(md.observed_temperature);
+md.n=3*ones(md.numberofelements,1);
+
+%Deal with boundary conditions:
+md=SetIceShelfBC(md,'Front.exp');
+
+%Parallel options
+md.np=3;
+md.time=50;
+md.waitonlock=30;
Index: /issm/trunk/test/Verification/test01_IceShelfIceFrontM2d/configuration.m
===================================================================
--- /issm/trunk/test/Verification/test01_IceShelfIceFrontM2d/configuration.m	(revision 2237)
+++ /issm/trunk/test/Verification/test01_IceShelfIceFrontM2d/configuration.m	(revision 2237)
@@ -0,0 +1,31 @@
+%name
+testname='test1_IceShelfIceFrontM2d';
+
+%meshing
+resolution=50000;
+domainname='DomainOutline.exp';
+riftname='';
+rifts=0;
+
+%geography
+iceshelfname='all';
+icesheetname='';
+
+%parameterization
+parametername='Square.par';
+			
+%elements type
+elementstype='''macayeal'',''all''';
+
+%extrusion
+numlayers=0;
+extrusionexponent=NaN;
+
+%sequences     analysis    sub_analysis  qmu   control  control_fit parallel
+sequences={  {'diagnostic',  'none',      0 ,     0,      'nan'        0    };...
+             {'prognostic',  'none',      0 ,     0,      'nan'        0    };...
+             {'transient',   'none',      0 ,     0,      'nan'        0    };...
+				 {'diagnostic',  'none',      0 ,     0,      'nan'        1    };...
+             {'prognostic',  'none',      0 ,     0,      'nan'        1    };...
+             {'transient',   'none',      0 ,     0,      'nan'        1    };...
+	};
Index: /issm/trunk/test/Verification/test01_IceShelfIceFrontM2d/testpresolve.m
===================================================================
--- /issm/trunk/test/Verification/test01_IceShelfIceFrontM2d/testpresolve.m	(revision 2237)
+++ /issm/trunk/test/Verification/test01_IceShelfIceFrontM2d/testpresolve.m	(revision 2237)
@@ -0,0 +1,3 @@
+if md.numberofgrids==340
+	load Velocities; md.vx=0.5*vx; md.vy=0.5*vy;
+end
Index: /issm/trunk/test/Verification/test02_IceShelfIceFrontM3d/DomainOutline.exp
===================================================================
--- /issm/trunk/test/Verification/test02_IceShelfIceFrontM3d/DomainOutline.exp	(revision 2237)
+++ /issm/trunk/test/Verification/test02_IceShelfIceFrontM3d/DomainOutline.exp	(revision 2237)
@@ -0,0 +1,10 @@
+## Name:domainoutline
+## Icon:0
+# Points Count  Value
+5 1.
+# X pos Y pos
+0 0
+1000000 0
+1000000 1000000
+0 1000000
+0 0
Index: /issm/trunk/test/Verification/test02_IceShelfIceFrontM3d/Front.exp
===================================================================
--- /issm/trunk/test/Verification/test02_IceShelfIceFrontM3d/Front.exp	(revision 2237)
+++ /issm/trunk/test/Verification/test02_IceShelfIceFrontM3d/Front.exp	(revision 2237)
@@ -0,0 +1,10 @@
+## Name:icefront
+## Icon:0
+# Points Count  Value
+5 1.
+# X pos Y pos
+-1000 900000
+-1000 1100000
+1100000 1100000
+1100000 900000
+-1000 900000
Index: /issm/trunk/test/Verification/test02_IceShelfIceFrontM3d/Square.par
===================================================================
--- /issm/trunk/test/Verification/test02_IceShelfIceFrontM3d/Square.par	(revision 2237)
+++ /issm/trunk/test/Verification/test02_IceShelfIceFrontM3d/Square.par	(revision 2237)
@@ -0,0 +1,38 @@
+%Start defining model parameters here
+
+%dynamics
+md.debug=0;
+md.dt=0.1; %1 year
+md.ndt=md.dt*3; 
+md.artificial_diffusivity=1;
+md.eps_abs=5000;
+
+hmin=300;
+hmax=1000;
+ymin=min(md.y);
+ymax=max(md.y);
+md.thickness=hmax+(hmin-hmax)*(md.y-ymin)/(ymax-ymin);
+md.bed=-md.rho_ice/md.rho_water*md.thickness;
+md.surface=md.bed+md.thickness;
+
+md.drag_type=2; %0 none 1 plastic 2 viscous
+md.drag=200*ones(md.numberofgrids,1); %q=1.
+%Take care of iceshelves: no basal drag
+pos=find(md.elementoniceshelf);
+md.drag(md.elements(pos,:))=0;
+md.p=ones(md.numberofelements,1);
+md.q=ones(md.numberofelements,1);
+md.viscosity_overshoot=0.0;
+
+md.observed_temperature=(273-20)*ones(md.numberofgrids,1);
+
+md.B=paterson(md.observed_temperature);
+md.n=3*ones(md.numberofelements,1);
+
+%Deal with boundary conditions:
+md=SetMarineIceSheetBC(md,'Front.exp');
+
+%Parallel options
+md.np=8;
+md.time=50;
+md.waitonlock=30;
Index: /issm/trunk/test/Verification/test02_IceShelfIceFrontM3d/configuration.m
===================================================================
--- /issm/trunk/test/Verification/test02_IceShelfIceFrontM3d/configuration.m	(revision 2237)
+++ /issm/trunk/test/Verification/test02_IceShelfIceFrontM3d/configuration.m	(revision 2237)
@@ -0,0 +1,37 @@
+%name
+testname='test2_IceShelfIceFrontM3d';
+
+%meshing
+resolution=100000;
+domainname='DomainOutline.exp';
+riftname='';
+rifts=0;
+
+%geography
+iceshelfname='all';
+icesheetname='';
+
+%parameterization
+parametername='Square.par';
+			
+%elements type
+elementstype='''macayeal'',''all''';
+
+%extrusion
+numlayers=8;
+extrusionexponent=4;
+
+%sequences     analysis    sub_analysis  qmu   control  control_fit parallel
+sequences={  {'diagnostic',  'none',      0 ,     0,      'nan'        0    };...
+             {'steadystate', 'none',      0 ,     0,      'nan'        0    };...
+             {'prognostic',  'none',      0 ,     0,      'nan'        0    };...
+             {'thermal',     'steady',    0 ,     0,      'nan'        0    };...
+             {'thermal',     'transient', 0 ,     0,      'nan'        0    };...
+             {'transient',   'none',      0 ,     0,      'nan'        0    };...
+				 {'diagnostic',  'none',      0 ,     0,      'nan'        1    };...
+             {'steadystate', 'none',      0 ,     0,      'nan'        1    };...
+             {'prognostic',  'none',      0 ,     0,      'nan'        1    };...
+             {'thermal',     'steady',    0 ,     0,      'nan'        1    };...
+             {'thermal',     'transient', 0 ,     0,      'nan'        1    };...
+             {'transient',   'none',      0 ,     0,      'nan'        1    };...
+	};
Index: /issm/trunk/test/Verification/test02_IceShelfIceFrontM3d/testpresolve.m
===================================================================
--- /issm/trunk/test/Verification/test02_IceShelfIceFrontM3d/testpresolve.m	(revision 2237)
+++ /issm/trunk/test/Verification/test02_IceShelfIceFrontM3d/testpresolve.m	(revision 2237)
@@ -0,0 +1,3 @@
+if md.numberofgrids==728,
+	load Velocities; md.vx=0.8*vx; md.vy=0.8*vy; md.vz=0.8*vz; md.temperature=temperature-1; md.pressure=pressure;
+end
Index: /issm/trunk/test/Verification/test03_IceShelfIceFrontMP3d/DomainOutline.exp
===================================================================
--- /issm/trunk/test/Verification/test03_IceShelfIceFrontMP3d/DomainOutline.exp	(revision 2237)
+++ /issm/trunk/test/Verification/test03_IceShelfIceFrontMP3d/DomainOutline.exp	(revision 2237)
@@ -0,0 +1,10 @@
+## Name:domainoutline
+## Icon:0
+# Points Count  Value
+5 1.
+# X pos Y pos
+0 0
+1000000 0
+1000000 1000000
+0 1000000
+0 0
Index: /issm/trunk/test/Verification/test03_IceShelfIceFrontMP3d/Front.exp
===================================================================
--- /issm/trunk/test/Verification/test03_IceShelfIceFrontMP3d/Front.exp	(revision 2237)
+++ /issm/trunk/test/Verification/test03_IceShelfIceFrontMP3d/Front.exp	(revision 2237)
@@ -0,0 +1,10 @@
+## Name:icefront
+## Icon:0
+# Points Count  Value
+5 1.
+# X pos Y pos
+-1000 900000
+-1000 1100000
+1100000 1100000
+1100000 900000
+-1000 900000
Index: /issm/trunk/test/Verification/test03_IceShelfIceFrontMP3d/Pattyn.exp
===================================================================
--- /issm/trunk/test/Verification/test03_IceShelfIceFrontMP3d/Pattyn.exp	(revision 2237)
+++ /issm/trunk/test/Verification/test03_IceShelfIceFrontMP3d/Pattyn.exp	(revision 2237)
@@ -0,0 +1,10 @@
+## Name:pattyndomain
+## Icon:0
+# Points Count  Value
+5 1.
+# X pos Y pos
+0 0
+500001 0
+500001 1000001
+0 1000001
+0 0
Index: /issm/trunk/test/Verification/test03_IceShelfIceFrontMP3d/Square.par
===================================================================
--- /issm/trunk/test/Verification/test03_IceShelfIceFrontMP3d/Square.par	(revision 2237)
+++ /issm/trunk/test/Verification/test03_IceShelfIceFrontMP3d/Square.par	(revision 2237)
@@ -0,0 +1,37 @@
+%Start defining model parameters here
+
+%dynamics
+md.debug=0;
+md.dt=0.1; %1 year
+md.ndt=md.dt*3; 
+md.artificial_diffusivity=1;
+md.eps_abs=5000;
+
+hmin=300;
+hmax=1000;
+ymin=min(md.y);
+ymax=max(md.y);
+md.thickness=hmax+(hmin-hmax)*(md.y-ymin)/(ymax-ymin);
+md.bed=-md.rho_ice/md.rho_water*md.thickness;
+md.surface=md.bed+md.thickness;
+
+md.drag_type=2; %0 none 1 plastic 2 viscous
+md.drag=200*ones(md.numberofgrids,1); %q=1.
+%Take care of iceshelves: no basal drag
+pos=find(md.elementoniceshelf);
+md.drag(md.elements(pos,:))=0;
+md.p=ones(md.numberofelements,1);
+md.q=ones(md.numberofelements,1);
+
+md.observed_temperature=(273-20)*ones(md.numberofgrids,1);
+
+md.B=paterson(md.observed_temperature);
+md.n=3*ones(md.numberofelements,1);
+
+%Deal with boundary conditions:
+md=SetMarineIceSheetBC(md,'Front.exp');
+
+%Parallel options
+md.np=3;
+md.time=50;
+md.waitonlock=30;
Index: /issm/trunk/test/Verification/test03_IceShelfIceFrontMP3d/configuration.m
===================================================================
--- /issm/trunk/test/Verification/test03_IceShelfIceFrontMP3d/configuration.m	(revision 2237)
+++ /issm/trunk/test/Verification/test03_IceShelfIceFrontMP3d/configuration.m	(revision 2237)
@@ -0,0 +1,31 @@
+%name
+testname='test3_IceShelfIceFrontMP3d';
+
+%meshing
+resolution=100000;
+domainname='DomainOutline.exp';
+riftname='';
+rifts=0;
+
+%geography
+iceshelfname='all';
+icesheetname='';
+
+%parameterization
+parametername='Square.par';
+			
+%elements type
+elementstype='''pattyn'',''Pattyn.exp'',''fill'',''macayeal''';
+
+%extrusion
+numlayers=8;
+extrusionexponent=4;
+
+%sequences     analysis    sub_analysis  qmu   control  control_fit parallel
+sequences={  {'prognostic',  'none',      0 ,     0,      'nan'        0    };...
+             {'thermal',     'steady',    0 ,     0,      'nan'        0    };...
+             {'thermal',     'transient', 0 ,     0,      'nan'        0    };...
+             {'prognostic',  'none',      0 ,     0,      'nan'        1    };...
+             {'thermal',     'steady',    0 ,     0,      'nan'        1    };...
+             {'thermal',     'transient', 0 ,     0,      'nan'        1    };...
+	};
Index: /issm/trunk/test/Verification/test03_IceShelfIceFrontMP3d/testpresolve.m
===================================================================
--- /issm/trunk/test/Verification/test03_IceShelfIceFrontMP3d/testpresolve.m	(revision 2237)
+++ /issm/trunk/test/Verification/test03_IceShelfIceFrontMP3d/testpresolve.m	(revision 2237)
@@ -0,0 +1,3 @@
+if md.numberofgrids==728
+	load Velocities; md.vx=0.8*vx; md.vy=0.8*vy; md.vz=0.8*vz; md.temperature=temperature-1; md.pressure=pressure;
+end
Index: /issm/trunk/test/Verification/test04_IceShelfIceFrontMPS3d/DomainOutline.exp
===================================================================
--- /issm/trunk/test/Verification/test04_IceShelfIceFrontMPS3d/DomainOutline.exp	(revision 2237)
+++ /issm/trunk/test/Verification/test04_IceShelfIceFrontMPS3d/DomainOutline.exp	(revision 2237)
@@ -0,0 +1,10 @@
+## Name:domainoutline
+## Icon:0
+# Points Count  Value
+5 1.
+# X pos Y pos
+0 0
+1000000 0
+1000000 1000000
+0 1000000
+0 0
Index: /issm/trunk/test/Verification/test04_IceShelfIceFrontMPS3d/Front.exp
===================================================================
--- /issm/trunk/test/Verification/test04_IceShelfIceFrontMPS3d/Front.exp	(revision 2237)
+++ /issm/trunk/test/Verification/test04_IceShelfIceFrontMPS3d/Front.exp	(revision 2237)
@@ -0,0 +1,10 @@
+## Name:icefront
+## Icon:0
+# Points Count  Value
+5 1.
+# X pos Y pos
+-1000 900000
+-1000 1100000
+1100000 1100000
+1100000 900000
+-1000 900000
Index: /issm/trunk/test/Verification/test04_IceShelfIceFrontMPS3d/Pattyn.exp
===================================================================
--- /issm/trunk/test/Verification/test04_IceShelfIceFrontMPS3d/Pattyn.exp	(revision 2237)
+++ /issm/trunk/test/Verification/test04_IceShelfIceFrontMPS3d/Pattyn.exp	(revision 2237)
@@ -0,0 +1,10 @@
+## Name:pattyndomain
+## Icon:0
+# Points Count  Value
+5 1.
+# X pos Y pos
+0 0
+500001 0
+500001 1000001
+0 1000001
+0 0
Index: /issm/trunk/test/Verification/test04_IceShelfIceFrontMPS3d/Square.par
===================================================================
--- /issm/trunk/test/Verification/test04_IceShelfIceFrontMPS3d/Square.par	(revision 2237)
+++ /issm/trunk/test/Verification/test04_IceShelfIceFrontMPS3d/Square.par	(revision 2237)
@@ -0,0 +1,37 @@
+%Start defining model parameters here
+
+%dynamics
+md.debug=0;
+md.dt=0.1; %1 year
+md.ndt=md.dt*3; 
+md.artificial_diffusivity=1;
+md.eps_abs=5000;
+
+hmin=300;
+hmax=1000;
+ymin=min(md.y);
+ymax=max(md.y);
+md.thickness=hmax+(hmin-hmax)*(md.y-ymin)/(ymax-ymin);
+md.bed=-md.rho_ice/md.rho_water*md.thickness;
+md.surface=md.bed+md.thickness;
+
+md.drag_type=2; %0 none 1 plastic 2 viscous
+md.drag=200*ones(md.numberofgrids,1); %q=1.
+%Take care of iceshelves: no basal drag
+pos=find(md.elementoniceshelf);
+md.drag(md.elements(pos,:))=0;
+md.p=ones(md.numberofelements,1);
+md.q=ones(md.numberofelements,1);
+
+md.observed_temperature=(273-20)*ones(md.numberofgrids,1);
+
+md.B=paterson(md.observed_temperature);
+md.n=3*ones(md.numberofelements,1);
+
+%Deal with boundary conditions:
+md=SetMarineIceSheetBC(md,'Front.exp');
+
+%Parallel options
+md.np=3;
+md.time=50;
+md.waitonlock=30;
Index: /issm/trunk/test/Verification/test04_IceShelfIceFrontMPS3d/Stokes.exp
===================================================================
--- /issm/trunk/test/Verification/test04_IceShelfIceFrontMPS3d/Stokes.exp	(revision 2237)
+++ /issm/trunk/test/Verification/test04_IceShelfIceFrontMPS3d/Stokes.exp	(revision 2237)
@@ -0,0 +1,10 @@
+## Name:stokesdomain
+## Icon:0
+# Points Count  Value
+5 1.
+# X pos Y pos
+250000 500000
+750001 500000
+750001 1000001
+250000 1000001
+250000 500000
Index: /issm/trunk/test/Verification/test04_IceShelfIceFrontMPS3d/configuration.m
===================================================================
--- /issm/trunk/test/Verification/test04_IceShelfIceFrontMPS3d/configuration.m	(revision 2237)
+++ /issm/trunk/test/Verification/test04_IceShelfIceFrontMPS3d/configuration.m	(revision 2237)
@@ -0,0 +1,31 @@
+%name
+testname='test4_IceShelfIceFrontMPS3d';
+
+%meshing
+resolution=100000;
+domainname='DomainOutline.exp';
+riftname='';
+rifts=0;
+
+%geography
+iceshelfname='all';
+icesheetname='';
+
+%parameterization
+parametername='Square.par';
+			
+%elements type
+elementstype='''pattyn'',''Pattyn.exp'',''fill'',''macayeal'',''stokes'',''Stokes.exp''';
+
+%extrusion
+numlayers=8;
+extrusionexponent=4;
+
+%sequences     analysis    sub_analysis  qmu   control  control_fit parallel
+sequences={  {'prognostic',  'none',      0 ,     0,      'nan'        0    };...
+             {'thermal',     'steady',    0 ,     0,      'nan'        0    };...
+             {'thermal',     'transient', 0 ,     0,      'nan'        0    };...
+             {'prognostic',  'none',      0 ,     0,      'nan'        1    };...
+             {'thermal',     'steady',    0 ,     0,      'nan'        1    };...
+             {'thermal',     'transient', 0 ,     0,      'nan'        1    };...
+	};
Index: /issm/trunk/test/Verification/test04_IceShelfIceFrontMPS3d/testpresolve.m
===================================================================
--- /issm/trunk/test/Verification/test04_IceShelfIceFrontMPS3d/testpresolve.m	(revision 2237)
+++ /issm/trunk/test/Verification/test04_IceShelfIceFrontMPS3d/testpresolve.m	(revision 2237)
@@ -0,0 +1,3 @@
+if md.numberofgrids==728
+	load Velocities; md.vx=0.8*vx; md.vy=0.8*vy; md.vz=0.8*vz; md.temperature=temperature-1; md.pressure=pressure;
+end
Index: /issm/trunk/test/Verification/test05_IceShelfIceFrontP3d/DomainOutline.exp
===================================================================
--- /issm/trunk/test/Verification/test05_IceShelfIceFrontP3d/DomainOutline.exp	(revision 2237)
+++ /issm/trunk/test/Verification/test05_IceShelfIceFrontP3d/DomainOutline.exp	(revision 2237)
@@ -0,0 +1,10 @@
+## Name:domainoutline
+## Icon:0
+# Points Count  Value
+5 1.
+# X pos Y pos
+0 0
+1000000 0
+1000000 1000000
+0 1000000
+0 0
Index: /issm/trunk/test/Verification/test05_IceShelfIceFrontP3d/Front.exp
===================================================================
--- /issm/trunk/test/Verification/test05_IceShelfIceFrontP3d/Front.exp	(revision 2237)
+++ /issm/trunk/test/Verification/test05_IceShelfIceFrontP3d/Front.exp	(revision 2237)
@@ -0,0 +1,10 @@
+## Name:icefront
+## Icon:0
+# Points Count  Value
+5 1.
+# X pos Y pos
+-1000 900000
+-1000 1100000
+1100000 1100000
+1100000 900000
+-1000 900000
Index: /issm/trunk/test/Verification/test05_IceShelfIceFrontP3d/Square.par
===================================================================
--- /issm/trunk/test/Verification/test05_IceShelfIceFrontP3d/Square.par	(revision 2237)
+++ /issm/trunk/test/Verification/test05_IceShelfIceFrontP3d/Square.par	(revision 2237)
@@ -0,0 +1,36 @@
+%Start defining model parameters here
+
+%dynamics
+md.debug=0;
+md.dt=0.1; %1 year
+md.ndt=md.dt*3; 
+md.artificial_diffusivity=0;
+
+hmin=300;
+hmax=1000;
+ymin=min(md.y);
+ymax=max(md.y);
+md.thickness=hmax+(hmin-hmax)*(md.y-ymin)/(ymax-ymin);
+md.bed=-md.rho_ice/md.rho_water*md.thickness;
+md.surface=md.bed+md.thickness;
+
+md.drag_type=2; %0 none 1 plastic 2 viscous
+md.drag=200*ones(md.numberofgrids,1); %q=1.
+%Take care of iceshelves: no basal drag
+pos=find(md.elementoniceshelf);
+md.drag(md.elements(pos,:))=0;
+md.p=ones(md.numberofelements,1);
+md.q=ones(md.numberofelements,1);
+
+md.observed_temperature=(273-20)*ones(md.numberofgrids,1);
+
+md.B=paterson(md.observed_temperature);
+md.n=3*ones(md.numberofelements,1);
+
+%Deal with boundary conditions:
+md=SetMarineIceSheetBC(md,'Front.exp');
+
+%Parallel options
+md.np=3;
+md.time=50;
+md.waitonlock=30;
Index: /issm/trunk/test/Verification/test05_IceShelfIceFrontP3d/configuration.m
===================================================================
--- /issm/trunk/test/Verification/test05_IceShelfIceFrontP3d/configuration.m	(revision 2237)
+++ /issm/trunk/test/Verification/test05_IceShelfIceFrontP3d/configuration.m	(revision 2237)
@@ -0,0 +1,35 @@
+%name
+testname='test5_IceShelfIceFrontP3d';
+
+%meshing
+resolution=100000;
+domainname='DomainOutline.exp';
+riftname='';
+rifts=0;
+
+%geography
+iceshelfname='all';
+icesheetname='';
+
+%parameterization
+parametername='Square.par';
+			
+%elements type
+elementstype='''pattyn'',''all''';
+
+%extrusion
+numlayers=8;
+extrusionexponent=4;
+
+%sequences     analysis    sub_analysis  qmu   control  control_fit parallel
+sequences={  {'diagnostic',  'none',      0 ,     0,      'nan'        0    };...
+             {'prognostic',  'none',      0 ,     0,      'nan'        0    };...
+             {'thermal',     'steady',    0 ,     0,      'nan'        0    };...
+             {'thermal',     'transient', 0 ,     0,      'nan'        0    };...
+             {'transient',   'none',      0 ,     0,      'nan'        0    };...
+				 {'diagnostic',  'none',      0 ,     0,      'nan'        1    };...
+             {'prognostic',  'none',      0 ,     0,      'nan'        1    };...
+             {'thermal',     'steady',    0 ,     0,      'nan'        1    };...
+             {'thermal',     'transient', 0 ,     0,      'nan'        1    };...
+             {'transient',   'none',      0 ,     0,      'nan'        1    };...
+	};
Index: /issm/trunk/test/Verification/test05_IceShelfIceFrontP3d/testpresolve.m
===================================================================
--- /issm/trunk/test/Verification/test05_IceShelfIceFrontP3d/testpresolve.m	(revision 2237)
+++ /issm/trunk/test/Verification/test05_IceShelfIceFrontP3d/testpresolve.m	(revision 2237)
@@ -0,0 +1,3 @@
+if md.numberofgrids==728
+	load Velocities; md.vx=0.8*vx; md.vy=0.8*vy; md.vz=0.8*vz; md.temperature=temperature-1; md.pressure=pressure;
+end
Index: /issm/trunk/test/Verification/test06_IceShelfIceFrontS3d/DomainOutline.exp
===================================================================
--- /issm/trunk/test/Verification/test06_IceShelfIceFrontS3d/DomainOutline.exp	(revision 2237)
+++ /issm/trunk/test/Verification/test06_IceShelfIceFrontS3d/DomainOutline.exp	(revision 2237)
@@ -0,0 +1,10 @@
+## Name:domainoutline
+## Icon:0
+# Points Count  Value
+5 1.
+# X pos Y pos
+0 0
+1000000 0
+1000000 1000000
+0 1000000
+0 0
Index: /issm/trunk/test/Verification/test06_IceShelfIceFrontS3d/Front.exp
===================================================================
--- /issm/trunk/test/Verification/test06_IceShelfIceFrontS3d/Front.exp	(revision 2237)
+++ /issm/trunk/test/Verification/test06_IceShelfIceFrontS3d/Front.exp	(revision 2237)
@@ -0,0 +1,10 @@
+## Name:icefront
+## Icon:0
+# Points Count  Value
+5 1.
+# X pos Y pos
+-1000 900000
+-1000 1100000
+1100000 1100000
+1100000 900000
+-1000 900000
Index: /issm/trunk/test/Verification/test06_IceShelfIceFrontS3d/Square.par
===================================================================
--- /issm/trunk/test/Verification/test06_IceShelfIceFrontS3d/Square.par	(revision 2237)
+++ /issm/trunk/test/Verification/test06_IceShelfIceFrontS3d/Square.par	(revision 2237)
@@ -0,0 +1,36 @@
+%Start defining model parameters here
+
+%dynamics
+md.debug=0;
+md.dt=0.1; %1 year
+md.ndt=md.dt*3; 
+md.artificial_diffusivity=1;
+
+hmin=300;
+hmax=1000;
+ymin=min(md.y);
+ymax=max(md.y);
+md.thickness=hmax+(hmin-hmax)*(md.y-ymin)/(ymax-ymin);
+md.bed=-md.rho_ice/md.rho_water*md.thickness;
+md.surface=md.bed+md.thickness;
+
+md.drag_type=2; %0 none 1 plastic 2 viscous
+md.drag=200*ones(md.numberofgrids,1); %q=1.
+%Take care of iceshelves: no basal drag
+pos=find(md.elementoniceshelf);
+md.drag(md.elements(pos,:))=0;
+md.p=ones(md.numberofelements,1);
+md.q=ones(md.numberofelements,1);
+
+md.observed_temperature=(273-20)*ones(md.numberofgrids,1);
+
+md.B=paterson(md.observed_temperature);
+md.n=3*ones(md.numberofelements,1);
+
+%Deal with boundary conditions:
+md=SetMarineIceSheetBC(md,'Front.exp');
+
+%Parallel options
+md.np=3;
+md.time=50;
+md.waitonlock=30;
Index: /issm/trunk/test/Verification/test06_IceShelfIceFrontS3d/configuration.m
===================================================================
--- /issm/trunk/test/Verification/test06_IceShelfIceFrontS3d/configuration.m	(revision 2237)
+++ /issm/trunk/test/Verification/test06_IceShelfIceFrontS3d/configuration.m	(revision 2237)
@@ -0,0 +1,35 @@
+%name
+testname='test6_IceShelfIceFrontS3d';
+
+%meshing
+resolution=100000;
+domainname='DomainOutline.exp';
+riftname='';
+rifts=0;
+
+%geography
+iceshelfname='all';
+icesheetname='';
+
+%parameterization
+parametername='Square.par';
+			
+%elements type
+elementstype='''pattyn'',''all'',''stokes'',''all''';
+
+%extrusion
+numlayers=8;
+extrusionexponent=4;
+
+%sequences     analysis    sub_analysis  qmu   control  control_fit parallel
+sequences={  {'diagnostic',  'none',      0 ,     0,      'nan'        0    };...
+             {'prognostic',  'none',      0 ,     0,      'nan'        0    };...
+             {'steadystate', 'none',      0 ,     0,      'nan'        0    };...
+             {'thermal',     'steady',    0 ,     0,      'nan'        0    };...
+             {'thermal',     'transient', 0 ,     0,      'nan'        0    };...
+				 {'diagnostic',  'none',      0 ,     0,      'nan'        1    };...
+             {'prognostic',  'none',      0 ,     0,      'nan'        1    };...
+             {'steadystate', 'none',      0 ,     0,      'nan'        1    };...
+             {'thermal',     'steady',    0 ,     0,      'nan'        1    };...
+             {'thermal',     'transient', 0 ,     0,      'nan'        1    };...
+	};
Index: /issm/trunk/test/Verification/test06_IceShelfIceFrontS3d/testpresolve.m
===================================================================
--- /issm/trunk/test/Verification/test06_IceShelfIceFrontS3d/testpresolve.m	(revision 2237)
+++ /issm/trunk/test/Verification/test06_IceShelfIceFrontS3d/testpresolve.m	(revision 2237)
@@ -0,0 +1,3 @@
+if md.numberofgrids==728
+	load Velocities; md.vx=0.8*vx; md.vy=0.8*vy; md.vz=0.8*vz; md.temperature=temperature-1; md.pressure=pressure;
+end
Index: /issm/trunk/test/Verification/test07_IceSheetIceFrontM2d/DomainOutline.exp
===================================================================
--- /issm/trunk/test/Verification/test07_IceSheetIceFrontM2d/DomainOutline.exp	(revision 2237)
+++ /issm/trunk/test/Verification/test07_IceSheetIceFrontM2d/DomainOutline.exp	(revision 2237)
@@ -0,0 +1,10 @@
+## Name:domainoutline
+## Icon:0
+# Points Count  Value
+5 1.
+# X pos Y pos
+0 0
+1000000 0
+1000000 1000000
+0 1000000
+0 0
Index: /issm/trunk/test/Verification/test07_IceSheetIceFrontM2d/Front.exp
===================================================================
--- /issm/trunk/test/Verification/test07_IceSheetIceFrontM2d/Front.exp	(revision 2237)
+++ /issm/trunk/test/Verification/test07_IceSheetIceFrontM2d/Front.exp	(revision 2237)
@@ -0,0 +1,10 @@
+## Name:icefront
+## Icon:0
+# Points Count  Value
+5 1.
+# X pos Y pos
+-1000 900000
+-1000 1100000
+1100000 1100000
+1100000 900000
+-1000 900000
Index: /issm/trunk/test/Verification/test07_IceSheetIceFrontM2d/Square.par
===================================================================
--- /issm/trunk/test/Verification/test07_IceSheetIceFrontM2d/Square.par	(revision 2237)
+++ /issm/trunk/test/Verification/test07_IceSheetIceFrontM2d/Square.par	(revision 2237)
@@ -0,0 +1,42 @@
+%Start defining model parameters here
+
+%dynamics
+md.debug=0;
+md.dt=1; %1 year
+md.ndt=md.dt*10; 
+md.artificial_diffusivity=0;
+
+hmin=300;
+hmax=1000;
+ymin=min(md.y);
+ymax=max(md.y);
+md.thickness=hmax+(hmin-hmax)*(md.y-ymin)/(ymax-ymin);
+md.bed=-md.rho_ice/md.rho_water*md.thickness;
+md.surface=md.bed+md.thickness;
+
+md.drag_type=2; %0 none 1 plastic 2 viscous
+md.drag=200*ones(md.numberofgrids,1); %q=1.
+%Take care of iceshelves: no basal drag
+pos=find(md.elementoniceshelf);
+md.drag(md.elements(pos,:))=0;
+md.p=ones(md.numberofelements,1);
+md.q=ones(md.numberofelements,1);
+md.viscosity_overshoot=0.3;
+
+md.observed_temperature=(273-20)*ones(md.numberofgrids,1);
+
+md.B=paterson(md.observed_temperature);
+md.n=3*ones(md.numberofelements,1);
+
+%Deal with boundary conditions:
+md=SetMarineIceSheetBC(md,'Front.exp');
+
+%Parallel options
+md.np=3;
+md.time=50;
+md.waitonlock=30;
+
+%load initial velocities
+if md.numberofgrids==340
+	load Velocities; md.vx=0.5*vx; md.vy=0.5*vy;
+end
Index: /issm/trunk/test/Verification/test07_IceSheetIceFrontM2d/configuration.m
===================================================================
--- /issm/trunk/test/Verification/test07_IceSheetIceFrontM2d/configuration.m	(revision 2237)
+++ /issm/trunk/test/Verification/test07_IceSheetIceFrontM2d/configuration.m	(revision 2237)
@@ -0,0 +1,30 @@
+%name
+testname='test7_IceSheetIceFrontM2d';
+
+%meshing
+resolution=50000;
+domainname='DomainOutline.exp';
+riftname='';
+rifts=0;
+
+%geography
+iceshelfname='';
+icesheetname='';
+
+%parameterization
+parametername='Square.par';
+			
+%elements type
+elementstype='''macayeal'',''all''';
+
+%extrusion
+numlayers=0;
+extrusionexponent=NaN;
+
+%sequences     analysis    sub_analysis  qmu   control  control_fit parallel
+sequences={  {'diagnostic',  'none',      0 ,     0,      'nan'        0    };...
+             {'prognostic',  'none',      0 ,     0,      'nan'        0    };...
+             {'transient',   'none',      0 ,     0,      'nan'        0    };...
+				 {'diagnostic',  'none',      0 ,     0,      'nan'        1    };...
+             {'transient',   'none',      0 ,     0,      'nan'        1    };...
+	};
Index: /issm/trunk/test/Verification/test08_IceSheetIceFrontM3d/DomainOutline.exp
===================================================================
--- /issm/trunk/test/Verification/test08_IceSheetIceFrontM3d/DomainOutline.exp	(revision 2237)
+++ /issm/trunk/test/Verification/test08_IceSheetIceFrontM3d/DomainOutline.exp	(revision 2237)
@@ -0,0 +1,10 @@
+## Name:domainoutline
+## Icon:0
+# Points Count  Value
+5 1.
+# X pos Y pos
+0 0
+1000000 0
+1000000 1000000
+0 1000000
+0 0
Index: /issm/trunk/test/Verification/test08_IceSheetIceFrontM3d/Front.exp
===================================================================
--- /issm/trunk/test/Verification/test08_IceSheetIceFrontM3d/Front.exp	(revision 2237)
+++ /issm/trunk/test/Verification/test08_IceSheetIceFrontM3d/Front.exp	(revision 2237)
@@ -0,0 +1,10 @@
+## Name:icefront
+## Icon:0
+# Points Count  Value
+5 1.
+# X pos Y pos
+-1000 900000
+-1000 1100000
+1100000 1100000
+1100000 900000
+-1000 900000
Index: /issm/trunk/test/Verification/test08_IceSheetIceFrontM3d/Square.par
===================================================================
--- /issm/trunk/test/Verification/test08_IceSheetIceFrontM3d/Square.par	(revision 2237)
+++ /issm/trunk/test/Verification/test08_IceSheetIceFrontM3d/Square.par	(revision 2237)
@@ -0,0 +1,38 @@
+%Start defining model parameters here
+
+%dynamics
+md.debug=0;
+md.dt=0.1; %1 year
+md.ndt=md.dt*3; 
+md.artificial_diffusivity=1;
+md.eps_abs=5000;
+
+hmin=300;
+hmax=1000;
+ymin=min(md.y);
+ymax=max(md.y);
+md.thickness=hmax+(hmin-hmax)*(md.y-ymin)/(ymax-ymin);
+md.bed=-md.rho_ice/md.rho_water*md.thickness;
+md.surface=md.bed+md.thickness;
+
+md.drag_type=2; %0 none 1 plastic 2 viscous
+md.drag=200*ones(md.numberofgrids,1); %q=1.
+%Take care of iceshelves: no basal drag
+pos=find(md.elementoniceshelf);
+md.drag(md.elements(pos,:))=0;
+md.p=ones(md.numberofelements,1);
+md.q=ones(md.numberofelements,1);
+md.viscosity_overshoot=0.0;
+
+md.observed_temperature=(273-20)*ones(md.numberofgrids,1);
+
+md.B=paterson(md.observed_temperature);
+md.n=3*ones(md.numberofelements,1);
+
+%Deal with boundary conditions:
+md=SetMarineIceSheetBC(md,'Front.exp');
+
+%Parallel options
+md.np=3;
+md.time=50;
+md.waitonlock=30;
Index: /issm/trunk/test/Verification/test08_IceSheetIceFrontM3d/configuration.m
===================================================================
--- /issm/trunk/test/Verification/test08_IceSheetIceFrontM3d/configuration.m	(revision 2237)
+++ /issm/trunk/test/Verification/test08_IceSheetIceFrontM3d/configuration.m	(revision 2237)
@@ -0,0 +1,35 @@
+%name
+testname='test8_IceSheetIceFrontM3d';
+
+%meshing
+resolution=100000;
+domainname='DomainOutline.exp';
+riftname='';
+rifts=0;
+
+%geography
+iceshelfname='';
+icesheetname='';
+
+%parameterization
+parametername='Square.par';
+			
+%elements type
+elementstype='''macayeal'',''all''';
+
+%extrusion
+numlayers=8;
+extrusionexponent=4;
+
+%sequences     analysis    sub_analysis  qmu   control  control_fit parallel
+sequences={  {'diagnostic',  'none',      0 ,     0,      'nan'        0    };...
+             {'prognostic',  'none',      0 ,     0,      'nan'        0    };...
+             {'thermal',     'steady',    0 ,     0,      'nan'        0    };...
+             {'thermal',     'transient', 0 ,     0,      'nan'        0    };...
+             {'transient',   'none',      0 ,     0,      'nan'        0    };...
+				 {'diagnostic',  'none',      0 ,     0,      'nan'        1    };...
+             {'prognostic',  'none',      0 ,     0,      'nan'        1    };...
+             {'thermal',     'steady',    0 ,     0,      'nan'        1    };...
+             {'thermal',     'transient', 0 ,     0,      'nan'        1    };...
+             {'transient',   'none',      0 ,     0,      'nan'        1    };...
+	};
Index: /issm/trunk/test/Verification/test08_IceSheetIceFrontM3d/testpresolve.m
===================================================================
--- /issm/trunk/test/Verification/test08_IceSheetIceFrontM3d/testpresolve.m	(revision 2237)
+++ /issm/trunk/test/Verification/test08_IceSheetIceFrontM3d/testpresolve.m	(revision 2237)
@@ -0,0 +1,3 @@
+if md.numberofgrids==728,
+	load Velocities; md.vx=0.3*vx; md.vy=0.3*vy; md.vz=0.3*vz; md.temperature=temperature-1; md.pressure=pressure;
+end
Index: /issm/trunk/test/Verification/test09_IceSheetIceFrontMP3d/DomainOutline.exp
===================================================================
--- /issm/trunk/test/Verification/test09_IceSheetIceFrontMP3d/DomainOutline.exp	(revision 2237)
+++ /issm/trunk/test/Verification/test09_IceSheetIceFrontMP3d/DomainOutline.exp	(revision 2237)
@@ -0,0 +1,10 @@
+## Name:domainoutline
+## Icon:0
+# Points Count  Value
+5 1.
+# X pos Y pos
+0 0
+1000000 0
+1000000 1000000
+0 1000000
+0 0
Index: /issm/trunk/test/Verification/test09_IceSheetIceFrontMP3d/Front.exp
===================================================================
--- /issm/trunk/test/Verification/test09_IceSheetIceFrontMP3d/Front.exp	(revision 2237)
+++ /issm/trunk/test/Verification/test09_IceSheetIceFrontMP3d/Front.exp	(revision 2237)
@@ -0,0 +1,10 @@
+## Name:icefront
+## Icon:0
+# Points Count  Value
+5 1.
+# X pos Y pos
+-1000 900000
+-1000 1100000
+1100000 1100000
+1100000 900000
+-1000 900000
Index: /issm/trunk/test/Verification/test09_IceSheetIceFrontMP3d/Pattyn.exp
===================================================================
--- /issm/trunk/test/Verification/test09_IceSheetIceFrontMP3d/Pattyn.exp	(revision 2237)
+++ /issm/trunk/test/Verification/test09_IceSheetIceFrontMP3d/Pattyn.exp	(revision 2237)
@@ -0,0 +1,10 @@
+## Name:pattyndomain
+## Icon:0
+# Points Count  Value
+5 1.
+# X pos Y pos
+0 0
+500001 0
+500001 1000001
+0 1000001
+0 0
Index: /issm/trunk/test/Verification/test09_IceSheetIceFrontMP3d/Square.par
===================================================================
--- /issm/trunk/test/Verification/test09_IceSheetIceFrontMP3d/Square.par	(revision 2237)
+++ /issm/trunk/test/Verification/test09_IceSheetIceFrontMP3d/Square.par	(revision 2237)
@@ -0,0 +1,38 @@
+%Start defining model parameters here
+
+%dynamics
+md.debug=0;
+md.dt=0.1; %1 year
+md.ndt=md.dt*3; 
+md.artificial_diffusivity=1;
+md.eps_abs=5000;
+md.min_thermal_constraints=1;
+
+hmin=300;
+hmax=1000;
+ymin=min(md.y);
+ymax=max(md.y);
+md.thickness=hmax+(hmin-hmax)*(md.y-ymin)/(ymax-ymin);
+md.bed=-md.rho_ice/md.rho_water*md.thickness;
+md.surface=md.bed+md.thickness;
+
+md.drag_type=2; %0 none 1 plastic 2 viscous
+md.drag=200*ones(md.numberofgrids,1); %q=1.
+%Take care of iceshelves: no basal drag
+pos=find(md.elementoniceshelf);
+md.drag(md.elements(pos,:))=0;
+md.p=ones(md.numberofelements,1);
+md.q=ones(md.numberofelements,1);
+
+md.observed_temperature=(273-20)*ones(md.numberofgrids,1);
+
+md.B=paterson(md.observed_temperature);
+md.n=3*ones(md.numberofelements,1);
+
+%Deal with boundary conditions:
+md=SetMarineIceSheetBC(md,'Front.exp');
+
+%Parallel options
+md.np=3;
+md.time=50;
+md.waitonlock=30;
Index: /issm/trunk/test/Verification/test09_IceSheetIceFrontMP3d/configuration.m
===================================================================
--- /issm/trunk/test/Verification/test09_IceSheetIceFrontMP3d/configuration.m	(revision 2237)
+++ /issm/trunk/test/Verification/test09_IceSheetIceFrontMP3d/configuration.m	(revision 2237)
@@ -0,0 +1,31 @@
+%name
+testname='test9_IceSheetIceFrontMP3d';
+
+%meshing
+resolution=100000;
+domainname='DomainOutline.exp';
+riftname='';
+rifts=0;
+
+%geography
+iceshelfname='';
+icesheetname='';
+
+%parameterization
+parametername='Square.par';
+			
+%elements type
+elementstype='''pattyn'',''Pattyn.exp'',''fill'',''macayeal''';
+
+%extrusion
+numlayers=8;
+extrusionexponent=4;
+
+%sequences     analysis    sub_analysis  qmu   control  control_fit parallel
+sequences={  {'prognostic',  'none',      0 ,     0,      'nan'        0    };...
+             {'thermal',     'steady',    0 ,     0,      'nan'        0    };...
+             {'thermal',     'transient', 0 ,     0,      'nan'        0    };...
+             {'prognostic',  'none',      0 ,     0,      'nan'        1    };...
+             {'thermal',     'steady',    0 ,     0,      'nan'        1    };...
+             {'thermal',     'transient', 0 ,     0,      'nan'        1    };...
+	};
Index: /issm/trunk/test/Verification/test09_IceSheetIceFrontMP3d/testpresolve.m
===================================================================
--- /issm/trunk/test/Verification/test09_IceSheetIceFrontMP3d/testpresolve.m	(revision 2237)
+++ /issm/trunk/test/Verification/test09_IceSheetIceFrontMP3d/testpresolve.m	(revision 2237)
@@ -0,0 +1,3 @@
+if md.numberofgrids==728,
+	load Velocities; md.vx=0.3*vx; md.vy=0.3*vy; md.vz=0.3*vz; md.temperature=temperature-1; md.pressure=pressure;
+end