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

Commit e6719183 authored by Stephan Thober's avatar Stephan Thober
Browse files

BUG FIX in deallocation of array C in mo_edk.f90 and updated edk.nml to make...

BUG FIX in deallocation of array C in mo_edk.f90 and updated edk.nml to make edk work out of the box
parent d8725316
......@@ -35,7 +35,7 @@ DataConvertFactor = 1.0e-1
! InterMth = 2 -> EDK
! InterMth = 1 -> OK
! InterMth = 0 -> No interpolation
InterMth = 0
InterMth = 2
!
! PROCESSING PERIOD
yStart = 1990
......@@ -49,7 +49,7 @@ dEnd = 31
maxDist = 1.5e5
!
! estimate variogram
flagVario = .True.
flagVario = .False.
!
! number of variogram parameters
nParam = 3
......@@ -57,10 +57,10 @@ nParam = 3
vType = 2
!
! file name where to store the variogram parameters
fNameVario = "var_param_de_pre.txt"
fNameVario = "check/var_param_de_pre.txt"
dh = 1.0e3 ! binsize for variogram [m]
hMax = 1.5e5 ! max distance h for variogram estimation [m]
! nc output specification
author_name = 'Friedrich Boeing'
author_name = 'XXX'
variable_name = 'pre'
/ !END*******Main Config***********
......@@ -214,8 +214,8 @@ contains
! call D_LSASF (A, B, X)
! print *, '<<<<<<<<<<<<<<'
! print *, 'rhs = ', B(:10)
! C = A
! X = B
C = A
X = B
! print *, 'A = ', C(:10, 1)
! print *, '<<<<<<<<<<<<<<'
call dgesv(size(A, 1), 1, A, size(A, 1), ipvt, B, size(A, 1), info)
......@@ -223,7 +223,7 @@ contains
! print *, 'B = ', B(:10)
! print *, 'A = ', A(:10, 1)
! print *, '<<<<<<<<<<<<<<'
if (maxval(abs(matmul(C, B) - X)) .gt. 1e-14) print *, 'result: ', maxval(abs(matmul(C, B) - X))
if (maxval(abs(matmul(C, B) - X)) .gt. 1e-10) print *, 'maximum error: ', maxval(abs(matmul(C, B) - X))
X = B
if (info .ne. 0_i4) then
print *, '***WARNING: calculation of weights failed'
......@@ -239,7 +239,7 @@ contains
! ! print *, 'ipvt: ', ipvt
! print *, 'info: ', info
! stop 'testing'
deallocate (A,B,ipvt)
deallocate (A, B, C, ipvt)
end subroutine get_kriging_weights
end module mo_EDK
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