From 7c5c3d305de0d6cc56c811b6f4ce6fea9d7fbb24 Mon Sep 17 00:00:00 2001 From: erickjomp Date: Sun, 23 Mar 2025 12:04:17 -0500 Subject: [PATCH] =?UTF-8?q?=E2=80=9CFInterpolation=20subroutine=20fixed?= MIME-Version: 1.0 Content-Type: text/plain; charset=UTF-8 Content-Transfer-Encoding: 8bit --- fortran/wrf_user.f90 | 2 +- 1 file changed, 1 insertion(+), 1 deletion(-) diff --git a/fortran/wrf_user.f90 b/fortran/wrf_user.f90 index 48963b7..a56c316 100644 --- a/fortran/wrf_user.f90 +++ b/fortran/wrf_user.f90 @@ -105,7 +105,7 @@ SUBROUTINE DINTERP3DZ(data3d, out2d, zdata, levels, nx, ny, nz, nlev, missingval desiredloc = levels(lev) DO WHILE ((.NOT. dointerp) .AND. (kp >= 2)) - IF (((zdata(i,j,kp-im) < desiredloc) .AND. (zdata(i,j,kp-ip) > desiredloc))) THEN + IF (((zdata(i,j,kp-im) < desiredloc) .AND. (zdata(i,j,kp-ip) >= desiredloc))) THEN w2 = (desiredloc - zdata(i,j,kp-im))/(zdata(i,j,kp-ip) - zdata(i,j,kp-im)) w1 = 1.D0 - w2 out2d(i,j,lev) = w1*data3d(i,j,kp-im) + w2*data3d(i,j,kp-ip)