anchor.diff:      !-- End: Raise cloud top if avg RH>RHSHmax and CAPE>0
       !
       !---------------------------SHALLOW CLOUD TOP---------------------------
  +            ! BK 2022/06/28: Warning and adjustment in case of low model top:
  +            if (LTOP<2) then 
  +                write(*,*) "   CU_BMJ WARNING: model top likely too low for correct convection simulation."!, LTOP, LBOT, LTP1,"[", this_image(),"]"
  +                LTOP=max(LTOP, 2)
  +            endif
               LBM1=LBOT-1
               PTPK=PTOP
               LTP1=LTOP-1
               DEPTH=PBOT-PTOP
  +
  +            ! BK 2022/06/28: Prevent LTP1 from going to zero:
  +        !   LTP1=max(LTP1,1) ! No longer needed with the above 'if' statement.
       !-----------------------------------------------------------------------
       !***  Begin debugging convection
               IF(PRINT_DIAG)THEN
  
positive.code:     MODULE MODULE_CU_BMJ
      !
      !-----------------------------------------------------------------------
              USE MODULE_MODEL_CONSTANTS
      !-----------------------------------------------------------------------
      !
              REAL,PARAMETER ::                                                 &
  
negative.code:             IF(L<=LMID) TREFK(L)=MAX(TREFK(L), TK(L)+DTSHAL)
              APEKXX=APEKXY
              PKXXXX=PKXXXY
              APEKXY=APEK(L-1)  ! error in ICAR; index below lower bound of 1
              PKXXXY=PK(L-1)
              ENDDO
      !
  