fix: missing allocation for profiles given vs rhot

2016-11-29 08:45:18 +00:00 · 2016-11-29 08:45:18 +00:00 · 0f2ff9a5a7
commit 0f2ff9a5a7
parent 5156f94f73
1 changed files with 4 additions and 2 deletions
--- 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