WPCD Summer Code Camp ("DIFFER")

8 Jun 2020, 09:00 → 19 Jun 2020, 18:00 Europe/Berlin

https://chalmers.zoom.us/j/65529293386 (Virtual)

David Coster (MPG), Michele Romanelli (European Task Force on Integrated Tokamak Modelling), Pär Strand, Rui Coelho (IPFN-IST)

WPCD aims during 2020 are the finalization of ETS-6 (IMAS version of ETS), the support to the exploitation of the released workflows in various tokamaks inside and outside EUROfusion and the development of new models/workflows. 

 

 During the two weeks code-camp, advancements in all the aspects of WPCD activities are expected to take place along with discussions and cross dissemination.

WPCD_Summer_General_Code_Camp_2020

Monday, 8 June

09:00 → 10:00 Plenary: Welcome

This is the main room to be used for the daily check-in and check-out events.

09:00 Welcome from TFL

09:10 Welcome from LOC

09:20 → 10:00 WIMAS: Integration of the Fokker-Planck solver FoPla in ETS workflow

WIMAS break out

Convener: Thomas Johnson (Fusion Plasma Physics (KTH/VR))

09:20 Discuss the integration of the Fokker-Planck solver FoPla in ETS workflow

FoPla is a Fokker-Planck code that can handle both ICRH and NBI, including the synergy between the two heating schemes.

10:00 → 11:00 Plenary: Check-In

This is the main room to be used for the daily check-in and check-out events.

11:00 → 12:00 EWE-Training: WPCD Tools

Training and demonstrations

11:00 WP CD Tools

This a refresher training on some of the key WP CD tools FC2K, Kepler, etc.

It will be slightly slanted towards the IMAS environment but differences to the CPO infrastructure whcih still is our production environemnt will be provided if different

Speaker: Dmitriy Yadikin (Chalmers, EUROfusion-VR)

refresh_fc2k_kepler_jun_20.pdf

16:00 → 17:00 Plenary: Check-out

This is the main room to be used for the daily check-in and check-out events.

Tuesday, 9 June

10:00 → 11:00 Plenary: Check-in

This is the main room to be used for the daily check-in and check-out events.

11:00 → 12:00 EWE-Training: ImasGo Training session

Training and demonstrations

11:00 IMASGO

Demonstration using the IMASGo! tools

Speakers: Michele Romanelli (European Task Force on Integrated Tokamak Modelling), Stuart Henderson (UKAEA)

14:00 → 16:00 WIMAS: NTM control with ECCD

WIMAS break out

Convener: Thomas Johnson (Fusion Plasma Physics (KTH/VR))

14:00 Planning NTM control with ECCD

Speakers: Lorenzo Figini, Olivier Sauter (EPFL), Silvana Nowak, Thomas Johnson (Fusion Plasma Physics (KTH/VR))

16:00 → 17:00 Plenary: Check-out

This is the main room to be used for the daily check-in and check-out events.

17:00 → 18:40 WIMAS: Discussion new actor to update core_profiles after equilibrium update

WIMAS break out

Convener: Thomas Johnson (Fusion Plasma Physics (KTH/VR))

17:00 Discuss new actor to update core_profiles after equilibrium is updated

Speaker: Thomas Johnson (Fusion Plasma Physics (KTH/VR))

Discuss new actor to update core_profiles after equilibrium is updated

Participants:
Thomas, Jorge, David, Rui, Nathan, Ernesto, Sterling, Dmitriy, Samuli

Background

The discussion was motivated by an old suggestion from Pär Strand and recent issues observed in simulations by Jorge, where B0 was not consistent in equilibrium and core_profiles. The reason for this inconsistency was that CHEASE was updating R0 inside a time-loop, but this change was lost in the output IDSs as R0 is time-independent in the IDS definitions.

The key problem that we wish to solve is that in the ETS-6 workflow the equilibrium and core_profiles are updated in different actors. In between these actors there are inconsistencies between fields that are common to both IDSs, e.g. volume, area, B0, etc. We already have an actor, eqinput, that updates equilibrium following updates to core_profiles. This discussion was to design an actor to update core_profiles after updating equilibrium.

Discussion on the design of the "new actor"

Action to solve issue with R0 and B0:

• Thomas should contact Olivier to add an option to keep R0 constant.

The main fields to be updated by the new actor; copy data from equilibrium to core_profiles:

• B0 to be update
• rho_tor[edge]
  • phi
• volume
  • rho_vol
• area
• surface area (to be added to profiles_1d/grid/)

In addition there are a number of quantities that depend on the depend on the main fields, that should be updated for consistency, or in order to conserve the correct quantities.

• Kinetic profiles: ne, ni, Te, Ti, ...
  • Pär has  been proposing to modify density and temperature to conserve particles and energy inside rho. 
  • Jorge: add this step and have a flag to activate/deactivate
• Quantities related to the plasma current, e.g. psi, J_parallel, J_tor, q, shear...
  • At present we use NCSCAL=4 in CHEASE, so poloidal flux (psi) is conserved (or differential flux).
  • If the plasma geometry is changing, while psi is conserved, does this mean that we introduce an artificial current source (as the area integrated current may not be conserved)? Is there a preferred quantity to conserve?
  • Questions:
    • Can we use CHEASE with prescribed Ip?​​​​​​
      • Answer: We can! In theory we should also be able to prescribed q.
    • ​​​​​​If we enforce boundary conditions to keep e.g. Ip constant, then what do we do with psi? Should we change
  • Conclusion:
    • It is not clear how to change any current-related quantities.
    • In initial versions of the code we should not modify any current-related quantities.

Discussion on field outside core_profiles that are affected by updates to the equilibrium

Core_transport IDS:

• D, v, fluxes - should these be invariant when the plasma geometry evolves?
  • Thomas: I don't think we can conserve fluxes, D and v for an arbitrary modification of the plasma geometry.
  • The quantity we'd like to conserve are the area integrated fluxes (at least there seemed to be consensus on this.)
    • However, core_transport include only the flux surface averaged flux density. To find the area integrated fluxes we need the surface area, which is only available in the equilibrium IDS.
    • The problem is that in some applications the equilibrium is updated more frequently than the transport coefficients. In this case we lose the value of the flux surfaces that was used when calculating the transport coefficients, thus we cannot reproduce the area integrated flux.
    • SOLUTION: Thomas should add a tracker to request that the flux surface area is added to the ids_core_radial_grid (i.e. in core_transport). 
  • If the plasma geometry is changing we can keep the total flux invariant. What about D and v?
    • We could apply D-v-splitting!
    • Do we need to know that values of <grad(rho)> and <grad(rho^2)> when the D and v were originally calculated?
      • No-one could give a strong argument for adding these in core_transport.

Action:

• Jorge: check time derivatives of V’, <grad(rho)>, <grad(rho^2)> and surfaces area

core_sources IDS:

• It is important to conserve cumulative sources when fetching sources from core_source
• Note that source may come from an old time step with a different cumulative volume, V(rho). 

Wednesday, 10 June

09:00 → 10:00 EWE

09:00 ETS at K-star

A discussion on some recent issues discvered throughthe K-star modelling activities

Speaker: Pär Strand

09:01 → 10:00 WIMAS: FoPla discussions

WIMAS break out

Convener: Thomas Johnson (Fusion Plasma Physics (KTH/VR))

09:01 FoPla discussions

Speaker: Thomas Johnson (Fusion Plasma Physics (KTH/VR))

10:00 → 11:00 Plenary: Check-in

This is the main room to be used for the daily check-in and check-out events.

11:00 → 12:00 EWE-Training: Trview Demonstration

Training and demonstrations

11:00 TRview

Demonstration on the TRview tool to set up AUG data for ETS (and others)

Speaker: Giovanni Tardini (IPP)

14:30 → 15:30 External: WPCD Technical Development Meeting

1. Status of GFORGE-NEXT
   1. machine has been available for a while
   2. problem being looked at by Daniel
2. Python IMAS interface [KvdP]
   1. postpone discussion?
3. IMASGO (EWE/CPT)
   1. EWE to look at extending data that is loaded
   2. Also ICRH data from TJ
4. UDA status
   1. MAST
   2. JET
   3. ITER
   4. others
      1. Need status and plans
5. C++ problems resolved?
   1. Problem mainly on the EFIT++ side?
   2. C++ templates at compile/link expensive
6. Memory consumption of AL with GGD
7. ETSVIZ -- extensions being discussed

16:00 → 17:00 Plenary: Check-out

This is the main room to be used for the daily check-in and check-out events.

16:45 → 17:30 NTM control with ECCD: Debrief

Convener: Thomas Johnson (Fusion Plasma Physics (KTH/VR))

16:45 Debrief

Speakers: Lorenzo Figini, Olivier Sauter, Silvana Nowak, Thomas Johnson (Fusion Plasma Physics (KTH/VR))

Thursday, 11 June

09:00 → 10:00 WIMAS: FoPla discussion

WIMAS break out

Convener: Thomas Johnson (Fusion Plasma Physics (KTH/VR))

09:00 FoPla discussions

Speaker: Thomas Johnson (Fusion Plasma Physics (KTH/VR))

10:00 → 11:00 Plenary: Check-in

This is the main room to be used for the daily check-in and check-out events.

11:00 → 12:00 EWE-Training: Basic introduction to GIT

Training and demonstrations

11:00 Basic Introducion to GIT

Just as the title says... as we have the means to move to GIT within the WP CD now seem to be a very opportune time to get some basics

Speaker: Karel van de Plassche (DIFFER)

20200611_wpcd_git_intro_plassche.pdf

11:01 → 13:01 WIMAS: Discussion on the PION code in IMAS

WIMAS break out

Convener: Thomas Johnson (Fusion Plasma Physics (KTH/VR))

11:01 Discussion on the PION code in the ETS-6

Speakers: Jorgi Manyer Fuertes (BSC), Mervi Mantsinen (BSC), Mireille Schneider (ITER Organisation), Thomas Johnson (Fusion Plasma Physics (KTH/VR)), Xavier Saez (BSC)

15:30 → 16:00 WIMAS: Run-in actor in IMAS

WIMAS break out

Convener: Thomas Johnson (Fusion Plasma Physics (KTH/VR))

15:30 Runin actor in IMAS

Speaker: Thomas Johnson (Fusion Plasma Physics (KTH/VR))

16:00 → 17:00 Plenary: Check-out

This is the main room to be used for the daily check-in and check-out events.

16:45 → 17:45 NTM control with ECCD: Debrief

Convener: Thomas Johnson (Fusion Plasma Physics (KTH/VR))

16:45 Debrief

Speaker: Thomas Johnson (Fusion Plasma Physics (KTH/VR))

Friday, 12 June

09:00 → 10:00 WIMAS: Discussion on FoPla

WIMAS break out

09:00 Discussion on FoPla

10:00 → 11:00 Plenary: Check-in

This is the main room to be used for the daily check-in and check-out events.

11:00 → 12:00 EWE-Training: CPO2IDS, IDS2CPO and friends

Training and demonstrations

11:00 CPO2IDS

Demonstration of CPO2IDS

Speaker: Nathan Cummings (CCFE)

ids2cpo-cpo2ids training.pptx

11:20 IDS2CPO

Demonstration of IDS2CPO

Speaker: Dr Jorge Ferreira (Instituto Superior Técnico / Instituto de Plasmas e Fusão Nuclear)

11:30 → 13:00 WIMAS: Discussion on sanitising input data in ETS-init

WIMAS break out

11:30 Discussion on sanitising input data in ETS-init

13:00 → 14:00 EWE: Discussion on internal boundary conditions and hyper diffusion

13:00 Discussion on internal boundary conditions and hyper diffusion

14:00 → 15:00 EWE: QLK/QLKNN workplan discussion

14:00 QLK/QLKNN workplan discussion

Agenda:

1. Status of Implementation Pär/Karel

2. Functional update requirements 
  a) Control of recalculation frequency of the integrals -  QLK updates needed? (Pär)
  b) Codeparams updates (Pär)

3. Benchmarking (Pärs Proposal) TBD 
  a) Benchmarking and validation EWE responsibility: ETS5 (a select few cases)
  b) ETSv6 testing:  On the same discharges but focus on ETSv5/ETS6 comparions (when ETS6 ready for exploitation switch over)

Need cases and procedures set up.

Quick summary:

1. We will run test cases form JINTRAC runs stored in IDS (equilibrium, profiles and core_sources is available BUT coresources might need updates (based on quick view on potential utodated mapped data)
2. The mkimas script is presented here https://users.euro-fusion.org/pages/data-cmg/wiki/JINTRAC_mkimas.html
3. Frida to start look at the script together with Jorge (?)
4. We will need to ask Nathan to provide  a coresource mapping for the WPCD side
5. Karel will provide references to proposed test cases (3-5 cases in priority order)
6. Karel and Pär to review which version of QLK/QLK-NN is needed and ensure that this is installed
7. Pär to update codeparams to  make QLK-NN available to end users
8.  Karel to provide run data (so that JAMS setting for the modules can be picked up)
9. Frida Pär and Dmitriy to validate setup and do runs

Reference:

• mkimas WIKI page: https://users.euro-fusion.org/pages/data-cmg/wiki/JINTRAC_mkimas.html
• QLKNN paper with base cases: https://aip.scitation.org/doi/full/10.1063/1.5134126, preprint: Fast modeling of turbulent transport in fusion plasmas using neural networks

 

16:00 → 17:00 Plenary: Check-out

This is the main room to be used for the daily check-in and check-out events.

16:45 → 17:30 NTM control with ECCD: Debrief

Convener: Thomas Johnson (Fusion Plasma Physics (KTH/VR))

16:45 Debrief

Speaker: Thomas Johnson (Fusion Plasma Physics (KTH/VR))

Monday, 15 June

10:00 → 11:00 Plenary: Check-in

This is the main room to be used for the daily check-in and check-out events.

11:00 → 12:00 EWE-Training: Training on autoGui

Training and demonstrations

Convener: Dmitriy Yadikin (Chalmers, EUROfusion-VR)

11:00 Training on autoGui

Speaker: Dmitriy Yadikin (Chalmers, EUROfusion-VR)

autogui_cc_jun20.pdf

13:00 → 13:30 Discuss ETS-6 Verification on ITER Scenarios

Speaker: Thomas Johnson (Fusion Plasma Physics (KTH/VR))

13:30 → 15:30 WIMAS: AMNS

WIMAS break out

Conveners: David Coster (MPG), Martin O'Mullane, David Tskhakaya

15:00 → 16:00 Discussion on using NEMO input in the PION code

Speaker: Thomas Johnson (Fusion Plasma Physics (KTH/VR))

16:00 → 17:00 Plenary: Check-out

This is the main room to be used for the daily check-in and check-out events.

16:45 → 17:30 NTM control with ECCD: Debrief

Convener: Thomas Johnson (Fusion Plasma Physics (KTH/VR))

16:45 Debrief

Speaker: Thomas Johnson (Fusion Plasma Physics (KTH/VR))

Tuesday, 16 June

10:00 → 11:00 Plenary: Check-in

This is the main room to be used for the daily check-in and check-out events.

11:00 → 12:00 EWE-Training: Training on the ETS-6 in IMAS

Training and demonstrations

Convener: Thomas Johnson (Fusion Plasma Physics (KTH/VR))

11:00 Training on the ETS-6 in IMAS

Speaker: Thomas Johnson (Fusion Plasma Physics (KTH/VR))

Training June2020 - How to run the ETS-6.pptx

15:30 → 16:00 WIMAS: Discussion on sawtooth actor in the ETS-6

WIMAS break out

Convener: Thomas Johnson (Fusion Plasma Physics (KTH/VR))

15:30 Discuss sawtooth actor in ETS-6

Speaker: Thomas Johnson (Fusion Plasma Physics (KTH/VR))

16:00 → 17:00 Plenary: Check-out

This is the main room to be used for the daily check-in and check-out events.

16:45 → 17:30 NTM control with ECCD: Debrief

Convener: Thomas Johnson (Fusion Plasma Physics (KTH/VR))

16:45 Debrief

Speaker: Thomas Johnson (Fusion Plasma Physics (KTH/VR))

Wednesday, 17 June

10:00 → 11:00 Plenary: Check-in

This is the main room to be used for the daily check-in and check-out events.

11:00 → 12:00 EWE-Training: Training on ReadAUG

Training and demonstrations

11:00 Training on ReeadAUG

Speaker: Mike Dunne (IPP)

13:00 → 15:00 Discussion on sawtooth modelling in ETS-6

Speaker: Thomas Johnson (Fusion Plasma Physics (KTH/VR))

Discussion on sawteeth in the ETS-6

Present: Antoine, Jorge, Olivier, Thomas

AGENDA:

- Status and plan for adapting sawtooth actors to IMAS

- The need to updates to the core_instant_changes IDS.

- Workflow implementation for instant-changes - indeed it is still non-existent.

 

Status of implementation:

• most of the work
• remaining is to fill core-instant-changes

The need to updates to the core_instant_changes IDS.

• Most changes are already implemented. The main issue now is that change in psi has to be put into the grid/psi.

How do we want to combine events?

1. priority and do one at a time
2. in serial and do equilibrium update between each event

 

Proposal: recompute the equilibrium after every event - before saving data to the UAL.

Olivier: if sawtooth model is included, then we need to recompute the equilibrium with high frequency after the sawtooth crash - possibly medium to high frequency in between.

Sawtooth-update would output core_profiles, in addition to core_instant_changes.

How to get consistency inside core_profiles?

Workflow:

  sawtooth-update  -> equilibrium -> update-core-profiles

Change the time-step, dt -> dt_min (1e-5 s), after every (sawtooth) event. Afterwards, the WF has to start increasing dt.

To do:

• Antoine, refresh actor
• Thomas, start building composite actor
• Let’s check the status in a month time…

 

14:30 → 15:30 External: WPCD Technical Development Meeting

1. Status of GFORGE-NEXT
   1. machine has been available for a while
   2. problem being looked at by Daniel
2. Python IMAS interface [KvdP]
   1. postpone discussion?
3. IMASGO (EWE/CPT)
   1. EWE to look at extending data that is loaded
   2. Also ICRH data from TJ
4. UDA status
   1. MAST
   2. JET
   3. ITER
   4. others
      1. Need status and plans
5. C++ problems resolved?
   1. Problem mainly on the EFIT++ side?
   2. C++ templates at compile/link expensive
6. Memory consumption of AL with GGD
7. ETSVIZ -- extensions being discussed

16:00 → 17:00 Plenary: Check-out

This is the main room to be used for the daily check-in and check-out events.

16:45 → 17:30 NTM control with ECCD: Debrief

Convener: Thomas Johnson (Fusion Plasma Physics (KTH/VR))

16:45 Debrief

Speaker: Thomas Johnson (Fusion Plasma Physics (KTH/VR))

Thursday, 18 June

10:00 → 11:00 Plenary: Check-in

This is the main room to be used for the daily check-in and check-out events.

15:00 → 16:00 WPCD Seminar

15:00 Data Analysis Workflows

Speaker: Rainer Fischer

fischer_ida.pdf

15:30 Data Analysis Workflows

Speaker: Marco Sertoli

Sertoli-2020_06_18-composition_code_workflow.pdf

16:00 → 17:00 Plenary: Check-out

This is the main room to be used for the daily check-in and check-out events.

16:45 → 17:30 NTM control with ECCD: Debrief

Convener: Thomas Johnson (Fusion Plasma Physics (KTH/VR))

16:45 Debrief

Speaker: Thomas Johnson (Fusion Plasma Physics (KTH/VR))

Friday, 19 June

10:00 → 12:00 Plenary: Final Reporting

This is the main room to be used for the daily check-in and check-out events.

14:00 → 17:00 Plenary: Final reporting, session 2

This is the main room to be used for the daily check-in and check-out events.