anchor.diff:          densf = interpol_densf([lat,lon], alt_table, -90._r8, 90._r8, -180._r8, 180._r8, nlat_alt, nlon_alt)
       end function smooth_andean_mountain
   
  +
   end module refinter
  
positive.code:             dens_f=b**4/gammas**4
          end if
  
      end function dens_andes
  
      !--------------------------------------------------------------------------
      !--------------------------------------------------------------------------
  
negative.code:         pi, &
          datadir, &
          pio2, &
          pi2, &
          r8, &
          rad2deg, &
          deg2rad, &
  