anchor.diff:      v = vnew
       if (.not. RedGrav) then
         eta = etanew
  -    endif
  +    end if
   
       ! Tendencies (old -> very old)
       dhdtveryold = dhdtold
  
positive.code:             form='unformatted')
          write(10) etaav(1:nx,1:ny)
          close(10)
        endif
  
        ! Check if there are NaNs in the data
        call break_if_NaN(h,nx,ny,layers,n)
  
negative.code:       print *, 'inconsistency between h and eta (in %):', &
            maxval(abs(h_norming - 1d0))*100d0
      end if
    end if
  
    !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
    !!!  Initialisation of the model STARTS HERE                              !!!
  