Changeset 26286
- Timestamp:
- 05/21/21 15:41:39 (4 years ago)
- Location:
- issm/trunk-jpl/examples
- Files:
-
- 2 edited
Legend:
- Unmodified
- Added
- Removed
-
issm/trunk-jpl/examples/EsaGRACE/runme.m
r26258 r26286 2 2 addpath('../Data','../Functions'); 3 3 4 steps=[ 5];4 steps=[1]; 5 5 6 if any(steps==1) 6 if any(steps==1) 7 7 disp(' Step 1: Global mesh creation'); 8 8 … … 18 18 19 19 plotmodel (md,'data',md.mask.ocean_levelset,'edgecolor','k'); 20 end 20 end 21 21 22 if any(steps==2) 22 if any(steps==2) 23 23 disp(' Step 2: Define loads in meters of ice height equivalent'); 24 24 md = loadmodel('./Models/EsaGRACE_Mesh'); … … 36 36 'view',[90 -90],'caxis',[-.1 .1],... 37 37 'title','Ice height equivalent [m]'); 38 end 38 end 39 39 40 if any(steps==3) 40 if any(steps==3) 41 41 disp(' Step 3: Parameterization'); 42 42 md = loadmodel('./Models/EsaGRACE_Loads'); … … 65 65 66 66 save ./Models/EsaGRACE_Parameterization md; 67 end 67 end 68 68 69 if any(steps==4) 69 if any(steps==4) 70 70 disp(' Step 4: Solve Esa solver'); 71 71 md = loadmodel('./Models/EsaGRACE_Parameterization'); … … 79 79 80 80 save ./Models/EsaGRACE_Solution md; 81 end 81 end 82 82 83 if any(steps==5) 83 if any(steps==5) 84 84 disp(' Step 5: Plot solutions'); 85 85 md = loadmodel('./Models/EsaGRACE_Solution'); … … 140 140 %export_fig(fig_name{kk}); 141 141 end 142 end 142 end -
issm/trunk-jpl/examples/EsaWahr/runme.m
r26258 r26286 2 2 addpath('../Functions'); 3 3 4 steps=[ 2];4 steps=[1]; 5 5 6 if any(steps==1) 6 if any(steps==1) 7 7 disp(' Step 1: Mesh creation'); 8 8 … … 12 12 13 13 plotmodel(md,'data','mesh'); 14 end 14 end 15 15 16 if any(steps==2) 16 if any(steps==2) 17 17 disp(' Step 2: Anisotropic mesh creation'); 18 18 … … 28 28 29 29 plotmodel (md,'data','mesh'); 30 end 30 end 31 31 32 if any(steps==3) 32 if any(steps==3) 33 33 disp(' Step 3: Define loads'); 34 34 md = loadmodel('./Models/EsaWahr_Mesh'); … … 48 48 49 49 plotmodel (md,'data',md.esa.deltathickness,'title','Ice height equivalent [m]'); 50 end 50 end 51 51 52 if any(steps==4) 52 if any(steps==4) 53 53 disp(' Step 4: Parameterization'); 54 54 md = loadmodel('./Models/EsaWahr_Loads'); … … 74 74 75 75 save ./Models/EsaWahr_Parameterization md; 76 end 76 end 77 77 78 if any(steps==5) 78 if any(steps==5) 79 79 disp(' Step 5: Solve Esa solver'); 80 80 md = loadmodel('./Models/EsaWahr_Parameterization'); … … 88 88 89 89 save ./Models/EsaWahr_Solution md; 90 end 90 end 91 91 92 if any(steps==6) 92 if any(steps==6) 93 93 disp(' Step 6: Plot solutions'); 94 94 md = loadmodel('./Models/EsaWahr_Solution'); … … 97 97 horz_n = md.results.EsaSolution.EsaYmotion*1000; % [mm] 98 98 horz_e = md.results.EsaSolution.EsaXmotion*1000; % [mm] 99 horz = sqrt(horz_n.^2+horz_e.^2); % [mm] 99 horz = sqrt(horz_n.^2+horz_e.^2); % [mm] 100 100 101 101 set(0,'DefaultAxesFontSize',24,'DefaultAxesLineWidth',1,'DefaultTextFontSize',24,'DefaultLineMarkerSize',6); … … 121 121 'colorbarpos',[0.53,0.065,0.18,0.02],'colorbartitle#4','East-west [mm]'); 122 122 %export_fig('Fig5.pdf'); 123 end 123 end 124 124 125 if any(steps==7) 125 if any(steps==7) 126 126 disp(' Step 7: Compare results against Wahr semi-analytic solutions'); 127 127 md = loadmodel('./Models/EsaWahr_Solution'); … … 130 130 horz_n = md.results.EsaSolution.EsaYmotion*1000; % [mm] 131 131 horz_e = md.results.EsaSolution.EsaXmotion*1000; % [mm] 132 horz = sqrt(horz_n.^2+horz_e.^2); % [mm] 132 horz = sqrt(horz_n.^2+horz_e.^2); % [mm] 133 133 134 134 xi=[0:500:100000]; % grid points [m] … … 162 162 set(gcf,'color','w'); 163 163 %export_fig('Fig6.pdf'); 164 end 164 end
Note:
See TracChangeset
for help on using the changeset viewer.