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/NCRTRP.f | 11 +++++++++-- 1 file changed, 9 insertions(+), 2 deletions(-) (limited to 'Donjon/src/NCRTRP.f') diff --git a/Donjon/src/NCRTRP.f b/Donjon/src/NCRTRP.f index 9b5203d..2298662 100644 --- a/Donjon/src/NCRTRP.f +++ b/Donjon/src/NCRTRP.f @@ -65,6 +65,13 @@ INTEGER, ALLOCATABLE, DIMENSION(:) :: JDEBAR,JARBVA REAL, ALLOCATABLE, DIMENSION(:) :: TERPA *---- +* TRIVAL CASE WHERE NCAL=1 +*---- + IF(NCAL.EQ.1) THEN + TERP(1)=1.0 + GO TO 110 + ENDIF +*---- * RECOVER TREE INFORMATION *---- JPCPO=LCMGID(IPCPO,'GLOBAL') @@ -205,11 +212,11 @@ TERP(ICAL)=TERP(ICAL)+TERTMP 100 CONTINUE ENDIF - IF(IMPX.GT.3) THEN + DEALLOCATE(JARBVA,JDEBAR) + 110 IF(IMPX.GT.3) THEN WRITE(IOUT,'(35H NCRTRP: TERP PARAMETERS IN MIXTURE,I4,1H:/(1X, 1 1P,10E12.4))') IBMOLD,(TERP(I),I=1,NCAL) ENDIF - DEALLOCATE(JARBVA,JDEBAR) RETURN *---- * MISSING ELEMENTARY CALCULATION EXCEPTION. -- cgit v1.2.3