From 3dc7a7f85fba765895c56e9c2d4c426c57924bcc Mon Sep 17 00:00:00 2001 From: RussTreadon-NOAA Date: Wed, 18 Dec 2024 00:43:41 +0000 Subject: [PATCH] replace work3d_inc with work2d_inc in add_psfc_increment (#56) --- src/netcdf_io/calc_analysis.fd/inc2anl.f90 | 9 ++++----- 1 file changed, 4 insertions(+), 5 deletions(-) diff --git a/src/netcdf_io/calc_analysis.fd/inc2anl.f90 b/src/netcdf_io/calc_analysis.fd/inc2anl.f90 index 66bc2e90..ffb0ccfa 100644 --- a/src/netcdf_io/calc_analysis.fd/inc2anl.f90 +++ b/src/netcdf_io/calc_analysis.fd/inc2anl.f90 @@ -344,8 +344,7 @@ subroutine add_psfc_increment use nemsio_module implicit none ! local variables - real, allocatable, dimension(:,:,:) :: work3d_inc - real, allocatable, dimension(:,:) :: ps_inc, work2d + real, allocatable, dimension(:,:) :: ps_inc, work2d, work2d_inc real, allocatable, dimension(:) :: bk5, work1d integer :: iret, j, jj type(Dataset) :: incncfile @@ -354,10 +353,10 @@ subroutine add_psfc_increment call read_attribute(fcstncfile, 'bk', bk5) ! read in delp increment to get ps increment incncfile = open_dataset(incr_file) - call read_vardata(incncfile, 'delp_inc', work3d_inc) + call read_vardata(incncfile, 'delp_inc', work2d_inc, nslice=nlev, slicedim=3) ! get ps increment from delp increment and bk allocate(ps_inc(nlon,nlat)) - ps_inc(:,:) = work3d_inc(:,:,nlev) / (bk5(nlev) - bk5(nlev-1)) + ps_inc(:,:) = work2d_inc(:,:) / (bk5(nlev) - bk5(nlev-1)) ! read in psfc background call read_vardata(fcstncfile, 'pressfc', work2d) ! add increment to background @@ -378,7 +377,7 @@ subroutine add_psfc_increment end if ! deallocate and close call close_dataset(incncfile) - deallocate(work2d,work3d_inc,ps_inc,bk5) + deallocate(work2d,work2d_inc,ps_inc,bk5) end subroutine add_psfc_increment