grans_obj_PrepareSolverBlock Subroutine

    subroutine grans_obj_PrepareSolverBlock(this)
      !> Prepares block for run.
      ! This subroutine reads pre-existing block file on disk and
      ! partitions domain based on parameters in the input file.
      implicit none
      class(grans_obj), intent(inout) :: this                                  !! GRANS solver
      ! Work variables
      integer :: Nb(3)
      integer :: ngc
      character(len=str64) :: filename

      call this%parser%Get("Block file",  filename, default = 'block'  )
      call this%parser%Get("Partition",   Nb,       default = [1,1,1]  )
      call this%parser%Get("Ghost cells", ngc,      default = 1        )

      ! Initialize and partition blocks
      call this%block%Initialize(ngc,this%parallel)
      call this%block%Read(filename)
      call this%block%Partition(Nb)

      return
    end subroutine grans_obj_PrepareSolverBlock