anchor.diff:          contains
           procedure :: initialise => initialise
           procedure :: write_grids => write_grids
  -        procedure :: store_max_stage => store_max_stage
  +        procedure :: store_static_variable => store_static_variable
           procedure :: store_priority_domain_cells => store_priority_domain_cells
           procedure :: finalise => finalise
   
  
positive.code:     implicit none
  
      private
      public :: nc_grid_output_type, check
  
      ! Indices of stage, uh, vh, elevation in domain%U
      integer(ip), parameter :: STG=1, UH=2, VH=3, ELV=4
      ! Note -- the above definitions must agree with those in domain_mod.f90
  
      type :: nc_grid_output_type
          !!
          !! Type for saving domain data to netcdf. Holds unit numbers, etc
  
negative.code:         integer :: dim_time_id, dim_x_id, dim_y_id
          integer :: var_time_id, var_x_id, var_y_id
  
  #ifdef DEBUG_ARRAY
          integer :: var_debug_array_id
              !! Time-varying array for debugging
  #endif
          integer :: var_is_priority_domain_id, var_priority_ind, var_priority_img
              !! Variables describing the priority domain
  
          integer :: flush_every_n_output_steps = 25
              !! It can be slow to flush every time we write to files. But
  