Skip to content
Snippets Groups Projects
Commit 682d8745 authored by Lucie Bakels's avatar Lucie Bakels
Browse files

Bugfix in nested fields: replace nint with int to not get negative weights in interpolation

parent 4a136e22
No related branches found
No related tags found
No related merge requests found
Show whitespace changes
Inline Side-by-side
src/interpol_mod.f90
+ 3
− 3
View file @ 682d8745
...@@ -136,8 +136,8 @@ subroutine find_grid_indices(xt,yt) ...@@ -136,8 +136,8 @@ subroutine find_grid_indices(xt,yt)
! jy=int(ytn) ! jy=int(ytn)
! nix=nint(xtn) ! nix=nint(xtn)
! njy=nint(ytn) ! njy=nint(ytn)
nix=max(min(nint(xtn),nxn(ngrid)-1),0) nix=max(min(int(xtn),nxn(ngrid)-1),0)
njy=max(min(nint(ytn),nyn(ngrid)-1),0) njy=max(min(int(ytn),nyn(ngrid)-1),0)
ix=nix ix=nix
jy=njy jy=njy
ixp=ix+1 ixp=ix+1
......
src/outgrid_mod.f90
+ 2
− 2
View file @ 682d8745
...@@ -273,8 +273,8 @@ subroutine outgrid_init ...@@ -273,8 +273,8 @@ subroutine outgrid_init
if (ngrid.gt.0) then if (ngrid.gt.0) then
xtn=(xl-xln(ngrid))*xresoln(ngrid) xtn=(xl-xln(ngrid))*xresoln(ngrid)
ytn=(yl-yln(ngrid))*yresoln(ngrid) ytn=(yl-yln(ngrid))*yresoln(ngrid)
ix=max(min(nint(xtn),nxn(ngrid)-1),0) ix=max(min(int(xtn),nxn(ngrid)-1),0)
jy=max(min(nint(ytn),nyn(ngrid)-1),0) jy=max(min(int(ytn),nyn(ngrid)-1),0)
! ix=int(xtn) ! ix=int(xtn)
! jy=int(ytn) ! jy=int(ytn)
ddy=ytn-real(jy) ddy=ytn-real(jy)
......
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Please to comment