Merge pull request #1080 from bnmurphy/bc_aeronum

Fix Aerosol BC Processing for non-CMAQ Inputs
USEPA · Apr 12, 2024 · 5638ccd · 5638ccd
2 parents c58dbf7 + ee7831e
commit 5638ccd
Showing 1 changed file with 3 additions and 3 deletions.
diff --git a/CCTM/src/aero/aero6/AERO_DATA.F b/CCTM/src/aero/aero6/AERO_DATA.F
@@ -1905,7 +1905,7 @@ Subroutine CHECK_AERO_ICBC( IBCON0, LM2WET, USE_M2, IS_BC,
          AER_PAR ( 1, IMODE, 8 ) = 0.    
 
          ! Begin Checking Aerosol Parameters. 
-         IF ( M3DRY .LT. 1.0e-30 .OR. M3WET .LT. 1.0e-30 ) THEN
+         IF ( M3DRY .LT. 1.1e-30 .OR. M3WET .LT. 1.1e-30 ) THEN
              ! Dry or Wet Mass is below limit -> reset distribution
              LSTAT( IMODE ) = 11
 
@@ -1943,7 +1943,7 @@ Subroutine CHECK_AERO_ICBC( IBCON0, LM2WET, USE_M2, IS_BC,
              DGWET = ( M3WET / NUM * EXP( -4.5 * L2SG ) )  ** ( ONE_THIRD )
              M2WET = EXP( ONE_THIRD * LOG( NUM ) + TWO_THIRDS * LOG( M3WET ) - L2SG )
 
-         ELSE IF ( NUM .LT. 1.0e-30 ) THEN
+         ELSE IF ( NUM .LT. 1.1e-30 ) THEN
              ! Number is below limit -> reset Num and Surf with
              !  default dry diameter and default standard deviation
              LSTAT( IMODE ) = 12
@@ -1958,7 +1958,7 @@ Subroutine CHECK_AERO_ICBC( IBCON0, LM2WET, USE_M2, IS_BC,
              DGWET = ( M3WET / NUM * EXP( -4.5 * L2SG ) )  ** ( ONE_THIRD )
              M2WET = EXP( ONE_THIRD * LOG( NUM ) + TWO_THIRDS * LOG( M3WET ) - L2SG )
 
-         ELSE IF ( M2 .LT. 1.0e-30 .OR. .NOT. USE_M2 ) THEN                
+         ELSE IF ( M2 .LT. 1.1e-30 .OR. .NOT. USE_M2 ) THEN                
              ! Mass and Number are valid. Surface area is either invalid
              ! or the user has issued an override. The mass and number
              ! concentrations specify the dry and wet diameters, and the