anchor.diff:      md%domains(3)%timestepping_method = 'rk2'
   
       ! The monai domain 
  -    ! (Peak stage seems to be affected by the westward extent
  -    ! of this, possibly mesh design needs more tuning)
       call md%domains(4)%match_geometry_to_parent(&
           parent_domain=md%domains(3), &
           lower_left  = [139.38_dp, 42.0828_dp], &
  
positive.code: 
      if(very_high_res_monai) then
          ! An even more detailed Monai domain
          ! Hard to get this one stable, although it's ok with EULER
          call md%domains(7)%match_geometry_to_parent(&
              parent_domain=md%domains(5), &
              !lower_left  = [139.4176_dp, 42.0935_dp], &
  
negative.code:             dx_refinement_factor = nest_ratio, &
              timestepping_refinement_factor = 6_ip,&
              rounding_method='nearest')
          ! Older versions of the code needed euler for stability. Does it still?
          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
  