Calculate moments (density, parallel flow, and parallel and perpendicular temperatures) in non-guiding centre coordinates and write to netCDF
subroutine do_write_full_moments_notgc()
use diagnostics_moments, only: getmoms_notgc
use theta_grid, only: ntgrid
use kt_grids, only: naky, ntheta0
use species, only: nspec
use mp, only: proc0
use gs2_io, only: nc_loop_fullmom
use fields_arrays, only: phinew, bparnew
implicit none
complex, dimension (-ntgrid:ntgrid,ntheta0,naky,nspec) :: ntot, density, &
upar, tpar, tperp
call getmoms_notgc(phinew, bparnew, density,upar,tpar,tperp,ntot)
if(proc0) then
call nc_loop_fullmom(get_netcdf_file_id(), nout, &
& ntot(igomega,:,:,:),density(igomega,:,:,:), &
& upar(igomega,:,:,:), &
& tpar(igomega,:,:,:),tperp(igomega,:,:,:) )
endif
end subroutine do_write_full_moments_notgc