Changeset 9223

Timestamp:

08/09/11 14:06:28 (14 years ago)

Author:

Mathieu Morlighem

Message:

better test 1201 EISMINT to be debugged

File:

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

issm/trunk/test/NightlyRun/test1201.m

@@ -4 +4 @@
 results={};
 
-for type=1:3;
+for stabilization=1:3;
         %The goal is to test the prognostic model
         md=bamg(model,'domain','../Exp/SquareEISMINT.exp','hmax',3000);
@@ -16 +16 @@
         md.vy=-400*ones(md.numberofnodes,1);
         md.vz=zeros(md.numberofnodes,1);
+        md.vel= 400*ones(md.numberofnodes,1);
 
-        %analytical test
-        md.spcvx=zeros(md.numberofnodes,1);
-        md.spcvy=zeros(md.numberofnodes,1);
-        md.spcvz=zeros(md.numberofnodes,1);
-        md.spcvx(:)=0;
-        md.spcvy(:)=-400;
-        md.spcvx(1)=0;
-        if type==1,
+        %Stabilization
+        if stabilization==1,
                 md.artificial_diffusivity=1;
                 md.prognostic_DG=0;
-        elseif type==2,
+        elseif stabilization==2,
                 md.artificial_diffusivity=0;
                 md.prognostic_DG=0;
-        elseif type==3,
+        elseif stabilization==3,
                 md.prognostic_DG=1;
         end
 
-        %Launch many prognostic as we don't really want a transient
-        %FOR NOW:
-        i=0;
-        time=0;
-        md.dummy=struct();
+        %spc thickness
         pos=find(md.y>199999.9);
-        connectivity=full(sparse(md.elements(:),1,1));
-        while(time<500),
-                disp(['step ' num2str(i) '/' num2str(500/md.dt)]);
-                i=i+1;
-                time=time+md.dt;
-                md=solve(md,PrognosticSolutionEnum);
-                md.thickness=PatchToVec(md.results.PrognosticSolution.Thickness);
-                md.thickness(pos)=500+100*sin(2*pi*time/200)*ones(size(pos,1),1);
-                md.surface=md.bed+md.thickness;
-                md.dummy(i).thickness=md.thickness;
-        end
-        results{type}=md.thickness;
+        times=0:1:500;
+        md.spcthickness=NaN*ones(md.numberofnodes+1,length(times));
+        md.spcthickness(end,:)=times;
+        md.spcthickness(pos,:)=repmat(500+100*sin(2*pi*times/200),length(pos),1);
+
+        %solve
+        md.isdiagnostic=0;
+        md.output_frequency=500; %keep only last step
+        md=solve(md,TransientSolutionEnum);
+        results{stabilization}=PatchToVec(md.results.TransientSolution(end).Thickness);
 end