anchor.diff:    ! ====================================================================
     LOGICAL                                      :: yn_amf_wfmod
     INTEGER                                      :: amf_wfmod_idx
  -  REAL(KIND=r8)                                :: amf_alb_lnd, amf_alb_sno, amf_wvl, amf_wvl2, amf_alb_cld, amf_max_sza
  +  REAL(KIND=r8)                                :: amf_alb_lnd, amf_alb_sno, amf_wvl, &
  +       amf_wvl2, amf_alb_cld, amf_max_sza
   
     ! ---------------------------------------
     ! Data obtained from the climatology file
  
positive.code:     ! This should be moved to the program that creates the
      ! look up tables
      ! ------------------------------- ------------------------
      DO issza = 1, sza_dim(1)
         re_sza(issza) = cos(d2r*REAL(lut_sza(sza_dim(1)+1-issza), KIND = r8))  ! JED fix
      END DO
      DO isvza = 1, vza_dim(1)
         re_vza(isvza) = cos(d2r*REAL(lut_vza(vza_dim(1)+1-isvza), KIND = r8)) ! JED fix
      END DO
      DO ispre = 1, clp_dim(1)
         re_pre(ispre) = REAL(lut_clp(clp_dim(1)+1-ispre), KIND = r8)
      END DO
      DO isalt = 1, alt_lay_dim(1)
         re_alt(isalt) = REAL(lut_pre_lay(isalt), KIND=8)
      END DO
  
      ! ---------------
  
negative.code:           ! First compute back from the parametrization the scattering
            ! weights for the given wavelength and albedo.
            ! ----------------------------------------------------------          
            DO iozo = 1, 1
               DO iwav = iwavs(1), iwavf(1) !vl_nwav
  
                  ! ---------------------------
                  ! Initialize scattwe to zeros
  