anchor.diff:  ! Local variables
   !
      real(8), dimension(nlev+1) ::  ppack, tpack, hpack, qpack, theta, omega_pack
  -   real(8), dimension(nlev+1) ::  t_updated, q_updated
  +   real(8), dimension(nlev+1) ::  t_updated, q_updated, p_updated
      real(8)                    ::  pblp, pbl_theta, pblh, latent_heat, sensible_heat   
      real(8)                    ::  pblp_updated, pressure_deficit, pw
      integer                    ::  xx, yy, dd, tt, ee
  
positive.code: !----------------------------------------------------------------------------------------------
  subroutine is_there_moist_convection(temperature, height, pressure, qhum, pblp,          &
                                       sensible   , latent, omega   , dt  , nlev, missing, &
                                       pressure_deficit, t_updated, q_updated, pblp_updated )
  
        integer, intent(in )  ::  nlev                   !*** # of atmospheric levels
        real(8), intent(in )  ::  temperature(nlev)      !*** temperature profile [K]
  
negative.code:             !-----------------------------------------------------------------------------
              tpack  =  t_updated
              qpack  =  q_updated
              ppack  =  p_updated
              pblp   =  pblp_updated
  
           end do time_of_day
  