diff --git a/src/oslo_aero_aerocom.F90 b/src/oslo_aero_aerocom.F90
index 0c825fd..d200614 100644
--- a/src/oslo_aero_aerocom.F90
+++ b/src/oslo_aero_aerocom.F90
@@ -268,6 +268,7 @@ subroutine aerocom1(lchnk, ncol, Cam, Nnatk, deltah_km, &
        end do
     enddo
     call outfld('ASYMMDRY',asydry_aer,pcols,lchnk)
+    call outfld('DELTAH',deltah_km,pcols,lchnk)
 
   end subroutine aerocom1
 
diff --git a/src/oslo_aero_diagnostics.F90 b/src/oslo_aero_diagnostics.F90
index 0d3f5e5..7402a2e 100644
--- a/src/oslo_aero_diagnostics.F90
+++ b/src/oslo_aero_diagnostics.F90
@@ -86,6 +86,8 @@ subroutine oslo_aero_diagnostics_init()
          call addfld ('PMTOT   ',horiz_only, 'A','ug/m3   ','Aerosol PM, all sizes')
          call addfld ('PM25    ',horiz_only, 'A','ug/m3   ','Aerosol PM2.5')
 
+         call addfld ('DELTAH  ',(/'lev'/), 'A', 'km      ','Layer thickness in km')
+
          call addfld ('PM2P5   ',(/'lev'/), 'A','ug/m3   ','3D aerosol PM2.5')
          call addfld ('MMRPM2P5',(/'lev'/), 'A','kg/kg   ','3D aerosol PM2.5 mass mixing ratio')
          call addfld ('MMRPM1  ',(/'lev'/), 'A','kg/kg   ','3D aerosol PM1.0 mass mixing ratio')
@@ -227,7 +229,6 @@ subroutine oslo_aero_diagnostics_init()
                call addfld(varName, horiz_only, 'A', 'unitless', 'relative exessive added mass column for mode'//trim(varname))
             end if
          enddo
-
          call add_default ('AKCXS   ', 1, ' ')
          call add_default ('PMTOT   ', 1, ' ')
          call add_default ('PM25    ', 1, ' ')
@@ -235,6 +236,7 @@ subroutine oslo_aero_diagnostics_init()
          call add_default ('MMRPM2P5', 1, ' ')
          call add_default ('MMRPM1  ', 1, ' ')
          call add_default ('GRIDAREA', 1, ' ')
+         call add_default ('DELTAH  ', 1, ' ')
          call add_default ('DAERH2O ', 1, ' ')
          call add_default ('MMR_AH2O', 1, ' ')
          call add_default ('ECDRYAER', 1, ' ')