aboutsummaryrefslogtreecommitdiff
path: root/Source/Utils/NCIGodfreyTables.H
blob: 8cb105aa001381e2d6c845aa47131cdd2fbdb627 (plain) (blame)
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
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
169
170
171
172
173
174
175
176
177
178
179
180
181
182
183
184
185
186
187
188
189
190
191
192
193
194
195
196
197
198
199
200
201
202
203
204
205
206
207
208
209
210
211
212
213
214
215
216
217
218
219
220
221
222
223
224
#include <AMReX_AmrCore.H>

#ifndef WARPX_GODFREY_COEFF_TABLE_H_
#define WARPX_GODFREY_COEFF_TABLE_H_

// Table width. This is related to the stencil length
const int tab_width = 4;
// table length. Each line correspond to 1 value of cdtodz
// (here 101 values).
const int tab_length = 101;

// Table of coefficient for Ex, Ey abd Bz
// We typically interpolate between two lines
const amrex::Real table_nci_godfrey_Ex_Ey_Bz[tab_length][tab_width]{
        -2.47536,2.04288,-0.598163,0.0314711,  
        -2.47536,2.04288,-0.598163,0.0314711,  
        -2.47545,2.04309,-0.598307,0.0315029,  
        -2.4756,2.04342,-0.598549,0.0315558,  
        -2.47581,2.0439,-0.598886,0.0316298,  
        -2.47608,2.0445,-0.59932,0.031725,  
        -2.47641,2.04525,-0.59985,0.0318412,  
        -2.4768,2.04612,-0.600477,0.0319785,  
        -2.47725,2.04714,-0.6012,0.0321367,  
        -2.47776,2.04829,-0.602019,0.0323158,  
        -2.47833,2.04957,-0.602934,0.0325158,  
        -2.47896,2.05099,-0.603944,0.0327364,  
        -2.47965,2.05254,-0.605051,0.0329777,  
        -2.4804,2.05423,-0.606253,0.0332396,  
        -2.48121,2.05606,-0.60755,0.0335218,  
        -2.48208,2.05802,-0.608942,0.0338243,  
        -2.48301,2.06012,-0.610429,0.0341469,  
        -2.48401,2.06235,-0.61201,0.0344895,  
        -2.48506,2.06471,-0.613685,0.0348519,  
        -2.48618,2.06721,-0.615453,0.0352339,  
        -2.48735,2.06984,-0.617314,0.0356353,  
        -2.48859,2.07261,-0.619268,0.0360559,  
        -2.48988,2.0755,-0.621312,0.0364954,  
        -2.49123,2.07853,-0.623447,0.0369536,  
        -2.49265,2.08169,-0.625672,0.0374302,  
        -2.49412,2.08498,-0.627986,0.0379248,  
        -2.49565,2.0884,-0.630386,0.0384372,  
        -2.49724,2.09194,-0.632873,0.0389669,  
        -2.49888,2.09561,-0.635443,0.0395135,  
        -2.50058,2.09939,-0.638096,0.0400766,  
        -2.50234,2.1033,-0.640829,0.0406557,  
        -2.50415,2.10732,-0.64364,0.0412502,  
        -2.50601,2.11145,-0.646526,0.0418594,  
        -2.50791,2.1157,-0.649485,0.0424828,  
        -2.50987,2.12004,-0.652512,0.0431196,  
        -2.51187,2.12448,-0.655604,0.0437688,  
        -2.51392,2.12901,-0.658756,0.0444297,  
        -2.516,2.13363,-0.661964,0.0451011,  
        -2.51812,2.13832,-0.665221,0.0457818,  
        -2.52027,2.14308,-0.668521,0.0464705,  
        -2.52244,2.14789,-0.671856,0.0471658,  
        -2.52464,2.15274,-0.675218,0.0478658,  
        -2.52684,2.15762,-0.678596,0.0485687,  
        -2.52906,2.16251,-0.68198,0.0492723,  
        -2.53126,2.16738,-0.685355,0.049974,  
        -2.53345,2.17222,-0.688706,0.0506708,  
        -2.53561,2.177,-0.692015,0.0513594,  
        -2.53773,2.18168,-0.69526,0.0520359,  
        -2.53978,2.18623,-0.698416,0.0526955,  
        -2.54175,2.19059,-0.701452,0.053333,  
        -2.5436,2.19471,-0.704331,0.0539417,  
        -2.54531,2.19852,-0.70701,0.0545141,  
        -2.54683,2.20193,-0.709433,0.0550409,  
        -2.5481,2.20483,-0.711533,0.0555106,  
        -2.54906,2.20709,-0.713224,0.0559094,  
        -2.54963,2.20852,-0.714397,0.0562198,  
        -2.54968,2.20888,-0.714907,0.0564196,  
        -2.54905,2.20785,-0.714562,0.0564797,  
        -2.54751,2.20496,-0.713094,0.0563618,  
        -2.54472,2.19955,-0.710118,0.0560124,  
        -2.54014,2.19058,-0.705048,0.0553544,  
        -2.53286,2.1763,-0.69693,0.0542684,  
        -2.52115,2.15344,-0.684027,0.05255,  
        -2.50098,2.11466,-0.66255,0.0497817,  
        -2.45797,2.03459,-0.620099,0.0446889,  
        -2.28371,1.72254,-0.465905,0.0283268,  
        -2.4885,2.04899,-0.599292,0.0390466,  
        -2.1433,1.36735,-0.220924,-0.00215633,  
        -2.4943,2.07019,-0.610552,0.035166,  
        -2.84529,2.77303,-1.00018,0.0724884,  
        -2.72242,2.51888,-0.847226,0.0509964,  
        -2.65633,2.3744,-0.750392,0.0326366,  
        -2.59601,2.23412,-0.646421,0.00868027,  
        -2.51477,2.0369,-0.491066,-0.0306397,  
        -2.35935,1.65155,-0.178971,-0.112713,  
        -1.84315,0.361693,0.876104,-0.393844,  
        -2.65422,2.39262,-0.789663,0.0516265,  
        -3.46529,4.42354,-2.45543,0.497097,  
        -3.15747,3.65311,-1.824,0.328432,  
        -3.04694,3.37613,-1.59668,0.267631,  
        -2.99205,3.23814,-1.48302,0.237103,  
        -2.96075,3.15894,-1.41733,0.219317,  
        -2.94172,3.11028,-1.37649,0.20811,  
        -2.92994,3.07962,-1.35025,0.200755,  
        -2.92283,3.06054,-1.33338,0.195859,  
        -2.91894,3.04938,-1.3229,0.192637,  
        -2.91736,3.04394,-1.31702,0.190612,  
        -2.91753,3.04278,-1.31456,0.189477,  
        -2.91905,3.04494,-1.31475,0.189026,  
        -2.92165,3.04973,-1.31705,0.189117,  
        -2.92512,3.05667,-1.32105,0.189646,  
        -2.92933,3.06539,-1.32646,0.190538,  
        -2.93416,3.07562,-1.33308,0.191735,  
        -2.93952,3.08715,-1.34072,0.193194,  
        -2.94535,3.09982,-1.34925,0.194881,  
        -2.95159,3.11349,-1.35858,0.196769,  
        -2.9582,3.12805,-1.36861,0.198838,  
        -2.96514,3.14342,-1.37929,0.201068,  
        -2.97239,3.15953,-1.39055,0.203448,  
        -2.97991,3.17632,-1.40234,0.205964,  
        -2.98769,3.19374,-1.41463,0.208607
            };

// Table of coefficient for Bx, By and Ez
// We typically interpolate between two lines
const amrex::Real table_nci_godfrey_Bx_By_Ez[tab_length][tab_width]{
        -2.80862,2.80104,-1.14615,0.154077,  
        -2.80862,2.80104,-1.14615,0.154077,  
        -2.80851,2.80078,-1.14595,0.154027,  
        -2.80832,2.80034,-1.14561,0.153945,  
        -2.80807,2.79973,-1.14514,0.153829,  
        -2.80774,2.79894,-1.14454,0.15368,  
        -2.80733,2.79798,-1.1438,0.153498,  
        -2.80685,2.79685,-1.14292,0.153284,  
        -2.8063,2.79554,-1.14192,0.153036,  
        -2.80568,2.79405,-1.14077,0.152756,  
        -2.80498,2.79239,-1.1395,0.152443,  
        -2.80421,2.79056,-1.13809,0.152098,  
        -2.80337,2.78856,-1.13656,0.151721,  
        -2.80246,2.78638,-1.13488,0.151312,  
        -2.80147,2.78404,-1.13308,0.150871,  
        -2.80041,2.78152,-1.13115,0.150397,  
        -2.79927,2.77882,-1.12908,0.149893,  
        -2.79807,2.77596,-1.12689,0.149356,  
        -2.79679,2.77292,-1.12456,0.148789,  
        -2.79543,2.76972,-1.12211,0.14819,  
        -2.79401,2.76634,-1.11953,0.14756,  
        -2.79251,2.76279,-1.11681,0.1469,  
        -2.79094,2.75907,-1.11397,0.146208,  
        -2.78929,2.75517,-1.111,0.145486,  
        -2.78757,2.7511,-1.10789,0.144733,  
        -2.78578,2.74686,-1.10466,0.14395,  
        -2.78391,2.74245,-1.1013,0.143137,  
        -2.78196,2.73786,-1.09781,0.142293,  
        -2.77994,2.73309,-1.09419,0.141419,  
        -2.77784,2.72814,-1.09043,0.140514,  
        -2.77566,2.72301,-1.08654,0.139578,  
        -2.7734,2.7177,-1.08252,0.138612,  
        -2.77106,2.7122,-1.07836,0.137614,  
        -2.76864,2.70651,-1.07406,0.136586,  
        -2.76613,2.70062,-1.06962,0.135525,  
        -2.76353,2.69453,-1.06503,0.134432,  
        -2.76084,2.68824,-1.0603,0.133307,  
        -2.75806,2.68173,-1.05541,0.132148,  
        -2.75518,2.675,-1.05037,0.130954,  
        -2.75219,2.66804,-1.04516,0.129725,  
        -2.7491,2.66084,-1.03978,0.12846,  
        -2.7459,2.65339,-1.03423,0.127156,  
        -2.74257,2.64566,-1.02848,0.125813,  
        -2.73912,2.63765,-1.02254,0.124428,  
        -2.73552,2.62934,-1.01638,0.122999,  
        -2.73178,2.62069,-1.01,0.121523,  
        -2.72787,2.61169,-1.00337,0.119996,  
        -2.72379,2.6023,-0.996479,0.118417,  
        -2.71951,2.59248,-0.989294,0.116778,  
        -2.71501,2.58218,-0.981786,0.115076,  
        -2.71026,2.57135,-0.97392,0.113303,  
        -2.70524,2.55991,-0.965651,0.111453,  
        -2.69989,2.54778,-0.956922,0.109514,  
        -2.69416,2.53484,-0.947666,0.107476,  
        -2.68799,2.52096,-0.937795,0.105324,  
        -2.68129,2.50596,-0.927197,0.103039,  
        -2.67394,2.48959,-0.915724,0.100597,  
        -2.66578,2.47153,-0.903179,0.097968,  
        -2.65657,2.4513,-0.889283,0.0951084,  
        -2.64598,2.42824,-0.873638,0.0919592,  
        -2.63347,2.40127,-0.855632,0.0884325,  
        -2.61813,2.36864,-0.834261,0.0843898,  
        -2.59821,2.32701,-0.807691,0.0795876,  
        -2.56971,2.26887,-0.77188,0.0735132,  
        -2.51823,2.16823,-0.713448,0.0645399,  
        -2.33537,1.8294,-0.533852,0.0409941,  
        -2.53143,2.14818,-0.670502,0.053982,  
        -2.17737,1.43641,-0.259095,0.00101255,  
        -2.51929,2.12931,-0.654743,0.0452381,  
        -2.86122,2.82221,-1.05039,0.0894636,  
        -2.72908,2.54506,-0.87834,0.0626188,  
        -2.6536,2.37954,-0.7665,0.0409117,  
        -2.58374,2.21923,-0.649738,0.0146791,  
        -2.49284,2.00346,-0.48457,-0.0255348,  
        -2.32762,1.60337,-0.1698,-0.105287,  
        -1.80149,0.316787,0.855414,-0.369652,  
        -2.60242,2.28418,-0.721378,0.040091,  
        -3.40335,4.25157,-2.29817,0.449834,  
        -3.0852,3.47341,-1.67791,0.28982,  
        -2.9642,3.17856,-1.44399,0.229852,  
        -2.89872,3.01966,-1.31861,0.197945,  
        -2.85668,2.91811,-1.23894,0.17783,  
        -2.82679,2.84621,-1.18287,0.163785,  
        -2.80401,2.79167,-1.14058,0.153278,  
        -2.78577,2.74819,-1.10706,0.145015,  
        -2.77061,2.7122,-1.07946,0.138267,  
        -2.75764,2.68152,-1.05606,0.132589,  
        -2.74627,2.65475,-1.03575,0.127695,  
        -2.73612,2.63093,-1.01777,0.123395,  
        -2.72692,2.6094,-1.00159,0.119553,  
        -2.71846,2.58968,-0.986841,0.116074,  
        -2.71061,2.57142,-0.973239,0.112887,  
        -2.70323,2.55434,-0.960573,0.109937,  
        -2.69626,2.53824,-0.948678,0.107185,  
        -2.68962,2.52294,-0.937429,0.104598,  
        -2.68327,2.50833,-0.926722,0.102151,  
        -2.67714,2.4943,-0.916477,0.0998223,  
        -2.67122,2.48076,-0.906627,0.0975966,  
        -2.66546,2.46764,-0.897118,0.0954599,  
        -2.65985,2.45489,-0.887903,0.0934011,  
        -2.65437,2.44244,-0.878945,0.0914107
            };

#endif // #ifndef WARPX_GODFREY_COEFF_TABLE_H_