Stressdiag

Open Scottyujia opened this issue 9 months ago • 1 comments

Purpose: The zonal flow dynamics is controlled by various turbulent stresses and the new features added are meant to calculate them.

Modifications:

fields.fpp

Functions similar to get_dgdx are added. Flags like 'phi', 'apar' and 'bpar' can be passed into the functions to calculate d/dx and d/dy in Fourier space for phi, apar and bpar separately, instead of simply calculating derivatives of chi.

2.stella_diagnostics.f90

The calculations of stresses are added. In the part that calculates omega, I also added the calculation of dphidt which is needed in the post-processing.

3.stella_io.fpp

Write the stresses to the nedcdf file.

Mar 17 '25 19:03 Scottyujia

Please merge the latest master into this before opening a pull request. This is based off of an old version of Stella, and is no longer compatible with the master.

Mar 27 '25 13:03 GeorgiaActon