Skip to content

Commit

Permalink
Revert "Calc_analysis fix for delp increment read to allow C1152 cycl…
Browse files Browse the repository at this point in the history
…ing (#56…"

This reverts commit 7aed28e.
  • Loading branch information
CoryMartin-NOAA authored Dec 17, 2024
1 parent 7aed28e commit 0f5e4e8
Showing 1 changed file with 2 additions and 4 deletions.
6 changes: 2 additions & 4 deletions src/netcdf_io/calc_analysis.fd/inc2anl.f90
Original file line number Diff line number Diff line change
Expand Up @@ -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))
Expand Down

0 comments on commit 0f5e4e8

Please sign in to comment.