Added checks for rhill and radius during io operations

MintonGroup · Aug 6, 2021 · 730753e · 730753e
1 parent 883f152
commit 730753e
Showing 1 changed file with 3 additions and 3 deletions.
diff --git a/src/io/io.f90 b/src/io/io.f90
@@ -984,7 +984,7 @@ module subroutine io_read_frame_body(self, iu, param, form, ierr)
             read(iu, iostat=ierr, err=100) pl%Gmass(:)
             pl%mass(:) = pl%Gmass(:) / param%GU 
             if (param%lrhill_present) read(iu, iostat=ierr, err=100) pl%rhill(:)
-            read(iu, iostat=ierr, err=100) pl%radius(:)
+            if (param%lclose) read(iu, iostat=ierr, err=100) pl%radius(:)
             if (param%lrotation) then
                read(iu, iostat=ierr, err=100) pl%rot(1, :)
                read(iu, iostat=ierr, err=100) pl%rot(2, :)
@@ -1347,8 +1347,8 @@ module subroutine io_write_frame_body(self, iu, param)
          select type(pl => self)  
          class is (swiftest_pl)  ! Additional output if the passed polymorphic object is a massive body
             write(iu) pl%Gmass(1:n)
-            write(iu) pl%rhill(1:n)
-            write(iu) pl%radius(1:n)
+            if (param%lrhill_present) write(iu) pl%rhill(1:n)
+            if (param%lclose) write(iu) pl%radius(1:n)
             if (param%lrotation) then
                write(iu) pl%Ip(1, 1:n)
                write(iu) pl%Ip(2, 1:n)