Index: /issm/trunk/src/m/classes/public/plot/plot_overlay.m
===================================================================
--- /issm/trunk/src/m/classes/public/plot/plot_overlay.m	(revision 2724)
+++ /issm/trunk/src/m/classes/public/plot/plot_overlay.m	(revision 2725)
@@ -31,5 +31,5 @@
 if redo,
 	t1=clock; fprintf('%s','Extracting radar image...');
-	md=radarpower(md,xlim,ylim,getfieldvalue(options,'highres',0));
+	md=radarpower(md,getfieldvalue(options,'hem','s'),xlim,ylim,getfieldvalue(options,'highres',0));
 	t2=clock;fprintf('%s\n',[' done (' num2str(etime(t2,t1)) ' seconds)']);
 end
@@ -75,5 +75,5 @@
 
 %intensity
-v=radar/max(max(radar));
+v=radar/max(radar(:));
 
 %hue
@@ -87,5 +87,5 @@
 
 %Include border
-v((1:border),:)=0;  v((end-border+1:end),:)=0; v(:,1:border)=0;v(:,(end-border+1:end))=0;
+%v((1:border),:)=0;  v((end-border+1:end),:)=0; v(:,1:border)=0;v(:,(end-border+1:end))=0;
 
 %Transform hsv to rgb
Index: /issm/trunk/src/m/classes/public/plot/plotdoc.m
===================================================================
--- /issm/trunk/src/m/classes/public/plot/plotdoc.m	(revision 2724)
+++ /issm/trunk/src/m/classes/public/plot/plotdoc.m	(revision 2725)
@@ -85,4 +85,5 @@
 disp('       ''overlay'': yes or no. This will overlay a radar amplitude image behind');
 disp('       ''highres'': resolution of overlayed radar amplitude image (default is 0, high resolution is 1).');
+disp('       ''hem'': specify hemisphere ''n'' or ''s'' (default is ''s'').');
 disp('       ''alpha'': transparency coefficient (the higher, the more transparent). Default is 1.5');
 disp('       ''scaling'': scaling factor used by quiver plots. Default is 0.4');
Index: /issm/trunk/src/m/classes/public/radarpower.m
===================================================================
--- /issm/trunk/src/m/classes/public/radarpower.m	(revision 2724)
+++ /issm/trunk/src/m/classes/public/radarpower.m	(revision 2725)
@@ -1,3 +1,3 @@
-function md=radarpower(md,xlim,ylim,highres)
+function md=radarpower(md,hemisphere,xlim,ylim,highres)
 %RADARPOWER - overlay a power radar image on an existing mesh
 %
@@ -7,5 +7,5 @@
 %
 %   Usage:
-%      md=radarpower(md,xlim,ylim,highres)
+%      md=radarpower(md,hemisphere,xlim,ylim,highres)
 
 global ISSM_DIR
@@ -31,27 +31,51 @@
 
 %the geotiff image is either 200m or 1km accuracy. 
-if highres,
-	if ~exist([MODELDATA '/MosaicTiffRsat/amm125m_v2_200m.tif']),
-		error(['radarpower error message: file ' MODELDATA '/MosaicTiffRsat/amm125m_v2_200m.tif not found. Check MODELDATA variable..']);
+if strcmpi(hemisphere,'n'),
+	if ~exist([MODELDATA '/MOG/mog150_greenland_map.jpg']),
+		error(['radarpower error message: file ' MODELDATA '/MOG/mog150_greenland_map.jpg not found. Check MODELDATA variable..']);
 	end
-	geotiff_name=[MODELDATA '/MosaicTiffRsat/amm125m_v2_200m.tif'];
+	jpgim=[MODELDATA '/MOG/mog150_greenland_map.jpg'];
+	geom=load([MODELDATA '/MOG/mog150_greenland_map.jpgw'],'ascii');
+	sizex=floor((x1-x0)/geom(1)); % x posting
+	sizey=floor((y1-y0)/geom(4)); % y posting
+	topleftx=floor((x0-geom(5))/geom(1)); % x min
+	toplefty=floor((geom(6)-y1)/geom(4)); % y max
+
+	%Read and crop file
+	tic
+	im=imread(jpgim);
+	toc
+	tic
+	im=im(toplefty:toplefty+sizey,topleftx:topleftx+sizex);
+	toc
+	md.sarpwr=double(flipud(im));
+	md.sarxm=(x0:(x1-x0)/(size(md.sarpwr,2)-1):x1);
+	md.sarym=(y0:(y1-y0)/(size(md.sarpwr,1)-1):y1);
+
 else
-	if ~exist([MODELDATA '/MosaicTiffRsat/amm125m_v2_1km.tif']),
-		error(['radarpower error message: file ' MODELDATA '/MosaicTiffRsat/amm125m_v2_1km.tif not found. Check MODELDATA variable..']);
+	if highres,
+		if ~exist([MODELDATA '/MosaicTiffRsat/amm125m_v2_200m.tif']),
+			error(['radarpower error message: file ' MODELDATA '/MosaicTiffRsat/amm125m_v2_200m.tif not found. Check MODELDATA variable..']);
+		end
+		geotiff_name=[MODELDATA '/MosaicTiffRsat/amm125m_v2_200m.tif'];
+	else
+		if ~exist([MODELDATA '/MosaicTiffRsat/amm125m_v2_1km.tif']),
+			error(['radarpower error message: file ' MODELDATA '/MosaicTiffRsat/amm125m_v2_1km.tif not found. Check MODELDATA variable..']);
+		end
+		geotiff_name=[MODELDATA '/MosaicTiffRsat/amm125m_v2_1km.tif'];
 	end
-	geotiff_name=[MODELDATA '/MosaicTiffRsat/amm125m_v2_1km.tif'];
+
+	%Name of image
+	size(imread(geotiff_name,'TIFF'))
+	inputname='./temp.tif';
+	system([path_gdal 'gdal_translate -quiet -projwin ' num2str(x0) ' ' num2str(y1) ' ' num2str(x1) ' ' num2str(y0) ' ' geotiff_name ' ' inputname ]);
+
+	%Read in temp.tif:
+	md.sarpwr=double(flipud(imread('temp.tif','TIFF')));
+	md.sarxm=(x0:(x1-x0)/(size(md.sarpwr,2)-1):x1);
+	md.sarym=(y0:(y1-y0)/(size(md.sarpwr,1)-1):y1);
+
+	%Erase image
+	system('rm -rf ./temp.tif');
+
 end
-
-%Name of image
-inputname='./temp.tif';
-
-command=['system(''' path_gdal 'gdal_translate -quiet -projwin ' num2str(x0) ' ' num2str(y1) ' ' num2str(x1) ' ' num2str(y0) ' ' geotiff_name ' ' inputname ''');'];
-eval(command);
-
-%Read in temp.tif:
-md.sarpwr=double(flipud(imread('temp.tif','TIFF')));
-md.sarxm=(x0:(x1-x0)/(size(md.sarpwr,2)-1):x1);
-md.sarym=(y0:(y1-y0)/(size(md.sarpwr,1)-1):y1);
-
-%Erase image
-system('rm -rf ./temp.tif');