Skip to content

Conversation

@imbeauf
Copy link
Collaborator

@imbeauf imbeauf commented Sep 23, 2025

</xs:element>
<xs:element name="motion_invariants_volume" type="generic_grid_scalar">
<xs:annotation>
<xs:documentation>Volume of the invariants of motion space volume elements</xs:documentation>
Copy link
Collaborator

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

Is there a simpler way to say this?

Copy link

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

I would use "Normalised volume of the CoM (constants of motion) cell elements". Constant is a more widely used terminology in the EP community than invariant.

Copy link
Collaborator Author

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

If it's normalized, to what is it normalized ? And what should be the units of this ? For the moment I indicated T.m^2.J

Copy link

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

For the energy I use eV instead of J, so that the CoM volume unit is T.m^2.eV.

The normalised CoM volume is normalised by the cartesian CoM volume element defined by the CoM grid steps ,stored in the GGD, \Delta V_{CoM} = \Delta Pphi \Delta \lambda \Delta E. In some parts of CoM space, this volume element is reduced when approaching the topological boundaries of the domain, which is why a 3D array is necessary for storing the volume element around each CoM cell.

Copy link
Collaborator Author

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

Definition clarified after talking to Guillaume, see last commit.
The CoM volume is normalised, so its units is "1", but we have restored "eV" for the energy since it makes sense for a single particule or particules distributions. This is allowed by the DD rules"SI+eV" and is used in several places of the DD for particle energies.

@DavidPCoster
Copy link
Contributor

I asked some of the fast particle people about this proposed change, and received the following:

Dave,

is seems the removing E is not what Guillaume wants: E should be the particle energy and not the electric field! Phase is probably the gyrophase in case we have a fully kinetic particle orbit. (add to description for clarity).

Maybe E has already been fixed meanwhile?

Did Guillaume comment on this?

Best,

Philipp

Copy link
Contributor

@DavidPCoster DavidPCoster left a comment

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

Has the comment from Philipp Lauber been addressed?

@SimonPinches
Copy link
Collaborator

@brochag, do you have any comments on this PR?

imbeauf and others added 3 commits October 2, 2025 13:53
Co-authored-by: Simon Pinches <SimonPinches@users.noreply.github.com>
Co-authored-by: Simon Pinches <SimonPinches@users.noreply.github.com>
@github-actions
Copy link

github-actions bot commented Oct 2, 2025

@github-actions
Copy link

github-actions bot commented Oct 2, 2025

@imbeauf
Copy link
Collaborator Author

imbeauf commented Oct 2, 2025

Has the comment from Philipp Lauber been addressed?

I have added "gyro"phase to the description.
The removal of the energy node was requested by @brochag, he told me that the Energy is redundant with the midplane invariants structure ... it's better if he comments directly on this though

@brochag
Copy link

brochag commented Oct 13, 2025

Hi all,

Sorry for my late response, I'm in between institutions and I haven't been given a work email yet ... Hopefully soon.

Indeed when storing the midplane invariants on the CoM grid, the energy leaf is redundant when using a constants of motion (CoM) grid that has energy as one of the dimension. As of now the only CoM grid option is (Pphi : toroidal canonical momentum, \lambda = mu B_0/E : pitch angle, E : energy).

However I will add later different options for the CoM grid dimensions (magnetic moment mu instead of energy E for example), so actually in that case the energy leaf will not be redundant anymore. I didn't think of that.

So in conclusion, we should leave the energy leaf structure. Sorry for the inaccurate request Frédéric.

Indeed the phase is the gyro phase when considering full-orbit particles.

Guillaume

@imbeauf
Copy link
Collaborator Author

imbeauf commented Oct 17, 2025

Ok, energy added again

@github-actions
Copy link

@imbeauf imbeauf requested a review from DavidPCoster October 17, 2025 14:31
imbeauf and others added 2 commits November 3, 2025 15:56
Co-authored-by: Simon Pinches <SimonPinches@users.noreply.github.com>
Co-authored-by: Simon Pinches <SimonPinches@users.noreply.github.com>
@github-actions
Copy link

github-actions bot commented Nov 3, 2025

@github-actions
Copy link

github-actions bot commented Nov 3, 2025

@github-actions
Copy link

@imbeauf
Copy link
Collaborator Author

imbeauf commented Nov 13, 2025

What do you think about integrating the two remaining ongoing PRs ?

  • Issue-85: I talked to Guillaume yesterday and have updated the related PR, still needs review
  • IMAS-5328 (additions to gyrokinetics): I have the impression that the CI is happy, but I am not sure this includes the compilation of AL. Can this be tested ?
    Of course, those two PRs are not essential for tag 4.1, but if they can be concluded rapidly it's worth including them.

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment

Labels

None yet

Projects

None yet

Development

Successfully merging this pull request may close these issues.

Adding new structures in the IMAS dictionary / distributions IDS

8 participants