Changeset 21942

Timestamp:

08/11/17 16:15:14 (8 years ago)

Author:

adhikari

Message:

CHG: further renaming of the txt files

File:

• issm/trunk-jpl/src/m/coordsystems/gmtmask.m (modified) (4 diffs)

issm/trunk-jpl/src/m/coordsystems/gmtmask.m

@@ -31 +31 @@
         %First, write our lat,long file for gmt:
         nv=length(lat);
-        filename = ['all_vertices-' num2str(feature('GetPid')) '.txt']
-        dlmwrite(filename,[long lat (1:nv)'],'delimiter','\t','precision',10);
-
+        filename_all = ['all_vertices-' num2str(feature('GetPid')) '.txt']; 
+        filename_oce = ['oce_vertices-' num2str(feature('GetPid')) '.txt']; 
+        dlmwrite(filename_all,[long lat (1:nv)'],'delimiter','\t','precision',10);
 
         %Avoid bypassing of the ld library path by Matlab (:()
@@ -42 +42 @@
 
         %figure out which vertices are on the ocean, which one on the continent:
-        [status,result] = system(['gmt gmtselect ./' filename ' -h0 -Df -R0/360/-90/90  -A0 -JQ180/200 -Nk/s/s/k/s > ./oce_vertices.txt']);
+        [status,result] = system(['gmt gmtselect ./' filename_all ' -h0 -Df -R0/360/-90/90  -A0 -JQ180/200 -Nk/s/s/k/s > ./' filename_oce]);
         if status~=0,
                 error(result);
@@ -52 +52 @@
         end
         %read the con_vertices.txt file and flag our mesh vertices on the continent
-        fid=fopen('./oce_vertices.txt','r');
+        fid=fopen(['./' filename_oce],'r');
         line=fgets(fid); 
         line=fgets(fid);
@@ -66 +66 @@
         mask(oce_vertices)=1;
         
-        system(['rm -rf ./' filename ' ./oce_vertices.txt ./gmt.history']);
+        system(['rm -rf ./' filename_all ' ./' filename_oce ' ./gmt.history']);
         if ~recursive, disp(sprintf('gmtmask: done')); end;