Changeset 22815

Timestamp:

05/30/18 22:29:37 (7 years ago)

Author:

adhikari

Message:

CHG: final minor changes to all tutorials

Location:

issm/trunk-jpl/examples

Files:

• EsaGRACE/runme.m (modified) (4 diffs)
• EsaWahr/runme.m (modified) (1 diff)
• SlrFarrell/runme.m (modified) (5 diffs)

issm/trunk-jpl/examples/EsaGRACE/runme.m

@@ -1 +1 @@
 
 clear all;
-steps=[0:5]; % [1:5]; 
+addpath('../Data','../Functions');
+
+steps=[0]; % [0:5]; 
 
 if any(steps==0) % Download GRACE land_mass data {{{
@@ -12 +14 @@
    mget(f,'GRCTellus.JPL.200204_201701.LND.RL05_1.DSTvSCS1411.nc');
 
+        !mv *.nc '../Data/'
+        
         % display the content of the netcdf file. 
         %ncdisp('GRCTellus.JPL.200204_201701.LND.RL05_1.DSTvSCS1411.nc')
@@ -122 +126 @@
         sol3 = md.results.EsaSolution.EsaNmotion*1000; % [mm] 
         sol4 = md.results.EsaSolution.EsaEmotion*1000; % [mm] 
+
         sol_name={'Change in water equivalent height [cm]', 'Vertical displacement [mm]',...
                 'Horizontal (NS) displacement [mm]', 'Horizontal (EW) displacement [mm]'}; 
+        fig_name={'Fig_dH.pdf','Fig_vert.pdf','Fig_horzNS.pdf','Fig_horzEW.pdf'}; 
 
         res = 1.0; % degree 
@@ -160 +166 @@
                 set(0,'DefaultAxesFontSize',18,'DefaultAxesLineWidth',1,'DefaultTextFontSize',18,'DefaultLineMarkerSize',8)
                 figure1=figure('Position', [100, 100, 1000, 500]); 
-                gcf; 
-                load coast; 
-                cla; 
+                gcf; load coast; cla; 
                 pcolor(lon_grid,lat_grid,sol_grid); shading flat; hold on; 
+                % mask out oceans {{{
+                if (kk==1)
+                        geoshow(flipud(lat),flipud(long),'DisplayType','polygon','FaceColor','white'); 
+                end % }}}
                 plot(long,lat,'k'); hold off; 
+                % colormap, xlim etc {{{
                 c1=colorbar;
-                colormap(jet); 
+                colormap('haxby');
+                caxis([-min(abs(min(sol)),abs(max(sol))) min(abs(min(sol)),abs(max(sol)))]); 
                 xlim([-180 180]); 
                 ylim([-90 90]); 
+                % }}}
                 grid on; 
                 title(sol_name(kk)); 
                 set(gcf,'color','w');
-
-                %export_fig('Fig5.pdf'); 
+                
+                %export_fig(fig_name{kk}); 
         end
 

issm/trunk-jpl/examples/EsaWahr/runme.m

@@ -3 +3 @@
 addpath('../Functions');
 
-steps=[0:6]; % [0:6]; 
+steps=[0]; % [0:6]; 
 
 if any(steps==0) % Simple mesh creation {{{

issm/trunk-jpl/examples/SlrFarrell/runme.m

@@ -1 +1 @@
 
 clear all;
-steps=[5]; % 
+
+steps=[1]; % [1:5]; 
 
 if any(steps==1) % Global mesh creation {{{ 
@@ -128 +129 @@
         md.verbose=verbose('111111111');
 
+        % Choose different convergence threshold. [10% 1% 0.1%] to match Farrell 3 panels in Fig. 1
+        md.slr.reltol = 0.1/100; % per cent change in solution 
+
         % Solve 
         md=solve(md,'Slr');
@@ -141 +145 @@
         % solutions. 
         sol = md.results.SealevelriseSolution.Sealevel*100; % per cent normalized by GMSL (which 1 m)  
-        
-        res = 1.0; % degree 
+
+        res = 1; % degree 
 
         % Make a grid of lats and lons, based on the min and max of the original vectors
@@ -150 +154 @@
         % Make a interpolation object
         F = scatteredInterpolant(md.mesh.lat,md.mesh.long,sol); 
-        F.Method = 'linear';
-        F.ExtrapolationMethod = 'linear'; 
+        F.Method = 'natural'; % for smooth contour 
+        F.ExtrapolationMethod = 'none'; 
         
         % Do the interpolation to get gridded solutions... 
@@ -161 +165 @@
         figure1=figure('Position', [100, 100, 1000, 500]); 
         gcf; load coast; cla; 
-        %pcolor(lon_grid,lat_grid,sol_grid); shading flat; hold on;
-        contourf(lon_grid,lat_grid,sol_grid,[96 98 100 102 104 105]); shading flat; hold on;
-        geoshow(lat,long,'DisplayType','polygon','FaceColor','white'); 
+        pcolor(lon_grid,lat_grid,sol_grid); shading flat; hold on;
+        [C,h]=contour(lon_grid,lat_grid,sol_grid,[96 98 100 102 104 105],'-k','linewidth',2);
+        clabel(C,h,'FontSize',18,'Color','red','LabelSpacing',500); 
+        geoshow(lat,long,'DisplayType','polygon','FaceColor',[.82 .82 .82]); 
         plot(long,lat,'k'); hold off; 
         % define colormap, caxis, xlim etc {{{
         c1=colorbar;
-        %colormap('haxby'); 
-        %caxis([96 104]); 
-        xlim([-180 180]); 
-        ylim([-90 90]); 
+        colormap(flipud(haxby)); 
+        caxis([96 105]); 
+        xlim([-170 170]); 
+        ylim([-85 85]); 
         % }}} 
         grid on; 
-        title('Relative sea-level [mm]'); 
+        title('Relative sea-level [% of GMSL]'); 
         set(gcf,'color','w');