Index

Rundeck: E6TomaF40

Created: Fri May 17 03:30:39 EDT 2024

nitrate_thermo_drv

Module:

File: TRACER_drv_eqsam.f

Summary: This routine sets up for and calls the thermodynamic module for aerosol gas-particle partitioning. A version of EQSAM (eqsam_v03d) is the current thermodynamic model.
Author : Susanne Bauer
Version:

---

Declaration:

subroutine nitrate_thermo_drv(ltop)

Variables:

ntm : used from tracer_com

trm : used from tracer_com

n_clay : used from tracer_com

n_hno3 : used from tracer_com

n_nh3 : used from tracer_com

n_nh4 : used from tracer_com

n_no3p : used from tracer_com

n_seasalt1 : used from tracer_com

n_silt1 : used from tracer_com

n_silt2 : used from tracer_com

n_silt3 : used from tracer_com

n_so4 : used from tracer_com

n_so4_d1 : used from tracer_com

n_so4_d2 : used from tracer_com

n_so4_d3 : used from tracer_com

off_hno3 : used from aerosol_sources

im : used from resolution

jm : used from resolution

lm : used from resolution

t : used from atm_com

q : used from atm_com

dtsrc : used from model_com

axyp : used from geom

byaxyp : used from geom

mair : used from constant

gasc : used from constant

lhe : used from constant

tr3dsource : used from fluxes

pmid : used from atm_com

pk : used from atm_com

ma : used from atm_com

grid : used from domain_decomp_atm

getdomainbounds : used from domain_decomp_atm

taijls : used from trdiag_com

Initial Value = >taijls_loc

ijlt_ah2o : used from trdiag_com

ijlt_aph : used from trdiag_com

j : integer

l : integer

i : integer

j_0 : integer

j_1 : integer

n : integer

i_0 : integer

i_1 : integer

ltop : integer

nca : integer, parameter

Initial Value = 11

nco : integer, parameter

Initial Value = 37

iopt : integer, parameter

Initial Value = 1

loop : integer, parameter

Initial Value = 1

imax : integer, parameter

Initial Value = 1

aunit1 : integer, parameter

Initial Value = 66

qsat : real*8

avol : real*8

aso4 : real(8)

ano3 : real(8)

anh4 : real(8)

ah2o : real(8)

gnh3 : real(8)

ghno3 : real(8)

dust : real(8)

salt : real(8)

tk : real(8)

rh : real(8)

rhd : real(8)

rhc : real(8)

yi : real(4), dimension(IMAX, NCA)

yo : real(4), dimension(IMAX, NCO)

mw_anh4 : real(8), parameter

Initial Value = 18.03850D+00

mw_gnh3 : real(8), parameter

Initial Value = MW_ANH4

mw_ano3 : real(8), parameter

Initial Value = 62.00494D+00

mw_ghno3 : real(8), parameter

Initial Value = MW_ANO3

mw_aso4 : real(8), parameter

Initial Value = 96.0636D+00

mw_na : real(8), parameter

Initial Value = 22.989768D+00

mw_cl : real(8), parameter

Initial Value = 35.4527D+00

mw_nacl : real(8), parameter

Initial Value = 58.442468D+00

mw_h2o : real(8), parameter

Initial Value = 18.01528D+00

mw_k : real(4), parameter

Initial Value = 39.0983

mw_ca : real(4), parameter

Initial Value = 40.078

mw_mg : real(4), parameter

Initial Value = 24.3050

mass_frac_k : real(4), parameter

Initial Value = 0.0028

mass_frac_ca : real(4), parameter

Initial Value = 0.024

mass_frac_mg : real(4), parameter

Initial Value = 0.0038

mass_frac_na : real(4), parameter

Initial Value = 0.014

frac_dust : real(4), parameter

Initial Value = 1.d0

frac_salt : real(4), parameter

Initial Value = 0.001

conv_kion : real(4), parameter

Initial Value = FRAC_DUST * MASS_FRAC_K / MW_K

conv_caion : real(4), parameter

Initial Value = FRAC_DUST * MASS_FRAC_CA / MW_CA

conv_mgion : real(4), parameter

Initial Value = FRAC_DUST * MASS_FRAC_MG / MW_MG

conv_naion : real(4), parameter

Initial Value = FRAC_DUST * MASS_FRAC_NA / MW_NA

rat_na : real(8), parameter

Initial Value = MW_NA / ( MW_NA + MW_CL )

rat_cl : real(8), parameter

Initial Value = MW_CL / ( MW_NA + MW_CL )

rmw_gnh3 : real(4), parameter

Initial Value = 1.d0 / MW_GNH3

rmw_anh4 : real(4), parameter

Initial Value = 1.d0 / MW_ANH4

rmw_ghno3 : real(4), parameter

Initial Value = 1.d0 / MW_GHNO3

rmw_ano3 : real(4), parameter

Initial Value = 1.d0 / MW_ANO3

rmw_aso4 : real(4), parameter

Initial Value = 1.d0 / MW_ASO4

rmw_nacl : real(4), parameter

Initial Value = 1.d0 / MW_NACL

rmw_na : real(4), parameter

Initial Value = 1.d0 / MW_NA

rmw_cl : real(8), parameter

Initial Value = 1.d0 / MW_CL

rhmax : real(8), parameter

Initial Value = 0.995D+00

rhmin : real(8), parameter

Initial Value = 0.010D+00

small_so4 : real(8), parameter

Initial Value = 1.0D-05

h : real(8)