Changeset 3581

Timestamp:

04/20/10 16:31:54 (15 years ago)

Author:

seroussi

Message:

mass conservation 2

Location:

issm/trunk/test/Validation/EISMINT/MassConservation

Files:

• CrossLine2.exp (modified) (1 diff)
• runme2.m (modified) (6 diffs)

issm/trunk/test/Validation/EISMINT/MassConservation/CrossLine2.exp

@@ -4 +4 @@
 2 1.
 # X pos Y pos
-10000  199999
-110000 0
+0 0
+200000 100000

issm/trunk/test/Validation/EISMINT/MassConservation/runme2.m

@@ -4 +4 @@
 %The goal is to test the prognostic model
 %md=bamg(model,'domain','DomainOutline.exp','hmax',4550);
-md=bamg(model,'domain','DomainOutline.exp','hmax',3000);
+md=bamg(model,'domain','DomainOutline.exp','hmax',2000);
 md=geography(md,'all','');
 md=parameterize(md,'Square.par');
@@ -15 +15 @@
 md.vx=160*sqrt(5)*ones(md.numberofgrids,1);
 md.vy= 80*sqrt(5)*ones(md.numberofgrids,1);
+sinus=1/5*sqrt(5);
+cosinus=2/5*sqrt(5);
 
 %0launch transient solution
@@ -20 +22 @@
 %FOR NOW:
 md.cluster=oshostname();
-md.np=10;
+md.np=14;
 i=0;
 time=0;
@@ -27 +29 @@
 posLEFT=find(md.x<0.01);
 connectivity=full(sparse(md.elements(:),1,1));
-while(time<100),
-        disp(['step ' num2str(i) '/' num2str(500/md.dt)]);
+while(time<600),
+        disp(['step ' num2str(i) '/' num2str(600/md.dt)]);
         i=i+1;
         time=time+md.dt;
@@ -35 +37 @@
         %md.thickness=thickness./connectivity;
         md.thickness=md.results.prognostic2.thickness;
-        md.thickness(posLEFT)=500+100*sin((2*pi/200)*(time-md.y(posLEFT)/md.vy(1)));
-        md.thickness(posDOWN)=500+100*sin((2*pi/200)*(time-md.x(posDOWN)/md.vx(1)));
+        md.thickness(posLEFT)=500+100*sin((2*pi/200)*(time-sinus*md.y(posLEFT)/sqrt(md.vx(1)^2+md.vy(1)^2)));
+        md.thickness(posDOWN)=500+100*sin((2*pi/200)*(time-cosinus*md.x(posDOWN)/sqrt(md.vx(1)^2+md.vy(1)^2)));
         md.surface=md.bed+md.thickness;
         md.dummy(i).thickness=md.thickness;
@@ -43 +45 @@
 
 %plot results
-plotmodel(md,'data',[md.thickness 500+100*sin(2*pi/200*(500-md.y/400))],'sectionvalue','CrossLine2.exp')
+plotmodel(md,'data',[md.thickness 500+100*sin(2*pi/200*(600-(md.x*cosinus+md.y*sinus)/sqrt(md.vx(1)^2+md.vy(1)^2)))],'sectionvalue','CrossLine2.exp')
 
 %Don't forget to add these lines in icetransient2d.m, just before the computation of the thickness to change the thickness on the upper boundary condition