anchor.diff:           local_srf, local_cld, local_cfr
       REAL    (KIND=r8), DIMENSION(ngeos5)         :: local_lev_pre
       REAL    (KIND=r8), DIMENSION(ngeos5-1)       :: local_lay_pre
  -    REAL    (KIND=r8), DIMENSION(srf_dim(1),alt_lay_dim(2)) :: re_pre_2D
  -    REAL    (KIND=r8), DIMENSION(alt_lay_dim(2))            :: lay_pre_1D
  +    REAL    (KIND=r8), DIMENSION(srf_dim(1),alt_lay_dim(2)) :: re_pre_2D, re_pre_cld_2D
  +    REAL    (KIND=r8), DIMENSION(alt_lay_dim(2))            :: lay_pre_1D, lay_pre_cld_1D
       REAL    (KIND=r8), PARAMETER :: d2r = 3.141592653589793d0/180.0  !! JED fix
   
       ! To select TOMS profile
  
positive.code:          dI0_clr_dim, hdferr)
  
      CALL h5dread_f(dI0_cld_did, H5T_NATIVE_REAL,  &
           lut_dI0_cld(1:toz_dim(1),1:srf_dim(1),1:alt_lay_dim(2),1:clp_dim(1),1:vza_dim(1),1:sza_dim(1)), &
           dI0_cld_dim, hdferr)
      CALL h5dread_f(dI1_cld_did, H5T_NATIVE_REAL,  &
           lut_dI1_cld(1:toz_dim(1),1:srf_dim(1),1:alt_lay_dim(2),1:clp_dim(1),1:vza_dim(1),1:sza_dim(1)), &
           dI0_cld_dim, hdferr)
      CALL h5dread_f(dI2_cld_did, H5T_NATIVE_REAL,  &
           lut_dI2_cld(1:toz_dim(1),1:srf_dim(1),1:alt_lay_dim(2),1:clp_dim(1),1:vza_dim(1),1:sza_dim(1)), &
           dI0_cld_dim, hdferr)
      
      CALL h5dread_f(air_did,     H5T_NATIVE_REAL, lut_air(1:toz_dim(1),1:srf_dim(1),1:alt_lay_dim(2)), &
  
negative.code:          dI0_clr_did, dI1_clr_did, dI2_clr_did, & 
           dI0_cld_did, dI1_cld_did, dI2_cld_did, & 
           air_did, alt_lev_did, alt_lay_did, pre_lev_did, pre_lay_did, ozo_did, tem_did, &
           dspace, toms_datatype_id, pre_lev_cld_did, pre_lay_cld_did
  
      INTEGER(SIZE_T)                :: size
      LOGICAL, SAVE :: h5inited = .FALSE.
  