From 223f09255ac19cae2d1c8ee64c37c8c30ef449d0 Mon Sep 17 00:00:00 2001
From: Dom Heinzeller <dom.heinzeller@icloud.com>
Date: Wed, 2 Sep 2020 13:13:37 -0600
Subject: [PATCH] Add dependencies to CCPP metadata (#153)

This PR removes dependency information from the CCPP prebuild config. See NCAR/ccpp-framework#308 and NCAR/ccpp-framework#317 for details on the motivation for this change and the actual implementation. It also removes some legacy code used by the dynamic CCPP build in the past.

This PR also contains the changes in #156, i.e. the completion of adding the active attribute to GFS_typedefs.F90. On top of this PR, the missing active attribute for phy_fctd is added.

Additionally, gfortran-10 compiler flags are added to CCPP's CMakeLists.txt.
---
 ccpp/CMakeLists.txt                      |   7 +-
 ccpp/config/ccpp_prebuild_config.py      | 190 ++----------
 ccpp/framework                           |   2 +-
 ccpp/physics                             |   2 +-
 gfsphysics/CCPP_layer/CCPP_data.meta     |   5 +
 gfsphysics/CCPP_layer/CCPP_typedefs.meta |  10 +
 gfsphysics/GFS_layer/GFS_typedefs.F90    |   3 +
 gfsphysics/GFS_layer/GFS_typedefs.meta   | 374 ++++++++++++++++++++++-
 8 files changed, 416 insertions(+), 177 deletions(-)

diff --git a/ccpp/CMakeLists.txt b/ccpp/CMakeLists.txt
index dd1c317e42..b1395b23bb 100644
--- a/ccpp/CMakeLists.txt
+++ b/ccpp/CMakeLists.txt
@@ -82,11 +82,12 @@ set (CMAKE_Fortran_FLAGS_OPT "")
 if (${CMAKE_Fortran_COMPILER_ID} MATCHES "GNU")
   set (CMAKE_Fortran_FLAGS  "${CMAKE_Fortran_FLAGS} -fcray-pointer -ffree-line-length-none -fno-range-check")
   set (CMAKE_Fortran_FLAGS  "${CMAKE_Fortran_FLAGS} -fbacktrace -cpp")
+  if(${CMAKE_Fortran_COMPILER_VERSION} VERSION_GREATER_EQUAL 10)
+    set(CMAKE_Fortran_FLAGS "${CMAKE_Fortran_FLAGS} -fallow-argument-mismatch")
+  endif()
   if (${CMAKE_BUILD_TYPE} MATCHES "Debug")
-    set (CMAKE_C_FLAGS        "${CMAKE_C_FLAGS} -Wall")
-    set (CMAKE_CXX_FLAGS      "${CMAKE_CXX_FLAGS} -Wall")
     set (CMAKE_Fortran_FLAGS  "${CMAKE_Fortran_FLAGS} -fno-unsafe-math-optimizations -frounding-math -fsignaling-nans")
-    set (CMAKE_Fortran_FLAGS  "${CMAKE_Fortran_FLAGS} -ffpe-trap=invalid,zero,overflow -fcheck=bounds -fbacktrace -fno-range-check -Wall")
+    set (CMAKE_Fortran_FLAGS  "${CMAKE_Fortran_FLAGS} -ffpe-trap=invalid,zero,overflow -fcheck=bounds -fbacktrace -fno-range-check")
   elseif (${CMAKE_BUILD_TYPE} MATCHES "Bitforbit")
     set (CMAKE_Fortran_FLAGS  "${CMAKE_Fortran_FLAGS}")
   endif (${CMAKE_BUILD_TYPE} MATCHES "Debug")
diff --git a/ccpp/config/ccpp_prebuild_config.py b/ccpp/config/ccpp_prebuild_config.py
index 0cc2fb605e..4d61fbbaed 100755
--- a/ccpp/config/ccpp_prebuild_config.py
+++ b/ccpp/config/ccpp_prebuild_config.py
@@ -11,8 +11,10 @@
 
 # Add all files with metadata tables on the host model side and in CCPP,
 # relative to basedir = top-level directory of host model. This includes
-# kind and type definitions used in CCPP physics.
+# kind and type definitions used in CCPP physics. Also add any internal
+# dependencies of these files to the list.
 VARIABLE_DEFINITION_FILES = [
+    # actual variable definition files
     'FV3/ccpp/physics/physics/machine.F',
     'FV3/ccpp/physics/physics/radsw_param.f',
     'FV3/ccpp/physics/physics/radlw_param.f',
@@ -85,128 +87,6 @@
         },
     }
 
-# Add all physics scheme dependencies relative to basedir - note that the CCPP
-# rules stipulate that dependencies are not shared between the schemes!
-SCHEME_FILES_DEPENDENCIES = [
-    'FV3/ccpp/physics/physics/GFDL_parse_tracers.F90',
-    'FV3/ccpp/physics/physics/aer_cloud.F',
-    'FV3/ccpp/physics/physics/aerclm_def.F',
-    'FV3/ccpp/physics/physics/aerinterp.F90',
-    'FV3/ccpp/physics/physics/calpreciptype.f90',
-    'FV3/ccpp/physics/physics/cldwat2m_micro.F',
-    'FV3/ccpp/physics/physics/cldmacro.F',
-    'FV3/ccpp/physics/physics/date_def.f',
-    'FV3/ccpp/physics/physics/funcphys.f90',
-    'FV3/ccpp/physics/physics/gcycle.F90',
-    'FV3/ccpp/physics/physics/gfs_phy_tracer_config.F',
-    'FV3/ccpp/physics/physics/gocart_tracer_config_stub.f',
-    'FV3/ccpp/physics/physics/h2o_def.f',
-    'FV3/ccpp/physics/physics/h2ointerp.f90',
-    'FV3/ccpp/physics/physics/iccn_def.F',
-    'FV3/ccpp/physics/physics/iccninterp.F90',
-    'FV3/ccpp/physics/physics/iounitdef.f',
-    'FV3/ccpp/physics/physics/machine.F',
-    'FV3/ccpp/physics/physics/mersenne_twister.f',
-    'FV3/ccpp/physics/physics/mfpbl.f',
-    'FV3/ccpp/physics/physics/micro_mg_utils.F90',
-    'FV3/ccpp/physics/physics/micro_mg2_0.F90',
-    'FV3/ccpp/physics/physics/micro_mg3_0.F90',
-    'FV3/ccpp/physics/physics/module_bfmicrophysics.f',
-    'FV3/ccpp/physics/physics/multi_gases.F90',
-    'FV3/ccpp/physics/physics/module_gfdl_cloud_microphys.F90',
-    'FV3/ccpp/physics/physics/module_nst_model.f90',
-    'FV3/ccpp/physics/physics/module_nst_parameters.f90',
-    'FV3/ccpp/physics/physics/module_nst_water_prop.f90',
-    'FV3/ccpp/physics/physics/module_mp_radar.F90',
-    'FV3/ccpp/physics/physics/module_mp_thompson.F90',
-    'FV3/ccpp/physics/physics/module_mp_thompson_make_number_concentrations.F90',
-    'FV3/ccpp/physics/physics/module_MP_FER_HIRES.F90',
-    'FV3/ccpp/physics/physics/module_bl_mynn.F90',
-    'FV3/ccpp/physics/physics/module_sf_mynn.F90',
-    'FV3/ccpp/physics/physics/module_SF_JSFC.F90',
-    'FV3/ccpp/physics/physics/module_BL_MYJPBL.F90',
-    'FV3/ccpp/physics/physics/module_sf_noahmp_glacier.f90',
-    'FV3/ccpp/physics/physics/module_sf_noahmplsm.f90',
-    'FV3/ccpp/physics/physics/cires_ugwp_module.F90',
-    'FV3/ccpp/physics/physics/ugwp_driver_v0.F',
-    'FV3/ccpp/physics/physics/cires_ugwp_triggers.F90',
-    'FV3/ccpp/physics/physics/cires_ugwp_initialize.F90',
-    'FV3/ccpp/physics/physics/cires_ugwp_solvers.F90',
-    'FV3/ccpp/physics/physics/cires_ugwp_utils.F90',
-    'FV3/ccpp/physics/physics/cires_orowam2017.f',
-    'FV3/ccpp/physics/physics/cires_vert_lsatdis.F90',
-    'FV3/ccpp/physics/physics/cires_vert_orodis.F90',
-    'FV3/ccpp/physics/physics/cires_vert_wmsdis.F90',
-    'FV3/ccpp/physics/physics/namelist_soilveg.f',
-    'FV3/ccpp/physics/physics/mfpblt.f',
-    'FV3/ccpp/physics/physics/mfpbltq.f',
-    'FV3/ccpp/physics/physics/mfscu.f',
-    'FV3/ccpp/physics/physics/mfscuq.f',
-    'FV3/ccpp/physics/physics/noahmp_tables.f90',
-    'FV3/ccpp/physics/physics/num_parthds.F',
-    'FV3/ccpp/physics/physics/ozne_def.f',
-    'FV3/ccpp/physics/physics/ozinterp.f90',
-    'FV3/ccpp/physics/physics/physcons.F90',
-    'FV3/ccpp/physics/physics/physparam.f',
-    'FV3/ccpp/physics/physics/radcons.f90',
-    'FV3/ccpp/physics/physics/radiation_aerosols.f',
-    'FV3/ccpp/physics/physics/radiation_astronomy.f',
-    'FV3/ccpp/physics/physics/radiation_clouds.f',
-    'FV3/ccpp/physics/physics/radiation_gases.f',
-    'FV3/ccpp/physics/physics/radiation_surface.f',
-    'FV3/ccpp/physics/physics/radlw_datatb.f',
-    'FV3/ccpp/physics/physics/radlw_param.f',
-    'FV3/ccpp/physics/physics/radsw_datatb.f',
-    'FV3/ccpp/physics/physics/radsw_param.f',
-    'FV3/ccpp/physics/physics/samfaerosols.F',
-    'FV3/ccpp/physics/physics/sfcsub.F',
-    'FV3/ccpp/physics/physics/sflx.f',
-    'FV3/ccpp/physics/physics/set_soilveg.f',
-    'FV3/ccpp/physics/physics/flake.F90',
-    'FV3/ccpp/physics/physics/surface_perturbation.F90',
-    'FV3/ccpp/physics/physics/cu_gf_deep.F90',
-    'FV3/ccpp/physics/physics/cu_gf_sh.F90',
-    'FV3/ccpp/physics/physics/tridi.f',
-    'FV3/ccpp/physics/physics/wv_saturation.F',
-    'FV3/ccpp/physics/physics/module_sf_ruclsm.F90',
-    'FV3/ccpp/physics/physics/namelist_soilveg_ruc.F90',
-    'FV3/ccpp/physics/physics/set_soilveg_ruc.F90',
-    'FV3/ccpp/physics/physics/module_soil_pre.F90',
-    # RRTMGP
-    'FV3/ccpp/physics/physics/rte-rrtmgp/rrtmgp/mo_gas_concentrations.F90',
-    'FV3/ccpp/physics/physics/rte-rrtmgp/rrtmgp/mo_gas_optics.F90',
-    'FV3/ccpp/physics/physics/rte-rrtmgp/rrtmgp/mo_gas_optics_rrtmgp.F90',
-    'FV3/ccpp/physics/physics/rte-rrtmgp/rrtmgp/mo_rrtmgp_constants.F90',
-    'FV3/ccpp/physics/physics/rte-rrtmgp/rrtmgp/mo_rrtmgp_util_reorder.F90',
-    'FV3/ccpp/physics/physics/rte-rrtmgp/rrtmgp/mo_rrtmgp_util_string.F90',
-    'FV3/ccpp/physics/physics/rte-rrtmgp/rrtmgp/kernels/mo_gas_optics_kernels.F90',
-    'FV3/ccpp/physics/physics/rte-rrtmgp/rrtmgp/kernels/mo_rrtmgp_util_reorder_kernels.F90',
-    'FV3/ccpp/physics/physics/rte-rrtmgp/rte/mo_fluxes.F90',
-    'FV3/ccpp/physics/physics/rte-rrtmgp/rte/mo_rte_util_array.F90',
-    'FV3/ccpp/physics/physics/rte-rrtmgp/rte/mo_optical_props.F90',
-    'FV3/ccpp/physics/physics/rte-rrtmgp/rte/mo_rte_kind.F90',
-    'FV3/ccpp/physics/physics/rte-rrtmgp/rte/mo_rte_lw.F90',
-    'FV3/ccpp/physics/physics/rte-rrtmgp/rte/mo_rte_sw.F90',
-    'FV3/ccpp/physics/physics/rte-rrtmgp/rte/mo_rte_config.F90',
-    'FV3/ccpp/physics/physics/rte-rrtmgp/rte/mo_source_functions.F90',
-    'FV3/ccpp/physics/physics/rte-rrtmgp/rte/kernels/mo_fluxes_broadband_kernels.F90',
-    'FV3/ccpp/physics/physics/rte-rrtmgp/rte/kernels/mo_optical_props_kernels.F90',
-    'FV3/ccpp/physics/physics/rte-rrtmgp/rte/kernels/mo_rte_solver_kernels.F90',
-    'FV3/ccpp/physics/physics/rte-rrtmgp/extensions/mo_compute_bc.F90',
-    'FV3/ccpp/physics/physics/rte-rrtmgp/extensions/mo_fluxes_byband.F90',
-    'FV3/ccpp/physics/physics/rte-rrtmgp/extensions/mo_fluxes_byband_kernels.F90',
-    'FV3/ccpp/physics/physics/rte-rrtmgp/extensions/mo_fluxes_bygpoint.F90',
-    'FV3/ccpp/physics/physics/rte-rrtmgp/extensions/mo_heating_rates.F90',
-    'FV3/ccpp/physics/physics/rte-rrtmgp/extensions/mo_rrtmgp_clr_all_sky.F90',
-    'FV3/ccpp/physics/physics/rte-rrtmgp/extensions/cloud_optics/mo_cloud_optics.F90',
-    'FV3/ccpp/physics/physics/rrtmg_lw_cloud_optics.F90'        ,
-    'FV3/ccpp/physics/physics/rrtmg_sw_cloud_optics.F90'        ,
-    'FV3/ccpp/physics/physics/rrtmgp_aux.F90'                   ,
-    # derived data type definitions
-    'FV3/gfsphysics/GFS_layer/GFS_typedefs.F90',
-    'FV3/gfsphysics/CCPP_layer/CCPP_typedefs.F90',
-    ]
-
 # Add all physics scheme files relative to basedir
 SCHEME_FILES = [
     # Relative path to source (from where ccpp_prebuild.py is called) : [ list of physics sets in which scheme may be called ];
@@ -300,25 +180,24 @@
     # HAFSFER_HIRES
     'FV3/ccpp/physics/physics/mp_fer_hires.F90',
     # RRTMGP
-    'FV3/ccpp/physics/physics/rrtmgp_lw_gas_optics.F90'         ,
-    'FV3/ccpp/physics/physics/rrtmgp_lw_cloud_optics.F90'       ,
-    'FV3/ccpp/physics/physics/rrtmgp_sw_gas_optics.F90'         ,
-    'FV3/ccpp/physics/physics/rrtmgp_sw_cloud_optics.F90'       ,
-    'FV3/ccpp/physics/physics/rrtmgp_sw_aerosol_optics.F90'     ,
-    'FV3/ccpp/physics/physics/rrtmgp_lw_rte.F90'                ,
-    'FV3/ccpp/physics/physics/rrtmgp_sw_rte.F90'                ,
-    'FV3/ccpp/physics/physics/rrtmgp_lw_aerosol_optics.F90'     ,
-    'FV3/ccpp/physics/physics/GFS_rrtmgp_setup.F90'             ,
-    'FV3/ccpp/physics/physics/GFS_rrtmgp_pre.F90'               ,
-    'FV3/ccpp/physics/physics/rrtmgp_lw_pre.F90'                ,
-    'FV3/ccpp/physics/physics/GFS_rrtmgp_sw_pre.F90'            ,
-    'FV3/ccpp/physics/physics/GFS_rrtmgp_lw_post.F90'           ,
-    'FV3/ccpp/physics/physics/rrtmgp_lw_cloud_sampling.F90' ,
-    'FV3/ccpp/physics/physics/rrtmgp_sw_cloud_sampling.F90' ,
-    'FV3/ccpp/physics/physics/GFS_cloud_diagnostics.F90'           ,
-    'FV3/ccpp/physics/physics/mo_cloud_sampling.F90'               ,
-    'FV3/ccpp/physics/physics/GFS_rrtmgp_gfdlmp_pre.F90'           ,
-    'FV3/ccpp/physics/physics/GFS_rrtmgp_zhaocarr_pre.F90'         ,
+    'FV3/ccpp/physics/physics/rrtmgp_lw_gas_optics.F90',
+    'FV3/ccpp/physics/physics/rrtmgp_lw_cloud_optics.F90',
+    'FV3/ccpp/physics/physics/rrtmgp_sw_gas_optics.F90',
+    'FV3/ccpp/physics/physics/rrtmgp_sw_cloud_optics.F90',
+    'FV3/ccpp/physics/physics/rrtmgp_sw_aerosol_optics.F90',
+    'FV3/ccpp/physics/physics/rrtmgp_lw_rte.F90',
+    'FV3/ccpp/physics/physics/rrtmgp_sw_rte.F90',
+    'FV3/ccpp/physics/physics/rrtmgp_lw_aerosol_optics.F90',
+    'FV3/ccpp/physics/physics/GFS_rrtmgp_setup.F90',
+    'FV3/ccpp/physics/physics/GFS_rrtmgp_pre.F90',
+    'FV3/ccpp/physics/physics/rrtmgp_lw_pre.F90',
+    'FV3/ccpp/physics/physics/GFS_rrtmgp_sw_pre.F90',
+    'FV3/ccpp/physics/physics/GFS_rrtmgp_lw_post.F90',
+    'FV3/ccpp/physics/physics/rrtmgp_lw_cloud_sampling.F90',
+    'FV3/ccpp/physics/physics/rrtmgp_sw_cloud_sampling.F90',
+    'FV3/ccpp/physics/physics/GFS_cloud_diagnostics.F90',
+    'FV3/ccpp/physics/physics/GFS_rrtmgp_gfdlmp_pre.F90',
+    'FV3/ccpp/physics/physics/GFS_rrtmgp_zhaocarr_pre.F90',
     'FV3/ccpp/physics/physics/GFS_rrtmgp_sw_post.F90'
     ]
 
@@ -336,13 +215,6 @@
 SCHEMES_CMAKEFILE  = '{build_dir}/ccpp/physics/CCPP_SCHEMES.cmake'
 SCHEMES_SOURCEFILE = '{build_dir}/ccpp/physics/CCPP_SCHEMES.sh'
 
-# CCPP host cap in which to insert the ccpp_field_add statements;
-# determines the directory to place ccpp_{modules,fields}.inc
-TARGET_FILES = [
-    'FV3/atmos_cubed_sphere/driver/fvGFS/atmosphere.F90',
-    'FV3/ccpp/driver/CCPP_Driver.F90',
-    ]
-
 # Auto-generated makefile/cmakefile snippets that contain all caps
 CAPS_MAKEFILE   = '{build_dir}/ccpp/physics/CCPP_CAPS.mk'
 CAPS_CMAKEFILE  = '{build_dir}/ccpp/physics/CCPP_CAPS.cmake'
@@ -453,12 +325,6 @@
     #'subroutine_name_2' : [ 'var1', 'var3'],
     }
 
-# Names of Fortran include files in the host model cap (do not change);
-# both files will be written to the directory of each target file, only
-# used by the dynamic builds
-MODULE_INCLUDE_FILE = 'ccpp_modules_{set}.inc'
-FIELDS_INCLUDE_FILE = 'ccpp_fields_{set}.inc'
-
 # Directory where to write static API to
 STATIC_API_DIR = '{build_dir}/ccpp/physics'
 STATIC_API_SRCFILE = '{build_dir}/ccpp/physics/CCPP_STATIC_API.sh'
@@ -471,17 +337,3 @@
 
 # LaTeX document containing the provided vs requested CCPP variables
 LATEX_VARTABLE_FILE = '{build_dir}/ccpp/framework/doc/DevelopersGuide/CCPP_VARIABLES_FV3.tex'
-
-
-###############################################################################
-# Template code to generate include files                                     #
-###############################################################################
-
-# Name of the CCPP data structure in the host model cap;
-# in the case of FV3, this is a 2-dimensional array with
-# the number of blocks as the first and the number of
-# OpenMP threads as the second dimension; nb is the loop
-# index for the current block, nt for the current thread.
-# Internally, the model uses an associate construct to
-# reference cdata(nb,nt) with cdata (recommended).
-CCPP_DATA_STRUCTURE = 'cdata'
diff --git a/ccpp/framework b/ccpp/framework
index 209f1c92d9..8365587133 160000
--- a/ccpp/framework
+++ b/ccpp/framework
@@ -1 +1 @@
-Subproject commit 209f1c92d99b7d4cc63e0d41c652fcfd730bd9fa
+Subproject commit 8365587133b40b1faa2831ac395bf098c3dac4a6
diff --git a/ccpp/physics b/ccpp/physics
index 4c17ff716a..0808cc2e89 160000
--- a/ccpp/physics
+++ b/ccpp/physics
@@ -1 +1 @@
-Subproject commit 4c17ff716a92d8ac0261d6cc2365bbd7a752b74a
+Subproject commit 0808cc2e8938ba66003b46746858143a9d75addb
diff --git a/gfsphysics/CCPP_layer/CCPP_data.meta b/gfsphysics/CCPP_layer/CCPP_data.meta
index e3f1249e5d..70c783820c 100644
--- a/gfsphysics/CCPP_layer/CCPP_data.meta
+++ b/gfsphysics/CCPP_layer/CCPP_data.meta
@@ -1,3 +1,8 @@
+[ccpp-table-properties]
+  name = CCPP_data
+  type = module
+  dependencies = ../../ccpp/framework/src/ccpp_types.F90,CCPP_typedefs.F90,../GFS_layer/GFS_typedefs.F90
+
 [ccpp-arg-table]
   name = CCPP_data
   type = module
diff --git a/gfsphysics/CCPP_layer/CCPP_typedefs.meta b/gfsphysics/CCPP_layer/CCPP_typedefs.meta
index d38fb86318..868dccebd1 100644
--- a/gfsphysics/CCPP_layer/CCPP_typedefs.meta
+++ b/gfsphysics/CCPP_layer/CCPP_typedefs.meta
@@ -1,3 +1,8 @@
+[ccpp-table-properties]
+  name = CCPP_interstitial_type
+  type = ddt
+  dependencies =
+
 [ccpp-arg-table]
   name = CCPP_interstitial_type
   type = ddt
@@ -334,6 +339,11 @@
   type = integer
 
 ########################################################################
+[ccpp-table-properties]
+  name = CCPP_typedefs
+  type = module
+  dependencies = ../../ccpp/physics/physics/machine.F
+
 [ccpp-arg-table]
   name = CCPP_typedefs
   type = module
diff --git a/gfsphysics/GFS_layer/GFS_typedefs.F90 b/gfsphysics/GFS_layer/GFS_typedefs.F90
index 4d3140d244..4291700597 100644
--- a/gfsphysics/GFS_layer/GFS_typedefs.F90
+++ b/gfsphysics/GFS_layer/GFS_typedefs.F90
@@ -7101,6 +7101,9 @@ subroutine interstitial_phys_reset (Interstitial, Model)
     !
   end subroutine interstitial_phys_reset
 
+  ! DH* 20200901: this routine is no longer used by CCPP's GFS_debug.F90. When new variables are
+  ! added to the GFS_interstitial_type, it is best to add the variable to both interstitial_print
+  ! below and to GFS_interstitialtoscreen in ccpp/physics/physics/GFS_debug.F90
   subroutine interstitial_print(Interstitial, Model, mpirank, omprank, blkno)
     !
     implicit none
diff --git a/gfsphysics/GFS_layer/GFS_typedefs.meta b/gfsphysics/GFS_layer/GFS_typedefs.meta
index dab6a5c173..84ba83e868 100644
--- a/gfsphysics/GFS_layer/GFS_typedefs.meta
+++ b/gfsphysics/GFS_layer/GFS_typedefs.meta
@@ -1,8 +1,8 @@
-[ccpp-arg-table]
-  name = GFS_init_type
+[ccpp-table-properties]
+  name = GFS_statein_type
   type = ddt
+  dependencies = 
 
-########################################################################
 [ccpp-arg-table]
   name = GFS_statein_type
   type = ddt
@@ -263,6 +263,11 @@
   kind = kind_phys
 
 ########################################################################
+[ccpp-table-properties]
+  name = GFS_stateout_type
+  type = ddt
+  dependencies = 
+
 [ccpp-arg-table]
   name = GFS_stateout_type
   type = ddt
@@ -439,6 +444,11 @@
   kind = kind_phys
 
 ########################################################################
+[ccpp-table-properties]
+  name = GFS_sfcprop_type
+  type = ddt
+  dependencies = 
+
 [ccpp-arg-table]
   name = GFS_sfcprop_type
   type = ddt
@@ -680,6 +690,7 @@
   dimensions = (horizontal_dimension)
   type = real
   kind = kind_phys
+  active = (flag_for_mass_flux_deep_convection_scheme == flag_for_gf_deep_convection_scheme)
 [hice]
   standard_name = sea_ice_thickness
   long_name = sea ice thickness
@@ -785,6 +796,7 @@
   dimensions = (horizontal_dimension)
   type = real
   kind = kind_phys
+  active = (flag_for_nsstm_run > 0)
 [z_c]
   standard_name = sub_layer_cooling_thickness
   long_name = sub-layer cooling thickness
@@ -792,6 +804,7 @@
   dimensions = (horizontal_dimension)
   type = real
   kind = kind_phys
+  active = (flag_for_nsstm_run > 0)
 [c_0]
   standard_name = coefficient_c_0
   long_name = coefficient 1 to calculate d(Tz)/d(Ts)
@@ -799,6 +812,7 @@
   dimensions = (horizontal_dimension)
   type = real
   kind = kind_phys
+  active = (flag_for_nsstm_run > 0)
 [c_d]
   standard_name = coefficient_c_d
   long_name = coefficient 2 to calculate d(Tz)/d(Ts)
@@ -806,6 +820,7 @@
   dimensions = (horizontal_dimension)
   type = real
   kind = kind_phys
+  active = (flag_for_nsstm_run > 0)
 [w_0]
   standard_name = coefficient_w_0
   long_name = coefficient 3 to calculate d(Tz)/d(Ts)
@@ -813,6 +828,7 @@
   dimensions = (horizontal_dimension)
   type = real
   kind = kind_phys
+  active = (flag_for_nsstm_run > 0)
 [w_d]
   standard_name = coefficient_w_d
   long_name = coefficient 4 to calculate d(Tz)/d(Ts)
@@ -820,6 +836,7 @@
   dimensions = (horizontal_dimension)
   type = real
   kind = kind_phys
+  active = (flag_for_nsstm_run > 0)
 [xt]
   standard_name = diurnal_thermocline_layer_heat_content
   long_name = heat content in diurnal thermocline layer
@@ -827,6 +844,7 @@
   dimensions = (horizontal_dimension)
   type = real
   kind = kind_phys
+  active = (flag_for_nsstm_run > 0)
 [xs]
   standard_name = sea_water_salinity
   long_name = salinity  content in diurnal thermocline layer
@@ -834,6 +852,7 @@
   dimensions = (horizontal_dimension)
   type = real
   kind = kind_phys
+  active = (flag_for_nsstm_run > 0)
 [xu]
   standard_name = diurnal_thermocline_layer_x_current
   long_name = u-current content in diurnal thermocline layer
@@ -841,6 +860,7 @@
   dimensions = (horizontal_dimension)
   type = real
   kind = kind_phys
+  active = (flag_for_nsstm_run > 0)
 [xv]
   standard_name = diurnal_thermocline_layer_y_current
   long_name = v-current content in diurnal thermocline layer
@@ -848,6 +868,7 @@
   dimensions = (horizontal_dimension)
   type = real
   kind = kind_phys
+  active = (flag_for_nsstm_run > 0)
 [xz]
   standard_name = diurnal_thermocline_layer_thickness
   long_name = diurnal thermocline layer thickness
@@ -855,6 +876,7 @@
   dimensions = (horizontal_dimension)
   type = real
   kind = kind_phys
+  active = (flag_for_nsstm_run > 0)
 [zm]
   standard_name = ocean_mixed_layer_thickness
   long_name = mixed layer thickness
@@ -862,6 +884,7 @@
   dimensions = (horizontal_dimension)
   type = real
   kind = kind_phys
+  active = (flag_for_nsstm_run > 0)
 [xtts]
   standard_name = sensitivity_of_dtl_heat_content_to_surface_temperature
   long_name = d(xt)/d(ts)
@@ -869,6 +892,7 @@
   dimensions = (horizontal_dimension)
   type = real
   kind = kind_phys
+  active = (flag_for_nsstm_run > 0)
 [xzts]
   standard_name = sensitivity_of_dtl_thickness_to_surface_temperature
   long_name = d(xz)/d(ts)
@@ -876,6 +900,7 @@
   dimensions = (horizontal_dimension)
   type = real
   kind = kind_phys
+  active = (flag_for_nsstm_run > 0)
 [d_conv]
   standard_name = free_convection_layer_thickness
   long_name = thickness of free convection layer (FCL)
@@ -883,6 +908,7 @@
   dimensions = (horizontal_dimension)
   type = real
   kind = kind_phys
+  active = (flag_for_nsstm_run > 0)
 [ifd]
   standard_name = index_of_dtlm_start
   long_name = index to start dtlm run or not
@@ -890,6 +916,7 @@
   dimensions = (horizontal_dimension)
   type = real
   kind = kind_phys
+  active = (flag_for_nsstm_run > 0)
 [dt_cool]
   standard_name = sub_layer_cooling_amount
   long_name = sub-layer cooling amount
@@ -897,6 +924,7 @@
   dimensions = (horizontal_dimension)
   type = real
   kind = kind_phys
+  active = (flag_for_nsstm_run > 0)
 [qrain]
   standard_name = sensible_heat_flux_due_to_rainfall
   long_name = sensible heat flux due to rainfall
@@ -904,6 +932,7 @@
   dimensions = (horizontal_dimension)
   type = real
   kind = kind_phys
+  active = (flag_for_nsstm_run > 0)
 [snowxy]
   standard_name = number_of_snow_layers
   long_name = number of snow layers
@@ -911,6 +940,7 @@
   dimensions = (horizontal_dimension)
   type = real
   kind = kind_phys
+  active = (flag_for_land_surface_scheme == flag_for_noahmp_land_surface_scheme)  
 [tvxy]
   standard_name = vegetation_temperature
   long_name = vegetation temperature
@@ -918,6 +948,7 @@
   dimensions = (horizontal_dimension)
   type = real
   kind = kind_phys
+  active = (flag_for_land_surface_scheme == flag_for_noahmp_land_surface_scheme)  
 [tgxy]
   standard_name = ground_temperature_for_noahmp
   long_name = ground temperature for noahmp
@@ -925,6 +956,7 @@
   dimensions = (horizontal_dimension)
   type = real
   kind = kind_phys
+  active = (flag_for_land_surface_scheme == flag_for_noahmp_land_surface_scheme)  
 [canicexy]
   standard_name = canopy_intercepted_ice_mass
   long_name = canopy intercepted ice mass
@@ -932,6 +964,7 @@
   dimensions = (horizontal_dimension)
   type = real
   kind = kind_phys
+  active = (flag_for_land_surface_scheme == flag_for_noahmp_land_surface_scheme)  
 [canliqxy]
   standard_name = canopy_intercepted_liquid_water
   long_name = canopy intercepted liquid water
@@ -939,6 +972,7 @@
   dimensions = (horizontal_dimension)
   type = real
   kind = kind_phys
+  active = (flag_for_land_surface_scheme == flag_for_noahmp_land_surface_scheme)  
 [eahxy]
   standard_name = canopy_air_vapor_pressure
   long_name = canopy air vapor pressure
@@ -946,6 +980,7 @@
   dimensions = (horizontal_dimension)
   type = real
   kind = kind_phys
+  active = (flag_for_land_surface_scheme == flag_for_noahmp_land_surface_scheme)  
 [tahxy]
   standard_name = canopy_air_temperature
   long_name = canopy air temperature
@@ -953,6 +988,7 @@
   dimensions = (horizontal_dimension)
   type = real
   kind = kind_phys
+  active = (flag_for_land_surface_scheme == flag_for_noahmp_land_surface_scheme)  
 [cmxy]
   standard_name = surface_drag_coefficient_for_momentum_for_noahmp
   long_name = surface drag coefficient for momentum for noahmp
@@ -960,6 +996,7 @@
   dimensions = (horizontal_dimension)
   type = real
   kind = kind_phys
+  active = (flag_for_land_surface_scheme == flag_for_noahmp_land_surface_scheme)  
 [chxy]
   standard_name = surface_drag_coefficient_for_heat_and_moisture_for_noahmp
   long_name = surface exchange coeff heat & moisture for noahmp
@@ -967,6 +1004,7 @@
   dimensions = (horizontal_dimension)
   type = real
   kind = kind_phys
+  active = (flag_for_land_surface_scheme == flag_for_noahmp_land_surface_scheme)  
 [fwetxy]
   standard_name = area_fraction_of_wet_canopy
   long_name = area fraction of canopy that is wetted/snowed
@@ -974,6 +1012,7 @@
   dimensions = (horizontal_dimension)
   type = real
   kind = kind_phys
+  active = (flag_for_land_surface_scheme == flag_for_noahmp_land_surface_scheme)  
 [sneqvoxy]
   standard_name = snow_mass_at_previous_time_step
   long_name = snow mass at previous time step
@@ -981,6 +1020,7 @@
   dimensions = (horizontal_dimension)
   type = real
   kind = kind_phys
+  active = (flag_for_land_surface_scheme == flag_for_noahmp_land_surface_scheme)  
 [alboldxy]
   standard_name = snow_albedo_at_previous_time_step
   long_name = snow albedo at previous time step
@@ -988,6 +1028,7 @@
   dimensions = (horizontal_dimension)
   type = real
   kind = kind_phys
+  active = (flag_for_land_surface_scheme == flag_for_noahmp_land_surface_scheme)  
 [qsnowxy]
   standard_name = snow_precipitation_rate_at_surface
   long_name = snow precipitation rate at surface
@@ -995,6 +1036,7 @@
   dimensions = (horizontal_dimension)
   type = real
   kind = kind_phys
+  active = (flag_for_land_surface_scheme == flag_for_noahmp_land_surface_scheme)  
 [wslakexy]
   standard_name = lake_water_storage
   long_name = lake water storage
@@ -1002,6 +1044,7 @@
   dimensions = (horizontal_dimension)
   type = real
   kind = kind_phys
+  active = (flag_for_land_surface_scheme == flag_for_noahmp_land_surface_scheme)  
 [zwtxy]
   standard_name = water_table_depth
   long_name = water table depth
@@ -1009,6 +1052,7 @@
   dimensions = (horizontal_dimension)
   type = real
   kind = kind_phys
+  active = (flag_for_land_surface_scheme == flag_for_noahmp_land_surface_scheme)  
 [waxy]
   standard_name = water_storage_in_aquifer
   long_name = water storage in aquifer
@@ -1016,6 +1060,7 @@
   dimensions = (horizontal_dimension)
   type = real
   kind = kind_phys
+  active = (flag_for_land_surface_scheme == flag_for_noahmp_land_surface_scheme)  
 [wtxy]
   standard_name = water_storage_in_aquifer_and_saturated_soil
   long_name = water storage in aquifer and saturated soil
@@ -1023,6 +1068,7 @@
   dimensions = (horizontal_dimension)
   type = real
   kind = kind_phys
+  active = (flag_for_land_surface_scheme == flag_for_noahmp_land_surface_scheme)  
 [tsnoxy]
   standard_name = snow_temperature
   long_name = snow_temperature
@@ -1030,6 +1076,7 @@
   dimensions = (horizontal_dimension, lower_bound_of_snow_vertical_dimension_for_land_surface_model:0)
   type = real
   kind = kind_phys
+  active = (flag_for_land_surface_scheme == flag_for_noahmp_land_surface_scheme)  
 [zsnsoxy]
   standard_name = layer_bottom_depth_from_snow_surface
   long_name = depth from the top of the snow surface at the bottom of the layer
@@ -1037,6 +1084,7 @@
   dimensions = (horizontal_dimension, lower_bound_of_snow_vertical_dimension_for_land_surface_model:soil_vertical_dimension_for_land_surface_model)
   type = real
   kind = kind_phys
+  active = (flag_for_land_surface_scheme == flag_for_noahmp_land_surface_scheme)  
 [snicexy]
   standard_name = snow_layer_ice
   long_name = snow layer ice
@@ -1044,6 +1092,7 @@
   dimensions = (horizontal_dimension, lower_bound_of_snow_vertical_dimension_for_land_surface_model:0)
   type = real
   kind = kind_phys
+  active = (flag_for_land_surface_scheme == flag_for_noahmp_land_surface_scheme)  
 [snliqxy]
   standard_name = snow_layer_liquid_water
   long_name = snow layer liquid water
@@ -1051,6 +1100,7 @@
   dimensions = (horizontal_dimension, lower_bound_of_snow_vertical_dimension_for_land_surface_model:0)
   type = real
   kind = kind_phys
+  active = (flag_for_land_surface_scheme == flag_for_noahmp_land_surface_scheme)  
 [lfmassxy]
   standard_name = leaf_mass
   long_name = leaf mass
@@ -1058,6 +1108,7 @@
   dimensions = (horizontal_dimension)
   type = real
   kind = kind_phys
+  active = (flag_for_land_surface_scheme == flag_for_noahmp_land_surface_scheme)  
 [rtmassxy]
   standard_name = fine_root_mass
   long_name = fine root mass
@@ -1065,6 +1116,7 @@
   dimensions = (horizontal_dimension)
   type = real
   kind = kind_phys
+  active = (flag_for_land_surface_scheme == flag_for_noahmp_land_surface_scheme)  
 [stmassxy]
   standard_name = stem_mass
   long_name = stem mass
@@ -1072,6 +1124,7 @@
   dimensions = (horizontal_dimension)
   type = real
   kind = kind_phys
+  active = (flag_for_land_surface_scheme == flag_for_noahmp_land_surface_scheme)  
 [woodxy]
   standard_name = wood_mass
   long_name = wood mass including woody roots
@@ -1079,6 +1132,7 @@
   dimensions = (horizontal_dimension)
   type = real
   kind = kind_phys
+  active = (flag_for_land_surface_scheme == flag_for_noahmp_land_surface_scheme)  
 [stblcpxy]
   standard_name = slow_soil_pool_mass_content_of_carbon
   long_name = stable carbon in deep soil
@@ -1086,6 +1140,7 @@
   dimensions = (horizontal_dimension)
   type = real
   kind = kind_phys
+  active = (flag_for_land_surface_scheme == flag_for_noahmp_land_surface_scheme)  
 [fastcpxy]
   standard_name = fast_soil_pool_mass_content_of_carbon
   long_name = short-lived carbon in shallow soil
@@ -1093,6 +1148,7 @@
   dimensions = (horizontal_dimension)
   type = real
   kind = kind_phys
+  active = (flag_for_land_surface_scheme == flag_for_noahmp_land_surface_scheme)  
 [xlaixy]
   standard_name = leaf_area_index
   long_name = leaf area index 
@@ -1100,6 +1156,7 @@
   dimensions = (horizontal_dimension)
   type = real
   kind = kind_phys
+  active = (flag_for_land_surface_scheme == flag_for_noahmp_land_surface_scheme .or. (flag_for_land_surface_scheme == flag_for_ruc_land_surface_scheme .and. flag_for_reading_leaf_area_index_from_input))
 [xsaixy]
   standard_name = stem_area_index
   long_name = stem area index 
@@ -1107,6 +1164,7 @@
   dimensions = (horizontal_dimension)
   type = real
   kind = kind_phys
+  active = (flag_for_land_surface_scheme == flag_for_noahmp_land_surface_scheme)  
 [taussxy]
   standard_name = nondimensional_snow_age
   long_name = non-dimensional snow age
@@ -1114,6 +1172,7 @@
   dimensions = (horizontal_dimension)
   type = real
   kind = kind_phys
+  active = (flag_for_land_surface_scheme == flag_for_noahmp_land_surface_scheme)  
 [smoiseq]
   standard_name = equilibrium_soil_water_content
   long_name = equilibrium soil water content
@@ -1121,6 +1180,7 @@
   dimensions = (horizontal_dimension,soil_vertical_dimension_for_land_surface_model)
   type = real
   kind = kind_phys
+  active = (flag_for_land_surface_scheme == flag_for_noahmp_land_surface_scheme)  
 [smcwtdxy]
   standard_name = soil_water_content_between_soil_bottom_and_water_table
   long_name = soil water content between the bottom of the soil and the water table
@@ -1128,6 +1188,7 @@
   dimensions = (horizontal_dimension)
   type = real
   kind = kind_phys
+  active = (flag_for_land_surface_scheme == flag_for_noahmp_land_surface_scheme)  
 [deeprechxy]
   standard_name = water_table_recharge_when_deep
   long_name = recharge to or from the water table when deep
@@ -1135,6 +1196,7 @@
   dimensions = (horizontal_dimension)
   type = real
   kind = kind_phys
+  active = (flag_for_land_surface_scheme == flag_for_noahmp_land_surface_scheme)  
 [rechxy]
   standard_name = water_table_recharge_when_shallow
   long_name = recharge to or from the water table when shallow
@@ -1142,6 +1204,7 @@
   dimensions = (horizontal_dimension)
   type = real
   kind = kind_phys
+  active = (flag_for_land_surface_scheme == flag_for_noahmp_land_surface_scheme)  
 [wetness]
   standard_name = normalized_soil_wetness_for_land_surface_model
   long_name = normalized soil wetness for lsm
@@ -1149,6 +1212,7 @@
   dimensions = (horizontal_dimension)
   type = real
   kind = kind_phys
+  active = (flag_for_land_surface_scheme == flag_for_ruc_land_surface_scheme) 
 [sh2o]
   standard_name = volume_fraction_of_unfrozen_soil_moisture_for_land_surface_model
   long_name = volume fraction of unfrozen soil moisture for lsm
@@ -1156,6 +1220,7 @@
   dimensions = (horizontal_dimension,soil_vertical_dimension_for_land_surface_model)
   type = real
   kind = kind_phys
+  active = (flag_for_land_surface_scheme == flag_for_ruc_land_surface_scheme) 
 [keepsmfr]
   standard_name = volume_fraction_of_frozen_soil_moisture_for_land_surface_model
   long_name = volume fraction of frozen soil moisture for lsm
@@ -1163,6 +1228,7 @@
   dimensions = (horizontal_dimension,soil_vertical_dimension_for_land_surface_model)
   type = real
   kind = kind_phys
+  active = (flag_for_land_surface_scheme == flag_for_ruc_land_surface_scheme) 
 [smois]
   standard_name = volume_fraction_of_soil_moisture_for_land_surface_model
   long_name = volumetric fraction of soil moisture for lsm
@@ -1170,6 +1236,7 @@
   dimensions = (horizontal_dimension,soil_vertical_dimension_for_land_surface_model)
   type = real
   kind = kind_phys
+  active = (flag_for_land_surface_scheme == flag_for_ruc_land_surface_scheme) 
 [tslb]
   standard_name = soil_temperature_for_land_surface_model
   long_name = soil temperature for land surface model
@@ -1177,6 +1244,7 @@
   dimensions = (horizontal_dimension,soil_vertical_dimension_for_land_surface_model)
   type = real
   kind = kind_phys
+  active = (flag_for_land_surface_scheme == flag_for_ruc_land_surface_scheme) 
 [zs]
   standard_name = depth_of_soil_levels_for_land_surface_model
   long_name = depth of soil levels for land surface model
@@ -1184,6 +1252,7 @@
   dimensions = (soil_vertical_dimension_for_land_surface_model)
   type = real
   kind = kind_phys
+  active = (flag_for_land_surface_scheme == flag_for_ruc_land_surface_scheme) 
 [clw_surf]
   standard_name = cloud_condensed_water_mixing_ratio_at_surface
   long_name = moist cloud water mixing ratio at surface
@@ -1191,6 +1260,7 @@
   dimensions = (horizontal_dimension)
   type = real
   kind = kind_phys
+  active = (flag_for_land_surface_scheme == flag_for_ruc_land_surface_scheme) 
 [qwv_surf]
   standard_name = water_vapor_mixing_ratio_at_surface
   long_name = water vapor mixing ratio at surface
@@ -1198,6 +1268,7 @@
   dimensions = (horizontal_dimension)
   type = real
   kind = kind_phys
+  active = (flag_for_land_surface_scheme == flag_for_ruc_land_surface_scheme) 
 [cndm_surf]
   standard_name = surface_condensation_mass
   long_name = surface condensation mass
@@ -1205,6 +1276,7 @@
   dimensions = (horizontal_dimension)
   type = real
   kind = kind_phys
+  active = (flag_for_land_surface_scheme == flag_for_ruc_land_surface_scheme) 
 [flag_frsoil]
   standard_name = flag_for_frozen_soil_physics
   long_name = flag for frozen soil physics (RUC)
@@ -1212,6 +1284,7 @@
   dimensions = (horizontal_dimension,soil_vertical_dimension_for_land_surface_model)
   type = real
   kind = kind_phys
+  active = (flag_for_land_surface_scheme == flag_for_ruc_land_surface_scheme) 
 [rhofr]
   standard_name = density_of_frozen_precipitation
   long_name = density of frozen precipitation
@@ -1219,6 +1292,7 @@
   dimensions = (horizontal_dimension)
   type = real
   kind = kind_phys
+  active = (flag_for_land_surface_scheme == flag_for_ruc_land_surface_scheme) 
 [tsnow]
   standard_name = snow_temperature_bottom_first_layer
   long_name = snow temperature at the bottom of the first snow layer
@@ -1226,6 +1300,7 @@
   dimensions = (horizontal_dimension)
   type = real
   kind = kind_phys
+  active = (flag_for_land_surface_scheme == flag_for_ruc_land_surface_scheme) 
 [snowfallac]
   standard_name = total_accumulated_snowfall
   long_name = run-total snow accumulation on the ground
@@ -1233,6 +1308,7 @@
   dimensions = (horizontal_dimension)
   type = real
   kind = kind_phys
+  active = (flag_for_land_surface_scheme == flag_for_ruc_land_surface_scheme) 
 [acsnow]
   standard_name = accumulated_water_equivalent_of_frozen_precip
   long_name = snow water equivalent of run-total frozen precip
@@ -1240,6 +1316,7 @@
   dimensions = (horizontal_dimension)
   type = real
   kind = kind_phys
+  active = (flag_for_land_surface_scheme == flag_for_ruc_land_surface_scheme) 
 [ustm]
   standard_name = surface_friction_velocity_drag
   long_name = friction velocity isolated for momentum only
@@ -1247,6 +1324,7 @@
   dimensions = (horizontal_dimension)
   type = real
   kind = kind_phys
+  active = (do_mynnsfclay)
 [zol]
   standard_name = surface_stability_parameter
   long_name = monin obukhov surface stability parameter
@@ -1254,6 +1332,7 @@
   dimensions = (horizontal_dimension)
   type = real
   kind = kind_phys
+  active = (do_mynnsfclay)
 [mol]
   standard_name = theta_star
   long_name = temperature flux divided by ustar (temperature scale)
@@ -1261,6 +1340,7 @@
   dimensions = (horizontal_dimension)
   type = real
   kind = kind_phys
+  active = (do_mynnsfclay)
 [rmol]
   standard_name = reciprocal_of_obukhov_length
   long_name = one over obukhov length
@@ -1268,6 +1348,7 @@
   dimensions = (horizontal_dimension)
   type = real
   kind = kind_phys
+  active = (do_mynnsfclay)
 [flhc]
   standard_name = surface_exchange_coefficient_for_heat
   long_name = surface exchange coefficient for heat
@@ -1275,6 +1356,7 @@
   dimensions = (horizontal_dimension)
   type = real
   kind = kind_phys
+  active = (do_mynnsfclay)
 [flqc]
   standard_name = surface_exchange_coefficient_for_moisture
   long_name = surface exchange coefficient for moisture
@@ -1282,6 +1364,7 @@
   dimensions = (horizontal_dimension)
   type = real
   kind = kind_phys
+  active = (do_mynnsfclay)
 [chs2]
   standard_name = surface_exchange_coefficient_for_heat_at_2m
   long_name = exchange coefficient for heat at 2 meters
@@ -1289,6 +1372,7 @@
   dimensions = (horizontal_dimension)
   type = real
   kind = kind_phys
+  active = (do_mynnsfclay)
 [cqs2]
   standard_name = surface_exchange_coefficient_for_moisture_at_2m
   long_name = exchange coefficient for moisture at 2 meters
@@ -1296,6 +1380,7 @@
   dimensions = (horizontal_dimension)
   type = real
   kind = kind_phys
+  active = (do_mynnsfclay)
 [lh]
   standard_name = surface_latent_heat
   long_name = latent heating at the surface (pos = up)
@@ -1303,6 +1388,7 @@
   dimensions = (horizontal_dimension)
   type = real
   kind = kind_phys
+  active = (do_mynnsfclay)
 [evap]
   standard_name = kinematic_surface_upward_latent_heat_flux
   long_name = kinematic surface upward latent heat flux
@@ -1331,6 +1417,7 @@
   dimensions = (horizontal_dimension)
   type = real
   kind = kind_phys
+  active = (flag_for_land_surface_scheme == flag_for_ruc_land_surface_scheme .or. flag_for_land_surface_scheme == flag_for_noahmp_land_surface_scheme)
 [rainncprv]
   standard_name = lwe_thickness_of_explicit_rainfall_amount_from_previous_timestep
   long_name = explicit rainfall from previous timestep
@@ -1338,6 +1425,7 @@
   dimensions = (horizontal_dimension)
   type = real
   kind = kind_phys
+  active = (flag_for_land_surface_scheme == flag_for_ruc_land_surface_scheme .or. flag_for_land_surface_scheme == flag_for_noahmp_land_surface_scheme)
 [iceprv]
   standard_name = lwe_thickness_of_ice_amount_from_previous_timestep
   long_name = ice amount from previous timestep
@@ -1345,6 +1433,7 @@
   dimensions = (horizontal_dimension)
   type = real
   kind = kind_phys
+  active = (flag_for_land_surface_scheme == flag_for_ruc_land_surface_scheme .or. flag_for_land_surface_scheme == flag_for_noahmp_land_surface_scheme)
 [snowprv]
   standard_name = lwe_thickness_of_snow_amount_from_previous_timestep
   long_name = snow amount from previous timestep
@@ -1352,6 +1441,7 @@
   dimensions = (horizontal_dimension)
   type = real
   kind = kind_phys
+  active = (flag_for_land_surface_scheme == flag_for_ruc_land_surface_scheme .or. flag_for_land_surface_scheme == flag_for_noahmp_land_surface_scheme)
 [graupelprv]
   standard_name = lwe_thickness_of_graupel_amount_from_previous_timestep
   long_name = graupel amount from previous timestep
@@ -1359,6 +1449,7 @@
   dimensions = (horizontal_dimension)
   type = real
   kind = kind_phys
+  active = (flag_for_land_surface_scheme == flag_for_ruc_land_surface_scheme .or. flag_for_land_surface_scheme == flag_for_noahmp_land_surface_scheme)
 [draincprv]
   standard_name = convective_precipitation_rate_from_previous_timestep
   long_name = convective precipitation rate from previous timestep
@@ -1366,6 +1457,7 @@
   dimensions = (horizontal_dimension)
   type = real
   kind = kind_phys
+  active = (flag_for_land_surface_scheme == flag_for_noahmp_land_surface_scheme)  
 [drainncprv]
   standard_name = explicit_rainfall_rate_from_previous_timestep
   long_name = explicit rainfall rate previous timestep
@@ -1373,6 +1465,7 @@
   dimensions = (horizontal_dimension)
   type = real
   kind = kind_phys
+  active = (flag_for_land_surface_scheme == flag_for_noahmp_land_surface_scheme)  
 [diceprv]
   standard_name = ice_precipitation_rate_from_previous_timestep
   long_name = ice precipitation rate from previous timestep
@@ -1380,6 +1473,7 @@
   dimensions = (horizontal_dimension)
   type = real
   kind = kind_phys
+  active = (flag_for_land_surface_scheme == flag_for_noahmp_land_surface_scheme)  
 [dsnowprv]
   standard_name = snow_precipitation_rate_from_previous_timestep
   long_name = snow precipitation rate from previous timestep
@@ -1387,6 +1481,7 @@
   dimensions = (horizontal_dimension)
   type = real
   kind = kind_phys
+  active = (flag_for_land_surface_scheme == flag_for_noahmp_land_surface_scheme)  
 [dgraupelprv]
   standard_name = graupel_precipitation_rate_from_previous_timestep
   long_name = graupel precipitation rate from previous timestep
@@ -1394,6 +1489,7 @@
   dimensions = (horizontal_dimension)
   type = real
   kind = kind_phys
+  active = (flag_for_land_surface_scheme == flag_for_noahmp_land_surface_scheme)  
 [alvsf]
   standard_name = mean_vis_albedo_with_strong_cosz_dependency
   long_name = mean vis albedo with strong cosz dependency
@@ -1424,6 +1520,11 @@
   kind = kind_phys
 
 ########################################################################
+[ccpp-table-properties]
+  name = GFS_coupling_type
+  type = ddt
+  dependencies = 
+
 [ccpp-arg-table]
   name = GFS_coupling_type
   type = ddt
@@ -1511,6 +1612,7 @@
   dimensions = (horizontal_dimension)
   type = real
   kind = kind_phys
+  active = (flag_for_flux_coupling .or. flag_for_stochastic_physics_perturbations .or. flag_for_chemistry_coupling .or. flag_for_global_cellular_automata)
 [rainc_cpl]
   standard_name = lwe_thickness_of_convective_precipitation_amount_for_coupling
   long_name = total convective precipitation
@@ -1525,6 +1627,7 @@
   dimensions = (horizontal_dimension)
   type = real
   kind = kind_phys
+  active = (flag_for_flux_coupling .or. flag_for_stochastic_physics_perturbations .or. flag_for_chemistry_coupling .or. flag_for_global_cellular_automata)
 [dusfc_cpl]
   standard_name = cumulative_surface_x_momentum_flux_for_coupling_multiplied_by_timestep
   long_name = cumulative sfc x momentum flux multiplied by timestep
@@ -1532,6 +1635,7 @@
   dimensions = (horizontal_dimension)
   type = real
   kind = kind_phys
+  active = (flag_for_flux_coupling)
 [dvsfc_cpl]
   standard_name = cumulative_surface_y_momentum_flux_for_coupling_multiplied_by_timestep
   long_name = cumulative sfc y momentum flux multiplied by timestep
@@ -1539,6 +1643,7 @@
   dimensions = (horizontal_dimension)
   type = real
   kind = kind_phys
+  active = (flag_for_flux_coupling)
 [dtsfc_cpl]
   standard_name = cumulative_surface_upward_sensible_heat_flux_for_coupling_multiplied_by_timestep
   long_name = cumulative sfc sensible heat flux multiplied by timestep
@@ -1546,6 +1651,7 @@
   dimensions = (horizontal_dimension)
   type = real
   kind = kind_phys
+  active = (flag_for_flux_coupling)
 [dqsfc_cpl]
   standard_name = cumulative_surface_upward_latent_heat_flux_for_coupling_multiplied_by_timestep
   long_name = cumulative sfc latent heat flux multiplied by timestep
@@ -1553,6 +1659,7 @@
   dimensions = (horizontal_dimension)
   type = real
   kind = kind_phys
+  active = (flag_for_flux_coupling)
 [dlwsfc_cpl]
   standard_name = cumulative_surface_downwelling_longwave_flux_for_coupling_multiplied_by_timestep
   long_name = cumulative sfc downward lw flux mulitplied by timestep
@@ -1560,6 +1667,7 @@
   dimensions = (horizontal_dimension)
   type = real
   kind = kind_phys
+  active = (flag_for_flux_coupling)
 [dswsfc_cpl]
   standard_name = cumulative_surface_downwelling_shortwave_flux_for_coupling_multiplied_by_timestep
   long_name = cumulative sfc downward sw flux multiplied by timestep
@@ -1567,6 +1675,7 @@
   dimensions = (horizontal_dimension)
   type = real
   kind = kind_phys
+  active = (flag_for_flux_coupling)
 [dnirbm_cpl]
   standard_name = cumulative_surface_downwelling_direct_near_infrared_shortwave_flux_for_coupling_multiplied_by_timestep
   long_name = cumulative sfc nir beam downward sw flux multiplied by timestep
@@ -1574,6 +1683,7 @@
   dimensions = (horizontal_dimension)
   type = real
   kind = kind_phys
+  active = (flag_for_flux_coupling)
 [dnirdf_cpl]
   standard_name = cumulative_surface_downwelling_diffuse_near_infrared_shortwave_flux_for_coupling_multiplied_by_timestep
   long_name = cumulative sfc nir diff downward sw flux multiplied by timestep
@@ -1581,6 +1691,7 @@
   dimensions = (horizontal_dimension)
   type = real
   kind = kind_phys
+  active = (flag_for_flux_coupling)
 [dvisbm_cpl]
   standard_name = cumulative_surface_downwelling_direct_ultraviolet_and_visible_shortwave_flux_for_coupling_multiplied_by_timestep
   long_name = cumulative sfc uv+vis beam dnwd sw flux multiplied by timestep
@@ -1588,6 +1699,7 @@
   dimensions = (horizontal_dimension)
   type = real
   kind = kind_phys
+  active = (flag_for_flux_coupling)
 [dvisdf_cpl]
   standard_name = cumulative_surface_downwelling_diffuse_ultraviolet_and_visible_shortwave_flux_for_coupling_multiplied_by_timestep
   long_name = cumulative sfc uv+vis diff dnwd sw flux multiplied by timestep
@@ -1595,6 +1707,7 @@
   dimensions = (horizontal_dimension)
   type = real
   kind = kind_phys
+  active = (flag_for_flux_coupling)
 [nlwsfc_cpl]
   standard_name = cumulative_surface_net_downward_longwave_flux_for_coupling_multiplied_by_timestep
   long_name = cumulative net downward lw flux multiplied by timestep
@@ -1602,6 +1715,7 @@
   dimensions = (horizontal_dimension)
   type = real
   kind = kind_phys
+  active = (flag_for_flux_coupling)
 [nswsfc_cpl]
   standard_name = cumulative_surface_net_downward_shortwave_flux_for_coupling_multiplied_by_timestep
   long_name = cumulative net downward sw flux multiplied by timestep
@@ -1609,6 +1723,7 @@
   dimensions = (horizontal_dimension)
   type = real
   kind = kind_phys
+  active = (flag_for_flux_coupling)
 [nnirbm_cpl]
   standard_name = cumulative_surface_net_downward_direct_near_infrared_shortwave_flux_for_coupling_multiplied_by_timestep
   long_name = cumulative net nir beam downward sw flux multiplied by timestep
@@ -1616,6 +1731,7 @@
   dimensions = (horizontal_dimension)
   type = real
   kind = kind_phys
+  active = (flag_for_flux_coupling)
 [nnirdf_cpl]
   standard_name = cumulative_surface_net_downward_diffuse_near_infrared_shortwave_flux_for_coupling_multiplied_by_timestep
   long_name = cumulative net nir diff downward sw flux multiplied by timestep
@@ -1623,6 +1739,7 @@
   dimensions = (horizontal_dimension)
   type = real
   kind = kind_phys
+  active = (flag_for_flux_coupling)
 [nvisbm_cpl]
   standard_name = cumulative_surface_net_downward_direct_ultraviolet_and_visible_shortwave_flux_for_coupling_multiplied_by_timestep
   long_name = cumulative net uv+vis beam downward sw rad flux multiplied by timestep
@@ -1630,6 +1747,7 @@
   dimensions = (horizontal_dimension)
   type = real
   kind = kind_phys
+  active = (flag_for_flux_coupling)
 [nvisdf_cpl]
   standard_name = cumulative_surface_net_downward_diffuse_ultraviolet_and_visible_shortwave_flux_for_coupling_multiplied_by_timestep
   long_name = cumulative net uv+vis diff downward sw rad flux multiplied by timestep
@@ -1637,6 +1755,7 @@
   dimensions = (horizontal_dimension)
   type = real
   kind = kind_phys
+  active = (flag_for_flux_coupling)
 [dusfci_cpl]
   standard_name = instantaneous_surface_x_momentum_flux_for_coupling
   long_name = instantaneous sfc x momentum flux
@@ -1644,6 +1763,7 @@
   dimensions = (horizontal_dimension)
   type = real
   kind = kind_phys
+  active = (flag_for_flux_coupling)
 [dvsfci_cpl]
   standard_name = instantaneous_surface_y_momentum_flux_for_coupling
   long_name = instantaneous sfc y momentum flux
@@ -1651,6 +1771,7 @@
   dimensions = (horizontal_dimension)
   type = real
   kind = kind_phys
+  active = (flag_for_flux_coupling)
 [dtsfci_cpl]
   standard_name = instantaneous_surface_upward_sensible_heat_flux_for_coupling
   long_name = instantaneous sfc sensible heat flux
@@ -1658,6 +1779,7 @@
   dimensions = (horizontal_dimension)
   type = real
   kind = kind_phys
+  active = (flag_for_flux_coupling)
 [dqsfci_cpl]
   standard_name = instantaneous_surface_upward_latent_heat_flux_for_coupling
   long_name = instantaneous sfc latent heat flux
@@ -1665,6 +1787,7 @@
   dimensions = (horizontal_dimension)
   type = real
   kind = kind_phys
+  active = (flag_for_flux_coupling)
 [dlwsfci_cpl]
   standard_name = instantaneous_surface_downwelling_longwave_flux_for_coupling
   long_name = instantaneous sfc downward lw flux
@@ -1672,6 +1795,7 @@
   dimensions = (horizontal_dimension)
   type = real
   kind = kind_phys
+  active = (flag_for_flux_coupling)
 [dswsfci_cpl]
   standard_name = instantaneous_surface_downwelling_shortwave_flux_for_coupling
   long_name = instantaneous sfc downward sw flux
@@ -1679,6 +1803,7 @@
   dimensions = (horizontal_dimension)
   type = real
   kind = kind_phys
+  active = (flag_for_flux_coupling)
 [dnirbmi_cpl]
   standard_name = instantaneous_surface_downwelling_direct_near_infrared_shortwave_flux_for_coupling
   long_name = instantaneous sfc nir beam downward sw flux
@@ -1686,6 +1811,7 @@
   dimensions = (horizontal_dimension)
   type = real
   kind = kind_phys
+  active = (flag_for_flux_coupling)
 [dnirdfi_cpl]
   standard_name = instantaneous_surface_downwelling_diffuse_near_infrared_shortwave_flux_for_coupling
   long_name = instantaneous sfc nir diff downward sw flux
@@ -1693,6 +1819,7 @@
   dimensions = (horizontal_dimension)
   type = real
   kind = kind_phys
+  active = (flag_for_flux_coupling)
 [dvisbmi_cpl]
   standard_name = instantaneous_surface_downwelling_direct_ultraviolet_and_visible_shortwave_flux_for_coupling
   long_name = instantaneous sfc uv+vis beam downward sw flux
@@ -1700,6 +1827,7 @@
   dimensions = (horizontal_dimension)
   type = real
   kind = kind_phys
+  active = (flag_for_flux_coupling)
 [dvisdfi_cpl]
   standard_name = instantaneous_surface_downwelling_diffuse_ultraviolet_and_visible_shortwave_flux_for_coupling
   long_name = instantaneous sfc uv+vis diff downward sw flux
@@ -1707,6 +1835,7 @@
   dimensions = (horizontal_dimension)
   type = real
   kind = kind_phys
+  active = (flag_for_flux_coupling)
 [nlwsfci_cpl]
   standard_name = instantaneous_surface_net_downward_longwave_flux_for_coupling
   long_name = instantaneous net sfc downward lw flux
@@ -1714,6 +1843,7 @@
   dimensions = (horizontal_dimension)
   type = real
   kind = kind_phys
+  active = (flag_for_flux_coupling)
 [nswsfci_cpl]
   standard_name = instantaneous_surface_net_downward_shortwave_flux_for_coupling
   long_name = instantaneous net sfc downward sw flux
@@ -1721,6 +1851,7 @@
   dimensions = (horizontal_dimension)
   type = real
   kind = kind_phys
+  active = (flag_for_flux_coupling)
 [nnirbmi_cpl]
   standard_name = instantaneous_surface_net_downward_direct_near_infrared_shortwave_flux_for_coupling
   long_name = instantaneous net nir beam sfc downward sw flux
@@ -1728,6 +1859,7 @@
   dimensions = (horizontal_dimension)
   type = real
   kind = kind_phys
+  active = (flag_for_flux_coupling)
 [nnirdfi_cpl]
   standard_name = instantaneous_surface_net_downward_diffuse_near_infrared_shortwave_flux_for_coupling
   long_name = instantaneous net nir diff sfc downward sw flux
@@ -1735,6 +1867,7 @@
   dimensions = (horizontal_dimension)
   type = real
   kind = kind_phys
+  active = (flag_for_flux_coupling)
 [nvisbmi_cpl]
   standard_name = instantaneous_surface_net_downward_direct_ultraviolet_and_visible_shortwave_flux_for_coupling
   long_name = instantaneous net uv+vis beam downward sw flux
@@ -1742,6 +1875,7 @@
   dimensions = (horizontal_dimension)
   type = real
   kind = kind_phys
+  active = (flag_for_flux_coupling)
 [nvisdfi_cpl]
   standard_name = instantaneous_surface_net_downward_diffuse_ultraviolet_and_visible_shortwave_flux_for_coupling
   long_name = instantaneous net uv+vis diff downward sw flux
@@ -1749,6 +1883,7 @@
   dimensions = (horizontal_dimension)
   type = real
   kind = kind_phys
+  active = (flag_for_flux_coupling)
 [t2mi_cpl]
   standard_name = instantaneous_temperature_at_2m_for_coupling
   long_name = instantaneous T2m
@@ -1756,6 +1891,7 @@
   dimensions = (horizontal_dimension)
   type = real
   kind = kind_phys
+  active = (flag_for_flux_coupling)
 [q2mi_cpl]
   standard_name = instantaneous_specific_humidity_at_2m_for_coupling
   long_name = instantaneous Q2m
@@ -1763,6 +1899,7 @@
   dimensions = (horizontal_dimension)
   type = real
   kind = kind_phys
+  active = (flag_for_flux_coupling)
 [u10mi_cpl]
   standard_name = instantaneous_x_wind_at_10m_for_coupling
   long_name = instantaneous U10m
@@ -1770,6 +1907,7 @@
   dimensions = (horizontal_dimension)
   type = real
   kind = kind_phys
+  active = (flag_for_flux_coupling .or. flag_for_wave_coupling)
 [v10mi_cpl]
   standard_name = instantaneous_y_wind_at_10m_for_coupling
   long_name = instantaneous V10m
@@ -1777,6 +1915,7 @@
   dimensions = (horizontal_dimension)
   type = real
   kind = kind_phys
+  active = (flag_for_flux_coupling .or. flag_for_wave_coupling)
 [tsfci_cpl]
   standard_name = instantaneous_surface_skin_temperature_for_coupling
   long_name = instantaneous sfc temperature
@@ -1784,6 +1923,7 @@
   dimensions = (horizontal_dimension)
   type = real
   kind = kind_phys
+  active = (flag_for_flux_coupling)
 [psurfi_cpl]
   standard_name = instantaneous_surface_air_pressure_for_coupling
   long_name = instantaneous sfc pressure
@@ -1791,6 +1931,7 @@
   dimensions = (horizontal_dimension)
   type = real
   kind = kind_phys
+  active = (flag_for_flux_coupling)
 [ulwsfcin_cpl]
   standard_name = surface_upwelling_longwave_flux_for_coupling
   long_name = surface upwelling LW flux for coupling
@@ -1798,6 +1939,7 @@
   dimensions = (horizontal_dimension)
   type = real
   kind = kind_phys
+  active = (flag_for_flux_coupling)
 [dusfcin_cpl]
   standard_name = surface_x_momentum_flux_for_coupling
   long_name = sfc x momentum flux for coupling
@@ -1805,6 +1947,7 @@
   dimensions = (horizontal_dimension)
   type = real
   kind = kind_phys
+  active = (flag_for_flux_coupling)
 [dvsfcin_cpl]
   standard_name = surface_y_momentum_flux_for_coupling
   long_name = sfc y momentum flux for coupling
@@ -1812,6 +1955,7 @@
   dimensions = (horizontal_dimension)
   type = real
   kind = kind_phys
+  active = (flag_for_flux_coupling)
 [dtsfcin_cpl]
   standard_name = surface_upward_sensible_heat_flux_for_coupling
   long_name = sfc sensible heat flux input
@@ -1819,6 +1963,7 @@
   dimensions = (horizontal_dimension)
   type = real
   kind = kind_phys
+  active = (flag_for_flux_coupling)
 [dqsfcin_cpl]
   standard_name = surface_upward_latent_heat_flux_for_coupling
   long_name = sfc latent heat flux input for coupling
@@ -1826,6 +1971,7 @@
   dimensions = (horizontal_dimension)
   type = real
   kind = kind_phys
+  active = (flag_for_flux_coupling)
 [hsnoin_cpl]
   standard_name = surface_snow_thickness_for_coupling
   long_name = sfc snow depth in meters over sea ice for coupling
@@ -1840,6 +1986,7 @@
   dimensions = (horizontal_dimension)
   type = real
   kind = kind_phys
+  active = (flag_for_flux_coupling)
 [ca_deep]
   standard_name = fraction_of_cellular_automata_for_deep_convection
   long_name = fraction of cellular automata for deep convection
@@ -1847,6 +1994,7 @@
   dimensions = (horizontal_dimension)
   type = real
   kind = kind_phys
+  active = (flag_for_cellular_automata)
 [vfact_ca]
   standard_name = vertical_weight_for_ca
   long_name = vertical weight for ca
@@ -1861,6 +2009,7 @@
   dimensions = (horizontal_dimension)
   type = real
   kind = kind_phys
+  active = (flag_for_cellular_automata)
 [condition]
   standard_name = physics_field_for_coupling
   long_name = physics_field_for_coupling
@@ -1875,6 +2024,7 @@
   dimensions = (horizontal_dimension,vertical_dimension)
   type = real
   kind = kind_phys
+  active = (flag_for_stochastic_shum_option)
 [sppt_wts]
   standard_name = weights_for_stochastic_sppt_perturbation
   long_name = weights for stochastic sppt perturbation
@@ -1882,6 +2032,7 @@
   dimensions = (horizontal_dimension,vertical_dimension)
   type = real
   kind = kind_phys
+  active = (flag_for_stochastic_physics_perturbations .or. flag_for_global_cellular_automata)
 [skebu_wts]
   standard_name = weights_for_stochastic_skeb_perturbation_of_x_wind
   long_name = weights for stochastic skeb perturbation of x wind
@@ -1889,6 +2040,7 @@
   dimensions = (horizontal_dimension,vertical_dimension)
   type = real
   kind = kind_phys
+  active = (flag_for_stochastic_skeb_option)
 [skebv_wts]
   standard_name = weights_for_stochastic_skeb_perturbation_of_y_wind
   long_name = weights for stochastic skeb perturbation of y wind
@@ -1896,6 +2048,7 @@
   dimensions = (horizontal_dimension,vertical_dimension)
   type = real
   kind = kind_phys
+  active = (flag_for_stochastic_skeb_option)
 [sfc_wts]
   standard_name = weights_for_stochastic_surface_physics_perturbation
   long_name = weights for stochastic surface physics perturbation
@@ -1903,6 +2056,7 @@
   dimensions = (horizontal_dimension,number_of_land_surface_variables_perturbed)
   type = real
   kind = kind_phys
+  active = (index_for_stochastic_land_surface_perturbation_type .ne. 0)
 [dqdti]
   standard_name = instantaneous_water_vapor_specific_humidity_tendency_due_to_convection
   long_name = instantaneous moisture tendency due to convection
@@ -1910,6 +2064,7 @@
   dimensions = (horizontal_dimension,vertical_dimension)
   type = real
   kind = kind_phys
+  active = (flag_for_chemistry_coupling)
 [nwfa2d]
   standard_name = tendency_of_water_friendly_aerosols_at_surface
   long_name = instantaneous water-friendly sfc aerosol source
@@ -1933,6 +2088,7 @@
   dimensions = (horizontal_dimension)
   type = real
   kind = kind_phys
+  active = (flag_for_chemistry_coupling)
 [dkt]
   standard_name = instantaneous_atmosphere_heat_diffusivity
   long_name = instantaneous atmospheric heat diffusivity
@@ -1940,6 +2096,7 @@
   dimensions = (horizontal_dimension,vertical_dimension)
   type = real
   kind = kind_phys
+  active = (flag_for_chemistry_coupling)
 [qci_conv]
   standard_name = convective_cloud_condesate_after_rainout
   long_name = convective cloud condesate after rainout
@@ -1947,7 +2104,13 @@
   dimensions = (horizontal_dimension,vertical_dimension)
   type = real
   kind = kind_phys
+  active = (flag_for_mass_flux_deep_convection_scheme == flag_for_gf_deep_convection_scheme)
 ########################################################################
+[ccpp-table-properties]
+  name = GFS_control_type
+  type = ddt
+  dependencies = 
+
 [ccpp-arg-table]
   name = GFS_control_type
   type = ddt
@@ -4329,6 +4492,11 @@
   type = logical
 
 ########################################################################
+[ccpp-table-properties]
+  name = GFS_grid_type
+  type = ddt
+  dependencies = 
+
 [ccpp-arg-table]
   name = GFS_grid_type
   type = ddt
@@ -4381,8 +4549,20 @@
   dimensions = (horizontal_dimension)
   type = real
   kind = kind_phys
+[xlon_d]
+  standard_name = longitude_in_degree
+  long_name = longitude in degree east
+  units = degree_east
+  dimensions = (horizontal_dimension)
+  type = real
+  kind = kind_phys
 
 ########################################################################
+[ccpp-table-properties]
+  name = GFS_tbd_type
+  type = ddt
+  dependencies = 
+
 [ccpp-arg-table]
   name = GFS_tbd_type
   type = ddt
@@ -4392,12 +4572,14 @@
   units = none
   dimensions = (horizontal_dimension)
   type = integer
+  active = (flag_for_lw_clouds_without_sub_grid_approximation == 2 .or. flag_for_sw_clouds_without_sub_grid_approximation == 2)
 [icsdlw]
   standard_name = seed_random_numbers_lw
   long_name = random seeds for sub-column cloud generators lw
   units = none
   dimensions = (horizontal_dimension)
   type = integer
+  active = (flag_for_lw_clouds_without_sub_grid_approximation == 2 .or. flag_for_sw_clouds_without_sub_grid_approximation == 2)
 [ozpl]
   standard_name = ozone_forcing
   long_name = ozone forcing data
@@ -4487,6 +4669,7 @@
   dimensions = (horizontal_dimension,vertical_dimension)
   type = real
   kind = kind_phys
+  active = (flag_for_stochastic_physics_perturbations .or. flag_for_global_cellular_automata)
 [drain_cpl]
   standard_name = tendency_of_lwe_thickness_of_precipitation_amount_for_coupling
   long_name = change in rain_cpl (coupling_type)
@@ -4494,6 +4677,7 @@
   dimensions = (horizontal_dimension)
   type = real
   kind = kind_phys
+  active = (flag_for_flux_coupling .or. flag_for_chemistry_coupling)
 [dsnow_cpl]
   standard_name = tendency_of_lwe_thickness_of_snow_amount_for_coupling
   long_name = change in show_cpl (coupling_type)
@@ -4501,6 +4685,7 @@
   dimensions = (horizontal_dimension)
   type = real
   kind = kind_phys
+  active = (flag_for_flux_coupling .or. flag_for_chemistry_coupling)
 [phy_fctd]
   standard_name = cloud_base_mass_flux
   long_name = cloud base mass flux for CS convection
@@ -4508,6 +4693,7 @@
   dimensions = (horizontal_dimension,number_of_cloud_types_CS)
   type = real
   kind = kind_phys
+  active = (number_of_cloud_types_CS > 0 .and. flag_for_Chikira_Sugiyama_deep_convection)
 [phy_f2d(:,1)]
   standard_name = surface_air_pressure_two_time_steps_back
   long_name = surface air pressure two time steps back
@@ -4641,6 +4827,7 @@
   dimensions = (horizontal_dimension,vertical_dimension)
   type = real
   kind = kind_phys
+  active = (flag_for_mass_flux_deep_convection_scheme == flag_for_gf_deep_convection_scheme .or. flag_for_mass_flux_deep_convection_scheme == flag_for_ntiedtke_deep_convection_scheme)
 [forceq]
   standard_name = moisture_tendency_due_to_dynamics
   long_name = moisture tendency due to dynamics only
@@ -4648,6 +4835,7 @@
   dimensions = (horizontal_dimension,vertical_dimension)
   type = real
   kind = kind_phys
+  active = (flag_for_mass_flux_deep_convection_scheme == flag_for_gf_deep_convection_scheme .or. flag_for_mass_flux_deep_convection_scheme == flag_for_ntiedtke_deep_convection_scheme)
 [prevst]
   standard_name = temperature_from_previous_timestep
   long_name = temperature from previous time step
@@ -4655,6 +4843,7 @@
   dimensions = (horizontal_dimension,vertical_dimension)
   type = real
   kind = kind_phys
+  active = (flag_for_mass_flux_deep_convection_scheme == flag_for_gf_deep_convection_scheme .or. flag_for_mass_flux_deep_convection_scheme == flag_for_ntiedtke_deep_convection_scheme)
 [prevsq]
   standard_name = moisture_from_previous_timestep
   long_name = moisture from previous time step
@@ -4662,12 +4851,14 @@
   dimensions = (horizontal_dimension,vertical_dimension)
   type = real
   kind = kind_phys
+  active = (flag_for_mass_flux_deep_convection_scheme == flag_for_gf_deep_convection_scheme .or. flag_for_mass_flux_deep_convection_scheme == flag_for_ntiedtke_deep_convection_scheme)
 [cactiv]
   standard_name = conv_activity_counter
   long_name = convective activity memory
   units = none
   dimensions = (horizontal_dimension)
   type = integer
+  active = (flag_for_mass_flux_deep_convection_scheme == flag_for_gf_deep_convection_scheme)
 [CLDFRA_BL]
   standard_name = subgrid_cloud_fraction_pbl
   long_name = subgrid cloud fraction from PBL scheme
@@ -4675,6 +4866,7 @@
   dimensions = (horizontal_dimension,vertical_dimension)
   type = real
   kind = kind_phys
+  active = (do_mynnedmf)
 [QC_BL]
   standard_name = subgrid_cloud_water_mixing_ratio_pbl
   long_name = subgrid cloud water mixing ratio from PBL scheme
@@ -4682,6 +4874,7 @@
   dimensions = (horizontal_dimension,vertical_dimension)
   type = real
   kind = kind_phys
+  active = (do_mynnedmf)
 [QI_BL]
   standard_name = subgrid_cloud_ice_mixing_ratio_pbl
   long_name = subgrid cloud ice mixing ratio from PBL scheme
@@ -4689,6 +4882,7 @@
   dimensions = (horizontal_dimension,vertical_dimension)
   type = real
   kind = kind_phys
+  active = (do_mynnedmf)
 [el_pbl]
   standard_name = mixing_length
   long_name = mixing length in meters
@@ -4696,6 +4890,7 @@
   dimensions = (horizontal_dimension,vertical_dimension)
   type = real
   kind = kind_phys
+  active = (do_mynnedmf)
 [Sh3D]
   standard_name = stability_function_for_heat
   long_name = stability function for heat
@@ -4703,6 +4898,7 @@
   dimensions = (horizontal_dimension,vertical_dimension)
   type = real
   kind = kind_phys
+  active = (do_mynnedmf)
 [qke]
   standard_name = tke_at_mass_points
   long_name = 2 x tke at mass points
@@ -4710,6 +4906,7 @@
   dimensions = (horizontal_dimension,vertical_dimension)
   type = real
   kind = kind_phys
+  active = (do_mynnedmf)
 [tsq]
   standard_name = t_prime_squared
   long_name = temperature fluctuation squared
@@ -4717,6 +4914,7 @@
   dimensions = (horizontal_dimension,vertical_dimension)
   type = real
   kind = kind_phys
+  active = (do_mynnedmf)
 [qsq]
   standard_name = q_prime_squared
   long_name = water vapor fluctuation squared
@@ -4724,6 +4922,7 @@
   dimensions = (horizontal_dimension,vertical_dimension)
   type = real
   kind = kind_phys
+  active = (do_mynnedmf)
 [cov]
   standard_name = t_prime_q_prime
   long_name = covariance of temperature and moisture
@@ -4731,6 +4930,7 @@
   dimensions = (horizontal_dimension,vertical_dimension)
   type = real
   kind = kind_phys
+  active = (do_mynnedmf)
 [phy_myj_qsfc]
   standard_name = surface_specific_humidity_for_MYJ_schemes
   long_name = surface air saturation specific humidity for MYJ schemes
@@ -4738,6 +4938,7 @@
   dimensions = (horizontal_dimension)
   type = real
   kind = kind_phys
+  active = (do_myjsfc .or. do_myjpbl)
 [phy_myj_thz0]
   standard_name = potential_temperature_at_viscous_sublayer_top
   long_name = potential temperature at viscous sublayer top over water
@@ -4745,6 +4946,7 @@
   dimensions = (horizontal_dimension)
   type = real
   kind = kind_phys
+  active = (do_myjsfc .or. do_myjpbl)
 [phy_myj_qz0]
   standard_name = specific_humidity_at_viscous_sublayer_top
   long_name = specific humidity at_viscous sublayer top over water
@@ -4752,6 +4954,7 @@
   dimensions = (horizontal_dimension)
   type = real
   kind = kind_phys
+  active = (do_myjsfc .or. do_myjpbl)
 [phy_myj_uz0]
   standard_name = u_wind_component_at_viscous_sublayer_top
   long_name = u wind component at viscous sublayer top over water
@@ -4759,6 +4962,7 @@
   dimensions = (horizontal_dimension)
   type = real
   kind = kind_phys
+  active = (do_myjsfc .or. do_myjpbl)
 [phy_myj_vz0]
   standard_name = v_wind_component_at_viscous_sublayer_top
   long_name = v wind component at viscous sublayer top over water
@@ -4766,6 +4970,7 @@
   dimensions = (horizontal_dimension)
   type = real
   kind = kind_phys
+  active = (do_myjsfc .or. do_myjpbl)
 [phy_myj_z0base]
   standard_name = baseline_surface_roughness_length
   long_name = baseline surface roughness length for momentum in meter
@@ -4773,6 +4978,7 @@
   dimensions = (horizontal_dimension)
   type = real
   kind = kind_phys
+  active = (do_myjsfc .or. do_myjpbl)
 [phy_myj_akhs]
   standard_name = heat_exchange_coefficient_for_MYJ_schemes
   long_name = surface heat exchange_coefficient for MYJ schemes
@@ -4780,6 +4986,7 @@
   dimensions = (horizontal_dimension)
   type = real
   kind = kind_phys
+  active = (do_myjsfc .or. do_myjpbl)
 [phy_myj_akms]
   standard_name = momentum_exchange_coefficient_for_MYJ_schemes
   long_name = surface momentum exchange_coefficient for MYJ schemes
@@ -4787,6 +4994,7 @@
   dimensions = (horizontal_dimension)
   type = real
   kind = kind_phys
+  active = (do_myjsfc .or. do_myjpbl)
 [phy_myj_chkqlm]
   standard_name = surface_layer_evaporation_switch
   long_name = surface layer evaporation switch
@@ -4794,6 +5002,7 @@
   dimensions = (horizontal_dimension)
   type = real
   kind = kind_phys
+  active = (do_myjsfc .or. do_myjpbl)
 [phy_myj_elflx]
   standard_name = kinematic_surface_latent_heat_flux
   long_name = kinematic surface latent heat flux
@@ -4801,6 +5010,7 @@
   dimensions = (horizontal_dimension)
   type = real
   kind = kind_phys
+  active = (do_myjsfc .or. do_myjpbl)
 [phy_myj_a1u]
   standard_name = weight_for_momentum_at_viscous_sublayer_top
   long_name = weight for momentum at viscous layer top
@@ -4808,6 +5018,7 @@
   dimensions = (horizontal_dimension)
   type = real
   kind = kind_phys
+  active = (do_myjsfc .or. do_myjpbl)
 [phy_myj_a1t]
   standard_name = weight_for_potental_temperature_at_viscous_sublayer_top
   long_name = weight for potental temperature at viscous layer top
@@ -4815,6 +5026,7 @@
   dimensions = (horizontal_dimension)
   type = real
   kind = kind_phys
+  active = (do_myjsfc .or. do_myjpbl)
 [phy_myj_a1q]
   standard_name = weight_for_specific_humidity_at_viscous_sublayer_top
   long_name = weight for Specfic Humidity at viscous layer top
@@ -4822,8 +5034,14 @@
   dimensions = (horizontal_dimension)
   type = real
   kind = kind_phys
+  active = (do_myjsfc .or. do_myjpbl)
 
 ########################################################################
+[ccpp-table-properties]
+  name = GFS_cldprop_type
+  type = ddt
+  dependencies = 
+
 [ccpp-arg-table]
   name = GFS_cldprop_type
   type = ddt
@@ -4850,6 +5068,11 @@
   kind = kind_phys
 
 ########################################################################
+[ccpp-table-properties]
+  name = GFS_radtend_type
+  type = ddt
+  dependencies = 
+
 [ccpp-arg-table]
   name = GFS_radtend_type
   type = ddt
@@ -4937,6 +5160,11 @@
   kind = kind_phys
 
 ########################################################################
+[ccpp-table-properties]
+  name = GFS_diag_type
+  type = ddt
+  dependencies = 
+
 [ccpp-arg-table]
   name = GFS_diag_type
   type = ddt
@@ -5057,6 +5285,7 @@
   dimensions = (horizontal_dimension)
   type = real 
   kind = kind_phys
+  active = (gwd_opt == 33)
 [dvsfc_ls]
   standard_name = integrated_y_momentum_flux_from_large_scale_gwd
   long_name = integrated y momentum flux from large scale gwd
@@ -5064,6 +5293,7 @@
   dimensions = (horizontal_dimension)
   type = real
   kind = kind_phys
+  active = (gwd_opt == 33)
 [dusfc_bl]
   standard_name = integrated_x_momentum_flux_from_blocking_drag
   long_name = integrated x momentum flux from blocking drag
@@ -5071,6 +5301,7 @@
   dimensions = (horizontal_dimension)
   type = real
   kind = kind_phys
+  active = (gwd_opt == 33)
 [dvsfc_bl]
   standard_name = integrated_y_momentum_flux_from_blocking_drag
   long_name = integrated y momentum flux from blocking drag
@@ -5078,6 +5309,7 @@
   dimensions = (horizontal_dimension)
   type = real
   kind = kind_phys
+  active = (gwd_opt == 33)
 [dusfc_ss]
   standard_name = integrated_x_momentum_flux_from_small_scale_gwd
   long_name = integrated x momentum flux from small scale gwd
@@ -5085,6 +5317,7 @@
   dimensions = (horizontal_dimension)
   type = real
   kind = kind_phys
+  active = (gwd_opt == 33)
 [dvsfc_ss]
   standard_name = integrated_y_momentum_flux_from_small_scale_gwd
   long_name = integrated y momentum flux from small scale gwd
@@ -5092,6 +5325,7 @@
   dimensions = (horizontal_dimension)
   type = real
   kind = kind_phys
+  active = (gwd_opt == 33)
 [dusfc_fd]
   standard_name = integrated_x_momentum_flux_from_form_drag
   long_name = integrated x momentum flux from form drag
@@ -5099,6 +5333,7 @@
   dimensions = (horizontal_dimension)
   type = real
   kind = kind_phys
+  active = (gwd_opt == 33)
 [dvsfc_fd]
   standard_name = integrated_y_momentum_flux_from_form_drag
   long_name = integrated y momentum flux from form drag
@@ -5106,6 +5341,7 @@
   dimensions = (horizontal_dimension)
   type = real
   kind = kind_phys
+  active = (gwd_opt == 33)
 [dtaux2d_ls]
   standard_name = x_momentum_tendency_from_large_scale_gwd
   long_name = x momentum tendency from large scale gwd
@@ -5113,6 +5349,7 @@
   dimensions = (horizontal_dimension,vertical_dimension)
   type = real
   kind = kind_phys
+  active = (gwd_opt == 33)
 [dtauy2d_ls]
   standard_name = y_momentum_tendency_from_large_scale_gwd
   long_name = y momentum tendency from large scale gwd
@@ -5120,6 +5357,7 @@
   dimensions = (horizontal_dimension,vertical_dimension)
   type = real
   kind = kind_phys
+  active = (gwd_opt == 33)
 [dtaux2d_bl]
   standard_name = x_momentum_tendency_from_blocking_drag
   long_name = x momentum tendency from blocking drag
@@ -5127,6 +5365,7 @@
   dimensions = (horizontal_dimension,vertical_dimension)
   type = real
   kind = kind_phys
+  active = (gwd_opt == 33)
 [dtauy2d_bl]
   standard_name = y_momentum_tendency_from_blocking_drag
   long_name = y momentum tendency from blocking drag
@@ -5134,6 +5373,7 @@
   dimensions = (horizontal_dimension,vertical_dimension)
   type = real
   kind = kind_phys
+  active = (gwd_opt == 33)
 [dtaux2d_ss]
   standard_name = x_momentum_tendency_from_small_scale_gwd
   long_name = x momentum tendency from small scale gwd
@@ -5141,12 +5381,14 @@
   dimensions = (horizontal_dimension,vertical_dimension)
   type = real
   kind = kind_phys
+  active = (gwd_opt == 33)
 [dtauy2d_ss]
   standard_name = y_momentum_tendency_from_small_scale_gwd
   long_name = y momentum tendency from small scale gwd
   units = m s-2
   dimensions = (horizontal_dimension,vertical_dimension)
   type = real
+  active = (gwd_opt == 33)
   kind = kind_phys
 [dtaux2d_fd]
   standard_name = x_momentum_tendency_from_form_drag
@@ -5155,6 +5397,7 @@
   dimensions = (horizontal_dimension,vertical_dimension)
   type = real
   kind = kind_phys
+  active = (gwd_opt == 33)
 [dtauy2d_fd]
   standard_name = y_momentum_tendency_from_form_drag
   long_name = y momentum tendency from form drag
@@ -5162,6 +5405,7 @@
   dimensions = (horizontal_dimension,vertical_dimension)
   type = real
   kind = kind_phys
+  active = (gwd_opt == 33)
 [totprcp]
   standard_name = accumulated_lwe_thickness_of_precipitation_amount
   long_name = accumulated total precipitation
@@ -5225,6 +5469,7 @@
   dimensions = (horizontal_dimension,vertical_dimension)
   type = real
   kind = kind_phys
+  active = (flag_for_microphysics_scheme == flag_for_fer_hires_microphysics_scheme)
 [gflux]
   standard_name = cumulative_surface_ground_heat_flux_multiplied_by_timestep
   long_name = cumulative groud conductive heat flux multiplied by timestep
@@ -5589,6 +5834,7 @@
   dimensions = (horizontal_dimension)
   type = real
   kind = kind_phys
+  active = (.not. flag_for_land_surface_scheme == flag_for_ruc_land_surface_scheme)
 [tdomr]
   standard_name = dominant_rain_type
   long_name = dominant rain type
@@ -5708,6 +5954,7 @@
   dimensions = (horizontal_dimension,vertical_dimension)
   type = real
   kind = kind_phys
+  active = (flag_diagnostics_3D)
 [dv3dt(:,:,1)]
   standard_name = cumulative_change_in_y_wind_due_to_PBL
   long_name = cumulative change in y wind due to PBL
@@ -5764,6 +6011,7 @@
   dimensions = (horizontal_dimension,vertical_dimension)
   type = real
   kind = kind_phys
+  active = (flag_diagnostics_3D)
 [dt3dt(:,:,1)]
   standard_name = cumulative_change_in_temperature_due_to_longwave_radiation
   long_name = cumulative change in temperature due to longwave radiation
@@ -5841,6 +6089,7 @@
   dimensions = (horizontal_dimension,vertical_dimension)
   type = real
   kind = kind_phys
+  active = (flag_diagnostics_3D)
 [dq3dt(:,:,1)]
   standard_name = cumulative_change_in_water_vapor_specific_humidity_due_to_PBL
   long_name = cumulative change in water vapor specific humidity due to PBL
@@ -5932,6 +6181,7 @@
   dimensions = (horizontal_dimension,vertical_dimension)
   type = real
   kind = kind_phys
+  active = (flag_diagnostics_3D .and. flag_tracer_diagnostics_3D)
 [refdmax]
   standard_name = maximum_reflectivity_at_1km_agl_over_maximum_hourly_time_interval
   long_name = maximum reflectivity at 1km agl over maximum hourly time interval
@@ -6034,6 +6284,7 @@
   dimensions = (horizonal_dimension,number_of_dust_bins_for_diagnostics)
   type = real
   kind = kind_phys
+  active = (number_of_dust_bins_for_diagnostics > 0)
 [ssem]
   standard_name = instantaneous_seasalt_emission_flux
   long_name = instantaneous sea salt emission flux
@@ -6041,6 +6292,7 @@
   dimensions = (horizonal_dimension,number_of_seasalt_bins_for_diagnostics)
   type = real
   kind = kind_phys
+  active = (number_of_seasalt_bins_for_diagnostics > 0)
 [sedim]
   standard_name = instantaneous_sedimentation
   long_name = instantaneous sedimentation
@@ -6048,6 +6300,7 @@
   dimensions = (horizonal_dimension,number_of_chemical_tracers_for_diagnostics)
   type = real
   kind = kind_phys
+  active = (number_of_chemical_tracers_for_diagnostics > 0)
 [drydep]
   standard_name = instantaneous_dry_deposition
   long_name = instantaneous dry deposition
@@ -6055,6 +6308,7 @@
   dimensions = (horizonal_dimension,number_of_chemical_tracers_for_diagnostics)
   type = real
   kind = kind_phys
+  active = (number_of_chemical_tracers_for_diagnostics > 0)
 [wetdpl]
   standard_name = instantaneous_large_scale_wet_deposition
   long_name = instantaneous large-scale wet deposition
@@ -6062,6 +6316,7 @@
   dimensions = (horizonal_dimension,number_of_chemical_tracers_for_diagnostics)
   type = real
   kind = kind_phys
+  active = (number_of_chemical_tracers_for_diagnostics > 0)
 [wetdpc]
   standard_name = instantaneous_convective_scale_wet_deposition
   long_name = instantaneous convective-scale wet deposition
@@ -6069,6 +6324,7 @@
   dimensions = (horizonal_dimension,number_of_chemical_tracers_for_diagnostics)
   type = real
   kind = kind_phys
+  active = (number_of_chemical_tracers_for_diagnostics > 0)
 [abem]
   standard_name = instantaneous_anthopogenic_and_biomass_burning_emissions
   long_name = instantaneous anthopogenic and biomass burning emissions for black carbon, organic carbon, and sulfur dioxide
@@ -6090,6 +6346,7 @@
   dimensions = (horizontal_dimension,vertical_dimension)
   type = real
   kind = kind_phys
+  active = (do_mynnedmf .and. (mynn_output_flag .ne. 0))
 [edmf_w]
   standard_name = emdf_updraft_vertical_velocity
   long_name = updraft vertical velocity from mass flux scheme
@@ -6097,6 +6354,7 @@
   dimensions = (horizontal_dimension,vertical_dimension)
   type = real
   kind = kind_phys
+  active = (do_mynnedmf .and. (mynn_output_flag .ne. 0))
 [edmf_qt]
   standard_name = emdf_updraft_total_water
   long_name = updraft total water from mass flux scheme
@@ -6104,6 +6362,7 @@
   dimensions = (horizontal_dimension,vertical_dimension)
   type = real
   kind = kind_phys
+  active = (do_mynnedmf .and. (mynn_output_flag .ne. 0))
 [edmf_thl]
   standard_name = emdf_updraft_theta_l
   long_name = updraft theta-l from mass flux scheme
@@ -6111,6 +6370,7 @@
   dimensions = (horizontal_dimension,vertical_dimension)
   type = real
   kind = kind_phys
+  active = (do_mynnedmf .and. (mynn_output_flag .ne. 0))
 [edmf_ent]
   standard_name = emdf_updraft_entrainment_rate
   long_name = updraft entranment rate from mass flux scheme
@@ -6118,6 +6378,7 @@
   dimensions = (horizontal_dimension,vertical_dimension)
   type = real
   kind = kind_phys
+  active = (do_mynnedmf .and. (mynn_output_flag .ne. 0))
 [edmf_qc]
   standard_name = emdf_updraft_cloud_water
   long_name = updraft cloud water from mass flux scheme
@@ -6125,6 +6386,7 @@
   dimensions = (horizontal_dimension,vertical_dimension)
   type = real
   kind = kind_phys
+  active = (do_mynnedmf .and. (mynn_output_flag .ne. 0))
 [sub_thl]
   standard_name = theta_subsidence_tendency
   long_name = updraft theta subsidence tendency
@@ -6132,6 +6394,7 @@
   dimensions = (horizontal_dimension,vertical_dimension)
   type = real
   kind = kind_phys
+  active = (do_mynnedmf .and. (mynn_output_flag .ne. 0))
 [sub_sqv]
   standard_name = water_vapor_subsidence_tendency
   long_name = updraft water vapor subsidence tendency
@@ -6139,6 +6402,7 @@
   dimensions = (horizontal_dimension,vertical_dimension)
   type = real
   kind = kind_phys
+  active = (do_mynnedmf .and. (mynn_output_flag .ne. 0))
 [det_thl]
   standard_name = theta_detrainment_tendency
   long_name = updraft theta detrainment tendency
@@ -6146,6 +6410,7 @@
   dimensions = (horizontal_dimension,vertical_dimension)
   type = real
   kind = kind_phys
+  active = (do_mynnedmf .and. (mynn_output_flag .ne. 0))
 [det_sqv]
   standard_name = water_vapor_detrainment_tendency
   long_name = updraft water vapor detrainment tendency
@@ -6153,12 +6418,14 @@
   dimensions = (horizontal_dimension,vertical_dimension)
   type = real
   kind = kind_phys
+  active = (do_mynnedmf .and. (mynn_output_flag .ne. 0))
 [nupdraft]
   standard_name = number_of_plumes
   long_name = number of plumes per grid column
   units = count
   dimensions = (horizontal_dimension)
   type = integer
+  active = (do_mynnedmf)
 [maxMF]
   standard_name = maximum_mass_flux
   long_name = maximum mass flux within a column
@@ -6166,6 +6433,7 @@
   dimensions = (horizontal_dimension)
   type = real
   kind = kind_phys
+  active = (do_mynnedmf)
 [ktop_shallow]
   standard_name = k_level_of_highest_reaching_plume
   long_name = k-level of highest reaching plume
@@ -6178,6 +6446,7 @@
   units = count
   dimensions = (horizontal_dimension)
   type = integer
+  active = (do_mynnedmf)
 [exch_h]
   standard_name = atmosphere_heat_diffusivity_for_mynnpbl
   long_name = diffusivity for heat for MYNN PBL (defined for all mass levels)
@@ -6185,6 +6454,7 @@
   dimensions = (horizontal_dimension,vertical_dimension)
   type = real
   kind = kind_phys
+  active = (do_mynnedmf)
 [exch_m]
   standard_name = atmosphere_momentum_diffusivity_for_mynnpbl
   long_name = diffusivity for momentum for MYNN PBL (defined for all mass levels)
@@ -6192,6 +6462,7 @@
   dimensions = (horizontal_dimension,vertical_dimension)
   type = real
   kind = kind_phys
+  active = (do_mynnedmf)
 [zmtb]
   standard_name = time_integral_of_height_of_mountain_blocking
   long_name = time integral of height of mountain blocking drag
@@ -6248,6 +6519,7 @@
   dimensions = (horizontal_dimension,vertical_dimension)
   type = real
   kind = kind_phys
+  active = (diag_ugwp_flag)
 [du3dt_ogw]
   standard_name = time_integral_of_change_in_x_wind_due_to_orographic_gravity_wave_drag
   long_name = time integral of change in x wind due to orographic gw drag
@@ -6255,6 +6527,7 @@
   dimensions = (horizontal_dimension,vertical_dimension)
   type = real
   kind = kind_phys
+  active = (diag_ugwp_flag)
 [du3dt_tms]
   standard_name = time_integral_of_change_in_x_wind_due_to_turbulent_orographic_form_drag
   long_name = time integral of change in x wind due to TOFD
@@ -6262,6 +6535,7 @@
   dimensions = (horizontal_dimension,vertical_dimension)
   type = real
   kind = kind_phys
+  active = (diag_ugwp_flag)
 [du3dt_ngw]
   standard_name = time_integral_of_change_in_x_wind_due_to_nonstationary_gravity_wave
   long_name = time integral of change in x wind due to NGW
@@ -6269,6 +6543,7 @@
   dimensions = (horizontal_dimension,vertical_dimension)
   type = real
   kind = kind_phys
+  active = (diag_ugwp_flag)
 [dv3dt_ngw]
   standard_name = time_integral_of_change_in_y_wind_due_to_nonstationary_gravity_wave
   long_name = time integral of change in y wind due to NGW
@@ -6276,6 +6551,7 @@
   dimensions = (horizontal_dimension,vertical_dimension)
   type = real
   kind = kind_phys
+  active = (diag_ugwp_flag)
 [aux2d]
   standard_name = auxiliary_2d_arrays
   long_name = auxiliary 2d arrays to output (for debugging)
@@ -6283,6 +6559,7 @@
   dimensions = (horizontal_dimension,number_of_3d_auxiliary_arrays)
   type = real
   kind = kind_phys
+  active = (number_of_2d_auxiliary_arrays > 0)
 [aux3d]
   standard_name = auxiliary_3d_arrays
   long_name = auxiliary 3d arrays to output (for debugging)
@@ -6290,9 +6567,15 @@
   dimensions = (horizontal_dimension,vertical_dimension,number_of_3d_auxiliary_arrays)
   type = real
   kind = kind_phys
+  active = (number_of_2d_auxiliary_arrays > 0)
 
 
 ########################################################################
+[ccpp-table-properties]
+  name = GFS_interstitial_type
+  type = ddt
+  dependencies = 
+
 [ccpp-arg-table]
   name = GFS_interstitial_type
   type = ddt
@@ -6303,6 +6586,7 @@
   dimensions = (horizontal_dimension,vertical_dimension)
   type = real
   kind = kind_phys
+  active = (flag_for_microphysics_scheme == flag_for_fer_hires_microphysics_scheme)
 [qc_r]
   standard_name = cloud_liquid_water_mixing_ratio
   long_name = the ratio of the mass of liquid water to the mass of dry air
@@ -6310,6 +6594,7 @@
   dimensions = (horizontal_dimension,vertical_dimension)
   type = real
   kind = kind_phys
+  active = (flag_for_microphysics_scheme == flag_for_fer_hires_microphysics_scheme)
 [qr_r]
   standard_name = cloud_rain_water_mixing_ratio
   long_name = the ratio of the mass rain water to the mass of dry air
@@ -6317,6 +6602,7 @@
   dimensions = (horizontal_dimension,vertical_dimension)
   type = real
   kind = kind_phys
+  active = (flag_for_microphysics_scheme == flag_for_fer_hires_microphysics_scheme)
 [qi_r]
   standard_name = cloud_ice_mixing_ratio
   long_name = the ratio of the mass of ice to the mass of dry air
@@ -6324,6 +6610,7 @@
   dimensions = (horizontal_dimension,vertical_dimension)
   type = real
   kind = kind_phys
+  active = (flag_for_microphysics_scheme == flag_for_fer_hires_microphysics_scheme)
 [qs_r]
   standard_name = cloud_snow_mixing_ratio
   long_name = the ratio of the mass of snow to mass of dry air
@@ -6331,6 +6618,7 @@
   dimensions = (horizontal_dimension,vertical_dimension)
   type = real
   kind = kind_phys
+  active = (flag_for_microphysics_scheme == flag_for_fer_hires_microphysics_scheme)
 [qg_r]
   standard_name = mass_weighted_rime_factor_mixing_ratio
   long_name = the ratio of the mass of rime factor to mass of dry air
@@ -6338,6 +6626,7 @@
   dimensions = (horizontal_dimension,vertical_dimension)
   type = real
   kind = kind_phys
+  active = (flag_for_microphysics_scheme == flag_for_fer_hires_microphysics_scheme)
 [f_ice]
   standard_name = fraction_of_ice_water_cloud
   long_name = fraction of ice water cloud
@@ -6345,6 +6634,7 @@
   dimensions = (horizontal_dimension,vertical_dimension)
   type = real
   kind = kind_phys
+  active = (flag_for_microphysics_scheme == flag_for_fer_hires_microphysics_scheme)
 [f_rain]
   standard_name = fraction_of_rain_water_cloud
   long_name = fraction of rain water cloud
@@ -6352,6 +6642,7 @@
   dimensions = (horizontal_dimension,vertical_dimension)
   type = real
   kind = kind_phys
+  active = (flag_for_microphysics_scheme == flag_for_fer_hires_microphysics_scheme)
 [f_rimef]
   standard_name = rime_factor
   long_name = rime factor
@@ -6359,6 +6650,7 @@
   dimensions = (horizontal_dimension,vertical_dimension)
   type = real
   kind = kind_phys
+  active = (flag_for_microphysics_scheme == flag_for_fer_hires_microphysics_scheme)
 [cwm]
   standard_name = total_cloud_condensate_mixing_ratio_updated_by_physics
   long_name = total cloud condensate mixing ratio (except water vapor) updated by physics
@@ -6366,6 +6658,7 @@
   dimensions = (horizontal_dimension,vertical_dimension)
   type = real
   kind = kind_phys
+  active = (flag_for_microphysics_scheme == flag_for_fer_hires_microphysics_scheme)
 [adjsfculw_ocean]
   standard_name = surface_upwelling_longwave_flux_over_ocean_interstitial
   long_name = surface upwelling longwave flux at current time over ocean (temporary use as interstitial)
@@ -6548,6 +6841,7 @@
   dimensions = (horizontal_dimension,vertical_dimension)
   type = real
   kind = kind_phys
+  active = (flag_for_microphysics_scheme == flag_for_morrison_gettelman_microphysics_scheme)
 [clcn]
   standard_name = convective_cloud_volume_fraction
   long_name = convective cloud volume fraction
@@ -6555,6 +6849,7 @@
   dimensions = (horizontal_dimension,vertical_dimension)
   type = real
   kind = kind_phys
+  active = (flag_for_microphysics_scheme == flag_for_morrison_gettelman_microphysics_scheme)
 [cldf]
   standard_name = cloud_area_fraction
   long_name = fraction of grid box area in which updrafts occur
@@ -6695,6 +6990,7 @@
   dimensions = (horizontal_dimension,4)
   type = real
   kind = kind_phys
+  active = (gwd_opt == 3 .or. gwd_opt == 33)
 [cmm_ocean]
   standard_name = surface_drag_wind_speed_for_momentum_in_air_over_ocean
   long_name = momentum exchange coefficient over ocean
@@ -6723,6 +7019,7 @@
   dimensions = (horizontal_dimension,vertical_dimension)
   type = real
   kind = kind_phys
+  active = (flag_for_microphysics_scheme == flag_for_morrison_gettelman_microphysics_scheme)
 [cnv_fice]
   standard_name = ice_fraction_in_convective_tower
   long_name = ice fraction in convective tower
@@ -6730,6 +7027,7 @@
   dimensions = (horizontal_dimension,vertical_dimension)
   type = real
   kind = kind_phys
+  active = (flag_for_microphysics_scheme == flag_for_morrison_gettelman_microphysics_scheme)
 [cnv_mfd]
   standard_name = detrained_mass_flux
   long_name = detrained mass flux
@@ -6737,6 +7035,7 @@
   dimensions = (horizontal_dimension,vertical_dimension)
   type = real
   kind = kind_phys
+  active = (flag_for_microphysics_scheme == flag_for_morrison_gettelman_microphysics_scheme)
 [cnv_ndrop]
   standard_name = number_concentration_of_cloud_liquid_water_particles_for_detrainment
   long_name = droplet number concentration in convective detrainment
@@ -6744,6 +7043,7 @@
   dimensions = (horizontal_dimension,vertical_dimension)
   type = real
   kind = kind_phys
+  active = (flag_for_microphysics_scheme == flag_for_morrison_gettelman_microphysics_scheme)
 [cnv_nice]
   standard_name = number_concentration_of_ice_crystals_for_detrainment
   long_name = crystal number concentration in convective detrainment
@@ -6751,6 +7051,7 @@
   dimensions = (horizontal_dimension,vertical_dimension)
   type = real
   kind = kind_phys
+  active = (flag_for_microphysics_scheme == flag_for_morrison_gettelman_microphysics_scheme)
 [cnvc]
   standard_name = convective_cloud_cover
   long_name = convective cloud cover
@@ -7486,6 +7787,7 @@
   dimensions = (horizontal_dimension)
   type = real
   kind = kind_phys
+  active = (flag_for_microphysics_scheme == flag_for_gfdl_microphysics_scheme .or. flag_for_microphysics_scheme == flag_for_thompson_microphysics_scheme)
 [gwdcu]
   standard_name = tendency_of_x_wind_due_to_convective_gravity_wave_drag
   long_name = zonal wind tendency due to convective gravity wave drag
@@ -7590,6 +7892,7 @@
   dimensions = (horizontal_dimension)
   type = real
   kind = kind_phys
+  active = (flag_for_microphysics_scheme == flag_for_gfdl_microphysics_scheme .or. flag_for_microphysics_scheme == flag_for_thompson_microphysics_scheme)
 [dry]
   standard_name = flag_nonzero_land_surface_fraction
   long_name = flag indicating presence of some land surface area fraction
@@ -7771,6 +8074,7 @@
   dimensions = (horizontal_dimension,vertical_dimension)
   type = real
   kind = kind_phys
+  active = (flag_for_microphysics_scheme == flag_for_morrison_gettelman_microphysics_scheme)
 [ncpi]
   standard_name = local_ice_number_concentration
   long_name = number concentration of ice local to physics
@@ -7778,6 +8082,7 @@
   dimensions = (horizontal_dimension,vertical_dimension)
   type = real
   kind = kind_phys
+  active = (flag_for_shoc)
 [ncpl]
   standard_name = local_condesed_water_number_concentration
   long_name = number concentration of condensed water local to physics
@@ -7785,6 +8090,7 @@
   dimensions = (horizontal_dimension,vertical_dimension)
   type = real
   kind = kind_phys
+  active = (flag_for_shoc)
 [ncpr]
   standard_name = local_rain_number_concentration
   long_name = number concentration of rain local to physics
@@ -7792,6 +8098,7 @@
   dimensions = (horizontal_dimension,vertical_dimension)
   type = real
   kind = kind_phys
+  active = (flag_for_microphysics_scheme == flag_for_morrison_gettelman_microphysics_scheme)
 [ncps]
   standard_name = local_snow_number_concentration
   long_name = number concentration of snow local to physics
@@ -7799,6 +8106,7 @@
   dimensions = (horizontal_dimension,vertical_dimension)
   type = real
   kind = kind_phys
+  active = (flag_for_microphysics_scheme == flag_for_morrison_gettelman_microphysics_scheme)
 [ncstrac]
   standard_name = number_of_tracers_for_CS
   long_name = number of convectively transported tracers in Chikira-Sugiyama deep convection scheme
@@ -7891,6 +8199,7 @@
   dimensions = (horizontal_dimension)
   type = real
   kind = kind_phys
+  active = (gwd_opt == 3 .or. gwd_opt == 33)
 [oa4ss]
   standard_name = asymmetry_of_subgrid_orography_small_scale
   long_name = asymmetry of subgrid orography small scale
@@ -7898,6 +8207,7 @@
   dimensions = (horizontal_dimension,4)
   type = real
   kind = kind_phys
+  active = (gwd_opt == 3 .or. gwd_opt == 33)
 [oc]
   standard_name = convexity_of_subgrid_orography
   long_name = convexity of subgrid orography
@@ -7912,6 +8222,7 @@
   dimensions = (horizontal_dimension)
   type = real
   kind = kind_phys
+  active = (gwd_opt == 3 .or. gwd_opt == 33)
 [olyr]
   standard_name = ozone_concentration_at_layer_for_radiation
   long_name = ozone concentration layer
@@ -7978,6 +8289,7 @@
   dimensions = (horizontal_dimension)
   type = real
   kind = kind_phys
+  active = (flag_for_land_surface_scheme == flag_for_noahmp_land_surface_scheme)
 [qgl]
   standard_name = local_graupel_mixing_ratio
   long_name = ratio of mass of graupel to mass of dry air plus vapor (without condensates) local to physics
@@ -7985,6 +8297,7 @@
   dimensions = (horizontal_dimension,vertical_dimension)
   type = real
   kind = kind_phys
+  active = (flag_for_microphysics_scheme == flag_for_morrison_gettelman_microphysics_scheme .or. flag_for_shoc)
 [qicn]
   standard_name = mass_fraction_of_convective_cloud_ice
   long_name = mass fraction of convective cloud ice water
@@ -7992,6 +8305,7 @@
   dimensions = (horizontal_dimension,vertical_dimension)
   type = real
   kind = kind_phys
+  active = (flag_for_microphysics_scheme == flag_for_morrison_gettelman_microphysics_scheme)
 [qlcn]
   standard_name = mass_fraction_of_convective_cloud_liquid_water
   long_name = mass fraction of convective cloud liquid water
@@ -7999,6 +8313,7 @@
   dimensions = (horizontal_dimension,vertical_dimension)
   type = real
   kind = kind_phys
+  active = (flag_for_microphysics_scheme == flag_for_morrison_gettelman_microphysics_scheme)
 [qlyr]
   standard_name = water_vapor_specific_humidity_at_layer_for_radiation
   long_name = specific humidity layer
@@ -8013,6 +8328,7 @@
   dimensions = (horizontal_dimension,vertical_dimension)
   type = real
   kind = kind_phys
+  active = (flag_for_microphysics_scheme == flag_for_morrison_gettelman_microphysics_scheme .or. flag_for_shoc)
 [qsnw]
   standard_name = local_snow_water_mixing_ratio
   long_name = ratio of mass of snow water to mass of dry air plus vapor (without condensates) local to physics
@@ -8020,6 +8336,7 @@
   dimensions = (horizontal_dimension,vertical_dimension)
   type = real
   kind = kind_phys
+  active = (flag_for_microphysics_scheme == flag_for_morrison_gettelman_microphysics_scheme .or. flag_for_shoc)
 [prcpmp]
   standard_name = lwe_thickness_of_explicit_precipitation_amount
   long_name = explicit precipitation (rain, ice, snow, graupel, ...) on physics timestep
@@ -8089,6 +8406,7 @@
   dimensions = (horizontal_dimension)
   type = real
   kind = kind_phys
+  active = (flag_for_microphysics_scheme == flag_for_gfdl_microphysics_scheme .or. flag_for_microphysics_scheme == flag_for_thompson_microphysics_scheme)
 [rainp]
   standard_name = tendency_of_rain_water_mixing_ratio_due_to_microphysics
   long_name = tendency of rain water mixing ratio due to microphysics
@@ -8330,6 +8648,7 @@
   dimensions = (horizontal_dimension)
   type = real
   kind = kind_phys
+  active = (flag_for_microphysics_scheme == flag_for_gfdl_microphysics_scheme .or. flag_for_microphysics_scheme == flag_for_thompson_microphysics_scheme)
 [snowmt]
   standard_name = surface_snow_melt
   long_name = snow melt during timestep
@@ -8378,6 +8697,7 @@
   dimensions = (horizontal_dimension)
   type = real
   kind = kind_phys
+  active = (flag_for_land_surface_scheme == flag_for_noahmp_land_surface_scheme)
 [theta]
   standard_name = angle_from_east_of_maximum_subgrid_orographic_variations
   long_name = angle with_respect to east of maximum subgrid orographic variations
@@ -8590,6 +8910,7 @@
   dimensions = (horizontal_dimension,vertical_dimension)
   type = real
   kind = kind_phys
+  active = (flag_for_microphysics_scheme == flag_for_morrison_gettelman_microphysics_scheme)
 [wcbmax]
   standard_name = maximum_updraft_velocity_at_cloud_base
   long_name = maximum updraft velocity at cloud base
@@ -8808,6 +9129,7 @@
   type = real
   kind = kind_phys
   optional = F
+  active = (flag_for_rrtmgp_radiation_scheme)
 [p_lev]
   standard_name = air_pressure_at_interface_for_RRTMGP_in_hPa
   long_name = air pressure level
@@ -8816,6 +9138,7 @@
   type = real
   kind = kind_phys
   optional = F
+  active = (flag_for_rrtmgp_radiation_scheme)
 [t_lay]
   standard_name = air_temperature_at_layer_for_RRTMGP
   long_name = air temperature layer
@@ -8824,6 +9147,7 @@
   type = real
   kind = kind_phys
   optional = F
+  active = (flag_for_rrtmgp_radiation_scheme)
 [t_lev]
   standard_name = air_temperature_at_interface_for_RRTMGP
   long_name = air temperature layer
@@ -8832,6 +9156,7 @@
   type = real
   kind = kind_phys
   optional = F
+  active = (flag_for_rrtmgp_radiation_scheme)
 [tv_lay]
   standard_name = virtual_temperature
   long_name = layer virtual temperature
@@ -8840,6 +9165,7 @@
   type = real
   kind = kind_phys
   optional = F
+  active = (flag_for_rrtmgp_radiation_scheme)
 [relhum]
   standard_name = relative_humidity
   long_name = layer relative humidity
@@ -8847,6 +9173,7 @@
   dimensions = (horizontal_dimension,vertical_dimension)
   type = real
   kind = kind_phys
+  active = (flag_for_rrtmgp_radiation_scheme)
 [deltaZ]
   standard_name = layer_thickness
   long_name = layer_thickness
@@ -8854,6 +9181,7 @@
   dimensions = (horizontal_dimension,vertical_dimension)
   type = real
   kind = kind_phys
+  active = (flag_for_rrtmgp_radiation_scheme)
 [tracer]
   standard_name = chemical_tracers
   long_name = chemical tracers
@@ -8861,6 +9189,7 @@
   dimensions = (horizontal_dimension,vertical_dimension,number_of_tracers)
   type = real
   kind = kind_phys
+  active = (flag_for_rrtmgp_radiation_scheme)
 [cloud_overlap_param]
   standard_name = cloud_overlap_param
   long_name = cloud overlap parameter
@@ -8868,6 +9197,7 @@
   dimensions = (horizontal_dimension,vertical_dimension)
   type = real
   kind = kind_phys
+  active = (flag_for_rrtmgp_radiation_scheme)
 [precip_overlap_param]
   standard_name = precip_overlap_param
   long_name = precipitation overlap parameter
@@ -8875,6 +9205,7 @@
   dimensions = (horizontal_dimension,vertical_dimension)
   type = real
   kind = kind_phys
+  active = (flag_for_rrtmgp_radiation_scheme)
 [ipsdsw0]
   standard_name = initial_permutation_seed_sw
   long_name = initial seed for McICA SW 
@@ -8910,6 +9241,7 @@
   type = real
   kind = kind_phys
   optional = F
+  active = (flag_for_rrtmgp_radiation_scheme)
 [fluxlwDOWN_allsky]
   standard_name = RRTMGP_lw_flux_profile_downward_allsky
   long_name = RRTMGP downward longwave all-sky flux profile
@@ -8918,6 +9250,7 @@
   type = real
   kind = kind_phys
   optional = F
+  active = (flag_for_rrtmgp_radiation_scheme)
 [fluxlwUP_clrsky]
   standard_name = RRTMGP_lw_flux_profile_upward_clrsky
   long_name = RRTMGP upward longwave clr-sky flux profile
@@ -8926,6 +9259,7 @@
   type = real
   kind = kind_phys
   optional = F
+  active = (flag_for_rrtmgp_radiation_scheme)
 [fluxlwDOWN_clrsky]
   standard_name = RRTMGP_lw_flux_profile_downward_clrsky
   long_name = RRTMGP downward longwave clr-sky flux profile
@@ -8933,6 +9267,7 @@
   dimensions = (horizontal_dimension,vertical_dimension_plus_one)
   type = real
   kind = kind_phys
+  active = (flag_for_rrtmgp_radiation_scheme)
 [sktp1r]
   standard_name = surface_skin_temperature_at_previous_time_step
   long_name = surface skin temperature at previous time step
@@ -8962,6 +9297,7 @@
   type = real
   kind = kind_phys
   optional = F
+  active = (flag_for_rrtmgp_radiation_scheme)
 [fluxswDOWN_allsky]
   standard_name = RRTMGP_sw_flux_profile_downward_allsky
   long_name = RRTMGP downward shortwave all-sky flux profile
@@ -8970,6 +9306,7 @@
   type = real
   kind = kind_phys
   optional = F
+  active = (flag_for_rrtmgp_radiation_scheme)
 [fluxswUP_clrsky]
   standard_name = RRTMGP_sw_flux_profile_upward_clrsky
   long_name = RRTMGP upward shortwave clr-sky flux profile
@@ -8978,6 +9315,7 @@
   type = real
   kind = kind_phys
   optional = F
+  active = (flag_for_rrtmgp_radiation_scheme)
 [fluxswDOWN_clrsky]
   standard_name = RRTMGP_sw_flux_profile_downward_clrsky
   long_name = RRTMGP downward shortwave clr-sky flux profile
@@ -8986,6 +9324,7 @@
   type = real
   kind = kind_phys
   optional = F
+  active = (flag_for_rrtmgp_radiation_scheme)
 [flxprf_lw]
   standard_name = RRTMGP_lw_fluxes
   long_name = lw fluxes total sky / csk and up / down at levels
@@ -8993,6 +9332,7 @@
   dimensions = (horizontal_dimension,vertical_dimension_plus_one)
   type = proflw_type
   optional = T
+  active = (flag_for_rrtmgp_radiation_scheme)
 [flxprf_sw]
   standard_name = RRTMGP_sw_fluxes
   long_name = sw fluxes total sky / csk and up / down at levels
@@ -9000,6 +9340,7 @@
   dimensions = (horizontal_dimension,vertical_dimension_plus_one)
   type = profsw_type
   optional = T
+  active = (flag_for_rrtmgp_radiation_scheme)
 [aerosolslw]
   standard_name = RRTMGP_aerosol_optical_properties_for_longwave_bands_01_16
   long_name = aerosol optical properties for longwave bands 01-16
@@ -9008,6 +9349,7 @@
   type = real
   kind = kind_phys
   optional = F
+  active = (flag_for_rrtmgp_radiation_scheme)
 [aerosolslw(:,:,:,1)]
   standard_name = RRTMGP_aerosol_optical_depth_for_longwave_bands_01_16
   long_name = aerosol optical depth for longwave bands 01-16
@@ -9036,6 +9378,7 @@
   dimensions = (horizontal_dimension,vertical_dimension, number_of_sw_bands_rrtmgp, number_of_aerosol_output_fields_for_shortwave_radiation)
   type = real
   kind = kind_phys
+  active = (flag_for_rrtmgp_radiation_scheme)
 [aerosolssw(:,:,:,1)]
   standard_name = RRTMGP_aerosol_optical_depth_for_shortwave_bands_01_16
   long_name = aerosol optical depth for shortwave bands 01-16
@@ -9063,12 +9406,14 @@
   units = none
   dimensions = (horizontal_dimension)
   type = integer
+  active = (flag_for_rrtmgp_radiation_scheme)
 [icseed_sw]
   standard_name = seed_random_numbers_sw_for_RRTMGP
   long_name = seed for random number generation for shortwave radiation
   units = none
   dimensions = (horizontal_dimension)
   type = integer
+  active = (flag_for_rrtmgp_radiation_scheme)
 [precip_frac]
   standard_name = precipitation_fraction_by_layer
   long_name = precipitation fraction in each layer
@@ -9076,6 +9421,7 @@
   dimensions = (horizontal_dimension,vertical_dimension)
   type = real
   kind = kind_phys
+  active = (flag_for_rrtmgp_radiation_scheme)
 [sw_gas_props]
   standard_name = coefficients_for_sw_gas_optics
   long_name = DDT containing spectral information for RRTMGP SW radiation scheme
@@ -9191,6 +9537,7 @@
   dimensions = (number_of_lw_bands_rrtmgp,horizontal_dimension)
   type = real
   kind = kind_phys
+  active = (flag_for_rrtmgp_radiation_scheme)
 [sec_diff_byband]
   standard_name = secant_of_diffusivity_angle_each_RRTMGP_LW_band
   long_name = secant of diffusivity angle in each RRTMGP LW band
@@ -9198,6 +9545,7 @@
   dimensions = (number_of_lw_bands_rrtmgp,horizontal_dimension)
   type = real
   kind = kind_phys
+  active = (flag_for_rrtmgp_radiation_scheme)
 [sfc_alb_nir_dir]
   standard_name = surface_albedo_nearIR_direct
   long_name = near-IR (direct) surface albedo (sfc_alb_nir_dir)
@@ -9205,6 +9553,7 @@
   dimensions = (number_of_sw_bands_rrtmgp,horizontal_dimension)
   type = real
   kind = kind_phys
+  active = (flag_for_rrtmgp_radiation_scheme)
 [sfc_alb_nir_dif]
   standard_name = surface_albedo_nearIR_diffuse
   long_name = near-IR (diffuse) surface albedo (sfc_alb_nir_dif)
@@ -9212,6 +9561,7 @@
   dimensions = (number_of_sw_bands_rrtmgp,horizontal_dimension)
   type = real
   kind = kind_phys
+  active = (flag_for_rrtmgp_radiation_scheme)
 [sfc_alb_uvvis_dir]
   standard_name =  surface_albedo_uvvis_dir
   long_name = UVVIS (direct) surface albedo (sfc_alb_uvvis_dir)
@@ -9219,6 +9569,7 @@
   dimensions = (number_of_sw_bands_rrtmgp,horizontal_dimension)
   type = real
   kind = kind_phys
+  active = (flag_for_rrtmgp_radiation_scheme)
 [sfc_alb_uvvis_dif]
   standard_name =  surface_albedo_uvvis_dif
   long_name = UVVIS (diffuse) surface albedo (sfc_alb_uvvis_dif)
@@ -9226,6 +9577,7 @@
   dimensions = (number_of_sw_bands_rrtmgp,horizontal_dimension)
   type = real
   kind = kind_phys
+  active = (flag_for_rrtmgp_radiation_scheme)
 [toa_src_lw]
   standard_name = toa_incident_lw_flux_by_spectral_point
   long_name = TOA longwave incident flux at each spectral points
@@ -9233,6 +9585,7 @@
   dimensions = (horizontal_dimension,number_of_lw_spectral_points_rrtmgp)
   type = real
   kind = kind_phys
+  active = (flag_for_rrtmgp_radiation_scheme)
 [toa_src_sw]
   standard_name = toa_incident_sw_flux_by_spectral_point
   long_name = TOA shortwave incident flux at each spectral points
@@ -9240,6 +9593,7 @@
   dimensions = (horizontal_dimension,number_of_sw_spectral_points_rrtmgp)
   type = real
   kind = kind_phys
+  active = (flag_for_rrtmgp_radiation_scheme)
 [active_gases_array]
   standard_name = list_of_active_gases_used_by_RRTMGP
   long_name = list of active gases used by RRTMGP
@@ -9247,8 +9601,14 @@
   dimensions =  (number_of_active_gases_used_by_RRTMGP)
   type = character
   kind = len=128
+  active = (flag_for_rrtmgp_radiation_scheme)
 
 ########################################################################
+[ccpp-table-properties]
+  name = GFS_data_type
+  type = ddt
+  dependencies = 
+
 [ccpp-arg-table]
   name = GFS_data_type
   type = ddt
@@ -9308,6 +9668,14 @@
   type = GFS_diag_type
 
 ########################################################################
+[ccpp-table-properties]
+  name = GFS_typedefs
+  type = module
+  relative_path = ../../ccpp/physics/physics
+  dependencies = machine.F,physcons.F90,radlw_param.f,radsw_param.f,GFDL_parse_tracers.F90
+  dependencies = rte-rrtmgp/rrtmgp/mo_gas_optics_rrtmgp.F90,rte-rrtmgp/rte/mo_optical_props.F90,rte-rrtmgp/extensions/cloud_optics/mo_cloud_optics.F90
+  dependencies = rte-rrtmgp/rrtmgp/mo_gas_concentrations.F90,rte-rrtmgp/rte/mo_rte_config.F90,rte-rrtmgp/rte/mo_source_functions.F90
+
 [ccpp-arg-table]
   name = GFS_typedefs
   type = module