anchor.diff:    !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
   
     do n = 1, nTimeSteps
  +    if (myid .eq. 0) then
   
       wind_x = base_wind_x*wind_mag_time_series(n)
       wind_y = base_wind_y*wind_mag_time_series(n)
  
positive.code:       etanew(i,j) =  values( ((i-1)*ny + j) )
      end do
    end do
  
  #endif
  
  
negative.code:     ! Apply the boundary conditions
      call apply_boundary_conditions(unew, hfacW, wetmask, nx, ny, layers)
      call apply_boundary_conditions(vnew, hfacS, wetmask, nx, ny, layers)
      end if
  
      ! Do the isopycnal layer physics
      if (.not. RedGrav) then
  