Basement debug output

NREL · Mar 29, 2022 · b7b82bc · b7b82bc · nrel-bot · Mar 29, 2022
1 parent 65f7235
commit b7b82bc
Showing 1 changed file with 9 additions and 3 deletions.
diff --git 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