Changeset 10129

Timestamp:

10/06/11 12:15:03 (13 years ago)

Author:

schlegel

Message:

Extract regional transient for forcing velocity at boundaries

Location:

issm/trunk/src/m/model

Files:

• plot/plot_BC.m (modified) (1 diff)
• plot/plot_transient_movie.m (modified) (1 diff)
• plot/plotmodel.m (modified) (1 diff)
• regionaltransient2d.m (added)

issm/trunk/src/m/model/plot/plot_BC.m

@@ -7 +7 @@
 
 %plot dirichlets
-h1=plot(md.mesh.x(find(~isnan(md.diagnostic.spcvx))),md.mesh.y(find(~isnan(md.diagnostic.spcvx))),'ro','MarkerSize',14,'MarkerFaceColor','r');
-h2=plot(md.mesh.x(find(~isnan(md.diagnostic.spcvy))),md.mesh.y(find(~isnan(md.diagnostic.spcvy))),'bo','MarkerSize',10,'MarkerFaceColor','b');
-h3=plot(md.mesh.x(find(~isnan(md.diagnostic.spcvz))),md.mesh.y(find(~isnan(md.diagnostic.spcvz))),'yo','MarkerSize',6 ,'MarkerFaceColor','y');
+h1=plot(md.mesh.x(find(~isnan(md.diagnostic.spcvx(1:md.mesh.numberofvertices,1)))),md.mesh.y(find(~isnan(md.diagnostic.spcvx(1:md.mesh.numberofvertices,1)))),'ro','MarkerSize',14,'MarkerFaceColor','r');
+h2=plot(md.mesh.x(find(~isnan(md.diagnostic.spcvy(1:md.mesh.numberofvertices,1)))),md.mesh.y(find(~isnan(md.diagnostic.spcvy(1:md.mesh.numberofvertices,1)))),'bo','MarkerSize',10,'MarkerFaceColor','b');
+h3=plot(md.mesh.x(find(~isnan(md.diagnostic.spcvz(1:md.mesh.numberofvertices,1)))),md.mesh.y(find(~isnan(md.diagnostic.spcvz(1:md.mesh.numberofvertices,1)))),'yo','MarkerSize',6 ,'MarkerFaceColor','y');
 
 %update legend

issm/trunk/src/m/model/plot/plot_transient_movie.m

@@ -16 +16 @@
         end
 
+        results=md.results.TransientSolution;
         %loop over the time steps
+        if exist(options,'transient_movie_limit'),
+                limit=getfieldvalue(options,'transient_movie_limit');
+                steps=[limit(1):limit(end)];
+        else
+                steps=1:length(results);
+        end
         results=md.results.TransientSolution;
-        for i=1:length(results)
+        nstep=1;
+        for i=steps
 
-                %process data
-                [x y z elements is2d isplanet]=processmesh(md,results(i).(field),options);
-                [data datatype]=processdata(md,results(i).(field),options);
+                if ~isempty(results(i).(field)),
+                        %process data
+                        [x y z elements is2d isplanet]=processmesh(md,results(i).(field),options);
+                        [data datatype]=processdata(md,results(i).(field),options);
 
-                titlestring=[field ' at time ' num2str(results(i).time/md.constants.yts) ' year'];
-                plot_unit(x,y,z,elements,data,is2d,isplanet,datatype,options)
-                apply_options_movie(md,options,titlestring);
-                
-                if exist(options,'transient_movie_output'),
-                        set(gcf,'Renderer','zbuffer','color','white'); %fixes a bug on Mac OS X (not needed in future Matlab version)
-                        if i==1,
-                                %initialize images and frame
-                                frame=getframe(gcf);
-                                [images,map]=rgb2ind(frame.cdata,256,'nodither');
-                                images(1,1,1,length(results))=0;
+                        titlestring=[field ' at time ' num2str(results(i).time/md.constants.yts) ' year'];
+                        plot_unit(x,y,z,elements,data,is2d,isplanet,datatype,options)
+                        apply_options_movie(md,options,titlestring);
+
+                        if exist(options,'transient_movie_output'),
+                                set(gcf,'Renderer','zbuffer','color','white'); %fixes a bug on Mac OS X (not needed in future Matlab version)
+                                if nstep==1,
+                                        %initialize images and frame
+                                        frame=getframe(gcf);
+                                        [images,map]=rgb2ind(frame.cdata,256,'nodither');
+                                        images(1,1,1,length(steps))=0;
+                                else
+                                        frame=getframe(gcf);
+                                        images(:,:,1,nstep) = rgb2ind(frame.cdata,map,'nodither');
+                                end
                         else
-                                frame=getframe(gcf);
-                                images(:,:,1,i) = rgb2ind(frame.cdata,map,'nodither');
+                                pause(0.1)
                         end
-                else
-                        pause(0.1)
+                        nstep=nstep+1;
                 end
         end

issm/trunk/src/m/model/plot/plotmodel.m

@@ -38 +38 @@
                 set(F,'Visible','Off');
         else
-                figure(figurenumber),clf;
+                figure(figurenumber);clf;
         end