add fwet for sulfate

Author: acollow

ESMF/GOCART2G_GridComp/SU2G_GridComp/SU2G_GridCompMod.F90

@@ -1232,7 +1232,7 @@ subroutine Run2 (GC, import, export, clock, RC)
                           delp, fMassSO4, fMassSO2, &
                           h2o2_init, ple, airdens, cn_prcp, ncn_prcp, pfl_lsan, pfi_lsan, t, &
                           nDMS, nSO2, nSO4, nMSA, DMS, SO2, SO4, dummyMSA, &
-                          SUWT, SUPSO4, SUPSO4WT, PSO4, PSO4WET, __RC__ )
+                          SUWT, SUPSO4, SUPSO4WT, PSO4, PSO4WET, self%fwet(3), __RC__ )
 
 !   Certain variables are multiplied by 1.0e-9 to convert from nanometers to meters
     call SU_Compute_Diags ( self%km, self%klid, self%rmed(nSO4), self%sigma(nSO4), self%rhop(nSO4), &

ESMF/GOCART2G_GridComp/SU2G_GridComp/SU2G_instance_SU.rc

@@ -23,6 +23,8 @@ aircraft_fuel_emission_factor: 1.0000
 # Scavenging efficiency per bin [km-1] (NOT USED UNLESS RAS IS CALLED)
 fscav: 0.0  0.0  0.4  0.4
 
+fwet: 1.0 1.0 1.0 1.0
+
 # Dry particle radius [um], used for settling
 particle_radius_microns: 0.0 0.0 0.35 0.0
 
@@ -55,4 +57,4 @@ using_GMI_H2O2: .false.
 
 # SettlingSolver options
 # Options: 'gocart' or 'ufs'
-settling_scheme: 'gocart'
+settling_scheme: 'gocart'

Process_Library/GOCART2G_Process.F90

@@ -7268,7 +7268,7 @@ subroutine SU_Wet_Removal ( km, nbins, klid, cdt, kin, grav, airMolWght, delp, f
                                h2o2_int, ple, rhoa, precc, precl, pfllsan, pfilsan, tmpu, &
                                nDMS, nSO2, nSO4, nMSA, DMS, SO2, SO4, MSA, &
                                fluxout, pSO4_colflux, pSO4wet_colflux, &
-                               pso4, pso4wet, rc )
+                               pso4, pso4wet, fwetSO4, rc )
 
 
 ! !USES:
@@ -7285,6 +7285,7 @@ subroutine SU_Wet_Removal ( km, nbins, klid, cdt, kin, grav, airMolWght, delp, f
    real, dimension(:,:,:), intent(in) :: delp   ! pressure thickness [Pa]
    real, intent(in) :: fMassSO4, fMassSO2
    real, dimension(:,:,:) :: h2o2_int
+   real, intent(in)    :: fwetSO4   ! large scale removal efficiency   
    real, pointer, dimension(:,:,:), intent(in) :: ple     ! level edge air pressure
    real, pointer, dimension(:,:,:), intent(in) :: rhoa    ! air density, [kg m-3]
    real, pointer, dimension(:,:), intent(in)   :: precc   ! total convective precip, [mm day-1]
@@ -7351,7 +7352,7 @@ subroutine SU_Wet_Removal ( km, nbins, klid, cdt, kin, grav, airMolWght, delp, f
 
 !  Rain parameters (from where?)
    real, parameter :: B0_ls = 1.0e-4
-   real, parameter :: F0_ls = 1.0
+   real :: F0_ls 
    real, parameter :: XL_ls = 5.0e-4
    real, parameter :: B0_cv = 1.5e-3
    real, parameter :: F0_cv = 0.3
@@ -7389,6 +7390,7 @@ subroutine SU_Wet_Removal ( km, nbins, klid, cdt, kin, grav, airMolWght, delp, f
 
    dims = shape(rhoa)
    i2 = dims(1); j2 = dims(2)
+   F0_ls = fwetSO4
 
 !  check if doing MSA and define nbins_ accordingly
 !   if (associated(MSA)) then