anchor.diff:    public
   
     integer, parameter :: dp = selected_real_kind(15,307)
  +  integer,parameter :: k15 = selected_int_kind(15)
     real(kind=dp), private, parameter::  fmwh2o = 18.01534d0
   
   contains
  
positive.code: !---bb(0-9,6)
                              9*0.d0        , 5.235718623d02/
      
      data xi1/4.260321148d0/
      data xl0,xl1,xl2/15.74373327d0,-34.17061978d0,19.31380707d0/
      
  !   save n,ll,x,xi1,xl0,xl1,xl2,z,b,bb,delt,delp
      
  
negative.code:     real(kind=dp), intent(out) :: hg,hgp,hgt
      integer, intent(out) :: ierr
    
      integer(kind = k15), save :: n(8),ll(8),x(8,2),z(8,3),xi1,xl0,xl1,xl2
      real(kind=dp), save :: b(8,2),bb(0:9,0:6)
      real(kind=dp) :: sumbx(8),sumbxt(8)
    
  