Skip to content

Commit d1e5989

Browse files
fix: outgoing vertical coordinate should be 'height' in meters using zGridCenter
fixes #524 Following the call tree, using oncenters always convert_vertical_obs -> convert_vert_distrib case (VERTISHEIGHT) -> find_vert_indices -> find_vert_level(oncenters = true) case (VERTEIGHT) if (oncenters) -> find_height_bounds(zGridCenter)
1 parent e218864 commit d1e5989

File tree

1 file changed

+2
-2
lines changed

1 file changed

+2
-2
lines changed

models/mpas_atm/model_mod.f90

Lines changed: 2 additions & 2 deletions
Original file line numberDiff line numberDiff line change
@@ -4686,8 +4686,8 @@ subroutine convert_vert_distrib(state_handle, ens_size, location, obs_kind, ztyp
46864686
fdata = 0.0_r8
46874687
do i = 1, n
46884688
where (istatus == 0)
4689-
fdata(i, :) = zGridFace(k_low(i, :),c(i))*(1.0_r8 - fract(i, :)) + &
4690-
zGridFace(k_up (i, :),c(i))*fract(i, :)
4689+
fdata(i, :) = zGridCenter(k_low(i, :),c(i))*(1.0_r8 - fract(i, :)) + &
4690+
zGridCenter(k_up (i, :),c(i))*fract(i, :)
46914691
end where
46924692
enddo
46934693

0 commit comments

Comments
 (0)