anchor.diff:  
           t1 = trim(domain%output_folder_name) // '/' // 'Grid_output_ID' // trim(t3) // '.nc'
           ! Character attributes
  -        natt = 7!8
  +        natt = 7
           allocate(attribute_names(natt), attribute_values(natt))
           !
           attribute_names(1) = 'timestepping_method'
  
positive.code:     use stop_mod, only: generic_stop
      use iso_fortran_env, only: output_unit, int32, int64
      use iso_c_binding, only: c_ptr, c_null_ptr
      use netcdf_util, only: nc_grid_output_type
      use logging_mod, only: log_output_unit
      use file_io_mod, only: mkdir_p
      use cliffs_tolkova_mod, only: cliffs, setSSLim
  
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
  
  