From 0f2ff9a5a755bb4a9c6aac4e5f3ee97f3d94b652 Mon Sep 17 00:00:00 2001
From: Lorenzo Figini <lorenzo.figini@istp.cnr.it>
Date: Tue, 29 Nov 2016 08:45:18 +0000
Subject: [PATCH] fix: missing allocation for profiles given vs rhot

---
 src/main.f90 | 6 ++++--
 1 file changed, 4 insertions(+), 2 deletions(-)

diff --git a/src/main.f90 b/src/main.f90
index ccfc85a..8e980f6 100644
--- a/src/main.f90
+++ b/src/main.f90
@@ -57,11 +57,13 @@ program main_std
   else
     call read_profiles(prfp%filenm, xrad, terad, derad, zfc)
     allocate(psrad(size(xrad)))
-    if(prfp%irho==0) then
+    if(prfp%irho==0) then ! xrad==rhot
+      allocate(rhot(size(psinr)))
       call setqphi_num(psinr,qpsi,psia,rhot)
       call set_rhospl(sqrt(psinr),rhot)
+      deallocate(rhot)
       psrad=frhopolv(xrad)
-    else if(prfp%irho == 1) then
+    else if(prfp%irho == 1) then ! xrad==rhop
       psrad=xrad**2
     else
       psrad=xrad