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

Index: ../trunk-jpl/src/m/contrib/larour/proj2shpprj.m
===================================================================
--- ../trunk-jpl/src/m/contrib/larour/proj2shpprj.m     (revision 24857)
+++ ../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);
Index: ../trunk-jpl/src/m/contrib/larour/applyqgisstyle.m
===================================================================
--- ../trunk-jpl/src/m/contrib/larour/applyqgisstyle.m  (nonexistent)
+++ ../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