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:     !***
      !***  CLOUD MUST BE AT LEAST TWO LAYERS THICK
      !***
      !        IF(LBOT-LTOP<2)LTOP=LBOT-2  (eliminate this criterion)
      !
      !-- End: Buoyancy check (24 Aug 2006)
      !
  
negative.code:             APEKXY=APEK(LBOT)
      !
              LMID=.5*(LBOT+LTOP)
              ! if((LTOP<2).OR.(LBOT<2)) write(*,*)"   LBOT,LTOP   ",LBOT,LTOP
              DO L=LBOT,LTOP,-1
              TREFKX=((PKXXXY-PKXXXX)*SMIX                                    &
              &          +TREFKX*APEKXX)/APEKXY
  