From 21253b193e759e843947a509d490c45823fa380e Mon Sep 17 00:00:00 2001 From: Alain Hebert Date: Fri, 12 Dec 2025 21:11:16 +0100 Subject: #18: Implement the trival one-calculation case in subroutines *CRTRP.f --- Donjon/src/ACRTRP.f | 11 +++++++++-- 1 file changed, 9 insertions(+), 2 deletions(-) (limited to 'Donjon/src/ACRTRP.f') diff --git a/Donjon/src/ACRTRP.f b/Donjon/src/ACRTRP.f index 6a3e875..1091bc4 100644 --- a/Donjon/src/ACRTRP.f +++ b/Donjon/src/ACRTRP.f @@ -62,6 +62,13 @@ REAL, ALLOCATABLE, DIMENSION(:) :: TERPA,VREAL CHARACTER(LEN=80), ALLOCATABLE, DIMENSION(:) :: PARNAM *---- +* TRIVAL CASE WHERE NCAL=1 +*---- + IF(NCAL.EQ.1) THEN + TERP(1)=1.0 + GO TO 110 + ENDIF +*---- * RECOVER TREE INFORMATION *---- CALL hdf5_read_data(IPAPX,"/paramtree/DEBTREE",JDEBAR) @@ -188,11 +195,11 @@ TERP(ICAL)=TERP(ICAL)+TERTMP 100 CONTINUE ENDIF - IF(IMPX.GT.3) THEN + DEALLOCATE(JARBVA,JDEBAR,NVALUE) + 110 IF(IMPX.GT.3) THEN WRITE(IOUT,'(25H ACRTRP: TERP PARAMETERS:/(1X,1P,10E12.4))') 1 (TERP(I),I=1,NCAL) ENDIF - DEALLOCATE(JARBVA,JDEBAR,NVALUE) RETURN *---- * MISSING ELEMENTARY CALCULATION EXCEPTION. -- cgit v1.2.3