anchor.diff:      md%domains(1)%nx = global_nx
       md%domains(1)%timestepping_method = timestepping_method
   
  -    ! Allocate domain -- must have set timestepping method BEFORE this
       call md%setup
   
       ! Call local routine to set initial conditions
  
positive.code:     tank_length = sea_length + land_length
   
      ! Resolution
      dx = 0.1_dp/mesh_refine !(0.5_dp/12.0_dp)
  
      ! Model parameters. Note we apply the boundary on the right. x = 0 is the shoreline
      global_lw = [tank_length, tank_width]
  
negative.code:     real(dp), parameter :: mesh_refine = 1.0_dp
  
      ! Domain info
      character(charlen) :: timestepping_method, tempchar
      
      real(dp) :: global_lw(2), global_ll(2)
      integer(ip) :: global_nx(2) 
  