revert 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, self%fwet(3), __RC__ )
+                          SUWT, SUPSO4, SUPSO4WT, PSO4, PSO4WET, __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), &

Process_Library/GOCART2G_Process.F90

@@ -3413,7 +3413,6 @@ subroutine WetRemovalGOCART2G ( km, klid, n1, n2, bin_ind, cdt, aero_type, kin,
         if (tmpu(i,j,k) < 258d0 .and. .not.snow_scavenging) then
             F = 0.d0
         endif
-
         effRemoval = fwet
         DC(n) = aerosol(i,j,k) * F * effRemoval *(1.-exp(-BT))
         if (DC(n).lt.0.) DC(n) = 0.
@@ -7268,7 +7267,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, fwetSO4, rc )
+                               pso4, pso4wet, rc )
 
 
 ! !USES:
@@ -7284,8 +7283,7 @@ subroutine SU_Wet_Removal ( km, nbins, klid, cdt, kin, grav, airMolWght, delp, f
    real, intent(in)    :: airMolWght ! air molecular weight [kg]
    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, dimension(:,:,:) :: h2o2_int   
    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]
@@ -7352,7 +7350,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 :: F0_ls 
+   real, parameter :: F0_ls = 1.0
    real, parameter :: XL_ls = 5.0e-4
    real, parameter :: B0_cv = 1.5e-3
    real, parameter :: F0_cv = 0.3
@@ -7390,7 +7388,6 @@ 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