anchor.diff:      
     END SUBROUTINE COMPUTE_SCATT
   
  -  SUBROUTINE compute_amf ( nt, nx, CmETA, climatology, &
  +  SUBROUTINE compute_amf ( nt, nx, CmETA, climatology, climatology_sd, &
                              scattw, saoamf, amfdiag, errstat)
   
       IMPLICIT NONE
  
positive.code:              aircolumn          = rho*lhgt*Navogadro*1.0E-4 ! # air/cm^2
               
               ! -------------------------------------------------------------
               climatology(ixtrack,itimes,n) = aircolumn * lgas(n) /1.0E9 ! [GAS]/cm^2
            END DO
            
            !  Set non-physical entries to zero.
            WHERE ( climatology(ixtrack,itimes,1:CmETA) < 0.0_r8 )
               climatology(ixtrack,itimes,1:CmETA) = 0.0_r8
            END WHERE
            
         END DO
      END DO
  
negative.code:     REAL    (KIND=r8), DIMENSION (1:nx,0:nt-1)       :: amfgeo
      REAL    (KIND=r8), DIMENSION (1:nx,0:nt-1)       :: l2cfr, l2ctp
      REAL    (KIND=r8), DIMENSION (1:nx,0:nt-1)       :: albedo, cli_psurface
      REAL    (KIND=r8), DIMENSION (1:nx,0:nt-1,CmETA) :: climatology, climatology_sd
      REAL    (KIND=r8), DIMENSION (1:nx,0:nt-1,CmETA) :: scattw !, akernels
  
  
  