From b7b82bcd38a9c596d3a57ecd41b9a94c831ae999 Mon Sep 17 00:00:00 2001 From: mjwitte Date: Mon, 28 Mar 2022 20:42:21 -0500 Subject: [PATCH] Basement debug output --- src/Basement/3DBasementHT.f90 | 12 +++++++++--- 1 file changed, 9 insertions(+), 3 deletions(-) diff --git a/src/Basement/3DBasementHT.f90 b/src/Basement/3DBasementHT.f90 index f890f99b0a4..c916a4cd1ca 100644 --- a/src/Basement/3DBasementHT.f90 +++ b/src/Basement/3DBasementHT.f90 @@ -9768,14 +9768,23 @@ SUBROUTINE CalcTearth(IEXT,JEXT,DZ,DZP,TG,CVG) DO COUNT1=1,NZBG-1 CONST(COUNT1,1)=TCOND*3600.d0/SoilDens/CG/DZ(COUNT1)/DZP(COUNT1-1) + IF (isnan(CONST(COUNT1,1))) THEN + WRITE (DebugOutFile,*) ' CalcTearth: COUNT1=', COUNT1,' CONST(COUNT1,1)=', CONST(COUNT1,1), ' TCOND=', TCOND, ' SoilDens=', SoilDens, ' CG=', CG, ' DZ(COUNT1)=', DZ(COUNT1), ' DZP(COUNT1-1)=', DZP(COUNT1-1) + ENDIF END DO !*** CONSTANTS IN POSITIVE CELL DIRECTION DO COUNT1=0,NZBG-2 CONST(COUNT1,2)=TCOND*3600.d0/SoilDens/CG/DZ(COUNT1)/DZP(COUNT1) + IF (isnan(CONST(COUNT1,2))) THEN + WRITE (DebugOutFile,*) ' CalcTearth: COUNT1=', COUNT1,' CONST(COUNT1,2)=', CONST(COUNT1,2), ' TCOND=', TCOND, ' SoilDens=', SoilDens, ' CG=', CG, ' DZ(COUNT1)=', DZ(COUNT1), ' DZP(COUNT1)=', DZP(COUNT1) + ENDIF END DO !*** CONSTANT IN POSITIVE CELL DIRECTION FOR DEEP GROUND CONDITION CONST(NZBG-1,2)=TCOND*7200.d0/SoilDens/CG/DZ(NZBG-1)/DZ(NZBG-1) + IF (isnan(CONST(COUNT1,2))) THEN + WRITE (DebugOutFile,*) ' CalcTearth: NZBG-1=', (NZBG-1),' CONST(NZBG-1,2)=', CONST(NZBG-1,2), ' TCOND=', TCOND, ' SoilDens=', SoilDens, ' CG=', CG, ' DZ(NZBG-1)=', DZ(NZBG-1) + ENDIF !*** FOR FIXED TEMPERATURE LOWER BOUNDARY CONDITION, SET !*** BOUNDARY VALUE IF (.not. SameString(FIXBC,'FALSE')) TDEEP=TG(NZBG) @@ -10436,9 +10445,6 @@ SUBROUTINE TRIDI1D (A,B,C,X,R,N) X(1)=R(1) DO COUNT1=2,N X(COUNT1)=R(COUNT1)-A(COUNT1)*X(COUNT1-1) - IF (isnan(X(COUNT1))) THEN - WRITE (DebugOutFile,*) ' TRIDI1D: X(COUNT1)=', COUNT1,' X(Count1)=', X(COUNT1), ' R(COUNT1)=', R(COUNT1), ' A(COUNT1)=', A(COUNT1) - ENDIF END DO END SUBROUTINE TRIDI1D