source: issm/oecreview/Archive/24684-25833/ISSM-24857-24858.diff

../trunk-jpl/src/m/contrib/larour/proj2shpprj.m

@@ -1 +1 @@
 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);

../trunk-jpl/src/m/contrib/larour/applyqgisstyle.m

@@ +1 @@
+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