Changeset 22701

Timestamp:

04/19/18 19:57:49 (7 years ago)

Author:

Mathieu Morlighem

Message:

CHG: added dem field if dem is not the same as field

File:

• issm/trunk-jpl/src/m/plot/plot_gridded.m (modified) (1 diff)

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

@@ -51 +51 @@
 image_rgb = ind2rgb(uint16((data_grid - data_min)*(length(map)/(data_max-data_min))),map);
 if exist(options,'shaded'),
+
+        if exist(options,'dem'),
+                dem_grid=InterpFromMeshToGrid(elements,x,y,getfieldvalue(options,'dem'),x_m,y_m,NaN);
+        else
+                dem_grid=data_grid_save;
+        end
         a    = -45;
         scut = 0.2;
         c    = 1;
         % computes lighting from elevation gradient
-        [fx,fy] = gradient(data_grid_save,x_m,y_m);
+        [fx,fy] = gradient(dem_grid,x_m,y_m);
         fxy = -fx*sind(a) - fy*cosd(a);
         clear fx fy % free some memory...