Update SRTM to 2.3 (#252)

* Update setup.sh

* Update version of SRTM

* changing SRTM15+ filename to remove version number so as to avoid needing to update. Also use wildcard in wget ftp download to automatically get the newest version

* removing unnecessary commited files

Co-authored-by: William Pringle <[email protected]>

Author: Keith Roberts

.gitignore

@@ -3,7 +3,7 @@ datasets/
 fort.*
 Examples/*
 !Examples/Example_*
-/*.mat
+*/*.mat
 *.15
 *.14
 *.13

Examples/Example_11_Remeshing_Patches.m

@@ -38,7 +38,7 @@
 % Get out the msh class and put on nodestrings
 m = mshopts.grd;
 % Interpolate topobathy data onto the vertices of the mesh.
-m = interp(m,'SRTM15+V2.1.nc');
+m = interp(m,'SRTM15+.nc');
 %% Extract the region which you want to remesh
 % For the purpose of this example, we have drawn a random polygon on top
 % of the mesh for which we would like to remesh.

Examples/Example_12_Riverine_flow.m

@@ -43,7 +43,7 @@
 
 %% STEP 2: specify geographical datasets and process the geographical data
 %% to be used later with other OceanMesh classes...
-dem       = 'SRTM15+V2.1.nc';
+dem       = 'SRTM15+.nc';
 coastline = 'GSHHS_f_L1';
 gdat = geodata('shp',coastline,'dem',dem,'h0',min_el,'bbox',bbox);
 %plot(gdat,'shp'); 

Examples/Example_3_ECGC.m

@@ -24,7 +24,7 @@
 
 %% STEP 2: specify geographical datasets and process the geographical data
 %% to be used later with other OceanMesh classes...
-dem       = 'SRTM15+V2.1.nc';
+dem       = 'SRTM15+.nc';
 coastline = 'GSHHS_f_L1';
 gdat1 = geodata('shp',coastline,'dem',dem,'h0',min_el,...
                 'bbox',bbox);

Examples/Example_4_PRVI.m

@@ -24,7 +24,7 @@
 max_el = 10e3;       	% maximum resolution in meters. 
 
 coastline = 'GSHHS_f_L1';
-dem       = 'SRTM15+V2.1.nc';
+dem       = 'SRTM15+.nc';
 gdat{1} = geodata('shp',coastline,...
                   'dem',dem,...
                   'bbox',bbox,...

Examples/Example_7_Global.m

@@ -22,7 +22,7 @@
 
 %% STEP 2: specify geographical datasets and process the geographical data
 %% to be used later with other OceanMesh classes...
-dem       = 'SRTM15+V2.1.nc';
+dem       = 'SRTM15+.nc';
 coastline = {'GSHHS_f_L1','GSHHS_f_L6'};
 gdat1 = geodata('shp',coastline,'dem',dem,...
                 'bbox',bbox,'h0',min_el);

Examples/Example_8_AK.m

@@ -25,7 +25,7 @@
 
 %% STEP 2: specify geographical datasets and process the geographical data 
 %%         to be used later with other OceanMesh classes...
-dem       = 'SRTM15+V2.1.nc';
+dem       = 'SRTM15+.nc';
 coastline = 'GSHHS_f_L1';
 gdat = geodata('shp',coastline,'bbox',bbox,'h0',min_el);
 % plotting the gdat to show it crosses the 180/-180 and 

Examples/Example_9_TBAY.m

@@ -31,7 +31,7 @@
 % Shoreline
 shoreline = 'GSHHS_f_L1';
 % Bathy data
-dem = 'SRTM15+V2.1.nc';
+dem = 'SRTM15+.nc';
 %% Define meshing domain and parameters
 bbox{1} = [-83, -82; 27 28.5];      % outermost bbox
 min_el = 100;                       % minimum edgelength

Tests/RunTests.m

@@ -18,13 +18,13 @@
 
    TestInterp
 
-   if exist('SRTM15+V2.1.nc')
+   if exist('SRTM15+.nc')
 
       TestECGC 
 
    else
 
-      warning('Need to download SRTM15+V2.1.nc to run TestECGC. Run setup.sh') 
+      warning('Need to download SRTM15+.nc to run TestECGC. Run setup.sh') 
 
    end
 

Tests/TestECGC.m

@@ -20,7 +20,7 @@
 R         = 3; 			        % Number of elements to resolve feature.
 
 
-dem       = 'SRTM15+V2.1.nc';
+dem       = 'SRTM15+.nc';
 coastline = 'GSHHS_f_L1';
 gdat1 = geodata('shp',coastline,'dem',dem,'h0',min_el,...
                 'bbox',bbox);