Commit 3e861c6e authored by Benoit Urruty's avatar Benoit Urruty

mesh_km

parent f16437f2
......@@ -241,7 +241,7 @@
},
{
"cell_type": "code",
"execution_count": 3,
"execution_count": 6,
"metadata": {},
"outputs": [
{
......@@ -381,24 +381,23 @@
"metadata": {},
"outputs": [],
"source": [
"# mask_out_error=np.loadtxt('../../../backup.txt',delimiter=',')\n",
"# # mask_out_error[index]=np.nan\n",
"# # mask_out_error[index2]=np.nan\n",
"# # np.savetxt('backup.txt',mask_out_error,delimiter=',')\n",
"\n",
"# fig=plt.figure(figsize=[12,10])\n",
"\n",
"# plt.rcParams.update({'font.size': 15})\n",
"# ax=plt.gca()\n",
"# sf.plot(ax=ax, color='k')\n",
"# im=plt.pcolormesh(x,y,mask_out_error)\n",
"# cbar=plt.colorbar(im)\n",
"# cbar.set_label('1 : overlap - 0 : outside')\n",
"\n",
"# # plt.title('link between the two beds +/- uncertainty')\n",
"# plt.tight_layout()\n",
"# plt.savefig(figure + 'uncertainty.png')\n",
"# del mask_out_error"
"mask_out_error=np.loadtxt('../../../backup.txt',delimiter=',')\n",
"# mask_out_error[index]=np.nan\n",
"# mask_out_error[index2]=np.nan\n",
"# np.savetxt('backup.txt',mask_out_error,delimiter=',')\n",
"\n",
"fig=plt.figure(figsize=[12,10])\n",
"\n",
"plt.rcParams.update({'font.size': 15})\n",
"ax=plt.gca()\n",
"sf.plot(ax=ax, color='k')\n",
"im=plt.pcolormesh(x,y,mask_out_error)\n",
"cbar=plt.colorbar(im)\n",
"cbar.set_label('1 : overlap - 0 : outside')\n",
"\n",
"# plt.title('link between the two beds +/- uncertainty')\n",
"plt.tight_layout()\n",
"plt.savefig(figure + 'uncertainty.svg')\n"
]
},
{
......@@ -406,7 +405,9 @@
"execution_count": 9,
"metadata": {},
"outputs": [],
"source": []
"source": [
"del mask_out_error"
]
},
{
"cell_type": "markdown",
FUNCTION HMin(Model,nodenumber,VarIn) RESULT(VarOut)
USE DefUtils
implicit none
!-----------------
TYPE(Model_t) :: Model
INTEGER :: nodenumber
REAL(kind=dp) :: VarIn,VarOut
REAL(kind=dp) :: Hc
LOGICAL :: Found
Hc = GetCReal( Model % Constants, 'HMin Critical Thickness', Found )
IF (.NOT.Found) &
CALL FATAL('HMin','<HMin Critical Thickness> not found')
VarOut = MAX(VarIn,Hc)
End FUNCTION HMin
146266 277434 15096
2
202 15096
303 277434
2164 2 1 2
2438 2 1 2
2885 2 1 7
5996 2 1 10
7384 2 1 5
7745 2 1 2
7830 2 1 3
9400 2 1 5
9676 2 1 5
9692 2 1 7
9959 2 1 5
11501 2 1 3
12034 2 1 5
12519 2 1 5
12689 2 1 2
13297 2 1 7
14260 2 1 7
14466 2 1 10
14521 2 1 10
14548 2 1 4
15669 3 1 3 2
16768 2 1 10
16831 2 1 5
16907 2 1 3
17064 2 1 10
17380 2 1 6
17735 2 1 10
18068 2 1 10
19603 2 1 6
19942 2 1 6
21465 2 1 7
21511 2 1 2
22181 2 1 7
22313 2 1 10
23019 2 1 7
23467 2 1 3
23644 2 1 2
26023 2 1 10
27365 2 1 7
28765 2 1 3
29532 2 1 2
31009 2 1 7
31108 2 1 5
31905 2 1 5
35005 2 1 3
35006 2 1 10
35706 2 1 3
36201 2 1 10
42869 2 1 6
44051 2 1 7
45149 2 1 2
46366 2 1 7
47488 2 1 2
51621 2 1 3
52831 2 1 2
52982 2 1 10
58623 2 1 5
59767 2 1 3
60090 2 1 5
61744 2 1 2
64445 2 1 7
65898 2 1 5
67687 2 1 2
70001 2 1 7
70541 2 1 7
71305 2 1 6
74435 2 1 7
74750 2 1 7
75843 2 1 5
79646 3 1 7 6
81526 2 1 10
82127 2 1 10
83106 2 1 6
83394 2 1 7
83762 2 1 5
83982 2 1 5
85075 2 1 10
85760 2 1 7
86485 2 1 7
86640 2 1 5
87185 2 1 7
87985 2 1 3
90498 2 1 7
91681 2 1 5
91885 3 1 4 2
92721 2 1 3
93446 2 1 3
93552 2 1 2
93632 2 1 5
95323 2 1 6
96067 2 1 2
96876 2 1 5
97444 2 1 6
102109 2 1 7
102614 2 1 5
103018 2 1 6
103026 2 1 7
103965 2 1 3
104240 2 1 5
104294 2 1 3
104354 2 1 5
105214 2 1 5
107934 2 1 5
108467 2 1 3
108760 2 1 7
109003 2 1 2
109535 2 1 3
109539 2 1 10
109766 2 1 3
110672 2 1 4
111057 2 1 5
111386 2 1 5
112575 2 1 5
113648 2 1 5
114644 2 1 6
114958 2 1 7
114984 2 1 10
115453 2 1 10
116140 2 1 6
116441 2 1 5
117327 2 1 10
117558 2 1 5
117565 3 1 7 10
119028 2 1 6
119257 2 1 7
119801 2 1 5
119819 2 1 2
120798 2 1 10
121256 2 1 7
121258 2 1 7
121540 2 1 2
121940 2 1 5
122503 2 1 10
123119 2 1 5
123559 3 1 4 5
124551 2 1 5
124674 2 1 7
125388 2 1 7
125887 2 1 7
127334 2 1 7
127516 2 1 10
128581 2 1 2
128730 2 1 5
128775 2 1 6
130443 2 1 5
131128 2 1 3
131545 3 1 6 5
132046 2 1 7
133144 2 1 10
133287 2 1 6
134777 2 1 3
135060 2 1 5
135143 2 1 5
135157 2 1 5
137380 2 1 7
138358 2 1 10
140569 2 1 5
140724 2 1 10
141282 2 1 4
142142 2 1 5
143994 2 1 6
144410 3 1 3 10
144469 2 1 5
144731 2 1 4
145085 2 1 5
145086 2 1 3
145138 2 1 5
146208 2 1 7
146253 2 1 4
3086 3 10 9 7
3802 2 10 3
4265 2 10 9
4700 2 10 9
5996 2 1 10
6122 2 10 9
7702 2 10 9
8480 2 12 10
8676 2 10 3
9680 2 10 3
9754 2 9 10
10047 2 10 3
10949 2 10 3
12344 2 12 10
12975 2 9 10
13728 2 10 9
13773 2 10 9
14466 2 1 10
14521 2 1 10
14608 2 11 10
16768 2 1 10
17064 2 1 10
17735 2 1 10
18068 2 1 10
18732 2 10 3
20123 2 9 10
21757 2 9 10
22313 2 1 10
22384 2 10 9
23549 2 3 10
24296 2 10 7
24757 2 10 3
25784 2 9 10
26023 2 1 10
27205 2 10 3
27647 2 10 3
29300 2 9 10
29336 2 10 7
30111 2 10 7
30124 2 10 9
30256 2 10 9
30327 2 10 3
30519 2 10 3
30615 2 9 10
30755 2 11 10
31112 2 10 9
32212 2 12 10
32387 2 10 9
32748 2 10 9
33304 3 12 10 9
34847 2 10 9
35006 2 1 10
35469 2 11 10
36201 2 1 10
38886 2 12 10
38964 2 10 9
40349 2 12 10
41722 2 10 9
41751 2 10 9
43849 2 12 10
44466 2 11 10
44734 2 10 9
45574 2 10 9
46020 2 11 10
46175 2 10 9
47228 2 9 10
47388 2 11 10
50478 2 9 10
51941 2 10 3
52982 2 1 10
54375 2 10 9
55926 2 12 10
56263 2 10 9
63498 3 11 10 12
64717 2 11 10
65602 2 12 10
65875 2 10 7
68564 2 10 9
72238 2 10 7
74095 2 10 3
74686 2 12 10
76642 2 12 10
76702 2 10 9
76708 2 11 10
79299 2 11 10
81526 2 1 10
81765 2 12 10
82127 2 1 10
82328 2 11 10
82391 2 10 9
83113 2 12 10
83436 2 12 10
83601 2 10 7
83808 2 10 3
83874 2 10 3
85075 2 1 10
86132 2 9 10
86511 2 10 9
87613 2 11 10
88572 2 10 9
88673 2 10 9
89547 2 10 7
91070 2 12 10
91167 2 10 7
91915 2 10 3
91985 2 12 10
94870 2 11 10
95456 2 10 7
96019 2 10 9
97470 2 12 10
99382 2 10 3
100798 2 12 10
101509 2 11 10
103012 2 11 10
104566 2 10 9
104840 2 10 9
105027 2 11 10
105030 2 10 9
105875 2 10 9
106330 2 10 3
109539 2 1 10
111318 2 10 3
111900 2 11 10
112151 2 9 10
112751 2 11 10
113028 2 10 7
113627 2 9 10
114866 2 10 9
114984 2 1 10
115453 2 1 10
115645 2 10 9
116322 2 11 10
117327 2 1 10
117565 3 1 7 10
120060 2 10 7
120798 2 1 10
121470 2 12 10
121556 2 10 3
122503 2 1 10
122848 2 10 7
124065 2 10 7
124181 3 11 10 3
126192 2 10 7
126995 2 12 10
127516 2 1 10
129588 2 10 3
132667 2 10 7
133144 2 1 10
133217 2 12 10
133766 2 12 10
134026 2 12 10
134477 2 10 9
135416 2 10 3
138358 2 1 10
138399 2 10 7
138572 2 10 7
138625 2 10 3
140719 2 10 7
140724 2 1 10
141743 2 9 10
141822 2 12 10
141961 2 12 10
143746 2 12 10
144410 3 1 3 10
145069 2 10 7
145494 2 10 9
2502 2 11 12
2888 2 11 15
3406 2 11 3
6344 2 11 12
6452 2 11 15
7368 2 11 15
7445 2 11 16
7565 2 11 16
8125 2 11 15
9339 2 11 12
10026 2 11 12
10487 2 11 15
10534 2 11 16
10783 2 11 15
14608 2 11 10
15059 2 11 12
15436 2 11 12
15988 2 11 12
18052 2 11 15
18993 2 11 12
19393 2 12 11
20477 2 11 16
21482 2 11 15
21738 2 11 16
22134 2 11 12
22895 2 11 12
23896 3 11 16 15
24968 2 11 15
26908 2 11 3
27038 2 11 16
27380 2 11 12
27493 2 11 16
29547 2 11 12
29956 2 11 16
30674 2 12 11
30755 2 11 10
31367 2 11 12
32318 2 11 16
33343 2 11 16
34675 2 11 16
34730 2 11 12
35469 2 11 10
36020 2 11 12
38568 2 11 12
38980 2 11 16
39753 2 11 12
40554 2 11 12
40878 2 11 15
41460 2 11 15
42995 2 11 15
43372 2 11 16
44126 2 11 15
44466 2 11 10
44633 2 11 16
44904 2 11 16
46020 2 11 10
46432 2 11 12
47388 2 11 10
52049 2 11 15
53976 2 11 12
54508 2 11 12
56745 2 11 12
56859 2 11 12
63498 3 11 10 12
63536 2 11 12
64717 2 11 10
65139 2 11 15
65697 2 12 11
66912 2 11 16
68277 2 11 16
68359 2 11 12
70028 2 11 16
71241 2 11 15
72436 2 11 15
74528 2 11 3
75205 2 11 3
75907 2 11 12
76345 2 12 11
76708 2 11 10
76955 2 11 15
77040 2 11 12