diff --git a/src/geometry_mod.F90 b/src/geometry_mod.F90 index 02a9edb5562bc4d220bb0f61fae9dab2bef497a6..4703bf1cfd12e70dc64dbbf6310e7fec846495ef 100644 --- a/src/geometry_mod.F90 +++ b/src/geometry_mod.F90 @@ -63,7 +63,15 @@ contains Ckxky(ikx,iky,iz) = -SIN(z)*kx -(COS(z)+shear*z*SIN(z))*ky ENDDO ENDDO - + IF (Nky .EQ. 1) THEN ! linear 1D run we switch kx and ky for parallel opt + DO iky = ikys, ikye + ky = kyarray(iky) + DO ikx= ikxs, ikxe + kx = kxarray(ikx) + Ckxky(ikx,iky,iz) = -SIN(z)*ky -(COS(z)+shear*z*SIN(z))*kx + ENDDO + ENDDO + ENDIF ENDDO zloop END SUBROUTINE eval_s_alpha_geometry