diff --git a/src/equilibrium.f90 b/src/equilibrium.f90 index d02f8cc..9a92494 100644 --- a/src/equilibrium.f90 +++ b/src/equilibrium.f90 @@ -192,8 +192,8 @@ contains q(1)=q0 q(2)=qa q(3)=alq - - if(ipass.ge.2) then + + if(ipass.ge.2) then ! get size of boundary and limiter arrays and allocate them read (u,*) nlim @@ -205,9 +205,9 @@ contains allocate(rlim(nlim),zlim(nlim)) read(u,*) (rlim(i),zlim(i),i=1,nlim) end if - - end if - + + end if + close(u) end subroutine read_equil_an