-
Notifications
You must be signed in to change notification settings - Fork 2
Expand file tree
/
Copy pathRunAllExamples.m
More file actions
122 lines (76 loc) · 3.98 KB
/
Copy pathRunAllExamples.m
File metadata and controls
122 lines (76 loc) · 3.98 KB
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
%% Master 17 March 2024, 01 August 2024, 23 Dec 2025
cd RadialIceCap\
Ua
cd ..
close all
cd 1dIceShelf
Ua
cd ..
close all
cd 1dIceStream
Ua
cd ..
close all
cd Crack
Ua
cd ..
cd GaussPeak
Ua
cd ..
close all
cd IceShelf
Ua
cd ..
close all
cd MismipPlus
Ua
cd ..
cd IceBerg
Ua
cd ..
cd PIG-TWG
UserVar.RunType='Inverse-MatOpt'; Ua(UserVar) ; % 01/03/2025, 19/07/2025
UserVar.RunType='Inverse-UaOpt' ; Ua(UserVar) ; % 01/03/2025, 19/07/2025
UserVar.RunType='Inverse-MatOpt' ; CtrlVar.Inverse.MinimisationMethod="MatlabOptimization-HessianBased"; Ua(UserVar,CtrlVar) ; % 01/03/2025, 19/07/2025
Klear
UserVar.RunType='Inverse-MatOpt' ; CtrlVar.Inverse.MinimisationMethod="MatlabOptimization-GradientBased"; Ua(UserVar,CtrlVar) ; % no-longer working as of Matlab 2021b...?
% The reasons are unclear, but for the time being gradient-based optimization
% with the Matlab toolbox can not be done with MATLAB R2022a.
% This is not too much of an issue as the default option is the Hessian-based approach anyhow,
% which is also the better option.
% But as of R2024b, and possibly earlier, this is again working... 01/03/2024
UserVar.RunType='TestingMeshOptions' ; Ua(UserVar) ; % 01/03/2025, 19/07/2025
cd ..
close all
cd MassBalanceFeedback
Ua
cd ..
%%
cd Calving
% A few examples:
% 1) Here calving is implemented by the user using a user-defined mass balance feedback as defined in DefineMassBalance.m
% This does not involve the level-set method implementation in Ua
UserVar.RunType="Test-1dAnalyticalIceShelf-CalvingThroughMassBalanceFeedback-"; Ua(UserVar) ;
% 2) Level-set is prescribed directly (in DefineCalving.m). This example shows how the user can
% define directly/manually the position of the calving front over time.
% No calving law is used and the calving rate is not specified.
UserVar.RunType="Test-1dAnalyticalIceShelf-CalvingThroughPrescribedLevelSet-" ; Ua(UserVar) ;
% 3) Similar to case 1) above but done for a MismipPlus style geometry.
UserVar.RunType="Test-CalvingThroughMassBalanceFeedback-"; Ua(UserVar) ;
% 4) Similar to case 2) above but done for a MismipPlus style geometry.
UserVar.RunType="Test-CalvingThroughPrescribedLevelSet-" ; Ua(UserVar) ;
% 5) Calving through element deactivation: Here calving is simulated by deactivating elements (done in
% DefineElementsToDeactivate.m). No level-set calculations/initialization are required.
%
UserVar.RunType="Test-ManuallyDeactivateElements-" ; Ua(UserVar) ;
cd ..
%%
cd Inverse
Ua
cd ..
close all
%%
cd Greenland\
Ua
cd ..
%%