anchor.diff:              !! 'nesting_boundary_flux_*')
   
           !
  -        ! Output file content
  +        ! Output file content -- bespoke binary format (that is no longer properly supported, consider removing)
           !
           character(len=charlen):: metadata_ascii_filename
               !! Name of ascii file where we output metadata
  
positive.code:         !
          ! Setup for Chezy or Manning friction
          !
          if(domain%friction_type == 'manning') then
              ! For efficiency we ensure the friction power is known at compile time
              ! Previously it was set here
  
          else if(domain%friction_type == 'chezy') then
              ! For efficiency we ensure the friction power is known at compile time
              ! Previously it was set here
  
              if(domain%timestepping_method == 'cliffs') then
                  write(log_output_unit, *) "ERROR: chezy_friction has not yet been implemented for our cliffs solver"
                  call generic_stop
              end if
  
          else
              write(log_output_unit, *) ' Error: Unrecognized value of domain%friction_type '
              call generic_stop
          end if
  
  #ifdef DEBUG_ARRAY
          ! Optional for debugging -- this will be written to the netcdf file
  
negative.code:             end do
  
              close(i)
          end if
  #else
  
              ! Save the nontemporal variables that require tracking the flow
              do j = 1, size(domain%max_U_variables)
                  call domain%nc_grid_output%store_static_variable(domain%max_U_variables(j), domain%max_U(:,:,j))
              end do
  
              ! Save manning^2
              if(allocated(domain%manning_squared) .and. &
                  any(domain%nontemporal_grids_to_store == 'manning_squared')) then
                  call domain%nc_grid_output%store_static_variable('manning_squared', domain%manning_squared)
              end if
  
              ! Save elevation
              if(any(domain%nontemporal_grids_to_store == 'elevation0')) then
                  call domain%nc_grid_output%store_static_variable('elevation0', domain%U(:,:,ELV))
              end if
  
  #endif
  
      end subroutine
  
  