diff --git a/parm/diag_table_aqm.FV3_GFS_v16 b/parm/diag_table_aqm.FV3_GFS_v16
index 1513e8ea74..b7bf4bfbd0 100644
--- a/parm/diag_table_aqm.FV3_GFS_v16
+++ b/parm/diag_table_aqm.FV3_GFS_v16
@@ -235,6 +235,10 @@
 "gfs_sfc",     "xzts",        "xzts",      "fv3_history2d",  "all",  .false.,  "none",  2
 "gfs_sfc",     "d_conv",      "dconv",     "fv3_history2d",  "all",  .false.,  "none",  2
 "gfs_sfc",     "qrain",       "qrain",     "fv3_history2d",  "all",  .false.,  "none",  2
+"gfs_sfc",     "canopylaixy",   "canopylaixy",   "fv3_history2d", "all",  .false.,  "none",  2
+"gfs_sfc",     "canopyfchxy",   "canopyfchxy",   "fv3_history2d", "all",  .false.,  "none",  2
+"gfs_sfc",     "canopyffracxy", "canopyffracxy", "fv3_history2d", "all",  .false.,  "none",  2
+"gfs_sfc",     "canopycluxy",   "canopycluxy",   "fv3_history2d", "all",  .false.,  "none",  2
 
 "gfs_phys",    "acond",        "acond",         "fv3_history2d",  "all",  .false.,  "none",  2
 "gfs_phys",    "cduvb_ave",    "cduvb_ave",     "fv3_history2d",  "all",  .false.,  "none",  2
diff --git a/ush/generate_FV3LAM_wflow.py b/ush/generate_FV3LAM_wflow.py
index f6311d7a9c..1a77e72d02 100755
--- a/ush/generate_FV3LAM_wflow.py
+++ b/ush/generate_FV3LAM_wflow.py
@@ -432,7 +432,11 @@ def generate_FV3LAM_wflow(
             if ( CCPP_PHYS_SUITE == "FV3_GFS_v16" or
                  CCPP_PHYS_SUITE == "FV3_GFS_v17_p8"):
                 gfs_physics_nml_dict.update({
-                "do_canopy": True         
+                "do_canopy": True,
+                "rdcanopylai": True,
+                "rdcanopyfch": True,
+                "rdcanopyffrac": True,
+                "rdcanopyclu": True
                 }) 
     settings["gfs_physics_nml"] = gfs_physics_nml_dict