Skip to content

Allow mkiau to handle case when analysis does not have delp #1488

New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

Sign up for GitHub

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Open
rtodling wants to merge 2 commits into
base: GEOS-FP-5_44
Choose a base branch
from
Open

Allow mkiau to handle case when analysis does not have delp #1488

Changes from all commits
Commits
Show all changes
2 commits
Select commit Hold shift + click to select a range
cc04a61
Allow mkiau to handle case when analysis does not have delp
rtodling Jul 24, 2026
a8f917c
Merge branch 'GEOS-FP-5_44' into feature/rtodling/mkiau_handle_absent…
rtodling Jul 24, 2026
File filter

Filter by extension

Filter by extension
Conversations
Failed to load comments.
Loading
Jump to
Jump to file
Failed to load files.
Loading
Diff view
Diff view
19 changes: 18 additions & 1 deletion GEOSmkiau_GridComp/GEOS_mkiauGridComp.F90
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -685,7 +685,9 @@ subroutine RUN ( GC, IMPORT, EXPORT, CLOCK, RC )
character(len=ESMF_MAXSTR) :: REPLAY_P, REPLAY_PS, REPLAY_DP
character(len=ESMF_MAXSTR) :: REPLAY_PHIS
character(len=ESMF_MAXSTR) :: REPLAY_T_TYPE
character(len=ESMF_MAXSTR) :: REPLAY_DP_FROM_PS
logical :: L_REPLAY_PHIS, L_REPLAY_P, L_REPLAY_U, L_REPLAY_V, L_REPLAY_T, L_REPLAY_QV, L_REPLAY_TS, L_REPLAY_O3
logical :: L_DP_FROM_PS
logical :: L_CUBE
logical :: do_transforms
logical :: USE_SPECFILT
Expand Down Expand Up @@ -857,6 +859,8 @@ subroutine RUN ( GC, IMPORT, EXPORT, CLOCK, RC )
VERIFY_(STATUS)
call MAPL_GetResource(MAPL, REPLAY_T_TYPE, Label="REPLAY_T_TYPE:", default='NULL', RC=STATUS)
VERIFY_(STATUS)
call MAPL_GetResource(MAPL, REPLAY_DP_FROM_PS, Label="REPLAY_DP_FROM_PS:", default='NULL', RC=STATUS)
VERIFY_(STATUS)

call MAPL_GetResource(MAPL, REPLAY_P_FACTOR, Label="REPLAY_P_FACTOR:", default=1.0 , RC=STATUS) ; VERIFY_(STATUS)
call MAPL_GetResource(MAPL, REPLAY_U_FACTOR, Label="REPLAY_U_FACTOR:", default=1.0 , RC=STATUS) ; VERIFY_(STATUS)
Expand Down Expand Up @@ -884,6 +888,7 @@ subroutine RUN ( GC, IMPORT, EXPORT, CLOCK, RC )
L_REPLAY_T = trim(REPLAY_T) .ne.'NO' .and. trim(REPLAY_T) .ne.'NULL'
L_REPLAY_QV = trim(REPLAY_QV) .ne.'NO' .and. trim(REPLAY_QV) .ne.'NULL'
L_REPLAY_O3 = trim(REPLAY_O3) .ne.'NO' .and. trim(REPLAY_O3) .ne.'NULL'
L_DP_FROM_PS = trim(REPLAY_DP_FROM_PS) .ne.'NO' .and. trim(REPLAY_DP_FROM_PS) .ne.'NULL'

call MAPL_GetResource(MAPL, DAMPBEG, LABEL="REPLAY_DAMPBEG:", default=1.0, RC=status)
VERIFY_(STATUS)
Expand Down Expand Up @@ -1789,7 +1794,9 @@ subroutine handleANA_
if( .not.FOUND .and. L_REPLAY_P ) then
write(STRING,'(A)') "ANA Variable: DP Not Found!"
call WRITE_PARALLEL( trim(STRING) )
RETURN_(ESMF_FAILURE)
if (.not.L_DP_FROM_PS) then
RETURN_(ESMF_FAILURE)
endif
endif

! ANA U-Wind
Expand Down Expand Up @@ -1938,6 +1945,16 @@ subroutine handleANA_
else
call set_eta ( LMana,ksdum,ptopdum,pintdum,ak_rep,bk_rep )
endif
! If so, create pressure thickness increment from surface pressure
! increment, regardless of whether thickness increment has been found or not
! --------------------------------------------------------------------------
if ( L_DP_FROM_PS ) then
write(STRING,'(A)') "ANA increment: DP derived from PS"
call WRITE_PARALLEL( trim(STRING) )
do L=1,LMana
dp_rep(:,:,L) = (ak_rep(L) - ak_rep(L-1)) + ps_ana(:,:)*(bk_rep(L) - bk_rep(L-1))
enddo
endif

ple_rep(:,:,0) = ak_rep(0)
do L=1,LMana
Expand Down
Loading