anchor.diff:              !-- Pack the omega variable to be same dimension as T,p,q, and h
               !-- Need to put this variable within the time loop to ensure that
               !-- omega varies throughout the day
  +            !-- Update the pressure field using the omega [Pa/s]
               !-----------------------------------------------------------------------------
               call packIt(omega(dd,(tt/18)+1,:,yy,xx) , 0.0 , nlev   , nlev1      , &
                           psfc (dd            ,yy,xx) , P   , missing, omega_pack   )
  
positive.code:             call is_there_moist_convection(tpack, hpack, ppack, qpack, pblp,          &
                                             sensible_heat, latent_heat,                &
                                             omega_pack, dt, nlev1, missing, &
                                             pressure_deficit, t_updated, q_updated, pblp_updated )
  
              !-----------------------------------------------------------------------------
              !----- Output the desired variables when moist convective initiation occurs
  
negative.code:       real(8), intent(out)  ::  pressure_deficit       !*** a measure saturation at the top of PBL [Pa]
        real(8), intent(out)  ::  t_updated  (nlev)      !*** updated temperature profile [K]
        real(8), intent(out)  ::  q_updated  (nlev)      !*** updated specific humidity profile [kg/kg]
        real(8), intent(out)  ::  p_updated  (nlev)      !*** updated pressure profile [Pa]
        real(8), intent(out)  ::  pblp_updated           !*** new pressure at the top of the boundary layer [Pa]
  
        real(8)               ::  updated_qhum(nlev), updated_temp(nlev), updated_pblp, updated_p(nlev)
        real(8)               ::  tbm, tdef, bclp, qbcl
  
  
        !--------------------------------------------------------------------------------------
        !----- New omega incorporation is to adjust the pressure levels by the omega profoile
        !----------------------------------------------------------------------------------------
        !call most_influential_layer( omega, updated_pblp, pressure, nlev, missing, strong_omega )
        call updated_pressure_with_omega(omega, pressure, dt, nlev,  missing, updated_p)
  
  
        !-----------------------------------------------------------------------------
        !----- Updated the profiles of temperature and humidity after injecting
        !----- the surface senisble and latent heat fluxes
        !-----------------------------------------------------------------------------
        call add_surface_fluxes (temperature, height, updated_p, qhum, pblp, sensible, latent, dt, nlev, missing, &
                                 updated_temp, updated_qhum, updated_pblp )
  
        !-----------------------------------------------------------------------------
        !----- Calculate HCF variables in new system
        !-----------------------------------------------------------------------------
        call hcfcalc ( nlev, missing, updated_temp, updated_p, updated_qhum, height, pblp, tbm, tdef, bclp, qbcl )
  
  
        !-----------------------------------------------------------------------------
        !----- Check for CI
        !-----------------------------------------------------------------------------
        pressure_deficit  =  updated_pblp - bclp
        t_updated         =  updated_temp
        q_updated         =  updated_qhum
        pblp_updated      =  updated_pblp
        p_updated         =  updated_p
  
  end subroutine is_there_moist_convection
  
  