diff --git a/src/incompressible/igrid.F90 b/src/incompressible/igrid.F90 index f07ce9ec..997b5bbb 100644 --- a/src/incompressible/igrid.F90 +++ b/src/incompressible/igrid.F90 @@ -542,6 +542,11 @@ subroutine init(this,inputfile, initialize2decomp) call decomp_info_init(nx, ny, nz, this%gpC) end if + if (any(this%gpC%xsz == 1) .or. any(this%gpC%ysz == 1) .or. any(this%gpC%zsz == 1))then + if(this%useWindTurbines)then + call gracefulExit("Pencil thickness = 1 detected in gpC. Wind turbine module may fail.", 901) + end if + end if call decomp_info_init(nx,ny,nz+1,this%gpE) if (this%useSystemInteractions) then