Index: /issm/trunk-jpl/test/NightlyRun/test2003.m
===================================================================
--- /issm/trunk-jpl/test/NightlyRun/test2003.m	(revision 26099)
+++ /issm/trunk-jpl/test/NightlyRun/test2003.m	(revision 26100)
@@ -5,11 +5,19 @@
 md.mesh=gmshplanet('radius',6.371012*10^3,'resolution',1000.); %1000 km resolution mesh
 
+%Geometry for the bed, arbitrary thickness of 1000: 
+md.geometry.bed=-ones(md.mesh.numberofvertices,1);
+md.geometry.base=md.geometry.bed;
+md.geometry.thickness=1000*ones(md.mesh.numberofvertices,1);
+md.geometry.surface=md.geometry.bed+md.geometry.thickness;
+
+
 %parameterize slc solution:
 %solidearth loading:  {{{
-md.masstransport.spcthickness=zeros(md.mesh.numberofelements,1);
+md.masstransport.spcthickness=[md.geometry.thickness;0];
 md.initialization.sealevel=zeros(md.mesh.numberofvertices,1);
 md.dsl.global_average_thermosteric_sea_level=[0;0];
 md.dsl.sea_surface_height_above_geoid=zeros(md.mesh.numberofvertices+1,1);
 md.dsl.sea_water_pressure_at_sea_floor=zeros(md.mesh.numberofvertices+1,1);
+md.smb.mass_balance=zeros(md.mesh.numberofvertices,1);
 
 %antarctica
@@ -17,5 +25,5 @@
 longe=sum(md.mesh.long(md.mesh.elements),2)/3;
 pos=find(late <-75 & longe <0);
-md.masstransport.spcthickness(pos)=-1;
+md.masstransport.spcthickness(md.mesh.elements(pos,:))= md.masstransport.spcthickness(md.mesh.elements(pos,:))-1;
 
 %elastic loading from love numbers:
@@ -45,15 +53,10 @@
 md.solidearth.settings.ocean_area_scaling = 0;
 
-%geometry
-di=md.materials.rho_ice/md.materials.rho_water;
-md.geometry.thickness=ones(md.mesh.numberofvertices,1);
-md.geometry.surface=(1-di)*zeros(md.mesh.numberofvertices,1);
-md.geometry.base=md.geometry.surface-md.geometry.thickness;
-md.geometry.bed=md.geometry.base;
-
 %materials
 md.initialization.temperature=273.25*ones(md.mesh.numberofvertices,1);
-md.materials.rheology_B=paterson(md.initialization.temperature);
-md.materials.rheology_n=3*ones(md.mesh.numberofelements,1);
+md.initialization.sealevel=zeros(md.mesh.numberofvertices,1);
+md.initialization.bottompressure=zeros(md.mesh.numberofvertices,1);
+md.initialization.dsl=zeros(md.mesh.numberofvertices,1);
+md.initialization.str=0;
 
 md.basalforcings.groundedice_melting_rate=zeros(md.mesh.numberofvertices,1);
@@ -78,5 +81,11 @@
 md.transient.isthermal=0;
 md.transient.ismasstransport=1;
+md.transient.isoceantransport=1;
 md.transient.isslc=1;
+
+ 
+md.timestepping.start_time=0;
+md.timestepping.time_step=1;
+md.timestepping.final_time=1;