Index: /issm/workshop/2012/Talks/09_PIG_Application/Exp_Par/Pig.par
===================================================================
--- /issm/workshop/2012/Talks/09_PIG_Application/Exp_Par/Pig.par	(revision 14113)
+++ /issm/workshop/2012/Talks/09_PIG_Application/Exp_Par/Pig.par	(revision 14114)
@@ -36,4 +36,5 @@
 md.geometry.bed     = 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);
+clear usrf, topg;
 
 disp('   Constructing thickness');
@@ -54,8 +55,10 @@
 disp('   Interpolating temperatures');
 md.initialization.temperature=InterpFromGridToMesh(x1,y1,temp,md.mesh.x,md.mesh.y,0)+273.15+Temp_change;
+clear temp;
 
 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;
+clear smb;
 
 disp('   Construct basal friction parameters');
@@ -75,4 +78,6 @@
 vx_obs=InterpFromGridToMesh(x2,y2,flipud(velx'),md.mesh.x,md.mesh.y,0);
 vy_obs=InterpFromGridToMesh(x2,y2,flipud(vely'),md.mesh.x,md.mesh.y,0);
+clear velx vely;
+
 vel_obs=sqrt(vx_obs.^2+vy_obs.^2);
 md.initialization.vx=vx_obs;
@@ -92,4 +97,5 @@
 disp('   Set geothermal heat flux');
 md.basalforcings.geothermalflux=InterpFromGridToMesh(x1,y1,gflux,md.mesh.x,md.mesh.y,0);
+clear gflux;
 
 disp('   Set Pressure');