Merge pull request #12482 from ericvmueller/master

FDS Source: prevent successive rescaling of LPC%LENGTH
firemodels · Feb 15, 2024 · b319a24 · b319a24
2 parents 39c2e1b + f221115
commit b319a24
Showing 1 changed file with 4 additions and 4 deletions.
diff --git a/Source/part.f90 b/Source/part.f90
@@ -1586,11 +1586,11 @@ SUBROUTINE INITIALIZE_SINGLE_PARTICLE
          SELECT CASE (LP_SF%GEOMETRY)
             CASE (SURF_CARTESIAN)
                B1%AREA = 2._EB*AREA
-               LPC%LENGTH = LPC%LENGTH*SQRT(SCALE_FACTOR)
+               LPC%LENGTH = SQRT(AREA)
                IF (LP_SF%THERMAL_BC_INDEX==THERMALLY_THICK) LP%MASS = LPC%INITIAL_MASS*SQRT(SCALE_FACTOR)
             CASE (SURF_CYLINDRICAL)
                B1%AREA = AREA*PI
-               LPC%LENGTH = LPC%LENGTH*SCALE_FACTOR
+               LPC%LENGTH = 0.5_EB*AREA/LP%RADIUS
                IF (LP_SF%THERMAL_BC_INDEX==THERMALLY_THICK) LP%MASS = LPC%INITIAL_MASS*SCALE_FACTOR
             CASE (SURF_SPHERICAL)
                B1%AREA = 4._EB*PI*LP%RADIUS**2._EB
@@ -1606,11 +1606,11 @@ SUBROUTINE INITIALIZE_SINGLE_PARTICLE
          SELECT CASE (LP_SF%GEOMETRY)
             CASE (SURF_CARTESIAN)
                B1%AREA = LP_VOLUME/LP%RADIUS
-               LPC%LENGTH = LPC%LENGTH*SQRT(SCALE_FACTOR)
+               LPC%LENGTH = SQRT(0.5_EB*B1%AREA)
                IF (LP_SF%THERMAL_BC_INDEX==THERMALLY_THICK) LP%MASS = LPC%INITIAL_MASS*SQRT(SCALE_FACTOR)
             CASE (SURF_CYLINDRICAL)
                B1%AREA = 2._EB*LP_VOLUME/LP%RADIUS
-               LPC%LENGTH = LPC%LENGTH*SCALE_FACTOR
+               LPC%LENGTH = LP_VOLUME/(PI*LP%RADIUS**2)
                IF (LP_SF%THERMAL_BC_INDEX==THERMALLY_THICK) LP%MASS = LPC%INITIAL_MASS*SCALE_FACTOR
             CASE (SURF_SPHERICAL)
                B1%AREA = 4._EB*PI*LP%RADIUS**2._EB