added missing allocation for simple limiter (when it is not read from eqdsk)

2016-12-19 16:38:07 +00:00 · 2016-12-19 16:38:07 +00:00 · 4a2b5fc638
commit 4a2b5fc638
parent 0f2ff9a5a7
1 changed files with 5 additions and 1 deletions
--- a/src/main.f90
+++ b/src/main.f90
@ -99,7 +99,11 @@ program main_std
  dzmx=rtrp%dst*rtrp%nstep*0.01_wp_
  z0m=z0*0.01_wp_
  if (.not.allocated(rlim).or.rtrp%ipass<0) then
-    rtrp%ipass=abs(rtrp%ipass)
+    if (allocated(rlim)) deallocate(rlim)
+    if (allocated(zlim)) deallocate(zlim)
+    allocate(rlim(5))
+    allocate(zlim(5))
+    if (rtrp%ipass<0) rtrp%ipass = -rtrp%ipass
    if(eqp%iequil<2) then
      rmxm=(rv(1)+rv(2))*0.01_wp_
    else