Thompson semi lag rain

.gitmodules

@@ -8,5 +8,5 @@
   branch = main
 [submodule "ccpp/physics"]
   path = ccpp/physics
-  url = https://github.com/NCAR/ccpp-physics
-  branch = main
+  url = https://github.com/RuiyuSun/ccpp-physics
+  branch = Thompson_semi_lag_rain

ccpp/data/GFS_typedefs.F90

@@ -825,6 +825,9 @@ module GFS_typedefs
     logical              :: ext_diag_thompson !< flag for extended diagnostic output from Thompson
     integer              :: thompson_ext_ndiag3d=37 !< number of 3d arrays for extended diagnostic output from Thompson
     real(kind=kind_phys) :: dt_inner        !< time step for the inner loop in s
+    logical              :: sedi_semi       !< flag for semi Lagrangian sedi of rain
+    logical              :: sedi_semi_update!< flag for v update in semi Lagrangian sedi of rain
+    logical              :: sedi_semi_decfl !< flag for interation with semi Lagrangian sedi of rain
 
     !--- GFDL microphysical paramters
     logical              :: lgfdlmprad      !< flag for GFDL mp scheme and radiation consistency
@@ -3257,6 +3260,9 @@ subroutine control_initialize (Model, nlunit, fn_nml, me, master, &
     real(kind=kind_phys) :: ttendlim       = -999.0             !< temperature tendency limiter, set to <0 to deactivate
     logical              :: ext_diag_thompson = .false.         !< flag for extended diagnostic output from Thompson
     real(kind=kind_phys) :: dt_inner       = -999.0             !< time step for the inner loop 
+    logical              :: sedi_semi      = .false.            !< flag for semi Lagrangian sedi of rain
+    logical              :: sedi_semi_update = .false.          !< flag for v update in semi Lagrangian sedi of rain
+    logical              :: sedi_semi_decfl = .false.           !< flag for interation with semi Lagrangian sedi of rain
 
     !--- GFDL microphysical parameters
     logical              :: lgfdlmprad     = .false.            !< flag for GFDLMP radiation interaction
@@ -3616,6 +3622,7 @@ subroutine control_initialize (Model, nlunit, fn_nml, me, master, &
                                mg_alf,   mg_qcmin, mg_do_ice_gmao, mg_do_liq_liu,           &
                                ltaerosol, lradar, nsradar_reset, lrefres, ttendlim,         &
                                ext_diag_thompson, dt_inner, lgfdlmprad,                     &
+                               sedi_semi, sedi_semi_update, sedi_semi_decfl,                &
                           !--- max hourly
                                avg_max_length,                                              &
                           !--- land/surface model control
@@ -4071,6 +4078,9 @@ subroutine control_initialize (Model, nlunit, fn_nml, me, master, &
     else
       Model%dt_inner       = Model%dtp
     endif
+    Model%sedi_semi        = sedi_semi
+    Model%sedi_semi_update = sedi_semi_update
+    Model%sedi_semi_decfl  = sedi_semi_decfl
 !--- F-A MP parameters
     Model%rhgrd            = rhgrd
     Model%spec_adv         = spec_adv
@@ -5149,6 +5159,9 @@ subroutine control_initialize (Model, nlunit, fn_nml, me, master, &
                                           ' ttendlim =',Model%ttendlim, &
                                           ' ext_diag_thompson =',Model%ext_diag_thompson, &
                                           ' dt_inner =',Model%dt_inner, &
+                                          ' sedi_semi=',Model%sedi_semi, & 
+                                          ' sedi_semi_update=',sedi_semi_update, & 
+                                          ' sedi_semi_decfl=',sedi_semi_decfl, &
                                           ' effr_in =',Model%effr_in, &
                                           ' lradar =',Model%lradar, &
                                           ' nsradar_reset =',Model%nsradar_reset, &
@@ -5567,6 +5580,9 @@ subroutine control_print(Model)
         print *, ' ttendlim          : ', Model%ttendlim
         print *, ' ext_diag_thompson : ', Model%ext_diag_thompson
         print *, ' dt_inner          : ', Model%dt_inner
+        print *, ' sedi_semi         : ', Model%sedi_semi
+        print *, ' sedi_semi_update  : ', Model%sedi_semi_update
+        print *, ' sedi_semi_decfl  : ', Model%sedi_semi_decfl
         print *, ' '
       endif
       if (Model%imp_physics == Model%imp_physics_mg) then

ccpp/data/GFS_typedefs.meta

@@ -3490,6 +3490,30 @@
   dimensions = ()
   type = real
   kind = kind_phys
+[sedi_semi]
+  standard_name = flag_for_semi_Lagrangian_sedi_rain
+  long_name = flag for semi Lagrangian sedi of rain
+  units = flag
+  dimensions = ()
+  type = logical
+  intent = in
+  optional = F
+[sedi_semi_update]
+  standard_name = flag_for_v_update_in_semi_Lagrangian_sedi
+  long_name = flag for v update in semi Lagrangian sedi of rain
+  units = flag
+  dimensions = ()
+  type = logical
+  intent = in
+  optional = F
+[sedi_semi_decfl]
+  standard_name = flag_for_iteration_with_semi_Lagrangian_sedi
+  long_name = flag for interation with semi Lagrangian sedi of rain
+  units = flag
+  dimensions = ()
+  type = logical
+  intent = in
+  optional = F
 [lgfdlmprad]
   standard_name = flag_for_GFDL_microphysics_radiation_interaction
   long_name = flag for GFDL microphysics-radiation interaction