Index: /issm/trunk-jpl/src/m/plot/plot_googlemaps.m
===================================================================
--- /issm/trunk-jpl/src/m/plot/plot_googlemaps.m	(revision 15168)
+++ /issm/trunk-jpl/src/m/plot/plot_googlemaps.m	(revision 15169)
@@ -46,6 +46,5 @@
 end
 
-%Retrieve image from md
-[X Y] = meshgrid(md.radaroverlay.x,md.radaroverlay.y);
+%Process image from model
 final = double(md.radaroverlay.pwr)/double(max(md.radaroverlay.pwr(:))); %rescale between 0 and 1
 
@@ -54,5 +53,18 @@
 
 %Prepare grid
-data_grid=InterpFromMeshToMesh2d(md.mesh.elements,md.mesh.x,md.mesh.y,data,X(:),Y(:),'default',NaN); data_grid=reshape(data_grid,size(X));
+if size(md.radaroverlay.x,1)==1 | size(md.radaroverlay.x,2)==1,
+	xmin=min(md.radaroverlay.x);
+	ymax=max(md.radaroverlay.y);
+	xspacing=md.radaroverlay.x(2)-md.radaroverlay.x(1);
+	yspacing=md.radaroverlay.y(2)-md.radaroverlay.y(1);
+	nlines=length(md.radaroverlay.y);
+	ncols =length(md.radaroverlay.x);
+		[x_m y_m data_grid]=InterpFromMeshToGrid(elements,x/getfieldvalue(options,'unit',1),y/getfieldvalue(options,'unit',1),...
+			data,xmin,ymax,xspacing,yspacing,nlines,ncols,NaN);
+else
+	[X Y] = meshgrid(md.radaroverlay.x,md.radaroverlay.y);
+	data_grid=InterpFromMeshToMesh2d(md.mesh.elements,md.mesh.x,md.mesh.y,data,X(:),Y(:),'default',NaN); data_grid=reshape(data_grid,size(X));
+	x_m=X(1,:); y_m=Y(:,1);
+end
 data_nan=isnan(data_grid);
 if exist(options,'caxis'),
@@ -69,5 +81,5 @@
 image_rgb = ind2rgb(uint16((data_grid - data_min)*(length(colorm)/(data_max-data_min))),colorm);
 
-alpha=ones(size(X));
+alpha=ones(size(data_grid));
 alpha(find(~data_nan))=transparency;
 alpha=repmat(alpha,[1 1 3]);
@@ -78,5 +90,5 @@
 subplotmodel(plotlines,plotcols,i,options);
 
-h=imagesc(X(1,:),Y(:,1),final);
+h=imagesc(x_m,y_m,final);
 
 %last step: mesh gridded?