anchor.diff:  
   end subroutine shcalc
   
  -!***************************************************************************************************
  +!*******************************************************************************
   
   module alfbasismodule
   
  
positive.code: 
      implicit none
  
      character(128)              :: datafilenew, datafilelast=''
      real(4)                     :: epochnew, epochlast=-999.0
      real(8)                     :: we0, we1
      integer(4)                  :: iepoch0, iepoch1, iterm, icoord
  
negative.code: 
  end subroutine alfbasis
  
  !*******************************************************************************
  