| 123456789101112131415161718192021222324252627282930313233343536373839404142434445464748495051525354555657585960616263646566676869707172737475767778798081828384858687888990919293949596979899100101102103104105106107108109110111112113114115116117118119120121122123124125126127128129130131132133134135136137138139140141142143144145146147148149150151152153154155156157158159160161162163164165166167168169170171172173174175176177178179180181182183184185186187188189190191192193194195196197198199200201202203204205206207208209210211212213214215216217218219220221222223224225226227228229 |
- !### macro's #####################################################
- !
- #define TRACEBACK write (gol,'("in ",a," (",a,", line",i5,")")') rname, __FILE__, __LINE__; call goErr
- #define IF_NOTOK_RETURN(action) if (status/=0) then; TRACEBACK; action; return; end if
- #define IF_ERROR_RETURN(action) if (status> 0) then; TRACEBACK; action; return; end if
- !
- #include "tm5.inc"
- !
- !#################################################################
- !------------------------------------------------------------------------------
- ! TM5 !
- !------------------------------------------------------------------------------
- !BOP
- !
- ! !MODULE: deposition_data
- !
- ! !DESCRIPTION: Dry Deposition Constants
- !\\
- !\\
- ! !INTERFACE:
- !
- module deposition_data
- !
- ! !USES:
- !
- use GO, only : gol, goPr, goErr
- implicit none
- public
- !
- ! !PUBLIC DATA MEMBERS:
- !
- !
- ! -- definition of terms which correct for the diffusivity
- ! for the computation of the boundary layer resistance, (v/D)**2/3
- ! v = 0.189 cm2 s-1 (heat) and D from DH2O/DX where DH2O=0.212
- !
- real, parameter :: diffrb_o3 = 1.2
- real, parameter :: diffrb_hno3 = 1.4
- real, parameter :: diffrb_no = 1.1
- real, parameter :: diffrb_no2 = 1.2
- real, parameter :: diffrb_so2 = 1.4
- real, parameter :: diffrb_so4 = 999.
- real, parameter :: diffrb_h2o2 = 1.2
- real, parameter :: diffrb_ald2 = 1.4
- real, parameter :: diffrb_mald = 1.4
- real, parameter :: diffrb_ch2o = 1.1
- real, parameter :: diffrb_ch3o2h = 1.3
- real, parameter :: diffrb_nh3 = 0.9
- real, parameter :: diffrb_pan = 1.7
- real, parameter :: diffrb_co = 1.2
- real, parameter :: diffrb_ch3coo2h = 1.5
- real, parameter :: diffrb_hcooh = 1.3
- real, parameter :: diffrb_hno2 = 1.3
- real, parameter :: diffrb_hno4 = 0.0 ! <--- dummy
- !
- ! -- parallel soil resistance, the resistances of the trace gases
- ! of the Wesely scheme are calculated assuming a O3 soil resistance
- ! of 400 s m-1, for SO2 rsoil=100 s m-1 and the wet skin resistance of
- ! O3 is 2000 s m-1 whereas rws SO2 is 100 s m-1.
- !
- real, parameter :: rsoil_o3 = 400.0
- real, parameter :: rsoil_hno3 = 1.0
- real, parameter :: rsoil_no = 1.0e+10
- real, parameter :: rsoil_no2 = 600.0
- real, parameter :: rsoil_so2 = 100.0
- real, parameter :: rsoil_so4 = 999.
- real, parameter :: rsoil_h2o2 = 80.0
- real, parameter :: rsoil_mald = 1.e+10
- real, parameter :: rsoil_ald2 = 1.0e+10
- real, parameter :: rsoil_ch2o = 1666.0
- real, parameter :: rsoil_ch3o2h = 3650.0
- real, parameter :: rsoil_nh3 = 100.0
- real, parameter :: rsoil_pan = 3994.0
- real, parameter :: rsoil_co = 5.0e+4
- real, parameter :: rsoil_ch3coo2h = 3290.0
- real, parameter :: rsoil_hcooh = 1.0
- real, parameter :: rsoil_hno2 = 97.0
- real, parameter :: rsoil_hno4 = 0.0 ! <--- dummy
- !
- ! -- sea water resistance, which is generally similar to the wet skin
- ! resistance except of that for SO2 due to a different pH of sea water
- ! compared to the pH of a wet canopy
- ! note for ammonia ocean emissions we use 2-way exchange approach FD
- !
- real, parameter :: rwat_o3 = 2000.0
- real, parameter :: rwat_hno3 = 1.0
- real, parameter :: rwat_no = 1.0e+10
- real, parameter :: rwat_no2 = 3000.0
- real, parameter :: rwat_so2 = 1.0
- real, parameter :: rwat_so4 = 999.0
- real, parameter :: rwat_h2o2 = 72.0
- real, parameter :: rwat_mald = 300.0
- real, parameter :: rwat_ald2 = 300.0
- real, parameter :: rwat_ch2o = 254.0
- real, parameter :: rwat_ch3o2h = 293.0
- real, parameter :: rwat_nh3 = 1.0
- real, parameter :: rwat_pan = 295.0
- real, parameter :: rwat_co = 1.0e+10
- real, parameter :: rwat_ch3coo2h = 291.0
- real, parameter :: rwat_hcooh = 1.0
- real, parameter :: rwat_hno2 = 75.0
- real, parameter :: rwat_hno4 = 0.0 ! <--- dummy
- !
- ! -- wet skin reservoir resistance
- !
- real, parameter :: rws_o3 = 2000.0
- real, parameter :: rws_hno3 = 1.0
- real, parameter :: rws_no = 1.0e+10
- real, parameter :: rws_no2 = 3000.0
- real, parameter :: rws_so2 = 100.0
- real, parameter :: rws_so4 = 999.0
- real, parameter :: rws_h2o2 = 72.0
- real, parameter :: rws_mald = 300.0
- real, parameter :: rws_ald2 = 300.0
- real, parameter :: rws_ch2o = 254.0
- real, parameter :: rws_ch3o2h = 293.0
- real, parameter :: rws_nh3 = 1.0
- real, parameter :: rws_pan = 295.0
- real, parameter :: rws_co = 1.0e+10
- real, parameter :: rws_ch3coo2h = 291.0
- real, parameter :: rws_hno2 = 75.0
- real, parameter :: rws_hcooh = 1.0
- real, parameter :: rws_hno4 = 0.0 ! <--- dummy
- !
- ! -- snow resistance, the snow/ice resistances of the trace gases
- ! of Wesely's scheme are taken same as that of the rsoil and corrected
- ! for temperatures smaller than 271 K
- !
- real, parameter :: rsnow_o3 = 2000.0
- real, parameter :: rsnow_hno3 = 1.0
- real, parameter :: rsnow_no = 1.0e+10
- real, parameter :: rsnow_no2 = 3000.0
- real, parameter :: rsnow_so2 = 1.0
- real, parameter :: rsnow_so4 = 999.0
- real, parameter :: rsnow_h2o2 = 80.0
- real, parameter :: rsnow_mald = 1.0e+10
- real, parameter :: rsnow_ald2 = 1.0e+10
- real, parameter :: rsnow_ch2o = 1666.0
- real, parameter :: rsnow_ch3o2h = 3650.0
- real, parameter :: rsnow_nh3 = 1.0e10
- real, parameter :: rsnow_pan = 3394.0
- real, parameter :: rsnow_co = 1.0e10
- real, parameter :: rsnow_ch3coo2h = 3290.0
- real, parameter :: rsnow_hcooh = 1.0e+10
- real, parameter :: rsnow_hno2 = 97.0
- real, parameter :: rsnow_hno4 = 0.0 ! <--- dummy
- !
- ! -- mesophyll resistance
- !
- real, parameter :: rmes_o3 = 1.0
- real, parameter :: rmes_hno3 = 1.0
- real, parameter :: rmes_no = 500.0
- real, parameter :: rmes_no2 = 1.0
- real, parameter :: rmes_so2 = 1.0
- real, parameter :: rmes_so4 = 999.0
- real, parameter :: rmes_h2o2 = 1.0
- real, parameter :: rmes_mald = 200.0
- real, parameter :: rmes_ald2 = 200.0
- real, parameter :: rmes_ch2o = 1.0
- real, parameter :: rmes_ch3o2h = 1.0
- real, parameter :: rmes_nh3 = 1.0
- real, parameter :: rmes_pan = 1.0
- real, parameter :: rmes_co = 1.0e10
- real, parameter :: rmes_ch3coo2h = 1.0
- real, parameter :: rmes_hcooh = 1.0
- real, parameter :: rmes_hno2 = 1.0
- real, parameter :: rmes_hno4 = 0.0 ! <--- dummy
- !
- ! -- cuticle resistance
- !
- real, parameter :: rcut_o3 = 1.0e+10
- real, parameter :: rcut_hno3 = 1.0
- real, parameter :: rcut_no = 1.0e+10
- real, parameter :: rcut_no2 = 1.0e+10
- real, parameter :: rcut_so2 = 1.0e+10
- real, parameter :: rcut_so4 = 999.
- real, parameter :: rcut_h2o2 = 1.0e+10
- real, parameter :: rcut_mald = 1.0e+10
- real, parameter :: rcut_ald2 = 1.0e+10
- real, parameter :: rcut_ch2o = 1.0e+10
- real, parameter :: rcut_ch3o2h = 1.0e+10
- real, parameter :: rcut_nh3 = 1.0e+10
- real, parameter :: rcut_pan = 1.0e+10
- real, parameter :: rcut_co = 1.0e+10
- real, parameter :: rcut_ch3coo2h = 1.0e+10
- real, parameter :: rcut_hcooh = 2500.0
- real, parameter :: rcut_hno2 = 1.0e+10
- real, parameter :: rcut_hno4 = 0.0 ! <--- dummy
- !
- ! -- Diffusivity coefficent, to correct stomatal resistance for differences
- ! in diffusivity between water vapour and the specific trace gas
- !
- real, parameter :: diffcf_o3 = 1.6
- real, parameter :: diffcf_hno3 = 1.9
- real, parameter :: diffcf_no = 1.3
- real, parameter :: diffcf_no2 = 1.6
- real, parameter :: diffcf_so2 = 1.9
- real, parameter :: diffcf_so4 = 999.
- real, parameter :: diffcf_h2o2 = 1.4
- real, parameter :: diffcf_mald = 1.6
- real, parameter :: diffcf_ald2 = 1.6
- real, parameter :: diffcf_ch2o = 1.3
- real, parameter :: diffcf_ch3o2h = 1.6
- real, parameter :: diffcf_nh3 = 1.0
- real, parameter :: diffcf_pan = 2.6
- real, parameter :: diffcf_co = 1.2
- real, parameter :: diffcf_ch3coo2h = 2.0
- real, parameter :: diffcf_hcooh = 1.6
- real, parameter :: diffcf_hno2 = 1.6
- real, parameter :: diffcf_hno4 = 0.0 ! <--- dummy
- !
- ! !REVISION HISTORY:
- ! 25 Mar 2010 - P. Le Sager - reintroduced SO4 parameters, added
- ! protex tags.
- ! 7 Feb 2013 - N Banda - changed high resistance values - from 1e+5 to 1e+10
- ! 3 Sep 2013 - P. Le Sager - declare all as real instead of integer
- !
- !
- ! !REMARKS:
- ! o extracted from dry_deposition.F90 codes;
- ! o dummy values for HNO4, actual deposition velocity is computed from:
- ! vd(HNO4) = ( vd(NO2) + vd(H2O2) )/2
- !
- !EOP
- !------------------------------------------------------------------------------
- end module deposition_data
|
