Index: /issm/trunk-jpl/src/m/contrib/larour/applyqgisstyle.m
===================================================================
--- /issm/trunk-jpl/src/m/contrib/larour/applyqgisstyle.m	(revision 24858)
+++ /issm/trunk-jpl/src/m/contrib/larour/applyqgisstyle.m	(revision 24858)
@@ -0,0 +1,17 @@
+function applyqgisstyle(filename,type); 
+
+	%Getting extension out: 
+	[path,name,ext]=fileparts(filename);
+
+	%New name: 
+	if isempty(path),
+		newname=[name '.qml'];
+	else
+		newname=[path '/' name '.qml'];
+	end
+
+	if strcmpi(type,'mesh'),
+		system(['cp /Users/larour/issm-jpl/usr/larour/Qgis/Mesh.qml ' newname]);
+	else
+		error(['applyqgisstyle error message: ' type ' style not supported yet']);
+	end
Index: /issm/trunk-jpl/src/m/contrib/larour/proj2shpprj.m
===================================================================
--- /issm/trunk-jpl/src/m/contrib/larour/proj2shpprj.m	(revision 24857)
+++ /issm/trunk-jpl/src/m/contrib/larour/proj2shpprj.m	(revision 24858)
@@ -1,7 +1,11 @@
 function proj2shpprj(filename,proj)
 
+	%filename: does it have an extension? remove it. 
+	[path,name,ext] = fileparts(filename);
+
+	%Figure out central meridian and latitude of origin: 
 	[lat,long]=projlatlong(proj);
-	
-	fid=fopen(filename,'w');
-	fprintf('PROJCS["Lambert_Azimuthal_Equal_Area",GEOGCS["GCS_WGS_1984",DATUM["D_unknown",SPHEROID["WGS84",6378137,298.257223563]],PRIMEM["Greenwich",0],UNIT["Degree",0.017453292519943295]],PROJECTION["Lambert_Azimuthal_Equal_Area"],PARAMETER["latitude_of_origin",%g],PARAMETER["central_meridian",%g],PARAMETER["false_easting",0],PARAMETER["false_northing",0],UNIT["Meter",1]]',lat,long);
+
+	fid=fopen([path '/' name '.prj'],'w');
+	fprintf(fid,'PROJCS["Lambert_Azimuthal_Equal_Area",GEOGCS["GCS_WGS_1984",DATUM["D_unknown",SPHEROID["WGS84",6378137,298.257223563]],PRIMEM["Greenwich",0],UNIT["Degree",0.017453292519943295]],PROJECTION["Lambert_Azimuthal_Equal_Area"],PARAMETER["latitude_of_origin",%g],PARAMETER["central_meridian",%g],PARAMETER["false_easting",0],PARAMETER["false_northing",0],UNIT["Meter",1]]',lat,long);
 	fclose(fid);