module seq_drydep_mod 2,6
!========================================================================
! Module for handling dry depostion of tracers.
! This module is shared by land and atmosphere models for the computations of
! dry deposition of tracers
!
! !REVISION HISTORY:
! 2008-Nov-12 - F. Vitt - creation.
! 2009-Feb-19 - E. Kluzek - merge shr_drydep_tables module in.
! 2009-Feb-20 - E. Kluzek - use shr_ coding standards, and check for namelist file.
! 2009-Feb-20 - E. Kluzek - Put _r8 on all constants, remove namelist read out.
! 2009-Mar-23 - F. Vitt - Some corrections/cleanup and addition of drydep_method.
! 2009-Mar-27 - E. Kluzek - Get description and units from J.F. Lamarque.
!========================================================================
! !USES:
use shr_sys_mod
, only : shr_sys_abort
use shr_log_mod, only : s_loglev => shr_log_Level
use shr_kind_mod, only : r8 => shr_kind_r8, CS => SHR_KIND_CS, CX => SHR_KIND_CX
use shr_const_mod, only : SHR_CONST_G, SHR_CONST_RDAIR, &
SHR_CONST_CPDAIR, SHR_CONST_MWWV
implicit none
save
private
! !PUBLIC MEMBER FUNCTIONS
public :: seq_drydep_read ! Read namelist
public :: seq_drydep_init ! Initialization of drydep data
public :: seq_drydep_setHCoeff ! Calculate Henry's law coefficients
! !PRIVATE ARRAY SIZES
integer, private, parameter :: maxspc = 100 ! Maximum number of species
integer, private, parameter :: n_species_table = 53 ! Number of species to work with
integer, private, parameter :: NSeas = 5 ! Number of seasons
integer, private, parameter :: NLUse = 11 ! Number of land-use types
! !PUBLIC DATA MEMBERS:
! method specification
character(16),public,parameter :: DD_XATM = 'xactive_atm'! dry-dep atmosphere
character(16),public,parameter :: DD_XLND = 'xactive_lnd'! dry-dep land
character(16),public,parameter :: DD_TABL = 'table' ! dry-dep table (atm and lnd)
character(16),public :: drydep_method = DD_XLND ! Which option choosen
real(r8), public, parameter :: ph = 1.e-5_r8 ! measure of the acidity (dimensionless)
logical, public :: lnd_drydep ! If dry-dep fields passed
integer, public :: n_drydep = 0 ! Number in drypdep list
character(len=32), public, dimension(maxspc) :: drydep_list = '' ! List of dry-dep species
character(len=CS), public :: drydep_fields_token = '' ! First drydep fields token
real(r8), public, allocatable, dimension(:) :: foxd ! reactivity factor for oxidation (dimensioness)
real(r8), public, allocatable, dimension(:) :: drat ! ratio of molecular diffusivity (D_H2O/D_species; dimensionless)
integer, public, allocatable, dimension(:) :: mapping ! mapping to species table
! --- Indices for each species ---
integer, public :: h2_ndx, ch4_ndx, co_ndx, pan_ndx, mpan_ndx, so2_ndx, o3_ndx, o3a_ndx
!---------------------------------------------------------------------------
! Table 1 from Wesely, Atmos. Environment, 1989, p1293
! Table 2 from Sheih, microfiche PB86-218104 and Walcek, Atmos. Environment, 1986, p949
! Table 3-5 compiled by P. Hess
!
! index #1 : season
! 1 -> midsummer with lush vegetation
! 2 -> autumn with unharvested cropland
! 3 -> late autumn after frost, no snow
! 4 -> winter, snow on ground, and subfreezing
! 5 -> transitional spring with partially green short annuals
!
! index #2 : landuse type
! 1 -> urban land
! 2 -> agricultural land
! 3 -> range land
! 4 -> deciduous forest
! 5 -> coniferous forest
! 6 -> mixed forest including wetland
! 7 -> water, both salt and fresh
! 8 -> barren land, mostly desert
! 9 -> nonforested wetland
! 10 -> mixed agricultural and range land
! 11 -> rocky open areas with low growing shrubs
!
! JFL August 2000
!---------------------------------------------------------------------------
!---------------------------------------------------------------------------
! table to parameterize the impact of soil moisture on the deposition of H2 and
! CO on soils (from Sanderson et al., J. Atmos. Chem., 46, 15-28, 2003).
!---------------------------------------------------------------------------
!--- deposition of h2 and CO on soils ---
real(r8), parameter, public :: h2_a(NLUse) = &
(/ 0.000_r8, 0.000_r8, 0.270_r8, 0.000_r8, 0.000_r8, &
0.000_r8, 0.000_r8, 0.000_r8, 0.000_r8, 0.000_r8, 0.000_r8/)
!--- deposition of h2 and CO on soils ---
real(r8), parameter, public :: h2_b(NLUse) = &
(/ 0.000_r8,-41.390_r8, -0.472_r8,-41.900_r8,-41.900_r8, &
-41.900_r8, 0.000_r8, 0.000_r8, 0.000_r8,-41.390_r8, 0.000_r8/)
!--- deposition of h2 and CO on soils ---
real(r8), parameter, public :: h2_c(NLUse) = &
(/ 0.000_r8, 16.850_r8, 1.235_r8, 19.700_r8, 19.700_r8, &
19.700_r8, 0.000_r8, 0.000_r8, 0.000_r8, 17.700_r8, 1.000_r8/)
!--- deposition of h2 and CO on soils
!
!--- ri: Richardson number (dimensionless)
!--- rlu: Resistance of leaves in upper canopy (s.m-1)
!--- rac: Aerodynamic resistance to lower canopy (s.m-1)
!--- rgss: Ground surface resistance for SO2 (s.m-1)
!--- rgso: Ground surface resistance for O3 (s.m-1)
!--- rcls: Lower canopy resistance for SO2 (s.m-1)
!--- rclo: Lower canopy resistance for O3 (s.m-1)
!
real(r8), public, dimension(NSeas,NLUse) :: ri, rlu, rac, rgss, rgso, rcls, rclo
data ri (1,1:NLUse) &
/1.e36_r8, 60._r8, 120._r8, 70._r8, 130._r8, 100._r8,1.e36_r8,1.e36_r8, 80._r8, 100._r8, 150._r8/
data rlu (1,1:NLUse) &
/1.e36_r8,2000._r8,2000._r8,2000._r8,2000._r8,2000._r8,1.e36_r8,1.e36_r8,2500._r8,2000._r8,4000._r8/
data rac (1,1:NLUse) &
/ 100._r8, 200._r8, 100._r8,2000._r8,2000._r8,2000._r8, 0._r8, 0._r8, 300._r8, 150._r8, 200._r8/
data rgss(1,1:NLUse) &
/ 400._r8, 150._r8, 350._r8, 500._r8, 500._r8, 100._r8, 0._r8,1000._r8, 0._r8, 220._r8, 400._r8/
data rgso(1,1:NLUse) &
/ 300._r8, 150._r8, 200._r8, 200._r8, 200._r8, 300._r8,2000._r8, 400._r8,1000._r8, 180._r8, 200._r8/
data rcls(1,1:NLUse) &
/1.e36_r8,2000._r8,2000._r8,2000._r8,2000._r8,2000._r8,1.e36_r8,1.e36_r8,2500._r8,2000._r8,4000._r8/
data rclo(1,1:NLUse) &
/1.e36_r8,1000._r8,1000._r8,1000._r8,1000._r8,1000._r8,1.e36_r8,1.e36_r8,1000._r8,1000._r8,1000._r8/
data ri (2,1:NLUse) &
/1.e36_r8,1.e36_r8,1.e36_r8,1.e36_r8, 250._r8, 500._r8,1.e36_r8,1.e36_r8,1.e36_r8,1.e36_r8,1.e36_r8/
data rlu (2,1:NLUse) &
/1.e36_r8,9000._r8,9000._r8,9000._r8,4000._r8,8000._r8,1.e36_r8,1.e36_r8,9000._r8,9000._r8,9000._r8/
data rac (2,1:NLUse) &
/ 100._r8, 150._r8, 100._r8,1500._r8,2000._r8,1700._r8, 0._r8, 0._r8, 200._r8, 120._r8, 140._r8/
data rgss(2,1:NLUse) &
/ 400._r8, 200._r8, 350._r8, 500._r8, 500._r8, 100._r8, 0._r8,1000._r8, 0._r8, 300._r8, 400._r8/
data rgso(2,1:NLUse) &
/ 300._r8, 150._r8, 200._r8, 200._r8, 200._r8, 300._r8,2000._r8, 400._r8, 800._r8, 180._r8, 200._r8/
data rcls(2,1:NLUse) &
/1.e36_r8,9000._r8,9000._r8,9000._r8,2000._r8,4000._r8,1.e36_r8,1.e36_r8,9000._r8,9000._r8,9000._r8/
data rclo(2,1:NLUse) &
/1.e36_r8, 400._r8, 400._r8, 400._r8,1000._r8, 600._r8,1.e36_r8,1.e36_r8, 400._r8, 400._r8, 400._r8/
data ri (3,1:NLUse) &
/1.e36_r8,1.e36_r8,1.e36_r8,1.e36_r8, 250._r8, 500._r8,1.e36_r8,1.e36_r8,1.e36_r8,1.e36_r8,1.e36_r8/
data rlu (3,1:NLUse) &
/1.e36_r8,1.e36_r8,9000._r8,9000._r8,4000._r8,8000._r8,1.e36_r8,1.e36_r8,9000._r8,9000._r8,9000._r8/
data rac (3,1:NLUse) &
/ 100._r8, 10._r8, 100._r8,1000._r8,2000._r8,1500._r8, 0._r8, 0._r8, 100._r8, 50._r8, 120._r8/
data rgss(3,1:NLUse) &
/ 400._r8, 150._r8, 350._r8, 500._r8, 500._r8, 200._r8, 0._r8,1000._r8, 0._r8, 200._r8, 400._r8/
data rgso(3,1:NLUse) &
/ 300._r8, 150._r8, 200._r8, 200._r8, 200._r8, 300._r8,2000._r8, 400._r8,1000._r8, 180._r8, 200._r8/
data rcls(3,1:NLUse) &
/1.e36_r8,1.e36_r8,9000._r8,9000._r8,3000._r8,6000._r8,1.e36_r8,1.e36_r8,9000._r8,9000._r8,9000._r8/
data rclo(3,1:NLUse) &
/1.e36_r8,1000._r8, 400._r8, 400._r8,1000._r8, 600._r8,1.e36_r8,1.e36_r8, 800._r8, 600._r8, 600._r8/
data ri (4,1:NLUse) &
/1.e36_r8,1.e36_r8,1.e36_r8,1.e36_r8, 400._r8, 800._r8,1.e36_r8,1.e36_r8,1.e36_r8,1.e36_r8,1.e36_r8/
data rlu (4,1:NLUse) &
/1.e36_r8,1.e36_r8,1.e36_r8,1.e36_r8,6000._r8,9000._r8,1.e36_r8,1.e36_r8,9000._r8,9000._r8,9000._r8/
data rac (4,1:NLUse) &
/ 100._r8, 10._r8, 10._r8,1000._r8,2000._r8,1500._r8, 0._r8, 0._r8, 50._r8, 10._r8, 50._r8/
data rgss(4,1:NLUse) &
/ 100._r8, 100._r8, 100._r8, 100._r8, 100._r8, 100._r8, 0._r8,1000._r8, 100._r8, 100._r8, 50._r8/
data rgso(4,1:NLUse) &
/ 600._r8,3500._r8,3500._r8,3500._r8,3500._r8,3500._r8,2000._r8, 400._r8,3500._r8,3500._r8,3500._r8/
data rcls(4,1:NLUse) &
/1.e36_r8,1.e36_r8,1.e36_r8,9000._r8, 200._r8, 400._r8,1.e36_r8,1.e36_r8,9000._r8,1.e36_r8,9000._r8/
data rclo(4,1:NLUse) &
/1.e36_r8,1000._r8,1000._r8, 400._r8,1500._r8, 600._r8,1.e36_r8,1.e36_r8, 800._r8,1000._r8, 800._r8/
data ri (5,1:NLUse) &
/1.e36_r8, 120._r8, 240._r8, 140._r8, 250._r8, 190._r8,1.e36_r8,1.e36_r8, 160._r8, 200._r8, 300._r8/
data rlu (5,1:NLUse) &
/1.e36_r8,4000._r8,4000._r8,4000._r8,2000._r8,3000._r8,1.e36_r8,1.e36_r8,4000._r8,4000._r8,8000._r8/
data rac (5,1:NLUse) &
/ 100._r8, 50._r8, 80._r8,1200._r8,2000._r8,1500._r8, 0._r8, 0._r8, 200._r8, 60._r8, 120._r8/
data rgss(5,1:NLUse) &
/ 500._r8, 150._r8, 350._r8, 500._r8, 500._r8, 200._r8, 0._r8,1000._r8, 0._r8, 250._r8, 400._r8/
data rgso(5,1:NLUse) &
/ 300._r8, 150._r8, 200._r8, 200._r8, 200._r8, 300._r8,2000._r8, 400._r8,1000._r8, 180._r8, 200._r8/
data rcls(5,1:NLUse) &
/1.e36_r8,4000._r8,4000._r8,4000._r8,2000._r8,3000._r8,1.e36_r8,1.e36_r8,4000._r8,4000._r8,8000._r8/
data rclo(5,1:NLUse) &
/1.e36_r8,1000._r8, 500._r8, 500._r8,1500._r8, 700._r8,1.e36_r8,1.e36_r8, 600._r8, 800._r8, 800._r8/
!---------------------------------------------------------------------------
! ... roughness length
!---------------------------------------------------------------------------
real(r8), public, dimension(NSeas,NLUse) :: z0
data z0 (1,1:NLUse) &
/1.000_r8,0.250_r8,0.050_r8,1.000_r8,1.000_r8,1.000_r8,0.0006_r8,0.002_r8,0.150_r8,0.100_r8,0.100_r8/
data z0 (2,1:NLUse) &
/1.000_r8,0.100_r8,0.050_r8,1.000_r8,1.000_r8,1.000_r8,0.0006_r8,0.002_r8,0.100_r8,0.080_r8,0.080_r8/
data z0 (3,1:NLUse) &
/1.000_r8,0.005_r8,0.050_r8,1.000_r8,1.000_r8,1.000_r8,0.0006_r8,0.002_r8,0.100_r8,0.020_r8,0.060_r8/
data z0 (4,1:NLUse) &
/1.000_r8,0.001_r8,0.001_r8,1.000_r8,1.000_r8,1.000_r8,0.0006_r8,0.002_r8,0.001_r8,0.001_r8,0.040_r8/
data z0 (5,1:NLUse) &
/1.000_r8,0.030_r8,0.020_r8,1.000_r8,1.000_r8,1.000_r8,0.0006_r8,0.002_r8,0.010_r8,0.030_r8,0.060_r8/
!real(r8), private, dimension(11,5), parameter :: z0xxx = reshape ( &
! (/ 1.000,0.250,0.050,1.000,1.000,1.000,0.0006,0.002,0.150,0.100,0.100 , &
! 1.000,0.100,0.050,1.000,1.000,1.000,0.0006,0.002,0.100,0.080,0.080 , &
! 1.000,0.005,0.050,1.000,1.000,1.000,0.0006,0.002,0.100,0.020,0.060 , &
! 1.000,0.001,0.001,1.000,1.000,1.000,0.0006,0.002,0.001,0.001,0.040 , &
! 1.000,0.030,0.020,1.000,1.000,1.000,0.0006,0.002,0.010,0.030,0.060 /), (/11,5/) )
!---------------------------------------------------------------------------
! public chemical data
!---------------------------------------------------------------------------
!--- data for foxd (reactivity factor for oxidation) ----
real(r8), public, parameter :: dfoxd(n_species_table) = &
(/ 1._r8 &
,1._r8 &
,1._r8 &
,.1_r8 &
,1.e-36_r8 &
,1.e-36_r8 &
,1._r8 &
,.1_r8 &
,1.e-36_r8 &
,0._r8 &
,0._r8 &
,.1_r8 &
,1.e-36_r8 &
,1.e-36_r8 &
,1.e-36_r8 &
,.1_r8 &
,1._r8 &
,1.e-36_r8 &
,.1_r8 &
,1._r8 &
,1.e-36_r8 &
,.1_r8 &
,.1_r8 &
,.1_r8 &
,.1_r8 &
,1.e-36_r8 &
,1.e-36_r8 &
,.1_r8 &
,1.e-36_r8 &
,.1_r8 &
,1.e-36_r8 &
,.1_r8 &
,.1_r8 &
,1.e-36_r8 &
,1.e-36_r8 &
,1.e-36_r8 &
,1.e-36_r8 &
,.1_r8 &
,1.e-36_r8 &
,.1_r8 &
,1.e-36_r8 &
,.1_r8 &
,.1_r8 &
,.1_r8 &
,1.e-36_r8 &
,1.e-36_r8 &
,1.e-36_r8 &
,1.e-36_r8 &
,1.e-36_r8 &
,.1_r8 &
,.1_r8 &
,.1_r8 &
,1.e-36_r8 &
/)
!EOP
! PRIVATE DATA:
Interface seq_drydep_setHCoeff ! overload subroutine
Module Procedure set_hcoeff_scalar
Module Procedure set_hcoeff_vector
End Interface
real(r8), private, parameter :: inf = O'0777600000000000000000' !--- infinity to initialize data to
real(r8), private, parameter :: small_value = 1.e-36_r8 !--- smallest value to use ---
!---------------------------------------------------------------------------
! private chemical data
!---------------------------------------------------------------------------
!--- Names of species that can work with ---
character(len=20), private, parameter :: species_name_table(n_species_table) = &
(/ 'OX ' &
,'H2O2 ' &
,'OH ' &
,'HO2 ' &
,'CO ' &
,'CH4 ' &
,'CH3O2 ' &
,'CH3OOH ' &
,'CH2O ' &
,'CHOOH ' &
,'NO ' &
,'NO2 ' &
,'HNO3 ' &
,'CO2 ' &
,'NH3 ' &
,'N2O5 ' &
,'NO3 ' &
,'CH3OH ' &
,'HO2NO2 ' &
,'O1D ' &
,'C2H6 ' &
,'C2H5O2 ' &
,'PO2 ' &
,'MACRO2 ' &
,'ISOPO2 ' &
,'C4H10 ' &
,'CH3CHO ' &
,'C2H5OOH ' &
,'C3H6 ' &
,'POOH ' &
,'C2H4 ' &
,'PAN ' &
,'CH3COOOH' &
,'C10H16 ' &
,'CHOCHO ' &
,'CH3COCHO' &
,'GLYALD ' &
,'CH3CO3 ' &
,'C3H8 ' &
,'C3H7O2 ' &
,'CH3COCH3' &
,'C3H7OOH ' &
,'RO2 ' &
,'ROOH ' &
,'Rn ' &
,'ISOP ' &
,'MVK ' &
,'MACR ' &
,'C2H5OH ' &
,'ONITR ' &
,'ONIT ' &
,'ISOPNO3 ' &
,'HYDRALD ' &
/)
!--- data for effective Henry's Law coefficient ---
real(r8), private, parameter :: dheff(n_species_table*6) = &
(/1.15e-02_r8, 2560._r8,0._r8 , 0._r8,0._r8 , 0._r8 &
,8.33e+04_r8, 7379._r8,2.2e-12_r8,-3730._r8,0._r8 , 0._r8 &
,3.00e+01_r8, 0._r8,0._r8 , 0._r8,0._r8 , 0._r8 &
,2.00e+03_r8, 6600._r8,3.5e-05_r8, 0._r8,0._r8 , 0._r8 &
,1.00e-03_r8, 0._r8,0._r8 , 0._r8,0._r8 , 0._r8 &
,1.70e-03_r8, 0._r8,0._r8 , 0._r8,0._r8 , 0._r8 &
,7.47e+00_r8, 5241._r8,0._r8 , 0._r8,0._r8 , 0._r8 &
,3.11e+02_r8, 5241._r8,0._r8 , 0._r8,0._r8 , 0._r8 &
,6.30e+03_r8, 6425._r8,0._r8 , 0._r8,0._r8 , 0._r8 &
,5.53e+03_r8, 5700._r8,1.8e-04_r8,-1510._r8,0._r8 , 0._r8 &
,1.90e-03_r8, 1480._r8,0._r8 , 0._r8,0._r8 , 0._r8 &
,6.40e-03_r8, 2500._r8,0._r8 , 0._r8,0._r8 , 0._r8 &
,0._r8 , 0._r8,2.6e+06_r8, 8700._r8,0._r8 , 0._r8 &
,3.40e-02_r8, 2420._r8,4.5e-07_r8,-1000._r8,3.6e-11_r8,-1760._r8 &
,7.40e+01_r8, 3400._r8,1.7e-05_r8, -450._r8,1.0e-14_r8,-6716._r8 &
,2.14e+00_r8, 3362._r8,0._r8 , 0._r8,0._r8 , 0._r8 &
,0.65e+00_r8, 0._r8,0._r8 , 0._r8,0._r8 , 0._r8 &
,2.20e+02_r8, 4934._r8,0._r8 , 0._r8,0._r8 , 0._r8 &
,0._r8 , 0._r8,3.2e+01_r8, 0._r8,0._r8 , 0._r8 &
,1.00e-16_r8, 0._r8,0._r8 , 0._r8,0._r8 , 0._r8 &
,1.70e-03_r8, 0._r8,0._r8 , 0._r8,0._r8 , 0._r8 &
,7.47e+00_r8, 5241._r8,0._r8 , 0._r8,0._r8 , 0._r8 &
,7.47e+00_r8, 5241._r8,0._r8 , 0._r8,0._r8 , 0._r8 &
,7.47e+00_r8, 5241._r8,0._r8 , 0._r8,0._r8 , 0._r8 &
,7.47e+00_r8, 5241._r8,0._r8 , 0._r8,0._r8 , 0._r8 &
,1.70e-03_r8, 0._r8,0._r8 , 0._r8,0._r8 , 0._r8 &
,1.14e+01_r8, 6267._r8,0._r8 , 0._r8,0._r8 , 0._r8 &
,3.36e+02_r8, 5995._r8,0._r8 , 0._r8,0._r8 , 0._r8 &
,1.70e-03_r8, 0._r8,0._r8 , 0._r8,0._r8 , 0._r8 &
,2.20e+02_r8, 5653._r8,0._r8 , 0._r8,0._r8 , 0._r8 &
,1.70e-03_r8, 0._r8,0._r8 , 0._r8,0._r8 , 0._r8 &
,5.00e+00_r8, 0._r8,0._r8 , 0._r8,0._r8 , 0._r8 &
,8.37e+02_r8, 5308._r8,1.8e-04_r8,-1510._r8,0._r8 , 0._r8 &
,1.70e-03_r8, 0._r8,0._r8 , 0._r8,0._r8 , 0._r8 &
,3.00e+05_r8, 0._r8,0._r8 , 0._r8,0._r8 , 0._r8 &
,3.71e+03_r8, 7541._r8,0._r8 , 0._r8,0._r8 , 0._r8 &
,4.14e+04_r8, 4630._r8,0._r8 , 0._r8,0._r8 , 0._r8 &
,7.47e+00_r8, 5241._r8,0._r8 , 0._r8,0._r8 , 0._r8 &
,1.45e-03_r8, 2700._r8,0._r8 , 0._r8,0._r8 , 0._r8 &
,3.00e+06_r8, 0._r8,0._r8 , 0._r8,0._r8 , 0._r8 &
,2.70e+01_r8, 5300._r8,0._r8 , 0._r8,0._r8 , 0._r8 &
,3.36e+02_r8, 5995._r8,0._r8 , 0._r8,0._r8 , 0._r8 &
,7.47e+00_r8, 5241._r8,0._r8 , 0._r8,0._r8 , 0._r8 &
,3.36e+02_r8, 5995._r8,0._r8 , 0._r8,0._r8 , 0._r8 &
,0.00e+00_r8, 0._r8,0._r8 , 0._r8,0._r8 , 0._r8 &
,1.70e-03_r8, 0._r8,0._r8 , 0._r8,0._r8 , 0._r8 &
,1.70e-03_r8, 0._r8,0._r8 , 0._r8,0._r8 , 0._r8 &
,1.70e-03_r8, 0._r8,0._r8 , 0._r8,0._r8 , 0._r8 &
,2.00e+02_r8, 6500._r8,0._r8 , 0._r8,0._r8 , 0._r8 &
,7.51e+03_r8, 6485._r8,0._r8 , 0._r8,0._r8 , 0._r8 &
,1.00e+03_r8, 6000._r8,0._r8 , 0._r8,0._r8 , 0._r8 &
,1.00e+01_r8, 0._r8,0._r8 , 0._r8,0._r8 , 0._r8 &
,7.00e+01_r8, 6000._r8,0._r8 , 0._r8,0._r8 , 0._r8 &
/)
real(r8), private, parameter :: wh2o = SHR_CONST_MWWV
real(r8), private, parameter :: mol_wgts(n_species_table) = &
(/ 47.9981995_r8, 34.0135994_r8, 17.0067997_r8, 33.0061989_r8, 28.0104008_r8, &
16.0405998_r8, 47.0320015_r8, 48.0393982_r8, 30.0251999_r8, 46.0246010_r8, &
30.0061398_r8, 46.0055389_r8, 63.0123405_r8, 44.0098000_r8, 17.0289402_r8, &
108.010483_r8, 62.0049400_r8, 32.0400009_r8, 79.0117416_r8, 15.9994001_r8, &
30.0664005_r8, 61.0578003_r8, 91.0830002_r8, 119.093399_r8, 117.119797_r8, &
58.1180000_r8, 44.0509987_r8, 62.0652008_r8, 42.0774002_r8, 92.0904007_r8, &
28.0515995_r8, 121.047943_r8, 76.0497971_r8, 136.228394_r8, 58.0355988_r8, &
72.0614014_r8, 60.0503998_r8, 75.0423965_r8, 44.0922012_r8, 75.0836029_r8, &
58.0768013_r8, 76.0910034_r8, 31.9988003_r8, 33.0061989_r8, 222.000000_r8, &
68.1141968_r8, 70.0877991_r8, 70.0877991_r8, 46.0657997_r8, 147.125946_r8, &
119.074341_r8, 162.117935_r8, 100.112999_r8 /)
!===============================================================================
CONTAINS
!===============================================================================
!====================================================================================
subroutine seq_drydep_read(NLFilename, seq_drydep_fields) 1,7
!========================================================================
! reads drydep_inparm namelist and sets up CCSM driver list of fields for
! land-atmosphere communications.
!
! !REVISION HISTORY:
! 2009-Feb-20 - E. Kluzek - Separate out as subroutine from previous input_init
!========================================================================
use shr_file_mod,only : shr_file_getUnit, shr_file_freeUnit
use shr_log_mod, only : s_logunit => shr_log_Unit
implicit none
character(len=*), intent(in) :: NLFilename ! Namelist filename
character(len=*), intent(out) :: seq_drydep_fields
!----- local -----
integer :: i ! Indices
integer :: unitn ! namelist unit number
integer :: ierr ! error code
logical :: exists ! if file exists or not
character(len=8) :: token ! dry dep field name to add
!----- formats -----
character(*),parameter :: subName = '(seq_drydep_read) '
character(*),parameter :: F00 = "('(seq_drydep_read) ',8a)"
character(*),parameter :: FI1 = "('(seq_drydep_init) ',a,I2)"
namelist /drydep_inparm/ drydep_list, drydep_method
!-----------------------------------------------------------------------------
! Read namelist and figure out the drydep field list to pass
! First check if file exists and if not, n_drydep will be zero
!-----------------------------------------------------------------------------
!--- Open and read namelist ---
if ( len_trim(NLFilename) == 0 )then
call shr_sys_abort( subName//'ERROR: nlfilename not set' )
end if
inquire( file=trim(NLFileName), exist=exists)
if ( exists ) then
unitn = shr_file_getUnit()
open( unitn, file=trim(NLFilename), status='old' )
if ( s_loglev > 0 ) write(s_logunit,F00) &
'Read in drydep_inparm namelist from: ', trim(NLFilename)
ierr = 1
do while ( ierr /= 0 )
read(unitn, drydep_inparm, iostat=ierr)
if (ierr < 0) then
call shr_sys_abort( subName//'ERROR: encountered end-of-file on namelist read' )
endif
end do
close( unitn )
call shr_file_freeUnit( unitn )
end if
n_drydep = 0
!--- Loop over species to fill list of fields to communicate for drydep ---
seq_drydep_fields = ' '
do i=1,maxspc
if ( len_trim(drydep_list(i))==0 ) exit
write(token,333) i
seq_drydep_fields = trim(seq_drydep_fields)//':'//trim(token)
if ( i == 1 ) then
drydep_fields_token = token
endif
n_drydep = n_drydep+1
enddo
!--- Make sure method is valid and determine if land is passing drydep fields ---
lnd_drydep = n_drydep>0 .and. drydep_method == DD_XLND
if ( s_loglev > 0 ) then
write(s_logunit,*) 'seq_drydep_read: drydep_method: ', trim(drydep_method)
if ( n_drydep == 0 )then
write(s_logunit,F00) 'No dry deposition fields will be transfered'
else
write(s_logunit,FI1) 'Number of dry deposition fields transfered is ', &
n_drydep
end if
end if
if ( trim(drydep_method)/=trim(DD_XATM) .and. &
trim(drydep_method)/=trim(DD_XLND) .and. &
trim(drydep_method)/=trim(DD_TABL) ) then
if ( s_loglev > 0 ) then
write(s_logunit,*) 'seq_drydep_read: drydep_method : ', trim(drydep_method)
write(s_logunit,*) 'seq_drydep_read: drydep_method must be set to : ', &
DD_XATM,', ', DD_XLND,', or ', DD_TABL
end if
call shr_sys_abort('seq_drydep_read: incorrect dry deposition method specification')
endif
333 format ('dd',i3.3)
end subroutine seq_drydep_read
!====================================================================================
subroutine seq_drydep_init( ) 1,1
!========================================================================
! Initialization of dry deposition fields
! reads drydep_inparm namelist and sets up CCSM driver list of fields for
! land-atmosphere communications.
! !REVISION HISTORY:
! 2008-Nov-12 - F. Vitt - first version
! 2009-Feb-20 - E. Kluzek - Check for existance of file if not return, set n_drydep=0
! 2009-Feb-20 - E. Kluzek - Move namelist read to separate subroutine
!========================================================================
use shr_log_mod, only : s_logunit => shr_log_Unit
implicit none
!----- local -----
integer :: i, l ! Indices
character(len=32) :: test_name ! field test name
!----- formats -----
character(*),parameter :: subName = '(seq_drydep_init) '
character(*),parameter :: F00 = "('(seq_drydep_init) ',8a)"
!-----------------------------------------------------------------------------
! Allocate and fill foxd, drat and mapping as well as species indices
!-----------------------------------------------------------------------------
if ( n_drydep > 0 ) then
allocate( foxd(n_drydep) )
allocate( drat(n_drydep) )
allocate( mapping(n_drydep) )
foxd = inf
drat = inf
mapping(:) = 0
end if
h2_ndx=-1; ch4_ndx=-1; co_ndx=-1; mpan_ndx = -1; pan_ndx = -1; so2_ndx=-1; o3_ndx=-1
!--- Loop over drydep species that need to be worked with ---
do i=1,n_drydep
if ( len_trim(drydep_list(i))==0 ) exit
test_name = drydep_list(i)
if( trim(test_name) == 'O3' ) then
test_name = 'OX'
end if
!--- Figure out if species maps to a species in the species table ---
do l = 1,n_species_table
if( trim( test_name ) == trim( species_name_table(l) ) ) then
mapping(i) = l
exit
end if
end do
!--- If it doesn't map to a species in the species table find species close enough ---
if( mapping(i) < 1 ) then
select case( trim(test_name) )
case( 'H2' )
test_name = 'CO'
case( 'HYAC', 'CH3COOH' )
test_name = 'CH2O'
case( 'O3S', 'O3INERT', 'MPAN' )
test_name = 'OX'
case( 'ISOPOOH', 'MACROOH', 'Pb', 'XOOH', 'H2SO4' )
test_name = 'HNO3'
case( 'ALKOOH', 'MEKOOH', 'TOLOOH', 'BENOOH', 'XYLOOH', 'TERPOOH','SOGM','SOGI','SOGT','SOGB','SOGX' )
test_name = 'CH3OOH'
case( 'SOA', 'SO2', 'SO4', 'CB1', 'CB2', 'OC1', 'OC2', 'NH3', 'NH4', 'SA1', 'SA2', 'SA3', 'SA4' )
test_name = 'OX' ! this is just a place holder. values are explicitly set below
case( 'SOAM', 'SOAI', 'SOAT', 'SOAB', 'SOAX' )
test_name = 'OX' ! this is just a place holder. values are explicitly set below
case( 'O3A', 'XMPAN' )
test_name = 'OX'
case( 'XPAN' )
test_name = 'PAN'
case( 'XNO' )
test_name = 'NO'
case( 'XNO2' )
test_name = 'NO2'
case( 'XHNO3' )
test_name = 'HNO3'
case( 'XONIT' )
test_name = 'ONIT'
case( 'XONITR' )
test_name = 'ONITR'
case( 'XHO2NO2')
test_name = 'HO2NO2'
case( 'XNH4NO3' )
test_name = 'HNO3'
case( 'NH4NO3' )
test_name = 'HNO3'
case default
test_name = 'blank'
end select
!--- If found a match check the species table again ---
if( trim(test_name) /= 'blank' ) then
do l = 1,n_species_table
if( trim( test_name ) == trim( species_name_table(l) ) ) then
mapping(i) = l
exit
end if
end do
else
if ( s_loglev > 0 ) write(s_logunit,F00) trim(drydep_list(i)), &
' not in tables; will have dep vel = 0'
end if
end if
!--- Figure out the specific species indices ---
if ( trim(drydep_list(i)) == 'H2' ) h2_ndx = i
if ( trim(drydep_list(i)) == 'CO' ) co_ndx = i
if ( trim(drydep_list(i)) == 'CH4' ) ch4_ndx = i
if ( trim(drydep_list(i)) == 'MPAN' ) mpan_ndx = i
if ( trim(drydep_list(i)) == 'PAN' ) pan_ndx = i
if ( trim(drydep_list(i)) == 'SO2' ) so2_ndx = i
if ( trim(drydep_list(i)) == 'OX' .or. trim(drydep_list(i)) == 'O3' ) o3_ndx = i
if ( trim(drydep_list(i)) == 'O3A' ) o3a_ndx = i
if( mapping(i) > 0) then
l = mapping(i)
foxd(i) = dfoxd(l)
drat(i) = sqrt(mol_wgts(l)/wh2o)
endif
enddo
where( rgss < 1._r8 )
rgss = 1._r8
endwhere
where( rac < small_value)
rac = small_value
endwhere
end subroutine seq_drydep_init
!====================================================================================
subroutine set_hcoeff_scalar( sfc_temp, heff ) 1,1
!========================================================================
! Interface to seq_drydep_setHCoeff when input is scalar
! wrapper routine used when surface temperature is a scalar (single column) rather
! than an array (multiple columns).
!
! !REVISION HISTORY:
! 2008-Nov-12 - F. Vitt - first version
!========================================================================
implicit none
real(r8), intent(in) :: sfc_temp ! Input surface temperature
real(r8), intent(out) :: heff(n_drydep) ! Output Henry's law coefficients
!----- local -----
real(r8) :: sfc_temp_tmp(1) ! surface temp
sfc_temp_tmp(:) = sfc_temp
call set_hcoeff_vector( 1, sfc_temp_tmp, heff(:n_drydep) )
end subroutine set_hcoeff_scalar
!====================================================================================
subroutine set_hcoeff_vector( ncol, sfc_temp, heff ) 2,2
!========================================================================
! Interface to seq_drydep_setHCoeff when input is vector
! sets dry depositions coefficients -- used by both land and atmosphere models
! !REVISION HISTORY:
! 2008-Nov-12 - F. Vitt - first version
!========================================================================
use shr_log_mod, only : s_logunit => shr_log_Unit
implicit none
integer, intent(in) :: ncol ! Input size of surface-temp vector
real(r8), intent(in) :: sfc_temp(ncol) ! Surface temperature
real(r8), intent(out) :: heff(ncol,n_drydep) ! Henry's law coefficients
!----- local -----
real(r8), parameter :: t0 = 298._r8 ! Standard Temperature
real(r8), parameter :: ph_inv = 1._r8/ph ! Inverse of PH
integer :: m, l, id ! indices
real(r8) :: e298 ! Henry's law coefficient @ standard temperature (298K)
real(r8) :: dhr ! temperature dependence of Henry's law coefficient
real(r8) :: dk1s(ncol) ! DK Work array 1
real(r8) :: dk2s(ncol) ! DK Work array 2
real(r8) :: wrk(ncol) ! Work array
!----- formats -----
character(*),parameter :: subName = '(seq_drydep_set_hcoeff) '
character(*),parameter :: F00 = "('(seq_drydep_set_hcoeff) ',8a)"
!-------------------------------------------------------------------------------
! notes:
!-------------------------------------------------------------------------------
wrk(:) = (t0 - sfc_temp(:))/(t0*sfc_temp(:))
do m = 1,n_drydep
l = mapping(m)
id = 6*(l - 1)
e298 = dheff(id+1)
dhr = dheff(id+2)
heff(:,m) = e298*exp( dhr*wrk(:) )
!--- Calculate coefficients based on the drydep tables ---
if( dheff(id+3) /= 0._r8 .and. dheff(id+5) == 0._r8 ) then
e298 = dheff(id+3)
dhr = dheff(id+4)
dk1s(:) = e298*exp( dhr*wrk(:) )
where( heff(:,m) /= 0._r8 )
heff(:,m) = heff(:,m)*(1._r8 + dk1s(:)*ph_inv)
elsewhere
heff(:,m) = dk1s(:)*ph_inv
endwhere
end if
!--- For coefficients that are non-zero AND CO2 or NH3 handle things this way ---
if( dheff(id+5) /= 0._r8 ) then
if( trim( drydep_list(m) ) == 'CO2' .or. trim( drydep_list(m) ) == 'NH3' ) then
e298 = dheff(id+3)
dhr = dheff(id+4)
dk1s(:) = e298*exp( dhr*wrk(:) )
e298 = dheff(id+5)
dhr = dheff(id+6)
dk2s(:) = e298*exp( dhr*wrk(:) )
!--- For Carbon dioxide ---
if( trim(drydep_list(m)) == 'CO2' ) then
heff(:,m) = heff(:,m)*(1._r8 + dk1s(:)*ph_inv)*(1._r8 + dk2s(:)*ph_inv)
!--- For NH3 ---
else if( trim( drydep_list(m) ) == 'NH3' ) then
heff(:,m) = heff(:,m)*(1._r8 + dk1s(:)*ph/dk2s(:))
!--- This can't happen ---
else
write(s_logunit,F00) 'Bad species ',drydep_list(m)
call shr_sys_abort( subName//'ERROR: in assigning coefficients' )
end if
end if
end if
end do
end subroutine set_hcoeff_vector
!===============================================================================
end module seq_drydep_mod