Commit 72b9a18e authored by Gregor Michalicek's avatar Gregor Michalicek

Bugfix: Calculate stepfunction also if sliceplot%iplot is true

...it is needed for the qfix call in optional/pldngen.
parent 60c0f1c8
......@@ -178,15 +178,12 @@
CALL prp_xcfft(stars,input, cell, xcpot)
!
ENDIF ! (mpi%irank == 0)
CALL stepf(sym,stars,atoms,oneD, input,cell, vacuum,mpi)
IF (.NOT.sliceplot%iplot) THEN
!
CALL stepf(sym,stars,atoms,oneD, input,cell, vacuum,mpi)
!
IF ( mpi%irank == 0 ) THEN
CALL convn(DIMENSION,atoms,stars)
!
!---> set up electric field parameters (if needed)
!
! CALL e_field(atoms, DIMENSION, stars, sym, vacuum, cell, input,field)
ENDIF
ENDIF
......
......@@ -540,8 +540,8 @@ SUBROUTINE postprocessInput(mpi,input,field,sym,stars,atoms,vacuum,obsolete,kpts
CALL MPI_BCAST(sliceplot%iplot,1,MPI_LOGICAL,0,mpi%mpi_comm,ierr)
#endif
IF (.NOT.sliceplot%iplot) THEN
CALL stepf(sym,stars,atoms,oneD,input,cell,vacuum,mpi)
CALL stepf(sym,stars,atoms,oneD,input,cell,vacuum,mpi)
IF (.NOT.sliceplot%iplot) THEN
IF (mpi%irank.EQ.0) THEN
CALL convn(DIMENSION,atoms,stars)
CALL e_field(atoms,DIMENSION,stars,sym,vacuum,cell,input,field%efield)
......
Markdown is supported
0% or
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment