Changeset 14075

Timestamp:

11/30/12 15:38:54 (12 years ago)

Author:

schlegel

Message:

CHG: update to match talk pdf

Location:

issm/workshop/2012/Talks/12_SeaRISE

Files:

• Code/plotparam.m (modified) (1 diff)
• runme.m (modified) (5 diffs)

issm/workshop/2012/Talks/12_SeaRISE/Code/plotparam.m

@@ -1 @@
->> plotmodel(md,'data',md.initialization.thickness);
+>> plotmodel(md,'data',md.geometry.thickness);
 >> plotmodel(md,'data',md.initialization.vel,'caxis',[1e-1 1e4],'log',10);

issm/workshop/2012/Talks/12_SeaRISE/runme.m

@@ -22 +22 @@
         %Mesh Greenland
         md=bamg(md,'hmax',400000,'hmin',5000,'gradation',1.7,'field',vel,'err',8);
-        plotmodel (md,'data','mesh');
         
         %convert x,y coordinates (Polar stereo) to lat/lon
@@ -28 +27 @@
 
         save ./Models/Greenland.Mesh_generation md;
+
+        plotmodel (md,'data','mesh');
 end
 
@@ -50 +51 @@
         md.inversion.step_threshold=0.99*ones(md.inversion.nsteps,1);
         md.inversion.maxiter_per_step=5*ones(md.inversion.nsteps,1);
-        md.verbose=verbose('solution',true,'control',true);
 
         %Cost functions
@@ -66 +66 @@
 
         %Additional parameters
-        md.diagnostic.restol=0.01;
-        md.diagnostic.reltol=0.1;
+        md.diagnostic.restol=0.01; md.diagnostic.reltol=0.1; 
         md.diagnostic.abstol=NaN;
 
@@ -91 +90 @@
         y1  = ncread(ncdata,'y1');
         smb = ncread(ncdata,'smb');
-        smb = InterpFromGridToMesh(x1,y1,smb',md.mesh.x,md.mesh.y,0)*1000/md.materials.rho_ice;
+        smb = InterpFromGridToMesh(x1,y1,smb',md.mesh.x,md.mesh.y,0);
+        smb = smb*md.materials.rho_freshwater/md.materials.rho_ice;
         smb = [smb smb smb-1.0];
         md.surfaceforcings.mass_balance = [smb;1 10 20];