Index: /issm/trunk-jpl/src/m/plot/googlemaps.m
===================================================================
--- /issm/trunk-jpl/src/m/plot/googlemaps.m	(revision 15154)
+++ /issm/trunk-jpl/src/m/plot/googlemaps.m	(revision 15155)
@@ -26,8 +26,17 @@
 	xlim=getfieldvalue(options,'xlim',[min(md.mesh.x) max(md.mesh.x)]);
 	ylim=getfieldvalue(options,'ylim',[min(md.mesh.y) max(md.mesh.y)]);
-	[latlist lonlist]= xy2ll(...
-		[linspace(xlim(1),xlim(2),100) linspace(xlim(2),xlim(2),100) linspace(xlim(2),xlim(1),100) linspace(xlim(1),xlim(1),100)],...
-		[linspace(ylim(1),ylim(1),100) linspace(ylim(1),ylim(2),100) linspace(ylim(2),ylim(2),100) linspace(ylim(2),ylim(1),100)],...
-		+1,45,70);
+	if strcmpi(md.mesh.hemisphere,'n'),
+		[latlist lonlist]= xy2ll(...
+			[linspace(xlim(1),xlim(2),100) linspace(xlim(2),xlim(2),100) linspace(xlim(2),xlim(1),100) linspace(xlim(1),xlim(1),100)],...
+			[linspace(ylim(1),ylim(1),100) linspace(ylim(1),ylim(2),100) linspace(ylim(2),ylim(2),100) linspace(ylim(2),ylim(1),100)],...
+			+1,45,70);
+	elseif strcmpi(md.mesh.hemisphere,'s'),
+		[latlist lonlist]= xy2ll(...
+			[linspace(xlim(1),xlim(2),100) linspace(xlim(2),xlim(2),100) linspace(xlim(2),xlim(1),100) linspace(xlim(1),xlim(1),100)],...
+			[linspace(ylim(1),ylim(1),100) linspace(ylim(1),ylim(2),100) linspace(ylim(2),ylim(2),100) linspace(ylim(2),ylim(1),100)],...
+			-1,0,71);
+	else
+		error('field hemisphere should either be ''n'' or ''s''');
+	end
 
 	%Image corners in lat/long
@@ -76,4 +85,6 @@
 		position = [num2str(latn) ',' num2str(lonn)];
 		disp(['Google Earth tile: ' num2str(x) '/' num2str(cols-1) ' ' num2str(y) '/' num2str(rows-1) ' (center: ' position ')']);
+		%Google maps API
+		%http://developers.google.com/maps/documentation/staticmaps/
 		params = [...
 			'center=' position ...
@@ -97,5 +108,11 @@
 [gX gY]=meshgrid(ulx:ulx+size(final,2)-1,uly:-1:uly-size(final,1)+1);
 [LAT LON]=pixelstolatlon(gX,gY, zoom);
-[X Y]=ll2xy(LAT,LON,+1,45,70);
+if strcmpi(md.mesh.hemisphere,'n'),
+	[X Y]=ll2xy(LAT,LON,+1,45,70);
+elseif strcmpi(md.mesh.hemisphere,'s'),
+	[X Y]=ll2xy(LAT,LON,-1,0,71);
+else
+	error('field hemisphere should either be ''n'' or ''s''');
+end
 
 md.radaroverlay.pwr=final;
Index: /issm/trunk-jpl/src/m/plot/plot_googlemaps.m
===================================================================
--- /issm/trunk-jpl/src/m/plot/plot_googlemaps.m	(revision 15154)
+++ /issm/trunk-jpl/src/m/plot/plot_googlemaps.m	(revision 15155)
@@ -25,8 +25,17 @@
 	xlim=getfieldvalue(options,'xlim',[min(x) max(x)]);
 	ylim=getfieldvalue(options,'ylim',[min(y) max(y)]);
-	[latlist lonlist]= xy2ll(...
-		[linspace(xlim(1),xlim(2),100) linspace(xlim(2),xlim(2),100) linspace(xlim(2),xlim(1),100) linspace(xlim(1),xlim(1),100)],...
-		[linspace(ylim(1),ylim(1),100) linspace(ylim(1),ylim(2),100) linspace(ylim(2),ylim(2),100) linspace(ylim(2),ylim(1),100)],...
-		+1,45,70);
+	if strcmpi(md.mesh.hemisphere,'n'),
+		[latlist lonlist]= xy2ll(...
+			[linspace(xlim(1),xlim(2),100) linspace(xlim(2),xlim(2),100) linspace(xlim(2),xlim(1),100) linspace(xlim(1),xlim(1),100)],...
+			[linspace(ylim(1),ylim(1),100) linspace(ylim(1),ylim(2),100) linspace(ylim(2),ylim(2),100) linspace(ylim(2),ylim(1),100)],...
+			+1,45,70);
+	elseif strcmpi(md.mesh.hemisphere,'s'),
+		[latlist lonlist]= xy2ll(...
+			[linspace(xlim(1),xlim(2),100) linspace(xlim(2),xlim(2),100) linspace(xlim(2),xlim(1),100) linspace(xlim(1),xlim(1),100)],...
+			[linspace(ylim(1),ylim(1),100) linspace(ylim(1),ylim(2),100) linspace(ylim(2),ylim(2),100) linspace(ylim(2),ylim(1),100)],...
+			-1,0,71);
+	else
+		error('field hemisphere should either be ''n'' or ''s''');
+	end
 
 	%Image corners in lat/long