diff --git a/src/physics/cam/gw_drag.F90 b/src/physics/cam/gw_drag.F90
index ae7bbfca6f..cc39d3bad6 100644
--- a/src/physics/cam/gw_drag.F90
+++ b/src/physics/cam/gw_drag.F90
@@ -1063,6 +1063,15 @@ subroutine gw_init()
      call add_default('EKGW', 1, ' ')
   end if
 
+  call addfld ('RHOI', (/ 'ilev' /), 'A', 'kg/m3', &
+       'density at interfaces')
+
+  call addfld ('DSE', (/ 'lev' /), 'I', 'J/kg', &
+       'dry static energy')
+
+  call addfld ('NMBV', (/ 'lev' /), 'I', 'J/kg', &
+       'Brunt Vaisala Frequency')
+
   call addfld ('UTGW_TOTAL',    (/ 'lev' /), 'A','m/s2', &
        'Total U tendency due to gravity wave drag')
   call addfld ('VTGW_TOTAL',    (/ 'lev' /), 'A','m/s2', &
@@ -2096,6 +2105,10 @@ subroutine gw_tend(state, pbuf, dt, ptend, cam_in, flx_heat)
   call outfld('CLDLIQTGW', ptend%q(:,:,ixcldliq), pcols, lchnk)
   call outfld('CLDICETGW', ptend%q(:,:,ixcldice), pcols, lchnk)
 
+  call outfld('RHOI', rhoi, pcols, lchnk)
+  call outfld('DSE', dse, ncol, lchnk)
+  call outfld('NMBV', nm, ncol, lchnk)
+
   ! Destroy objects.
   call p%finalize()