gkowal/amun-code
1
0
Fork 0
Code Issues Pull Requests Actions Packages Projects Releases Wiki Activity

INTERPOLATIONS: Return error in prepare_mgp() or prepare_gp() fail.

Browse Source 
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
This commit is contained in:
Grzegorz Kowal 2018-08-27 19:15:21 -03:00
parent b8b9f6feac
commit 32822544b0
1 changed files with 18 additions and 7 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

25
src/interpolations.F90
View File

@ -319,7 +319,8 @@ module interpolations
! prepare matrix coefficients
!
call prepare_gp()
call prepare_gp(iret)
if (iret > 0) return
interfaces => interfaces_dir
reconstruct_states => reconstruct_gp
@ -341,7 +342,8 @@ module interpolations
! prepare matrix coefficients
!
call prepare_mgp()
call prepare_mgp(iret)
if (iret > 0) return
interfaces => interfaces_mgp
if (verbose .and. 2 * ng <= ngp - 1) &
@ -355,7 +357,8 @@ module interpolations
write(error_unit,"('[',a,']: ',a)") trim(loc) &
, "The selected reconstruction method is not " // &
"implemented: " // trim(sreconstruction)
stop
iret = 200
return
end if
end select
@ -530,7 +533,7 @@ module interpolations
!
!===============================================================================
!
subroutine prepare_mgp()
subroutine prepare_mgp(iret)
! include external procedures
!
@ -542,6 +545,10 @@ module interpolations
!
implicit none
! subroutine arguments
!
integer, intent(inout) :: iret
! local variables
!
logical :: flag
@ -674,7 +681,7 @@ module interpolations
if (.not. flag) then
write(error_unit,"('[',a,']: ',a)") trim(loc) &
, "Could not invert the covariance matrix!"
stop
iret = 201
end if
!-------------------------------------------------------------------------------
@ -4406,7 +4413,7 @@ module interpolations
!
!===============================================================================
!
subroutine prepare_gp()
subroutine prepare_gp(iret)
! include external procedures
!
@ -4418,6 +4425,10 @@ module interpolations
!
implicit none
! subroutine arguments
!
integer, intent(inout) :: iret
! local variables
!
logical :: flag
@ -4506,7 +4517,7 @@ module interpolations
if (.not. flag) then
write(error_unit,"('[',a,']: ',a)") trim(loc) &
, "Could not invert the covariance matrix!"
stop
iret = 202
end if
!-------------------------------------------------------------------------------