! $Id: emissions_mod.f,v 1.20 2007/11/05 16:16:17 bmy Exp $
MODULE EMISSIONS_MOD 1
!
!******************************************************************************
! Module EMISSIONS_MOD is used to call the proper emissions subroutine
! for the various GEOS-CHEM simulations. (bmy, 2/11/03, 9/18/07)
!
! Module Routines:
! ============================================================================
! (1 ) DO_EMISSIONS : Driver which calls various emissions routines
!
! GEOS-CHEM modules referenced by emissions_mod.f
! ============================================================================
! (1 ) bravo_mod.f :
! (2 ) c2h6_mod.f : Module w/ routines for C2H6 chemistry
! (3 ) carbon_mod.f : Module w/ routines for carbon arsl emissions
! (4 ) ch3i_mod.f : Module w/ routines for CH3I chemistry
! (5 ) co2_mod.f : Module w/ routines for CO2 chemistry
! (6 ) dust_mod.f : Module w/ routines for dust aerosol emissions
! (7 ) emep_mod.f : Module w/ routines to read EMEP (Europe) emissions
! (8 ) epa_nei_mod.f : Module w/ routines to read EPA/NEI99 (USA) emissions
! (9 ) error_mod.f : Module w/ NaN and other error checks
! (10) global_ch4_mod.f : Module w/ routines for CH4 emissions
! (11) h2_hd_mod.f : Module w/ routines for H2 and HD chemistry
! (12) hcn_ch3cn_mod.f : Module w/ routines for HCN and CH3CN emissions
! (13) Kr85_mod.f : Module w/ routines for Kr85 emissions
! (14) logical_mod.f : Module w/ GEOS-CHEM logical switches
! (15) mercury_mod.f : Module w/ routines for mercury chemistry
! (16) RnPbBe_mod.f : Module w/ routines for Rn-Pb-Be emissions
! (17) tagged_co_mod.f : Module w/ routines for Tagged CO emissions
! (18) time_mod.f : Module w/ routines to compute date & time
! (19) tracer_mod.f : Module w/ GEOS-CHEM tracer array STT etc.
! (20) seasalt_mod.f : Module w/ routines for seasalt emissions
! (21) sulfate_mod.f : Module w/ routines for sulfate emissions
!
! NOTES:
! (1 ) Now references DEBUG_MSG from "error_mod.f"
! (2 ) Now references "Kr85_mod.f" (jsw, bmy, 8/20/03)
! (3 ) Now references "carbon_mod.f" and "dust_mod.f" (rjp, tdf, bmy, 4/2/04)
! (4 ) Now references "seasalt_mod.f" (rjp, bmy, bec, 4/20/04)
! (5 ) Now references "logical_mod" & "tracer_mod.f" (bmy, 7/20/04)
! (6 ) Now references "epa_nei_mod.f" and "time_mod.f" (bmy, 11/5/04)
! (7 ) Now references "emissions_mod.f" (bmy, 12/7/04)
! (8 ) Now calls EMISSSULFATE if LCRYST=T. Also read EPA/NEI emissions for
! the offline aerosol simulation. (bmy, 1/11/05)
! (9 ) Remove code for the obsolete CO-OH param simulation (bmy, 6/24/05)
! (10) Now references "co2_mod.f" (pns, bmy, 7/25/05)
! (11) Now references "emep_mod.f" (bdf, bmy, 10/1/05)
! (12) Now references "gfed2_biomass_mod.f" (bmy, 3/30/06)
! (13) Now references "bravo_mod.f" (rjp, kfb, bmy, 6/26/06)
! (14) Now references "edgar_mod.f" (avd, bmy, 7/6/06)
! (15) Now references "streets_anthro_mod.f" (yxw, bmy, 8/18/06)
! (16) Now references "h2_hd_mod.f" (lyj, phs, 9/18/07)
!******************************************************************************
!
IMPLICIT NONE
!=================================================================
! MODULE ROUTINES -- follow below the "CONTAINS" statement
!=================================================================
CONTAINS
!------------------------------------------------------------------------------
SUBROUTINE DO_EMISSIONS 1,81
!
!******************************************************************************
! Subroutine DO_EMISSIONS is the driver routine which calls the appropriate
! emissions subroutine for the various GEOS-CHEM simulations.
! (bmy, 2/11/03, 8/17/06)
!
! NOTES:
! (1 ) Now references DEBUG_MSG from "error_mod.f" (bmy, 8/7/03)
! (2 ) Now calls Kr85 emissions if NSRCX == 12 (jsw, bmy, 8/20/03)
! (3 ) Now calls EMISSCARBON and EMISSDUST for carbon aerosol and dust
! aerosol chemistry (rjp, tdf, bmy, 4/2/04)
! (4 ) Now calls EMISSSEASALT for seasalt emissions (rjp, bec, bmy, 4/20/04)
! (5 ) Now use inquiry functions from "tracer_mod.f". Now references
! "logical_mod.f" (bmy, 7/20/04)
! (6 ) Now references ITS_A_NEW_MONTH from "time_mod.f". Now references
! EMISS_EPA_NEI from "epa_nei_mod.f" (bmy, 11/5/04)
! (7 ) Now calls EMISSMERCURY from "mercury_mod.f" (eck, bmy, 12/7/04)
! (8 ) Now calls EMISSSULFATE if LCRYST=T. Also read EPA/NEI emissions for
! the offline sulfate simulation. Also call EMISS_EPA_NEI for the
! tagged CO simulation. (cas, bmy, stu, 1/10/05).
! (9 ) Now call EMISSSEASALT before EMISSSULFATE (bec, bmy, 4/13/05)
! (10) Now call EMISS_HCN_CH3CN from "hcn_ch3cn_mod.f". Also remove all
! references to the obsolete CO-OH param simulation. (xyp, bmy, 6/23/05)
! (11) Now call EMISSCO2 from "co2_mod.f" (pns, bmy, 7/25/05)
! (12) Now references EMISS_EMEP from "emep_mod.f" (bdf, bmy, 11/1/05)
! (13) Now call GFED2_COMPUTE_BIOMASS to read 1x1 biomass emissions and
! regrid to the model resolution once per month. (bmy, 3/30/06)
! (14) Now references EMISS_BRAVO from "bravo_mod.f" (rjp, kfb, bmy, 6/26/06)
! (15) Now references EMISS_EDGAR from "edgar_mod.f" (avd, bmy, 7/6/06)
! (16) Now references EMISS_STREETS_ANTHRO from "streets_anthro_mod.f"
! (yxw, bmy, 8/17/06)
!******************************************************************************
!
! References to F90 modules
USE BIOMASS_MOD
, ONLY : NBIOMAX
USE BIOMASS_MOD, ONLY : COMPUTE_BIOMASS_EMISSIONS
USE BRAVO_MOD, ONLY : EMISS_BRAVO
USE C2H6_MOD, ONLY : EMISSC2H6
USE CARBON_MOD, ONLY : EMISSCARBON
USE CH3I_MOD, ONLY : EMISSCH3I
USE CO2_MOD, ONLY : EMISSCO2
USE DUST_MOD, ONLY : EMISSDUST
USE EDGAR_MOD, ONLY : EMISS_EDGAR
USE EMEP_MOD, ONLY : EMISS_EMEP
USE EPA_NEI_MOD, ONLY : EMISS_EPA_NEI
USE ERROR_MOD, ONLY : DEBUG_MSG
USE GLOBAL_CH4_MOD, ONLY : EMISSCH4
USE H2_HD_MOD, ONLY : EMISS_H2_HD
USE HCN_CH3CN_MOD, ONLY : EMISS_HCN_CH3CN
USE Kr85_MOD, ONLY : EMISSKr85
USE LOGICAL_MOD
USE MERCURY_MOD, ONLY : EMISSMERCURY
USE RnPbBe_MOD, ONLY : EMISSRnPbBe
USE SEASALT_MOD, ONLY : EMISSSEASALT
USE STREETS_ANTHRO_MOD, ONLY : EMISS_STREETS_ANTHRO
USE SULFATE_MOD, ONLY : EMISSSULFATE
USE TIME_MOD, ONLY : GET_MONTH, GET_YEAR
USE TIME_MOD, ONLY : ITS_A_NEW_MONTH, ITS_A_NEW_YEAR
USE TRACER_MOD
USE TAGGED_CO_MOD, ONLY : EMISS_TAGGED_CO
# include "CMN_SIZE" ! Size parameters
! Local Variables
INTEGER :: MONTH, YEAR
REAL*8 :: BIOMASS(IIPAR,JJPAR,NBIOMAX)
!=================================================================
! DO_EMISSIONS begins here!
!=================================================================
! Get year and month
YEAR = GET_YEAR()
MONTH = GET_MONTH()
! Get biomass burning emissions for use below
IF ( LBIOMASS ) THEN
CALL COMPUTE_BIOMASS_EMISSIONS( YEAR, MONTH )
ENDIF
! Test by simulation type
IF ( ITS_A_FULLCHEM_SIM() ) THEN
!--------------------
! NOx-Ox-HC-aerosol
!--------------------
! Read David Streets' emisisons over China / SE ASia
IF ( LSTREETS .and. ITS_A_NEW_YEAR() ) THEN
CALL EMISS_STREETS_ANTHRO
ENDIF
! Read EDGAR emissions once per month
IF ( LEDGAR .and. ITS_A_NEW_MONTH() ) THEN
CALL EMISS_EDGAR( YEAR, MONTH )
ENDIF
! Read EPA/NEI99 (USA) emissions once per month
IF ( LNEI99 .and. ITS_A_NEW_MONTH() ) CALL EMISS_EPA_NEI
! Read BRAVO (Mexico) emissions once per year
IF ( LBRAVO .and. ITS_A_NEW_YEAR() ) CALL EMISS_BRAVO
! Read EMEP (Europe) emissions once per year
IF ( LEMEP .and. ITS_A_NEW_YEAR() ) CALL EMISS_EMEP
! NOx-Ox-HC (w/ or w/o aerosols)
CALL EMISSDR
! Emissions for various aerosol types
IF ( LSSALT ) CALL EMISSSEASALT
IF ( LSULF .or. LCRYST ) CALL EMISSSULFATE
IF ( LCARB ) CALL EMISSCARBON
IF ( LDUST ) CALL EMISSDUST
ELSE IF ( ITS_AN_AEROSOL_SIM() ) THEN
!--------------------
! Offline aerosol
!--------------------
! Read David Streets' emisisons over China / SE ASia
IF ( LSTREETS .and. ITS_A_NEW_YEAR() ) THEN
CALL EMISS_STREETS_ANTHRO
ENDIF
! Read EDGAR emissions once per month
IF ( LEDGAR .and. ITS_A_NEW_MONTH() ) THEN
CALL EMISS_EDGAR( YEAR, MONTH )
ENDIF
! Read EPA/NEI99 emissions once per month
IF ( LNEI99 .and. ITS_A_NEW_MONTH() ) CALL EMISS_EPA_NEI
! Read BRAVO (Mexico) emissions once per year
IF ( LBRAVO .and. ITS_A_NEW_YEAR() ) CALL EMISS_BRAVO
! Read EMEP (Europe) emissions once per year
IF ( LEMEP .and. ITS_A_NEW_YEAR() ) CALL EMISS_EMEP
! Emissions for various aerosol types
IF ( LSSALT ) CALL EMISSSEASALT
IF ( LSULF .or. LCRYST ) CALL EMISSSULFATE
IF ( LCARB ) CALL EMISSCARBON
IF ( LDUST ) CALL EMISSDUST
ELSE IF ( ITS_A_RnPbBe_SIM() ) THEN
!--------------------
! Rn-Pb-Be
!--------------------
CALL EMISSRnPbBe
ELSE IF ( ITS_A_CH3I_SIM() ) THEN
!--------------------
! CH3I
!--------------------
! Emit CH3I
CALL EMISSCH3I
ELSE IF ( ITS_A_HCN_SIM() ) THEN
!--------------------
! HCN - CH3CN
!--------------------
CALL EMISS_HCN_CH3CN( N_TRACERS, STT )
ELSE IF ( ITS_A_TAGCO_SIM() ) THEN
!--------------------
! Tagged CO
!--------------------
! Read David Streets' emisisons over China / SE ASia
IF ( LSTREETS .and. ITS_A_NEW_YEAR() ) THEN
CALL EMISS_STREETS_ANTHRO
ENDIF
! Read EDGAR emissions once per month
IF ( LEDGAR .and. ITS_A_NEW_MONTH() ) THEN
CALL EMISS_EDGAR( YEAR, MONTH )
ENDIF
! Read EPA (USA) emissions once per month
IF ( LNEI99 .and. ITS_A_NEW_MONTH() ) CALL EMISS_EPA_NEI
! Read BRAVO (Mexico) emissions once per year
IF ( LBRAVO .and. ITS_A_NEW_YEAR() ) CALL EMISS_BRAVO
! Read EPA (Europe) emissions once per year
IF ( LEMEP .and. ITS_A_NEW_YEAR() ) CALL EMISS_EMEP
! Emit tagged CO
CALL EMISS_TAGGED_CO
ELSE IF ( ITS_A_C2H6_SIM() ) THEN
!--------------------
! C2H6
!--------------------
! Emit C2H6
CALL EMISSC2H6
ELSE IF ( ITS_A_CH4_SIM() ) THEN
!--------------------
! CH4
!--------------------
! Read David Streets' emisisons over China / SE ASia
IF ( LSTREETS .and. ITS_A_NEW_YEAR() ) THEN
CALL EMISS_STREETS_ANTHRO
ENDIF
! Emit CH4
CALL EMISSCH4
ELSE IF ( ITS_A_MERCURY_SIM() ) THEN
!--------------------
! Mercury
!--------------------
CALL EMISSMERCURY
ELSE IF ( ITS_A_CO2_SIM() ) THEN
!--------------------
! CO2
!--------------------
! Read David Streets' emisisons over China / SE ASia
IF ( LSTREETS .and. ITS_A_NEW_YEAR() ) THEN
CALL EMISS_STREETS_ANTHRO
ENDIF
! Emit CO2
CALL EMISSCO2
ELSE IF ( ITS_A_H2HD_SIM() ) THEN
!--------------------
! Offline H2/HD
!--------------------
! Read David Streets' emisisons over China / SE ASia
IF ( LSTREETS .and. ITS_A_NEW_YEAR() ) THEN
CALL EMISS_STREETS_ANTHRO
ENDIF
! Read EDGAR emissions once per month
IF ( LEDGAR .and. ITS_A_NEW_MONTH() ) THEN
CALL EMISS_EDGAR( YEAR, MONTH )
ENDIF
! Read EPA (USA) emissions once per month
IF ( LNEI99 .and. ITS_A_NEW_MONTH() ) CALL EMISS_EPA_NEI
! Read BRAVO (Mexico) emissions once per year
IF ( LBRAVO .and. ITS_A_NEW_YEAR() ) CALL EMISS_BRAVO
! Read EPA (Europe) emissions once per year
IF ( LEMEP .and. ITS_A_NEW_YEAR() ) CALL EMISS_EMEP
! Emit H2/HD
CALL EMISS_H2_HD
ENDIF
!### Debug
IF ( LPRT ) CALL DEBUG_MSG ( '### DO_EMISSIONS: a EMISSIONS' )
! Return to calling program
END SUBROUTINE DO_EMISSIONS
!------------------------------------------------------------------------------
! End of module
END MODULE EMISSIONS_MOD