anchor.diff:                           PF_POST_BLOCK,       &
                            fecho_output_invariants)
           
  +        call pf_print_options(pf,un_opt=6)
  +        
           ! advance in time with libpfasst
           level = nlevs
   
  
positive.code:         pf%nlevels           = nlevs                         ! number of SDC levels
          pf%niters            = niters                        ! number of SDC iterations
          pf%save_timings      = 1                             ! output the timings in fort.601 file
          qtype                = translate_qtype(qtype_name, & ! select the type of nodes
                                              qnl)
  
          if (nlevs == 1) then
              nvars = [nfields*nvars_per_field(1)]    ! number of degrees of freedom for the levels
          else 
              stop 'This number of levels is not supported yet'
          end if
  
          ! initialize level-specific data structures
  
negative.code:         allocate(sweet_level_t::pf%levels(level)%ulevel)
          allocate(sweet_data_factory_t::pf%levels(level)%ulevel%factory)
          allocate(sweet_sweeper_t::pf%levels(level)%ulevel%sweeper)
  
          ! cast the object into sweet data objects
          sd_factory_ptr    => as_sweet_data_factory(pf%levels(level)%ulevel%factory)
  