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/PCRTRP.f | 9 ++++++++- 1 file changed, 8 insertions(+), 1 deletion(-) (limited to 'Donjon/src/PCRTRP.f') diff --git a/Donjon/src/PCRTRP.f b/Donjon/src/PCRTRP.f index a19ba85..5e32eb7 100644 --- a/Donjon/src/PCRTRP.f +++ b/Donjon/src/PCRTRP.f @@ -61,6 +61,13 @@ *---- REAL, ALLOCATABLE, DIMENSION(:) :: TERPA *---- +* TRIVAL CASE WHERE NCAL=1 +*---- + IF(NCAL.EQ.1) THEN + TERP(1)=1.0 + GO TO 110 + ENDIF +*---- * COMPUTE TERP FACTORS *---- TERP(:NCAL)=0.0 @@ -172,7 +179,7 @@ TERP(ICAL)=TERP(ICAL)+TERTMP 100 CONTINUE ENDIF - IF(IMPX.GT.3) THEN + 110 IF(IMPX.GT.3) THEN WRITE(IOUT,'(25H PCRTRP: TERP PARAMETERS:/(1X,1P,10E12.4))') 1 (TERP(I),I=1,NCAL) ENDIF -- cgit v1.2.3