Update example use of IMAS heat flux model

ProjectTorreyPines · Aug 1, 2024 · 17d4598 · 17d4598
1 parent 2512248
commit 17d4598
Show file tree

Hide file tree

Showing 2 changed files with 7 additions and 3 deletions.
diff --git a/Project.toml b/Project.toml
@@ -7,9 +7,11 @@ version = "1.0.1"
 ArgParse = "c7e460c6-2fb9-53a9-8c5b-16f535851c63"
 DSP = "717857b8-e6f2-59f4-9121-6e50c889abd2"
 DelimitedFiles = "8bb1440f-4735-579b-a4ab-409b98df4dab"
+EFIT = "cda752c5-6b03-55a3-9e33-132a441b0c17"
 GGDUtils = "b7b5e640-9b39-4803-84eb-376048795def"
-IMASDD = "06b86afa-9f21-11ec-2ef8-e51b8960cfc5"
 IMAS = "13ead8c1-b7d1-41bb-a6d0-5b8b65ed587a"
+IMASDD = "06b86afa-9f21-11ec-2ef8-e51b8960cfc5"
+IMASdd = "c5a45a97-b3f9-491c-b9a7-aa88c3bc0067"
 Interpolations = "a98d9a8b-a2ab-59e6-89dd-64a1c18fca59"
 LsqFit = "2fda8390-95c7-5789-9bda-21331edee243"
 PhysicalConstants = "5ad8b20f-a522-5ce9-bfc9-ddf1d5bda6ab"

diff --git a/examples/wall_heat_flux_demo.jl b/examples/wall_heat_flux_demo.jl
@@ -18,7 +18,9 @@ dd = SD4SOLPS.preparation(
     output_format=output_format,
     allow_boundary_flux_correction=true,
 )
-
+for eqt in dd.equilibrium.time_slice
+    IMAS.flux_surfaces(eqt)
+end
 Rwall, Zwall, rwall, zwall, s, SOL, r, q = IMAS.mesher_HF(dd)
-IMAS.particle_HF(dd.equilibrium.time_slice[1], SOL, rwall, zwall, r, q)
+Qwall, Qpara = IMAS.particle_HF(dd.equilibrium.time_slice[1], SOL, rwall, zwall, r, q)
 Plots.plot(s, Qwall)