anchor.diff:          ! Finish looping at some point
           if (md%domains(1)%time > final_time) exit
   
  -        if(very_high_res_monai) then
  -            ! Take a different time step once the high res domain starts evolving.
  -            ! This will introduce a (formal) first-order error into the linear leap-frog scheme,
  -            ! because that requires a fixed time step, but only at the change in time point.
  -            ! But practically no problem (?).
  -            if(md%domains(1)%time < md%domains(7)%static_before_time) then
  -                call md%evolve_one_step(global_dt)
  -            else
  -                call md%evolve_one_step(global_dt/very_high_res_timestep_reduction)
  -            end if
  -        else
  +        !if(very_high_res_monai) then
  +        !    ! Take a different time step once the high res domain starts evolving.
  +        !    ! This will introduce a (formal) first-order error into the linear leap-frog scheme,
  +        !    ! because that requires a fixed time step, but only at the change in time point.
  +        !    ! But practically no problem (?).
  +        !    if(md%domains(1)%time < md%domains(7)%static_before_time) then
  +        !        call md%evolve_one_step(global_dt)
  +        !    else
  +        !        call md%evolve_one_step(global_dt/very_high_res_timestep_reduction)
  +        !    end if
  +        !else
               ! Regular case
               call md%evolve_one_step(global_dt)
  -        end if
  +        !end if
   
       end do
   
  
positive.code:     md%domains(3)%timestepping_method = default_nonlinear_timestepping_method
      md%domains(3)%cliffs_minimum_allowed_depth = 1.0_dp
  
      ! The monai domain 
      call md%domains(4)%match_geometry_to_parent(&
          parent_domain=md%domains(3), &
          lower_left  = [139.38_dp, 42.0828_dp], &
  
negative.code:     md%domains(6)%timestepping_method = default_nonlinear_timestepping_method  
      md%domains(6)%cliffs_minimum_allowed_depth = 0.2_dp
  
      !if(very_high_res_monai) then
      !    ! An even more detailed Monai domain
      !    call md%domains(7)%match_geometry_to_parent(&
      !        parent_domain=md%domains(5), &
      !        !lower_left  = [139.4176_dp, 42.0935_dp], &
      !        lower_left  = [139.4230_dp, 42.09896_dp], &
      !        !upper_right = [139.4299_dp, 42.1046_dp], &
      !        upper_right = [139.4254_dp, 42.09995_dp], &
      !        dx_refinement_factor = nest_ratio, &
      !        timestepping_refinement_factor = 6_ip,&
      !        rounding_method='nearest')
      !    ! Older versions of the code needed euler for stability.
      !    md%domains(7)%timestepping_method = 'euler'  
      !    ! Should be dry before this time -- note also we change the time-step 
      !    ! in the evolve loop on the basis of this value
      !    md%domains(7)%static_before_time = very_high_res_static_before_time
      !end if
  
      ! Allocate domains and prepare comms
      call md%setup()
  