Changeset 1238 for issm/trunk/test/Miscellaneous/Bump/runme.m

Timestamp:

07/06/09 11:22:46 (16 years ago)

Author:

seroussi

Message:

better bump + results in models

File:

• issm/trunk/test/Miscellaneous/Bump/runme.m (modified) (3 diffs)

issm/trunk/test/Miscellaneous/Bump/runme.m

@@ -3 +3 @@
 %initialize model with Stokes model
 md=model;
-md=mesh(md,'DomainOutline.exp',20000);
+md=mesh(md,'DomainOutline.exp',2000);
+md=meshexprefine(md,'ContourBump.exp',500);
 md=geography(md,'','');
-md=parameterize(md,'Square.par');
+md=parameterize(md,'Bump.par');
+md.drag=50*ones(md.numberofgrids,1);
 md=extrude(md,8,4);
 md.pressure=zeros(md.numberofgrids,1);
@@ -22 +24 @@
 %now run a control method on macayeal
 md=model;
-md=mesh(md,'DomainOutline.exp',20000);
+md=mesh(md,'DomainOutline.exp',2000);
+md=meshexprefine(md,'ContourBump.exp',500);
 md=geography(md,'','');
-md=parameterize(md,'Square.par');
+md=parameterize(md,'Bump.par');
 md.pressure=zeros(md.numberofgrids,1);
 md.vx=zeros(md.numberofgrids,1);
@@ -40 +43 @@
 md=solve(md,'analysis_type','control')
 
-%%now run a control method on pattyn
-%md=model;
-%md=mesh(md,'DomainOutline.exp',20000);
-%md=geography(md,'','');
-%md=parameterize(md,'Square.par');
-%md=extrude(md,8,4);
-%md.pressure=zeros(md.numberofgrids,1);
-%md.vx=zeros(md.numberofgrids,1);
-%md.vy=zeros(md.numberofgrids,1);
-%md.vz=zeros(md.numberofgrids,1);
-%md=setelementstype(md,'pattyn','all');
-%
-%%Initialize the CM
-%md.vx_obs=mds.results.diagnostic.vx;
-%md.vy_obs=mds.results.diagnostic.vy;
-%md.vel_obs=mds.results.diagnostic.vel;
-%
-%%Compute the CM
-%md=solve(md,'analysis_type','control')
+%save model macayeal
+mdm=md;
+save modelmacayeal mdm
+
+%now run a control method on pattyn
+md=model;
+md=mesh(md,'DomainOutline.exp',2000);
+md=meshexprefine(md,'ContourBump.exp',500);
+md=geography(md,'','');
+md=parameterize(md,'Bump.par');
+md=extrude(md,8,4);
+md.pressure=zeros(md.numberofgrids,1);
+md.vx=zeros(md.numberofgrids,1);
+md.vy=zeros(md.numberofgrids,1);
+md.vz=zeros(md.numberofgrids,1);
+md=setelementstype(md,'pattyn','all');
+
+%Initialize the CM
+md.vx_obs=mds.results.diagnostic.vx;
+md.vy_obs=mds.results.diagnostic.vy;
+md.vel_obs=mds.results.diagnostic.vel;
+
+%Compute the CM
+md=solve(md,'analysis_type','control')
+
+%save model pattyn
+mdp=md;
+save modelpattyn mdp