diff --git a/src/geometry_mod.F90 b/src/geometry_mod.F90 index 7105088dadc3279d17df84ea96554b9d4312b4af..41c42936a76e2fc50a7511000bd1f22f19f0fada 100644 --- a/src/geometry_mod.F90 +++ b/src/geometry_mod.F90 @@ -369,7 +369,9 @@ CONTAINS SELECT CASE(parallel_bc) CASE ('dirichlet')! connected to 0 ikx_zBC_L(iky,ikx) = -99 - CASE ('periodic') !reroute it by cycling through modes + CASE ('periodic') + ikx_zBC_L(iky,ikx) = ikx + CASE ('cyclic')! reroute it by cycling through modes ikx_zBC_L(iky,ikx) = MODULO(ikx_zBC_L(iky,ikx)-1,Nkx)+1 END SELECT ENDIF @@ -395,8 +397,10 @@ CONTAINS SELECT CASE(parallel_bc) CASE ('dirichlet') ! connected to 0 ikx_zBC_R(iky,ikx) = -99 - CASE ('periodic') !reroute it by cycling through modes - write(*,*) 'check',ikx,iky, kxarray(ikx) + shift, '>', kx_max + CASE ('periodic') ! connected to itself as for shearless + ikx_zBC_R(iky,ikx) = ikx + CASE ('cyclic') + ! write(*,*) 'check',ikx,iky, kxarray(ikx) + shift, '>', kx_max ikx_zBC_R(iky,ikx) = MODULO(ikx_zBC_R(iky,ikx)-1,Nkx)+1 END SELECT ENDIF