Should not set Scenario["plasma"]["2D_prof"] = []

[WolfgangSuttrop]

Hi Severin et al,

ECRadScenario.set_up_profiles_from_imas() is a bit too eager to reset everything to an (empty) list:

ECRad_PyLib/src/ecrad_pylib/ECRad_Scenario.py

Lines 232 to 233 in f6322e9

	for key in self["plasma"]:
	self["plasma"][key] = []

ECRad_PyLib/src/ecrad_pylib/ECRad_Scenario.py

Lines 240 to 241 in f6322e9

	for key in self["plasma"]:
	self["plasma"][key] = []

E.g. Scenario['plasma']['2D_prof'] should remain a scalar boolean in order to be able to write it to netcdf output later on:

ECRad_PyLib/src/ecrad_pylib/ECRad_Scenario.py

Lines 672 to 673 in f6322e9

	var = rootgrp["Scenario"].createVariable("plasma" + "_" + "2D_prof", "b")
	var[...] = int(self["plasma"]["2D_prof"])

I got it to work after rewriting the two loops above as:

for key in ["rhot_prof", "rhop_prof", "Te", "ne"]:
    self["plasma"][key] = []

Greetings,
Wolfgang

[AreWeDreaming]

Hi Wolfgang, sorry for ignoring this. I think I have created a quite general fix in 47e0910. Do you use the conda build? I can create a new release for you to test it.

[AreWeDreaming]

Also new mamba release for Pylib now,