anchor.diff:      end subroutine moving_average2
   
   
  -
  -
  -
       !--------------------------------------------------------------------------
       ! Andes density data for interpolation
       !--------------------------------------------------------------------------
  -
  -
       subroutine andes_density_table(alt_table)
           real (kind=8), dimension(nlat_alt*nlon_alt,3), intent(inout) :: alt_table
           integer :: i, j
  
positive.code:         real (kind=8), dimension(1:n_lon)  :: y
          real (kind=8), dimension(n_lat*n_lon, 3)  :: table
  
  
          xmin = latmin
          xmax = latmax
          ymin = lonmin
  
negative.code:         real (kind=8), dimension(nlat_alt*nlon_alt,3), intent(in) :: dens_table
          real (kind=8) :: lat, lon,latc, lonc, coslat
          real (kind=8) :: epsilons, gammas, sx, maxdist, dists, a, b
  
          latc=-20.0*pi/180._r8
          lonc=-60.0*pi/180._r8
  