|
115 | 115 | "metadata": {}, |
116 | 116 | "outputs": [], |
117 | 117 | "source": [ |
118 | | - "# initialise the profiles\n", |
119 | | - "from freegsnke.jtor_update import ConstrainBetapIp\n", |
120 | | - "profiles = ConstrainBetapIp(\n", |
121 | | - " eq=eq, # equilibrium object\n", |
122 | | - " betap=0.15, # poloidal beta\n", |
123 | | - " Ip=11e6, # plasma current\n", |
124 | | - " fvac=0.5, # fvac = rB_{tor}\n", |
125 | | - " alpha_m=2.0, # profile function parameter\n", |
126 | | - " alpha_n=1.0 # profile function parameter\n", |
| 118 | + "from freegsnke.jtor_update import Fiesta_Topeol\n", |
| 119 | + "profiles = Fiesta_Topeol(\n", |
| 120 | + " eq=eq,\n", |
| 121 | + " Beta0=0.5978,\n", |
| 122 | + " Ip=15e6,\n", |
| 123 | + " fvac=0.5,\n", |
| 124 | + " Raxis=6.2,\n", |
| 125 | + " alpha_m=2.0,\n", |
| 126 | + " alpha_n=1.395\n", |
127 | 127 | ")" |
128 | 128 | ] |
129 | 129 | }, |
|
257 | 257 | " # max_rel_psit=.05,\n", |
258 | 258 | " verbose=True, # print output\n", |
259 | 259 | " l2_reg=1e-14,\n", |
| 260 | + " Picard_handover=1e-4,\n", |
| 261 | + " full_jacobian_handover=[1e-3, 1e-2]\n", |
260 | 262 | " )\n" |
261 | 263 | ] |
262 | 264 | }, |
|
291 | 293 | "# with open('simple_diverted_currents_PaxisIp.pk', 'wb') as f:\n", |
292 | 294 | "# pickle.dump(obj=inverse_current_values, file=f)" |
293 | 295 | ] |
294 | | - }, |
295 | | - { |
296 | | - "cell_type": "code", |
297 | | - "execution_count": null, |
298 | | - "metadata": {}, |
299 | | - "outputs": [], |
300 | | - "source": [] |
301 | | - }, |
302 | | - { |
303 | | - "cell_type": "code", |
304 | | - "execution_count": null, |
305 | | - "metadata": {}, |
306 | | - "outputs": [], |
307 | | - "source": [] |
308 | 296 | } |
309 | 297 | ], |
310 | 298 | "metadata": { |
311 | 299 | "kernelspec": { |
312 | | - "display_name": "freegsnke4e", |
| 300 | + "display_name": "freegsnke-dev", |
313 | 301 | "language": "python", |
314 | 302 | "name": "python3" |
315 | 303 | }, |
|
323 | 311 | "name": "python", |
324 | 312 | "nbconvert_exporter": "python", |
325 | 313 | "pygments_lexer": "ipython3", |
326 | | - "version": "3.10.16" |
| 314 | + "version": "3.13.12" |
327 | 315 | } |
328 | 316 | }, |
329 | 317 | "nbformat": 4, |
|
0 commit comments