anchor.diff:      character(30) wetMaskFile
   
   !   Numerics
  -    double precision pi,au,ah(layers),ar,dt
  +    double precision pi, dt
  +    double precision au,ah(layers),ar, botDrag
       double precision slip
       double precision hmin
       integer nTimeSteps
  
positive.code: !   Set default values here
  !   Possibly wait until the model is split into multiple files, then hide the long unsightly code there.
  
      NAMELIST /NUMERICS/ au,ah,ar,dt,slip,nTimeSteps,dumpFreq, & 
                          avFreq,hmin, maxits, freesurfFac, eps
  
      NAMELIST /MODEL/ hmean, depthFile, H0, RedGrav
  
  
negative.code: !     Calculate dhdt, dudt, dvdt at current time step
      call evaluate_dhdt(dhdtold, h,u,v,ah,dx,dy,nx,ny,layers, spongeHTimeScale,spongeH,wetmask)
  
      call evaluate_dudt(dudtold, h,u,v,b,zeta,wind_x,fu, & 
          au,ar,slip,dx,dy,hfacN,hfacS,nx,ny,layers,rho0, &
          spongeUTimeScale,spongeU,RedGrav,botDrag)
  
      !if the wind is changed to be meridional this code will need changing
  
      call evaluate_dvdt(dvdtold, h,u,v,b,zeta,wind_y,fv, & 
          au,ar,slip, dx,dy,hfacW,hfacE,nx,ny,layers,rho0, &
          spongeVTimeScale,spongeV,RedGrav,botDrag)
  
  !    Calculate the values at half the time interval with Forward Euler
      hhalf = h+0.5d0*dt*dhdtold
  