anchor.diff:          real(dp) :: dx_refinement_factor_of_parent_domain = ONE_dp
               !! This information is useful in the context of nesting, where interior
               !! domains have cell sizes that are an integer divisor of the
  -            !! coarse-domain dx. 
  +            !! coarse-domain dx.
           integer(ip) :: timestepping_refinement_factor = 1_ip
               !! Number of timesteps taken by this domain inside each multidomain timestep.
               !! Only used in conjunction with the multidomain class
  
positive.code:             !END DO
              !$OMP END DO
              !$OMP END PARALLEL
              
  EVOLVE_TIMER_STOP('update_max_quantities')
          end if
  
  
negative.code:         integer(ip):: i, metadata_unit, natt
          character(len=charlen), allocatable :: attribute_names(:), attribute_values(:)
  
  
          call domain%create_output_folders()
  
          ! Get domain id as a character
  