anchor.diff:  
     ! Numerics
     double precision :: pi
  -  integer :: nwrite, avwrite
     double precision :: rjac
  +
  +  ! dumping output
  +  integer :: nwrite, avwrite, checkpointwrite
   
     ! External solver variables
     integer          :: offsets(2,5)
  
positive.code:         dudt, dvdt, dhdt, &
          wind_x, wind_y, nx, ny, layers, &
          n, nwrite, avwrite, RedGrav, DumpWind, debug_level)
  
      cur_time = time()
      if (cur_time - last_report_time > 3) then
  
negative.code: 
  subroutine model_run(h, u, v, eta, depth, dx, dy, wetmask, fu, fv, &
      dt, au, ar, botDrag, ah, slip, hmin, nTimeSteps, dumpFreq, avFreq, &
      checkpointFreq, maxits, eps, freesurfFac, thickness_error, &
      debug_level, g_vec, rho0, &
      base_wind_x, base_wind_y, wind_mag_time_series, &
      spongeHTimeScale, spongeUTimeScale, spongeVTimeScale, &
      spongeH, spongeU, spongeV, &
  