-
Notifications
You must be signed in to change notification settings - Fork 10
/
BoundaryConditions.m
executable file
·74 lines (59 loc) · 1.97 KB
/
BoundaryConditions.m
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
classdef BoundaryConditions
%%
%
% Very simple class specifying the different types of Dirichlet boundary conditions.
%
%
%
%%
properties
% basal velocities
ubFixedNode=[];
ubFixedValue=[];
vbFixedNode=[];
vbFixedValue=[];
ubTiedNodeA=[];
ubTiedNodeB=[];
vbTiedNodeA=[];
vbTiedNodeB=[];
ubvbFixedNormalNode=[];
ubvbFixedNormalValue=[];
% deformational velocity, only used when solving SIA or hybrid
udFixedNode=[];
udFixedValue=[];
vdFixedNode=[];
vdFixedValue=[];
udTiedNodeA=[];
udTiedNodeB=[];
vdTiedNodeA=[];
vdTiedNodeB=[];
udvdFixedNormalNode=[];
udvdFixedNormalValue=[];
% thickness - only used if solving for thickness, for example in a uvh step
hFixedNode=[];
hFixedValue=[];
hTiedNodeA=[];
hTiedNodeB=[];
% positive thickness constraints - these are introduced automatically internally, do not
% prescribe directly!
hPosNode=[];
hPosValue=[];
hPosNodeActivated=[];
hPosNodeDeActivated=[] ;
% Boundary conditions for the level set field
% Only need when using the level-set method (currenlty under development, so don't
% use!)
LSFFixedNode=[];
LSFFixedValue=[];
LSFTiedNodeA=[];
LSFTiedNodeB=[];
LSFL=[] ;
LSFrhs=[] ;
% rate of thickness change - only used when using calculating dh/dt in combination with the ajoint methods
% for example when using measurements of dh/dt in an inversion
dhdtFixedNode=[];
dhdtFixedValue=[];
dhdtTiedNodeA=[];
dhdtTiedNodeB=[];
end
end