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