anchor.diff:            ! local_heights are pressures in hPa at layer centers (I should get rid of this variable)
             DO ilevel = 1, CmETA ! for layer center
                layer_press(ilevel) = (lpre(ilevel-1) + lpre(ilevel))/2.0 ! Pa
  -             local_heights(ixtrack,itimes,ilevel)  = layer_press(ilevel)/1.E2 ! hPa
  +             local_heights(ilevel)  = layer_press(ilevel)/1.E2 ! hPa
                clima_layer_press(ilevel) = (clima_lpre(ilevel-1) + clima_lpre(ilevel))/2.0 ! Pa
                clima_local_heights(ilevel)  = clima_layer_press(ilevel)/1.E2 ! hPa
                re_tmp(ilevel) = REAL(Temperature(idx_lon,idx_lat,CmETA+1-ilevel),KIND=r8)
  
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:        ! Now it is only needed to interpolate to the pixels of the granule.
         ! It was read there to obtain the dimensions of the number of levels.
         ! ---------------------------------------------------------------------
         CALL omi_climatology (climatology, cli_psurface, cli_temperature, terrain_height, &
              lat, lon, nt, nx, xtrange, locerrstat)
         
         ! -------------------------------------
  