Changeset 4850 for issm/trunk/test/Validation/ControlMethods/runmeB.m

Timestamp:

07/28/10 16:55:17 (15 years ago)

Author:

seroussi

Message:

updated Validation/ControlMethods for today's fields

File:

• issm/trunk/test/Validation/ControlMethods/runmeB.m (modified) (4 diffs)

issm/trunk/test/Validation/ControlMethods/runmeB.m

@@ -24 +24 @@
                 for j=1:division
                         md.observed_temperature(find(md.x>=(minx+(i-1)*Lx/division) & md.x<=(minx+i*Lx/division) & md.y>=(miny+(j-1)*Ly/division) & md.y<=(miny+j*Ly/division)))=temperature(count);
-                        md.B=paterson(md.observed_temperature);
+                        md.rheology_B=paterson(md.observed_temperature);
                         count=count+1;
                 end
@@ -31 +31 @@
         %diagnostic with this B
         md.cluster='wilkes'; 
-        md=solve(md,'analysis_type','diagnostic');
+        md=solve(md,'analysis_type',DiagnosticSolutionEnum);
 
         %save observed velocities and exact B
-        vx_obs=md.vx; vy_obs=md.vy; vel_obs=md.vel;
-        B=md.B;
-        save ObservedVelocities vx_obs vy_obs vel_obs
+        vx_obs=zeros(md.numberofgrids,1);
+        vx_obs(md.results.DiagnosticSolution.Vx.index)=md.results.DiagnosticSolution.Vx.value;
+        vy_obs=zeros(md.numberofgrids,1);
+        vy_obs(md.results.DiagnosticSolution.Vy.index)=md.results.DiagnosticSolution.Vy.value;
+        B=md.rheology_B;
+        save ObservedVelocities vx_obs vy_obs
         save ExactB B
         clear md
@@ -54 +57 @@
         %plug computed velocities
         load ObservedVelocities
-        md.vx_obs=vx_obs; md.vy_obs=vy_obs; md.vel_obs=vel_obs;
+        md.vx_obs=vx_obs; md.vy_obs=vy_obs;
 
         %control method
         md.nsteps=5;
-        md.control_type='B';
-        md.mincontrolconstraint=10^7;
-        md.maxcontrolconstraint=1.2*10^9;
+        md.control_type='rheology_B';
+        md.cm_min=10^7;
+        md.cm_max=1.2*10^9;
+        md.cm_jump=0.99*ones(md.nsteps,1);
         md.plot=1;
-        md.debug=0;
-        md.optscal=10^8*ones(md.nsteps,1); md.optscal(1)=2*10^8;
+        md.verbose=0;
+        md.optscal=2*10^8*ones(md.nsteps,1); md.optscal(1)=5*10^8;
         md.fit=0*ones(md.nsteps,1);
+        md.weights=1*ones(md.numberofgrids,1);
         md.maxiter=10*ones(md.nsteps,1);
         md.eps_rel=0.0001;
@@ -71 +76 @@
 
         md.np=8;
-        %md.cluster='wilkes';
-        %md=solve(md,'analysis_type','control');
-        md=solve(md,'analysis_type','control');
-        %md=solve(md,'analysis_type','control','package','macayeal');
+        md.control_analysis=1;
+        md=solve(md,'analysis_type',DiagnosticSolutionEnum);
 end