Changeset 19197 for issm/trunk-jpl/src/m/contrib/paraview/exportVTK.m

Timestamp:

03/16/15 10:23:06 (10 years ago)

Author:

bdef

Message:

NEW: adding support for forcings

File:

• issm/trunk-jpl/src/m/contrib/paraview/exportVTK.m (modified) (2 diffs)

issm/trunk-jpl/src/m/contrib/paraview/exportVTK.m

@@ -59 +59 @@
                 if(size(sol_struct{i},2)>num_of_timesteps);
                         num_of_timesteps=size(sol_struct{i},2);
-    end
+          end
+                outstep=model.timestepping.time_step*model.settings.output_frequency;
   end
 else
@@ -145 +146 @@
                                 s='%e\n';
                                 fprintf(fid,s,res_struct.(fieldnames{k}));
-            end         
+                                %check for forcings     
+                        elseif (size(res_struct.(fieldnames{k}),1)==num_of_points+1);
+                                %paraview does not like NaN, replacing
+                                nanval=find(isnan(res_struct.(fieldnames{k})));
+                                res_struct.(fieldnames{k})(nanval)=-9999;
+                                %also checking for verry small value that mess up
+                                smallval=(abs(res_struct.(fieldnames{k}))<1.0e-20);
+                                res_struct.(fieldnames{k})(smallval)=0.0;
+                                if (size(res_struct.(fieldnames{k}),2)==num_of_timesteps),
+                                        fprintf(fid,'SCALARS %s float 1 \n',fieldnames{k});
+                                        fprintf(fid,'LOOKUP_TABLE default\n');
+                                        s='%e\n';
+                                        fprintf(fid,s,res_struct.(fieldnames{k})(1:end-1,timestep));
+                                else,
+                                        %forcing and results not on the same timestep,need some treatment
+                                        fprintf(fid,'SCALARS %s float 1 \n',fieldnames{k});
+                                        fprintf(fid,'LOOKUP_TABLE default\n');
+                                        index=1;
+                                        currenttime=((timestep-1)*outstep)+model.timestepping.start_time;
+                                        while (res_struct.(fieldnames{k})(end,index)<=currenttime);
+                                                if index==size(res_struct.(fieldnames{k}),2)
+                                                        break
+                                                end     
+                                                index=index+1;
+                      end
+                                        uptime=res_struct.(fieldnames{k})(end,index);
+                                        uplim=res_struct.(fieldnames{k})(1:end-1,index);
+                                        while (res_struct.(fieldnames{k})(end,index)>=currenttime);
+                                                if index==1
+                                                        break
+                              end
+                                                index=index-1;
+                      end
+                                        lowtime=res_struct.(fieldnames{k})(end,index);
+                                        lowlim=res_struct.(fieldnames{k})(1:end-1,index);
+                                        if uptime==currenttime,
+                                                interp=uplim;
+                                        elseif lowtime==currenttime,
+                                                interp=lowlim;
+                                        else
+                                                interp=lowlim+(uplim-lowlim)*((currenttime-lowtime)/(uptime-lowtime));
+                                        end
+                                        s='%e\n';
+                                        fprintf(fid,s,interp);
+                                end
+                  end           
                 end 
         end