From ee7831e6cdeaeefad5cd3e36929dc0ec14500a4d Mon Sep 17 00:00:00 2001 From: Ben Murphy Date: Fri, 2 Feb 2024 08:45:48 -0500 Subject: [PATCH] Implements stricter tolerance for detecting lower bound number, third moment and surface area in aerosol boundary condition processing. On branch bc_aeronum modified: CCTM/src/aero/aero6/AERO_DATA.F --- CCTM/src/aero/aero6/AERO_DATA.F | 6 +++--- 1 file changed, 3 insertions(+), 3 deletions(-) diff --git a/CCTM/src/aero/aero6/AERO_DATA.F b/CCTM/src/aero/aero6/AERO_DATA.F index 8ea8d1ab2..05da16897 100644 --- a/CCTM/src/aero/aero6/AERO_DATA.F +++ b/CCTM/src/aero/aero6/AERO_DATA.F @@ -1897,7 +1897,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 @@ -1935,7 +1935,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 @@ -1950,7 +1950,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