Step 1 (forward): normalise + Z periodic post-process Z is dim1 (serial j loop), periodic R2C — no sign flip needed since j only goes to nz/2+1. nx, ny are the local loop bounds (equal to nx_glob, ny_glob on a single rank); the normalisation divisor always needs the globals.
| Type | Intent | Optional | Attributes | Name | ||
|---|---|---|---|---|---|---|
| complex(kind=dp), | intent(inout), | device, dimension(:, :, :) | :: | div_u | ||
| integer, | intent(in), | value | :: | nz_h | ||
| integer, | intent(in), | value | :: | nx | ||
| integer, | intent(in), | value | :: | ny | ||
| integer, | intent(in), | value | :: | nz | ||
| integer, | intent(in), | value | :: | nx_glob | ||
| integer, | intent(in), | value | :: | ny_glob | ||
| real(kind=dp), | intent(in), | device, dimension(:) | :: | az | ||
| real(kind=dp), | intent(in), | device, dimension(:) | :: | bz |