Changeset 18205

Timestamp:

06/30/14 18:34:01 (11 years ago)

Author:

schlegel

Message:

CHG: Greenland and melting rates

Location:

issm/trunk-jpl/examples

Files:

• Greenland/Greenland.par (modified) (1 diff)
• Greenland/runme.m (modified) (3 diffs)
• Jakobshavn/Jks.par (modified) (1 diff)
• Pig/Pig.par (modified) (1 diff)

issm/trunk-jpl/examples/Greenland/Greenland.par

@@ -58 +58 @@
 disp('   Set other boundary conditions');
 md.mask.ice_levelset(md.mesh.vertexonboundary==1)=0;
-md.basalforcings.melting_rate = zeros(md.mesh.numberofvertices,1);
+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);

issm/trunk-jpl/examples/Greenland/runme.m

@@ -1 +1 @@
 clear all;
-steps=[ 1 ];
+steps=[1];
 
 %Location of SeaRISE dataset
@@ -9 +9 @@
 
         %Generate initial uniform mesh (resolution = 20000 m)
-        md=triangle(model,'./Exp_Par/DomainOutline.exp',20000);
+        md=triangle(model,'./DomainOutline.exp',20000);
 
         % Get velocities (Note: You can use ncdisp('file') to see an ncdump)
@@ -36 +36 @@
 
         md = setmask(md,'','');
-        md = parameterize(md,'./Exp_Par/Greenland.par');
+        md = parameterize(md,'./Greenland.par');
         md = setflowequation(md,'SSA','all');
 

issm/trunk-jpl/examples/Jakobshavn/Jks.par

@@ -59 +59 @@
 md=SetMarineIceSheetBC(md,'./Front.exp');
 md.basalforcings.floatingice_melting_rate=zeros(md.mesh.numberofvertices,1);
+md.basalforcings.groundedice_melting_rate=zeros(md.mesh.numberofvertices,1);

issm/trunk-jpl/examples/Pig/Pig.par

@@ -130 +130 @@
 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);