From 682d87450fa4b6a99b98b12592cebc1d60a4d22e Mon Sep 17 00:00:00 2001
From: Lucie Bakels <lbakels@jet01.img.univie.ac.at>
Date: Tue, 2 Jul 2024 10:56:37 +0200
Subject: [PATCH] Bugfix in nested fields: replace nint with int to not get
negative weights in interpolation
---
src/interpol_mod.f90 | 6 +++---
src/outgrid_mod.f90 | 4 ++--
2 files changed, 5 insertions(+), 5 deletions(-)
diff --git a/src/interpol_mod.f90 b/src/interpol_mod.f90
index 8cd29e3f..8620892b 100644
--- a/src/interpol_mod.f90
+++ b/src/interpol_mod.f90
@@ -136,8 +136,8 @@ subroutine find_grid_indices(xt,yt)
! jy=int(ytn)
! nix=nint(xtn)
! njy=nint(ytn)
- nix=max(min(nint(xtn),nxn(ngrid)-1),0)
- njy=max(min(nint(ytn),nyn(ngrid)-1),0)
+ nix=max(min(int(xtn),nxn(ngrid)-1),0)
+ njy=max(min(int(ytn),nyn(ngrid)-1),0)
ix=nix
jy=njy
ixp=ix+1
@@ -1987,4 +1987,4 @@ subroutine stdev_meter(iw)
end subroutine stdev_meter
#endif
-end module interpol_mod
\ No newline at end of file
+end module interpol_mod
diff --git a/src/outgrid_mod.f90 b/src/outgrid_mod.f90
index 97138884..42da3f16 100644
--- a/src/outgrid_mod.f90
+++ b/src/outgrid_mod.f90
@@ -273,8 +273,8 @@ subroutine outgrid_init
if (ngrid.gt.0) then
xtn=(xl-xln(ngrid))*xresoln(ngrid)
ytn=(yl-yln(ngrid))*yresoln(ngrid)
- ix=max(min(nint(xtn),nxn(ngrid)-1),0)
- jy=max(min(nint(ytn),nyn(ngrid)-1),0)
+ ix=max(min(int(xtn),nxn(ngrid)-1),0)
+ jy=max(min(int(ytn),nyn(ngrid)-1),0)
! ix=int(xtn)
! jy=int(ytn)
ddy=ytn-real(jy)
--
GitLab