Changeset 26865

Timestamp:

02/09/22 05:37:44 (3 years ago)

Author:

Mathieu Morlighem

Message:

NEW: added grounding line and icefront option to plot_transient_movie, and better handle output

File:

• issm/trunk-jpl/src/m/plot/plot_transient_movie.m (modified) (4 diffs)

issm/trunk-jpl/src/m/plot/plot_transient_movie.m

@@ -42 +42 @@
         end
 
+        %Do we have an output?
+        isavi = 0;
+        isgif = 0;
+        ismp4 = 0;
+        if exist(options,'transient_movie_output'),
+                filename=getfieldvalue(options,'transient_movie_output');
+                [pathstr,name,ext] = fileparts(filename);
+                if strcmp(ext,'.gif')
+                        isgif = 1;
+                elseif strcmp(ext,'.mp4')
+                        ismp4 = 1;
+                elseif strcmp(ext,'.avi')
+                        isavi = 1;
+                end
+        end
+        if isavi || ismp4
+                vid=VideoWriter([filename(1:end-4),'.avi'],'Motion JPEG AVI');
+                vid.FrameRate = 2; 
+                open(vid); 
+        end
+
         %calculate caxis
         if ~exist(options,'caxis'),
@@ -53 +74 @@
         end
 
+
+        %Process mesh once for all
+        [x y z elements is2d isplanet]=processmesh(md,results(i).(field),options);
+
         %display movie
         nstep=1;
@@ -60 +85 @@
                 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);
 
                         clf;
-                        titlestring=[field ' at time ' num2str(results(i).time,'%7.2f') ' year'];
+                        titlestring=[field ' (time ' num2str(results(i).time,'%7.2f') ' yr)'];
                         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)
+                        %Add grounding line
+                        if exist(options,'groundingline')
+                                contours=isoline(md, results(i).MaskOceanLevelset,'output','matrix');
+                                hold on
+                                plot(contours(:,1),contours(:,2),getfieldvalue(options,'groundingline'));
+                        end
+
+                        %Add ice front
+                        if exist(options,'icefront')
+                                contours=isoline(md, results(i).MaskIceLevelset,'output','matrix');
+                                hold on
+                                plot(contours(:,1),contours(:,2),getfieldvalue(options,'icefront'));
+                        end
+
+                        if isgif
+                                frame=getframe(gcf);
+                                im = frame2im(frame);
+                                [imind,cmap] = rgb2ind(im,256);
+                                if i==1
+                                        imwrite(imind, cmap, filename, 'DelayTime',getfieldvalue(options,'transient_movie_time',.5), 'LoopCount',inf)
+                                else
+                                        imwrite(imind, cmap, filename, 'WriteMode','append');
+                                end
+                        elseif isavi || ismp4
+                                F=getframe(gcf);
+                                writeVideo(vid, F);
+                        end
+
+                        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
@@ -87 +139 @@
 
         %output movie if requested.
-        if exist(options,'transient_movie_output'),
-                filename=getfieldvalue(options,'transient_movie_output');
-                imwrite(images,map,filename,'DelayTime',getfieldvalue(options,'transient_movie_time',2),'LoopCount',inf)
+        if isavi || ismp4
+                close(vid);
+        end
+        if ismp4
+                while(~exist([filename]))
+                        disp(['Waiting for ' filename ]);
+                        pause(1)
+                end
+                filename = filename(1:end-4);
+                command=sprintf('ffmpeg -y -i %s.avi -c:v libx264 -crf 19 -preset slow -c:a libfaac -b:a 192k -ac 2 %s.mp4',filename,filename);
+                system(command);
+                filename = [filename '.mp4'];
+                disp(['Movie saved also as ' filename]);
+        end
+        if isgif
+                imwrite(images,map,filename,'DelayTime',getfieldvalue(options,'transient_movie_time',.5),'LoopCount',inf)
         end