initial implementation of new tendency application control options

Author: grantfirl

ccpp/physics

@@ -282,6 +282,10 @@ module CCPP_typedefs
     real (kind=kind_phys), pointer      :: stress_land(:)     => null()  !<
     real (kind=kind_phys), pointer      :: stress_water(:)    => null()  !<
     real (kind=kind_phys), pointer      :: t2mmp(:)           => null()  !<
+    real (kind=kind_phys), pointer      :: ten_q(:,:,:)       => null()
+    real (kind=kind_phys), pointer      :: ten_t(:,:)         => null()
+    real (kind=kind_phys), pointer      :: ten_u(:,:)         => null()
+    real (kind=kind_phys), pointer      :: ten_v(:,:)         => null()
     real (kind=kind_phys), pointer      :: theta(:)           => null()  !<
     real (kind=kind_phys), pointer      :: tlvl(:,:)          => null()  !<
     real (kind=kind_phys), pointer      :: tlyr(:,:)          => null()  !<
@@ -630,6 +634,10 @@ subroutine gfs_interstitial_create (Interstitial, IM, Model)
     allocate (Interstitial%stress_ice      (IM))
     allocate (Interstitial%stress_land     (IM))
     allocate (Interstitial%stress_water    (IM))
+    allocate (Interstitial%ten_q           (IM,Model%levs,Model%ntrac))
+    allocate (Interstitial%ten_t           (IM,Model%levs))
+    allocate (Interstitial%ten_u           (IM,Model%levs))
+    allocate (Interstitial%ten_v           (IM,Model%levs))
     allocate (Interstitial%theta           (IM))
     allocate (Interstitial%tlvl            (IM,Model%levr+1+LTP))
     allocate (Interstitial%tlyr            (IM,Model%levr+LTP))
@@ -1312,6 +1320,10 @@ subroutine gfs_interstitial_phys_reset (Interstitial, Model)
     Interstitial%stress_ice      = Model%huge
     Interstitial%stress_land     = Model%huge
     Interstitial%stress_water    = Model%huge
+    Interstitial%ten_q           = clear_val
+    Interstitial%ten_t           = clear_val
+    Interstitial%ten_u           = clear_val
+    Interstitial%ten_v           = clear_val
     Interstitial%theta           = clear_val
     Interstitial%tprcp_ice       = Model%huge
     Interstitial%tprcp_land      = Model%huge

scm/src/CCPP_typedefs.F90

@@ -2026,6 +2026,34 @@
   type = real
   kind = kind_phys
   active = (control_for_land_surface_scheme == identifier_for_noahmp_land_surface_scheme)
+[ten_q]
+  standard_name = tendency_of_tracer_concentration
+  long_name = tendency of tracer concentration calculated by one physics scheme
+  units = kg kg-1 s-1
+  dimensions = (horizontal_loop_extent,vertical_layer_dimension,number_of_tracers)
+  type = real
+  kind = kind_phys
+[ten_t]
+  standard_name = tendency_of_air_temperature
+  long_name = tendency of air temperature calculated by one physics scheme
+  units = K s-1
+  dimensions = (horizontal_loop_extent,vertical_layer_dimension)
+  type = real
+  kind = kind_phys
+[ten_u]
+  standard_name = tendency_of_x_wind
+  long_name = tendency of x wind calculated by one physics scheme
+  units = m s-2
+  dimensions = (horizontal_loop_extent,vertical_layer_dimension)
+  type = real
+  kind = kind_phys
+[ten_v]
+  standard_name = tendency_of_y_wind
+  long_name = tendency of y wind calculated by one physics scheme
+  units = m s-2
+  dimensions = (horizontal_loop_extent,vertical_layer_dimension)
+  type = real
+  kind = kind_phys
 [theta]
   standard_name = angle_from_east_of_maximum_subgrid_orographic_variations
   long_name = angle with_respect to east of maximum subgrid orographic variations

scm/src/CCPP_typedefs.meta

@@ -792,6 +792,18 @@ module GFS_typedefs
                                             !< (yr, mon, day, t-zone, hr, min, sec, mil-sec)
     integer              :: idate(4)        !< initial date with different size and ordering
                                             !< (hr, mon, day, yr)
+!--- tendency control
+    integer              :: tend_opt_swrad
+    integer              :: tend_opt_lwrad
+    integer              :: tend_opt_rad_scaler
+    integer              :: tend_opt_surface
+    integer              :: tend_opt_pbl
+    integer              :: tend_opt_gwd
+    integer              :: tend_opt_photochem
+    integer              :: tend_opt_deep_conv
+    integer              :: tend_opt_shal_conv
+    integer              :: tend_opt_mp
+    integer              :: tend_opt_stoch
 !--- radiation control parameters
     real(kind=kind_phys) :: fhswr           !< frequency for shortwave radiation (secs)
     real(kind=kind_phys) :: fhlwr           !< frequency for longwave radiation (secs)
@@ -3375,6 +3387,23 @@ subroutine control_initialize (Model, nlunit, fn_nml, me, master, &
     integer              :: thermodyn_id   =  1              !< valid for GFS only for get_prs/phi
     integer              :: sfcpress_id    =  1              !< valid for GFS only for get_prs/phi
 
+    !--- time-coupling options after a scheme completes
+    ! 1 = immediately apply tendencies
+    ! 2 = add tendencies to a sum to be applied later
+    ! 3 = add tendencies to a sum and apply the accumulated sum to the state
+    ! 4 = ignore output tendencies (e.g. some other scheme may use/apply them)
+    integer              :: tend_opt_swrad      = 4
+    integer              :: tend_opt_lwrad      = 4
+    integer              :: tend_opt_rad_scaler = 2
+    integer              :: tend_opt_surface    = 2
+    integer              :: tend_opt_pbl        = 2
+    integer              :: tend_opt_gwd        = 3
+    integer              :: tend_opt_photochem  = 1
+    integer              :: tend_opt_deep_conv  = 1
+    integer              :: tend_opt_shal_conv  = 1
+    integer              :: tend_opt_mp         = 1
+    integer              :: tend_opt_stoch      = 1
+
     !--- coupling parameters
     logical              :: cplflx         = .false.         !< default no cplflx collection
     logical              :: cplice         = .false.         !< default no cplice collection (used together with cplflx)
@@ -4045,6 +4074,11 @@ subroutine control_initialize (Model, nlunit, fn_nml, me, master, &
                                fhzero, fhzero_array, fhzero_fhour, ldiag3d, qdiag3d, lssav, &
                                naux2d, dtend_select, naux3d, aux2d_time_avg,                &
                                aux3d_time_avg, fhcyc, thermodyn_id, sfcpress_id,            &
+                          !--- tendency application controls
+                               tend_opt_swrad, tend_opt_lwrad, tend_opt_rad_scaler,         &
+                               tend_opt_surface, tend_opt_pbl, tend_opt_gwd,                &
+                               tend_opt_photochem, tend_opt_deep_conv, tend_opt_shal_conv,  &
+                               tend_opt_mp, tend_opt_stoch,                                 &
                           !--- coupling parameters
                                cplflx, cplice, cplocn2atm, cplwav, cplwav2atm, cplaqm,      &
                                cplchm, cpllnd, cpllnd2atm, cpl_imp_mrg, cpl_imp_dbg,        &
@@ -4444,7 +4478,18 @@ subroutine control_initialize (Model, nlunit, fn_nml, me, master, &
         Model%chunk_begin(i) = Model%chunk_end(i-1) + 1
         Model%chunk_end(i) = Model%chunk_begin(i) + blksz(i) - 1
     end do
-    
+!--- tendency controls
+    Model%tend_opt_swrad      = tend_opt_swrad
+    Model%tend_opt_lwrad      = tend_opt_lwrad
+    Model%tend_opt_rad_scaler = tend_opt_rad_scaler
+    Model%tend_opt_surface    = tend_opt_surface
+    Model%tend_opt_pbl        = tend_opt_pbl
+    Model%tend_opt_gwd        = tend_opt_gwd
+    Model%tend_opt_photochem  = tend_opt_photochem
+    Model%tend_opt_deep_conv  = tend_opt_deep_conv
+    Model%tend_opt_shal_conv  = tend_opt_shal_conv
+    Model%tend_opt_mp         = tend_opt_mp
+    Model%tend_opt_stoch      = tend_opt_stoch  
 !--- coupling parameters
     Model%cplflx           = cplflx
     Model%cplice           = cplice