-
Notifications
You must be signed in to change notification settings - Fork 0
/
Copy pathtemplates.py
257 lines (251 loc) · 14 KB
/
templates.py
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
225
226
227
228
229
230
231
232
233
234
235
236
237
238
239
240
241
242
243
244
245
246
247
248
249
250
251
252
253
254
255
256
257
# '<SURFACE NUM> <MATERIAL NUM> -<DENSITY> -<SURFACE NUM> $ imp=1 <DESCRIPTION>'
surface_cell_planet = '%s 1 -%s -%s IMP:N=1 $ Planet\n'
# '<SURFACE NUM> <MATERIAL NUM> -<SURFACE NUM> <PREVIOUS SURFACE> $ imp=1 <DESCRIPTION>'
surface_cell_alt = '%s 0 -%s %s IMP:N=1 $ %s\n'
# '99 0 <FINAL SURFACE NUM> $ imp=0 outside'
surface_cell_outside = '99 0 %s IMP:N=0 $ outside\n'
# '<SURFACE NUM> so <SURFACE ALTITUDE IN CM>'
surface_card = '%s so %f\n'
# 'm1 $ Geochemistry with <WT % HYDROGEN> wt\% H\n'
material_header = 'm1 $ Geochemistry with %f wt. pct H\n'
# TODO: Ask Lucy: Should I use elemental (ZZZ000) unless specified atomic mass?
h = ' 1001 -%f $ H\n'
c = ' 6012 -%s $ C\n'
n = ' 7014 -%s $ N\n'
o16 = ' 8016 -%f $ O-16\n'
o17 = ' 8017 -%s $ O-17\n'
na = ' 11023 -%s $ Na\n'
mg = ' 12024 -%s $ Mg\n'
al = ' 13027 -%s $ Al\n'
si28 = ' 14028 -%s $ Si-28\n'
si29 = ' 14029 -%s $ Si-29\n'
si30 = ' 14030 -%s $ Si-30\n'
p = ' 15031 -%s $ P\n'
s = ' 16032 -%s $ S\n'
cl = ' 17000 -%s $ Cl\n' # Not sure what mass to use (35.453)
k = ' 19039 -%s $ K\n'
ca = ' 20040 -%s $ Ca\n'
sc = ' 21045 -%s $ Sc\n'
ti = ' 22048 -%s $ Ti\n'
cr50 = ' 24050 -%s $ Cr-50\n'
cr52 = ' 24052 -%s $ Cr-52\n'
cr53 = ' 24053 -%s $ Cr-53\n'
cr54 = ' 24054 -%s $ Cr-54\n'
mn55 = ' 25055 -%s $ Mn-55\n'
fe54 = ' 26054 -%s $ Fe-54\n'
fe56 = ' 26056 -%s $ Fe-56\n'
fe57 = ' 26057 -%s $ Fe-57\n'
fe58 = ' 26058 -%s $ Fe-58\n'
co = ' 27059 -%s $ Co\n'
ni = ' 28000 -%s $ Ni\n' # Not sure what mass to use (58.6934)
cu = ' 29000 -%s $ Cu\n' # Not sure what mass to use (63.546)
zn = ' 30000 -%s $ Zn\n' # Not sure what mass to use (65.38)
ga = ' 31000 -%s $ Ga\n' # Not sure what mass to use (69.723)
ge = ' 32073 -%s $ Ge\n'
ars = ' 33075 -%s $ As\n'
se = ' 34079 -%s $ Se\n'
w = ' 74184 -%s $ W\n'
re = ' 75187 -%s $ Re\n'
ir = ' 77193 -%s $ Ir\n'
pt = ' 78000 -%s $ Pt\n'
au = ' 79197 -%s $ Au\n'
sm = ' 62150 -%s $ Sm\n'
gd = ' 64157 -%s $ Gd\n'
th232 = ' 90232 -%s $ Th\n'
u234 = ' 92234 -%s $ U\n'
sdef = '''sdef sur=11 nrm=-1 erg=d1 par=9
si1 h 0.00000
10.2329 10.4713 10.7152 10.9648 11.2202 11.4815
11.749 12.0226 12.3027 12.5893 12.8825 13.1826
13.4896 13.8038 14.1254 14.4544 14.7911 15.1356
15.4882 15.8489 16.2181 16.5959 16.9824 17.378
17.7828 18.197 18.6209 19.0546 19.4984 19.9526
20.4174 20.893 21.3796 21.8776 22.3872 22.9087
23.4423 23.9883 24.5471 25.1189 25.704 26.3027
26.9153 27.5423 28.1838 28.8403 29.5121 30.1995
30.903 31.6228 32.3594 33.1131 33.8844 34.6737
35.4813 36.3078 37.1535 38.0189 38.9045 39.8107
40.738 41.6869 42.658 43.6516 44.6684 45.7088
46.7735 47.863 48.9779 50.1187 51.2861 52.4807
53.7032 54.9541 56.2341 57.544 58.8844 60.256
61.6595 63.0957 64.5654 66.0693 67.6083 69.1831
70.7946 72.4436 74.131 75.8578 77.6247 79.4328
81.2831 83.1764 85.1138 87.0964 89.1251 91.2011
93.3254 95.4993 97.7237 100.000 102.329 104.713
107.152 109.648 112.202 114.815 117.49 120.226
123.027 125.893 128.825 131.826 134.896 138.038
141.254 144.544 147.911 151.356 154.882 158.489
162.181 165.959 169.824 173.78 177.828 181.97
186.209 190.546 194.984 199.526 204.174 208.93
213.796 218.776 223.872 229.087 234.423 239.883
245.471 251.189 257.04 263.027 269.153 275.423
281.838 288.403 295.121 301.995 309.03 316.228
323.594 331.131 338.844 346.737 354.813 363.078
371.535 380.189 389.045 398.107 407.38 416.869
426.58 436.516 446.684 457.088 467.735 478.63
489.779 501.187 512.861 524.807 537.032 549.541
562.341 575.44 588.844 602.56 616.595 630.957
645.654 660.693 676.083 691.831 707.946 724.436
741.31 758.578 776.247 794.328 812.831 831.764
851.138 870.964 891.251 912.011 933.254 954.993
977.237 1000.00 1023.29 1047.13 1071.52 1096.48
1122.02 1148.15 1174.9 1202.26 1230.27 1258.93
1288.25 1318.26 1348.96 1380.38 1412.54 1445.44
1479.11 1513.56 1548.82 1584.89 1621.81 1659.59
1698.24 1737.8 1778.28 1819.7 1862.09 1905.46
1949.84 1995.26 2041.74 2089.3 2137.96 2187.76
2238.72 2290.87 2344.23 2398.83 2454.71 2511.89
2570.4 2630.27 2691.53 2754.23 2818.38 2884.03
2951.21 3019.95 3090.3 3162.28 3235.94 3311.31
3388.44 3467.37 3548.13 3630.78 3715.35 3801.89
3890.45 3981.07 4073.8 4168.69 4265.8 4365.16
4466.84 4570.88 4677.35 4786.3 4897.79 5011.87
5128.61 5248.07 5370.32 5495.41 5623.41 5754.4
5888.44 6025.6 6165.95 6309.57 6456.54 6606.93
6760.83 6918.31 7079.46 7244.36 7413.1 7585.78
7762.47 7943.28 8128.31 8317.64 8511.38 8709.64
8912.51 9120.11 9332.54 9549.93 9772.37 10000.0
10232.9 10471.3 10715.2 10964.8 11220.2 11481.5
11749.0 12022.6 12302.7 12589.3 12882.5 13182.6
13489.6 13803.8 14125.4 14454.4 14791.1 15135.6
15488.2 15848.9 16218.1 16595.9 16982.4 17378.0
17782.8 18197.0 18620.9 19054.6 19498.4 19952.6
20417.4 20893.0 21379.6 21877.6 22387.2 22908.7
23442.3 23988.3 24547.1 25118.9 25704.0 26302.7
26915.3 27542.3 28183.8 28840.3 29512.1 30199.5
30903.0 31622.8 32359.4 33113.1 33884.4 34673.7
35481.3 36307.8 37153.5 38018.9 38904.5 39810.7
40738.0 41686.9 42658.0 43651.6 44668.4 45708.8
46773.5 47863.0 48977.9 50118.7 51286.1 52480.7
53703.2 54954.1 56234.1 57544.0 58884.4 60256.0
61659.5 63095.7 64565.4 66069.3 67608.3 69183.1
70794.6 72443.6 74131.0 75857.8 77624.7 79432.8
81283.1 83176.4 85113.8 87096.4 89125.1 91201.1
93325.4 95499.3 97723.7 100000.0
sp1 0.0 7.29193E-05 7.6264E-05 7.97598E-05
8.34136E-05 8.72323E-05 9.12232E-05 9.53937E-05
9.97519E-05 1.04306E-04 1.09064E-04 1.14036E-04
1.19231E-04 1.24658E-04 1.30327E-04 1.36249E-04
1.42435E-04 1.48897E-04 1.55646E-04 1.62694E-04
1.70055E-04 1.77742E-04 1.85768E-04 1.94149E-04
2.02899E-04 2.12034E-04 2.21571E-04 2.31525E-04
2.41916E-04 2.5276E-04 2.64078E-04 2.75889E-04
2.88213E-04 3.01072E-04 3.14489E-04 3.28485E-04
3.43084E-04 3.58313E-04 3.74196E-04 3.90759E-04
4.08031E-04 4.26041E-04 4.44817E-04 4.64391E-04
4.84794E-04 5.0606E-04 5.28223E-04 5.51318E-04
5.75381E-04 6.00452E-04 6.26568E-04 6.53771E-04
6.82103E-04 7.11606E-04 7.42326E-04 7.74309E-04
8.07602E-04 8.42255E-04 8.78318E-04 9.15845E-04
9.54889E-04 9.95506E-04 1.03775E-03 1.08169E-03
1.12738E-03 1.17488E-03 1.22425E-03 1.27557E-03
1.32891E-03 1.38432E-03 1.44189E-03 1.50168E-03
1.56377E-03 1.62825E-03 1.69518E-03 1.76466E-03
1.83675E-03 1.91156E-03 1.98916E-03 2.06964E-03
2.1531E-03 2.23962E-03 2.32931E-03 2.42225E-03
2.51855E-03 2.61829E-03 2.7216E-03 2.82855E-03
2.93927E-03 3.05386E-03 3.17241E-03 3.29505E-03
3.42187E-03 3.553E-03 3.68853E-03 3.82858E-03
3.97327E-03 4.12271E-03 4.277E-03 4.43626E-03
4.60062E-03 4.77017E-03 4.94503E-03 5.12531E-03
5.31114E-03 5.5026E-03 5.69982E-03 5.90291E-03
6.11196E-03 6.32708E-03 6.54837E-03 6.77593E-03
7.00986E-03 7.25025E-03 7.49719E-03 7.75076E-03
8.01103E-03 8.2781E-03 8.55203E-03 8.83287E-03
9.1207E-03 9.41557E-03 9.71751E-03 1.00266E-02
1.03428E-02 1.06662E-02 1.09968E-02 1.13346E-02
1.16796E-02 1.20317E-02 1.23911E-02 1.27575E-02
1.31311E-02 1.35117E-02 1.38992E-02 1.42937E-02
1.46949E-02 1.51029E-02 1.55174E-02 1.59384E-02
1.63657E-02 1.67991E-02 1.72385E-02 1.76837E-02
1.81345E-02 1.85907E-02 1.9052E-02 1.95183E-02
1.99891E-02 2.04643E-02 2.09436E-02 2.14267E-02
2.19131E-02 2.24027E-02 2.2895E-02 2.33897E-02
2.38863E-02 2.43845E-02 2.48839E-02 2.53841E-02
2.58845E-02 2.63848E-02 2.68845E-02 2.73831E-02
2.78801E-02 2.8375E-02 2.88673E-02 2.93565E-02
2.98421E-02 3.03235E-02 3.08001E-02 3.12715E-02
3.17371E-02 3.21963E-02 3.26486E-02 3.30933E-02
3.35301E-02 3.39582E-02 3.43772E-02 3.47864E-02
3.51854E-02 3.55736E-02 3.59505E-02 3.63155E-02
3.66682E-02 3.70079E-02 3.73343E-02 3.76469E-02
3.79451E-02 3.82286E-02 3.84969E-02 3.87496E-02
3.89864E-02 3.92068E-02 3.94104E-02 3.95971E-02
3.97664E-02 3.99182E-02 4.00521E-02 4.01679E-02
4.02655E-02 4.03446E-02 4.04051E-02 4.0447E-02
4.047E-02 4.04743E-02 4.04597E-02 4.04263E-02
4.0374E-02 4.0303E-02 4.02134E-02 4.01053E-02
3.99788E-02 3.98341E-02 3.96715E-02 3.94911E-02
3.92933E-02 3.90783E-02 3.88465E-02 3.85982E-02
3.83337E-02 3.80535E-02 3.7758E-02 3.74476E-02
3.71227E-02 3.67838E-02 3.64315E-02 3.60661E-02
3.56883E-02 3.52986E-02 3.48975E-02 3.44855E-02
3.40632E-02 3.36313E-02 3.31902E-02 3.27405E-02
3.22829E-02 3.18179E-02 3.13461E-02 3.08681E-02
3.03845E-02 2.98958E-02 2.94027E-02 2.89056E-02
2.84052E-02 2.79021E-02 2.73967E-02 2.68896E-02
2.63814E-02 2.58725E-02 2.53634E-02 2.48547E-02
2.43469E-02 2.38403E-02 2.33354E-02 2.28328E-02
2.23326E-02 2.18355E-02 2.13417E-02 2.08517E-02
2.03657E-02 1.98842E-02 1.94073E-02 1.89355E-02
1.8469E-02 1.8008E-02 1.75529E-02 1.71037E-02
1.66608E-02 1.62244E-02 1.57945E-02 1.53714E-02
1.49553E-02 1.45461E-02 1.41442E-02 1.37495E-02
1.33621E-02 1.29821E-02 1.26096E-02 1.22446E-02
1.18872E-02 1.15373E-02 1.1195E-02 1.08603E-02
1.05331E-02 1.02134E-02 9.90124E-03 9.59653E-03
9.29923E-03 9.00928E-03 8.72661E-03 8.45115E-03
8.18282E-03 7.92153E-03 7.66721E-03 7.41974E-03
7.17903E-03 6.94499E-03 6.71749E-03 6.49643E-03
6.2817E-03 6.07318E-03 5.87074E-03 5.67428E-03
5.48366E-03 5.29877E-03 5.11948E-03 4.94567E-03
4.77721E-03 4.61397E-03 4.45584E-03 4.30268E-03
4.15438E-03 4.0108E-03 3.87183E-03 3.73735E-03
3.60723E-03 3.48136E-03 3.35962E-03 3.24189E-03
3.12806E-03 3.01802E-03 2.91166E-03 2.80887E-03
2.70954E-03 2.61357E-03 2.52086E-03 2.4313E-03
2.34481E-03 2.26128E-03 2.18062E-03 2.10274E-03
2.02756E-03 1.95498E-03 1.88492E-03 1.8173E-03
1.75204E-03 1.68907E-03 1.62829E-03 1.56966E-03
1.51308E-03 1.4585E-03 1.40585E-03 1.35505E-03
1.30606E-03 1.2588E-03 1.21321E-03 1.16925E-03
1.12686E-03 1.08597E-03 1.04654E-03 1.00852E-03
9.71862E-04 9.36514E-04 9.02432E-04 8.69572E-04
8.37892E-04 8.0735E-04 7.77907E-04 7.49523E-04
7.22161E-04 6.95786E-04 6.70362E-04 6.45856E-04
6.22236E-04 5.99469E-04 5.77527E-04 5.56378E-04
5.35996E-04 5.16353E-04 4.97423E-04 4.79179E-04
4.61599E-04 4.44657E-04 4.28331E-04 4.12599E-04
3.9744E-04 3.82833E-04 3.68759E-04 3.55197E-04
3.42131E-04 3.29541E-04 3.17411E-04 3.05724E-04
2.94464E-04 2.83616E-04 2.73165E-04 2.63097E-04
2.53396E-04 2.44052E-04 2.35049E-04 2.26377E-04
2.18023E-04 2.09975E-04 2.02222E-04 1.94754E-04
1.87561E-04 1.80631E-04 1.73957E-04 1.67527E-04
1.61334E-04
'''
physics = '''c =========Physics============================================
c
mode n p h
phys:p 1e5 2j -1 $ j -102
phys:n 1e5 5j 2
phys:h 100000.5
c mphys
cut:n 2j 0 0
lca 8j 1 1
lcb 1000 1000 1000 1000 1000
totnu
c'''
tally_header = '''
c =========Tally============================================
'''
# <8.617E-09 for # of surfaces, last (99) is 0>
tally_tmp = '''tmp1 %s\n'''
tally_surface = '''f%d2:n %s
c fm%d2 3.14159E+14
fc%d2 surface
e%d2 0.0 1e-10 74log 1e5 $ energy bins
c *c%d2 180 90 0 $ angle bins\n'''
nps = '''print
NPS %s
c ctme 1\n'''