anchor.diff:            DO ilevel = 1, CmETA
                CALL ezspline_1d_interpolation (INT(CmETA,KIND=i4), clima_local_heights, &
                     re_tmp, &
  -                  one, local_heights(ixtrack,itimes,ilevel), &
  +                  one, local_heights(ilevel), &
                     ltmp(ilevel), status)
                CALL ezspline_1d_interpolation (INT(CmETA,KIND=i4), clima_local_heights, &
                     re_gas, &
  -                  one, local_heights(ixtrack,itimes,ilevel), &
  +                  one, local_heights(ilevel), &
                     lgas(ilevel), status)
                CALL ezspline_1d_interpolation (INT(CmETA,KIND=i4), clima_local_heights, &
                     re_h2o, &
  -                  one, local_heights(ixtrack,itimes,ilevel), &
  +                  one, local_heights(ilevel), &
                     lh2o(ilevel), status)
             END DO         
            
  
positive.code:           saoamf(ixtrack,itimes) = SUM(scattw(ixtrack, itimes, 1:CmETA) * &
                                          climatology(ixtrack,itimes,1:CmETA))     / &
                                     SUM(climatology(ixtrack,itimes,1:CmETA))
         END DO ! Finish xtrack pixel loop
      END DO ! Finish 
      
  
negative.code:     
    END SUBROUTINE COMPUTE_SCATT
  
    SUBROUTINE compute_amf ( nt, nx, CmETA, climatology, cli_temperature, cli_psurface, &
                             scattw, saoamf, amfdiag, errstat)
  
      IMPLICIT NONE
  