From 4a2b5fc6382154bede23c15e32787284f919b47f Mon Sep 17 00:00:00 2001
From: Lorenzo Figini <lorenzo.figini@istp.cnr.it>
Date: Mon, 19 Dec 2016 16:38:07 +0000
Subject: [PATCH] added missing allocation for simple limiter (when it is not
 read from eqdsk)

---
 src/main.f90 | 6 +++++-
 1 file changed, 5 insertions(+), 1 deletion(-)

diff --git a/src/main.f90 b/src/main.f90
index 8e980f6..c889d76 100644
--- 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