Feel free to join the next Helmholtz Hacky Hour #30 on Wednesday, June 16, 2021 from 2PM to 3PM!

Commit 7f5f6b25 authored by Stephan Thober's avatar Stephan Thober
Browse files

BUG FIX: Nk_old was always initialized with zero in subroutine EDK in mo_EDK.f90

parent b85f3810
File deleted
......@@ -209,7 +209,7 @@ program ED_Kriging
Nk_old = cell(iCell)%Nk_old
!write(*,*),"Flag 1.5"
!write(*,*),"X before call EDK: ",X
! interploation
select case (interMth)
case (1)
......@@ -221,8 +221,9 @@ program ED_Kriging
!call EDK(iCell, jStartTmp, jEndTmp, dCS, MetSta, dS, cell, tempX, tempNkOld)
call EDK(iCell, jStartTmp, jEndTmp, dCS, MetSta, dS, cell, X, Nk_old)
end select
cell(iCell)%W = X
cell(iCell)%Nk_old = Nk_old
cell(iCell)%W = X
cell(iCell)%Nk_old = Nk_old
!write(*,*),"X after call EDK = ",X
end do ncellsloop
!$OMP end do
......@@ -309,7 +310,7 @@ program ED_Kriging
do iCell = 1, nCell
! initialize cell
deallocate(cell(iCell)%W)
if (allocated(cell(iCell)%W)) deallocate(cell(iCell)%W)
!cell(iCell)%z = noDataValue
end do
......
......@@ -55,7 +55,6 @@ contains
l = 0
ll = 0
Nk = 0
Nk_old = 0
! switch ordinary kriging off if not explicitly given
doOK_loc = .False.
......@@ -64,7 +63,7 @@ contains
! IF NK changed -> re-estimate weights
timeloop: do jd = jStart, jEnd
Nk_old = Nk
if (jd > jStart) Nk_old = Nk
Nk = 0_i4
l = 0
ll = 0
......
Markdown is supported
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment