MODULE module_shallowcu_driver
CONTAINS
   SUBROUTINE shallowcu_driver( &
                      ids,ide, jds,jde, kds,kde &
                     ,ims,ime, jms,jme, kms,kme &
                     ,ips,ipe, jps,jpe, kps,kpe &
                     ,i_start,i_end,j_start,j_end,kts,kte,num_tiles &
                     ,u,v,th,t &
                     ,p,pi,rho,moist &
                     ,num_moist &
                     ,itimestep,dt,dx,cudt,curr_secs,adapt_step_flag &
                     ,rainsh,pratesh,nca,rainshv &
                     ,z,z_at_w,dz8w,mavail,pblh,p8w &
                     ,tke_pbl &
                     ,cldfra,cldfra_old,cldfrash &
                     ,htop,hbot &
                     ,shcu_physics &
                     ,qv_curr, qc_curr, qr_curr &
                     ,qi_curr, qs_curr, qg_curr &
                     ,dlf, rliq, rliq2 &
                     ,cmfmc, cmfmc2 &
                     ,cush, snowsh, icwmrsh, rprdsh, cbmf, cmfsl &
                     ,cmflq, evapcsh &
                     ,rqvshten,rqcshten,rqrshten &
                     ,rqishten,rqsshten,rqgshten &
                     ,rqvblten,rqvften &
                     ,rushten,rvshten &
                     ,rthshten,rthraten,rthblten,rthften &
                     ,f_qv,f_qc,f_qr &
                     ,f_qi,f_qs,f_qg &
                     ,ht &
                                                                      )
   USE module_model_constants
   USE module_state_description, ONLY: CAMUWSHCUSCHEME, G3SHCUSCHEME
   USE module_cu_camuwshcu_driver, ONLY : camuwshcu_driver
   USE module_dm
   USE module_domain, ONLY: domain
   IMPLICIT NONE
   INTEGER, INTENT(IN ) :: &
                                      ids,ide, jds,jde, kds,kde, &
                                      ims,ime, jms,jme, kms,kme, &
                                                        kts,kte, &
                                      itimestep, num_tiles
   INTEGER, DIMENSION(num_tiles), INTENT(IN) :: &
     & i_start,i_end,j_start,j_end
   INTEGER, INTENT(IN ) :: &
                  num_moist
   INTEGER, INTENT(IN ) :: shcu_physics
   REAL, DIMENSION( ims:ime, kms:kme, jms:jme, num_moist ), &
         INTENT(INOUT) :: &
                                                          moist
   REAL, DIMENSION( ims:ime, kms:kme, jms:jme ), &
         INTENT(IN ) :: &
                                                         cldfra &
                                                    ,cldfra_old &
                                                      , z &
                                                      , z_at_w &
                                                      , dz8w &
                                                      , p8w &
                                                      , p &
                                                      , pi &
                                                      , u &
                                                      , v &
                                                      , th &
                                                      , t &
                                                      , tke_pbl &
                                                      , rho
   REAL, DIMENSION( ims:ime , jms:jme ), INTENT(IN) :: &
                  MAVAIL,PBLH,ht
   REAL, DIMENSION( ims:ime , jms:jme ), &
          INTENT(INOUT) :: RAINSH &
                                                    , NCA &
                                                    , HTOP &
                                                    , HBOT
   REAL, DIMENSION( ims:ime , jms:jme ),INTENT(INOUT),OPTIONAL :: &
        PRATESH, RAINSHV
   REAL, DIMENSION( ims:ime , jms:jme ) :: tmppratesh
   REAL, INTENT(IN ) :: DT, DX
   INTEGER, INTENT(IN ),OPTIONAL :: &
                                      ips,ipe, jps,jpe, kps,kpe
   REAL, INTENT(IN ),OPTIONAL :: CUDT
   REAL, INTENT(IN ),OPTIONAL :: CURR_SECS
   LOGICAL,INTENT(IN ),OPTIONAL :: adapt_step_flag
   REAL :: cudt_pass, curr_secs_pass
   LOGICAL :: adapt_step_flag_pass
   REAL, DIMENSION( ims:ime, kms:kme, jms:jme ), &
         OPTIONAL, INTENT(INOUT) :: &
                      qv_curr, qc_curr, qr_curr &
                     ,qi_curr, qs_curr, qg_curr &
                     ,rqvshten,rqcshten,rqrshten &
                     ,rqishten,rqsshten,rqgshten &
                     ,rqvblten,rqvften &
                     ,rthraten,rthblten &
                     ,rthften,rushten,rvshten,rthshten
   REAL, DIMENSION( ims:ime , jms:jme ), &
                    OPTIONAL, INTENT(INOUT) :: &
                rliq, rliq2 &
               ,cbmf, cush, snowsh
   REAL, DIMENSION( ims:ime, kms:kme, jms:jme ), &
         OPTIONAL, INTENT(INOUT) :: &
                  cldfrash, cmfsl, cmflq, icwmrsh, &
                  dlf, evapcsh, &
                  cmfmc, cmfmc2, rprdsh
   LOGICAL, INTENT(IN), OPTIONAL :: &
                                                      f_qv &
                                                     ,f_qc &
                                                     ,f_qr &
                                                     ,f_qi &
                                                     ,f_qs &
                                                     ,f_qg
   INTEGER :: i,j,k,its,ite,jts,jte,ij
   CHARACTER(len=200) :: message
   if (.not. PRESENT(CURR_SECS)) then
      curr_secs_pass = -1
   else
      curr_secs_pass = curr_secs
   endif
   if (.not. PRESENT(CUDT)) then
      cudt_pass = -1
   else
      cudt_pass = cudt
   endif
   if (.not. PRESENT(adapt_step_flag)) then
      adapt_step_flag_pass = .false.
   else
      adapt_step_flag_pass = adapt_step_flag
   endif
   if ( PRESENT ( pratesh ) ) then
      tmppratesh(:,:) = pratesh(:,:)
   else
      tmppratesh(:,:) = 0.
   end if
   IF (shcu_physics .eq. 0) return
!$OMP PARALLEL DO &
!$OMP PRIVATE ( ij ,its,ite,jts,jte, i,j,k)
   DO ij = 1 , num_tiles
      its = i_start(ij)
      ite = i_end(ij)
      jts = j_start(ij)
      jte = j_end(ij)
   scps_select: SELECT CASE(shcu_physics)
   CASE (G3SHCUSCHEME)
   CASE (CAMUWSHCUSCHEME)
      CALL wrf_debug(100,'in camuw_scps')
      IF(.not.f_qi)THEN
         WRITE( message , * ) 'This shallow cumulus option requires ice microphysics option: f_qi = ', f_qi
         CALL wrf_error_fatal3("<stdin>",370,&
message )
      ENDIF
      CALL camuwshcu_driver( &
            IDS=ids,IDE=ide,JDS=jds,JDE=jde,KDS=kds,KDE=kde &
           ,IMS=ims,IME=ime,JMS=jms,JME=jme,KMS=kms,KME=kme &
           ,ITS=its,ITE=ite,JTS=jts,JTE=jte,KTS=kts,KTE=kte &
           ,NUM_MOIST=num_moist, DT=dt &
           ,P=p, P8W=p8w, PI_PHY=pi &
           ,Z=z, Z_AT_W=z_at_w, DZ8W=dz8w &
           ,T_PHY=t, U_PHY=u, V_PHY=v &
           ,MOIST=moist, QV=qv_curr, QC=qc_curr, QI=qi_curr &
           ,PBLH_IN=pblh, TKE_PBL=tke_pbl &
           ,CLDFRA=cldfra, CLDFRA_OLD=cldfra_old, CLDFRASH=cldfrash &
           ,CUSH_INOUT=cush, RAINSH=rainsh, PRATESH=tmppratesh &
           ,SNOWSH=snowsh &
           ,ICWMRSH=icwmrsh, CMFMC=cmfmc, CMFMC2_INOUT=cmfmc2 &
           ,RPRDSH_INOUT=rprdsh, CBMF_INOUT=cbmf &
           ,CMFSL=cmfsl, CMFLQ=cmflq, DLF=dlf, EVAPCSH_INOUT=evapcsh &
           ,RLIQ=rliq, RLIQ2_INOUT=rliq2, CUBOT=hbot, CUTOP=htop &
           ,RUSHTEN=rushten, RVSHTEN=rvshten, RTHSHTEN=rthshten &
           ,RQVSHTEN=rqvshten, RQCSHTEN=rqcshten, RQRSHTEN=rqrshten &
           ,RQISHTEN=rqishten, RQSSHTEN=rqsshten, RQGSHTEN=rqgshten &
           ,HT=ht &
                                                                         )
   CASE DEFAULT
      WRITE( message , * ) 'The shallow cumulus option does not exist: shcu_physics = ', shcu_physics
      CALL wrf_error_fatal3("<stdin>",398,&
message )
   END SELECT scps_select
   ENDDO
   !$OMP END PARALLEL DO
   if (PRESENT(PRATESH)) then
      pratesh(:,:) = tmppratesh(:,:)
      if (PRESENT(RAINSHV)) then
         rainshv(:,:) = pratesh(:,:)*dt
      endif
   endif
   END SUBROUTINE shallowcu_driver
END MODULE module_shallowcu_driver