anchor.diff:               
                ! -------------------------------------------------------------
                climatology(ixtrack,itimes,n) = aircolumn * lgas(n) / 1.0E9 ! [GAS]/cm^2
  +             climatology_sd(ixtrack,itimes,n) = aircolumn * lgas_sd(n) / 1.0E9 ! [GAS]/cm^2
             END DO
             
             !  Set non-physical entries to zero.
  
positive.code:              climatology_sd(ixtrack,itimes,1:CmETA) = 0.0_r8
            END WHERE
  
            
         END DO
      END DO
  
  
negative.code:              clima_layer_press(ilevel) = DLOG((clima_lpre(ilevel-1) + clima_lpre(ilevel))/2.0) ! Pa
               re_tmp(ilevel) = Temperature(idx_lon,idx_lat,CmETA+1-ilevel)
               re_gas(ilevel) = Gas_profiles(idx_lon,idx_lat,CmETA+1-ilevel)
               re_gas_sd(ilevel) = Gas_profiles_sd(idx_lon,idx_lat,CmETA+1-ilevel)
               re_h2o(ilevel) = H2O_profiles(idx_lon,idx_lat,CmETA+1-ilevel)
            END DO
  
            ! ----------------------------------------------------------
            ! Interpolate temperature, gas, gas standard deviation
            ! and h2o profiles from clima_local_heights to local_heights
            ! ----------------------------------------------------------
            DO ilevel = 1, CmETA
               ltmp(ilevel) = linInterpol( CmETA, clima_layer_press(1:CmETA), &
                    re_tmp(1:CmETA), layer_press(ilevel), status=status)
               lgas(ilevel) = linInterpol( CmETA, clima_layer_press(1:CmETA), &
                    re_gas(1:CmETA), layer_press(ilevel), status=status)
               lgas_sd(ilevel) = linInterpol( CmETA, clima_layer_press(1:CmETA), &
                    re_gas_sd(1:CmETA), layer_press(ilevel), status=status)
               lh2o(ilevel) = linInterpol( CmETA, clima_layer_press(1:CmETA), &
                    re_h2o(1:CmETA), layer_press(ilevel), status=status)
            END DO         
  