From e03ca414a841c9739cbdcd38e1ce43d79695acc1 Mon Sep 17 00:00:00 2001 From: Helen Kershaw Date: Wed, 4 Oct 2023 12:03:40 -0600 Subject: [PATCH] fix: use lattidude of observation location for geopot to geometric height This is quite different from main. I'm not sure if this is a bug on my end or the order of caluculation makes a big difference. See https://github.com/NCAR/DART/issues/404#issuecomment-1747401428 for details --- models/wrf/model_mod.f90 | 2 +- 1 file changed, 1 insertion(+), 1 deletion(-) diff --git a/models/wrf/model_mod.f90 b/models/wrf/model_mod.f90 index 7fc9592aba..4d72585162 100644 --- a/models/wrf/model_mod.f90 +++ b/models/wrf/model_mod.f90 @@ -2005,7 +2005,7 @@ subroutine convert_vertical_obs(state_handle, num, locs, loc_qtys, loc_types, & zk = geopotential_height_interpolate(ens_size, state_handle, QTY_GEOPOTENTIAL_HEIGHT, id, ll, ul, lr, ur, k, dxm, dx, dy, dym) zk1 = geopotential_height_interpolate(ens_size, state_handle, QTY_GEOPOTENTIAL_HEIGHT, id, ll, ul, lr, ur, k+1, dxm, dx, dy, dym) geop = vertical_interpolation(ens_size, zloc, zk, zk1) - zout = compute_geometric_height(geop(1), grid(id)%latitude(i, j)) + zout = compute_geometric_height(geop(1), lon_lat_vert(2)) print*, 'HK geop, zout', geop, zout endif