anchor.diff:  
   
     ! Physics
  -  double precision :: g_vec(layers)
  +  double precision :: g_vec(layerwise_input_length)
     double precision :: rho0
   
     ! Wind
  
positive.code:   ! Numerics
    double precision :: pi, dt
    double precision :: au, ar, botDrag
    double precision :: ah(layers)
    double precision :: slip
    double precision :: hmin
    integer nTimeSteps
  
negative.code:   double precision :: thickness_error
  
    ! Model
    double precision :: hmean(layerwise_input_length)
    ! Switch for using n + 1/2 layer physics, or using n layer physics
    logical :: RedGrav
  
  