Changeset 9233

Timestamp:

08/09/11 16:20:32 (14 years ago)

Author:

Mathieu Morlighem

Message:

Fixed DG transient

File:

• issm/trunk/test/NightlyRun/test1201.m (modified) (3 diffs)

issm/trunk/test/NightlyRun/test1201.m

@@ -9 +9 @@
         md=geography(md,'all','');
         md=parameterize(md,'../Par/SquareEISMINT.par');
+        md.surface_mass_balance(:)=0;
         md=setelementstype(md,'macayeal','all');
         md.cluster=generic('name',oshostname(),'np',8);
@@ -15 +16 @@
         md.vx=zeros(md.numberofnodes,1);
         md.vy=-400*ones(md.numberofnodes,1);
-        md.vz=zeros(md.numberofnodes,1);
-        md.vel= 400*ones(md.numberofnodes,1);
 
         %Stabilization
@@ -35 +34 @@
         md.spcthickness(end,:)=times;
         md.spcthickness(pos,:)=repmat(500+100*sin(2*pi*times/200),length(pos),1);
-        md.surface_mass_balance(:)=0;
+        if stabilization==3,
+                pos=find(isnan(md.spcthickness)); md.spcthickness(pos)=500; %No NaN for DG
+        end
 
         %solve