diff --git a/src/netcdf_io/calc_analysis.fd/inc2anl.f90 b/src/netcdf_io/calc_analysis.fd/inc2anl.f90
index e8d290eb..66bc2e90 100644
--- a/src/netcdf_io/calc_analysis.fd/inc2anl.f90
+++ b/src/netcdf_io/calc_analysis.fd/inc2anl.f90
@@ -347,16 +347,14 @@ subroutine add_psfc_increment
     real, allocatable, dimension(:,:,:) :: work3d_inc
     real, allocatable, dimension(:,:) :: ps_inc, work2d
     real, allocatable, dimension(:) :: bk5, work1d
-    integer :: iret, j, jj, k
+    integer :: iret, j, jj
     type(Dataset) :: incncfile
 
     ! get bk5 from attributes
     call read_attribute(fcstncfile, 'bk', bk5)
     ! read in delp increment to get ps increment
     incncfile = open_dataset(incr_file)
-    do k=1,nlev
-       call read_vardata(incncfile, 'delp_inc', work3d_inc, nslice=k, slicedim=3)
-    enddo
+    call read_vardata(incncfile, 'delp_inc', work3d_inc)
     ! get ps increment from delp increment and bk
     allocate(ps_inc(nlon,nlat))
     ps_inc(:,:) = work3d_inc(:,:,nlev) / (bk5(nlev) - bk5(nlev-1))