anchor.diff:  ! File: anaxv.f90
   !
   ! This file is part of the SpanLib library.
  -! Copyright (C) 2012  Stephane Raynaud
  +! Copyright (C) 2012-2015  Stephane Raynaud
   ! Contact: stephane dot raynaud at gmail dot com
   !
   ! This library is free software; you can redistribute it and/or
  
positive.code:         end if
      end do
      close(11)
      
  end subroutine read_anaxv_field
  
  
  
negative.code:     integer :: it, igap
      real :: zx(ns),zn
      logical :: bad(ns)
  
      ! Write
      zn = -1.
      open(11,file=trim(fname),form='unformatted')
      do it = 1, nt
          zx = x(:, it)
          bad = zx==mv
          where(bad)zx = sqrt(zn)
          if(all(bad))then
              igap = 0
          elseif(any(bad))then
  