Difference between revisions of "Tagged O3 simulation"

From Geos-chem
Jump to: navigation, search
(Now let user select the year of archived ozone data)
(Recent modifications to the tagged O3 simulation)
Line 58: Line 58:
 
--[[User:Bmy|Bob Y.]] 11:38, 11 April 2014 (EDT)
 
--[[User:Bmy|Bob Y.]] 11:38, 11 April 2014 (EDT)
  
== Recent modifications to the tagged O3 simulation ==
+
== Modifications to the tagged O3 simulation ==
  
 
=== Write initial and final Ox mass to file ===
 
=== Write initial and final Ox mass to file ===
Line 101: Line 101:
  
 
--[[User:Bmy|Bob Y.]] 17:04, 30 May 2014 (EDT)<br>--[[User:Bmy|Bob Yantosca]] ([[User talk:Bmy|talk]]) 17:51, 8 January 2016 (UTC)
 
--[[User:Bmy|Bob Y.]] 17:04, 30 May 2014 (EDT)<br>--[[User:Bmy|Bob Yantosca]] ([[User talk:Bmy|talk]]) 17:51, 8 January 2016 (UTC)
 
=== Reactivate stratospheric P(Ox) in tagged Ox simulation ===
 
 
<span style="color:green">'''''These issues were resolved during the [[GEOS-Chem_v9-02#Period_of_public_comment|public comment period]] following the provisional release of GEOS-Chem v9-02 and were included in the public release of [[GEOS-Chem v9-02]] (03 Mar 2014).  This update is included in Adjoint [[GEOS-Chem_Adjoint_v35 | v35j]].'''''</span>
 
 
[mailto:jlin5@pku.edu.cn Jintai Lin] noted that the stratospheric P(Ox) produced by the Synoz boundary condition should be passed to the Tagged Ox simulation source code in <tt>tagged_ox_mod.F</tt>.  This is only for runs which use the Synoz boundary condition.  If you are using Linoz, then this should be disabled.
 
 
The updates are as follows:
 
 
(1) In routine <tt>ADD_STRAT_POx</tt> (in <tt>tagged_ox_mod.F</tt>), we need to uncomment this line:
 
 
      STT(I,J,L,N_STRAT) = STT(I,J,L,N_STRAT) + POx
 
 
This will restore the production of Ox in the stratosphere.
 
 
(2) In routine <tt>UPBDFLX_O3</tt> (in module <tt>upbdflx_mod.F</tt>), we must first reference the <tt>LLINOZ</tt> flag from <tt>logical_mod.f</tt>. 
 
 
      USE LOGICAL_MOD, ONLY : LLINOZ
 
 
*If <tt>LLINOZ</tt> is TRUE, then we are using Linoz linear stratospheric chemistry scheme. 
 
*If <tt>LLINOZ</tt> is FALSE, then we default to the Synoz stratopsheric boundary condition.
 
 
(3) Also in <tt>UPBDFLX_O3</tt>, we need to define <tt>J30S</tt> and <tt>J30N</tt> parameters for nested grid regions.  This will allow Tagged Ox simulations for the nested grid regions to be conducted.
 
 
#elif defined( GRID05x0666 )
 
      !INTEGER, PARAMETER  :: J30S = 1, J30N = JJPAR
 
! Lin_20111026
 
#if  defined( NESTED_CH )
 
      INTEGER, PARAMETER  :: J30S = 1,  J30N = 83
 
#elif  defined( NESTED_NA )
 
      INTEGER, PARAMETER  :: J30S = 1,  J30N = 41
 
#elif  defined( NESTED_EU )
 
      INTEGER, PARAMETER  :: J30S = 1,  J30N = 1        ! add later
 
#else
 
      INTEGER, PARAMETER  :: J30S = 1, J30N = JJPAR
 
#endif
 
 
(4) Further down in routine <tt>UPBDFLX_O3</tt>, we must call <tt>ADD_STRAT_POx</tt>.  This will pass the P(Ox) produced by Synoz to the source code in <tt>tagged_ox_mod.F</tt>.  This call is placed immediately after updating the <tt>STT</tt> tracer array:
 
 
            ! Store O3 flux in the proper tracer number
 
            STT(I,J,L,NTRACER) = STT(I,J,L,NTRACER) + PO3
 
 
            ! Store O3 flux for strat Ox tracer (Tagged Ox only)
 
            IF ( .not. LLINOZ ) THEN
 
                IF ( ITS_A_TAGOX_SIM() ) CALL ADD_STRAT_POX( I, J, L, PO3 )
 
            ENDIF
 
 
NOTE: The call to <tt>ADD_STRAT_POx</tt> will be ignored if we are using the Linoz linear stratospheric O3 chemistry scheme.
 
 
--[[User:Bmy|Bob Y.]] 16:06, 1 November 2011 (EDT)<br>--[[User:Bmy|Bob Yantosca]] ([[User talk:Bmy|talk]]) 17:55, 8 January 2016 (UTC)
 
  
  
 
----
 
----
 
'''''[[Tagged CO simulation|Previous]] | [[TransportTracers simulation|Next]] | [[Guide to GEOS-Chem simulations]]'''''
 
'''''[[Tagged CO simulation|Previous]] | [[TransportTracers simulation|Next]] | [[Guide to GEOS-Chem simulations]]'''''

Revision as of 20:28, 19 September 2022

Previous | Next | Guide to GEOS-Chem simulations

  1. Simulations overview
  2. Mechanisms for full-chemistry simulations (e.g. Standard, Tropchem, etc.)
  3. [[Aerosol-only simulation]
  4. CH4 simulation
  5. CO2 simulation
  6. Hg simulation
  7. POPs simulation
  8. Tagged CO simulation
  9. Tagged O3 simulation
  10. TransportTracers simulation


This page contains information about the Tagged O3 simulation in GEOS-Chem (formerly known as the "Tagged Ox" simulation).

Overview

The tagged O3 simulation allows you to use archived ozone production and loss rates to perform a simulation for geographically tagged ozone tracers without having to run the full-chemistry simulation. The tagged O3 simulation is also used to spin up the ozone into steady-state when validating a new meterorolgical field product.

Nomenclature

We now prefer to refer to this simulation as the "Tagged O3 Simulation" instead of the "Tagged Ox Simulation". This reflects the removal of the Ox tracer from GEOS-Chem v9-02 and higher versions. Ox has now been replaced by O3.

--Bob Y. 11:08, 11 April 2014 (EDT)

Types of tagged O3 simulations

You may perform two different types of tagged O3 simulations. These simulations differ in the number of tagged O3 tracers that are available:

Simple tagged O3 simulation

The default setting is the SIMPLE TAGGED O3 SIMULATION, which carries 2 tracers:

(1 ) O3      : Total O3
(2 ) O3Strt  : O3 from the Stratosphere

Extended tagged O3 simulation

You may also choose to run the EXTENDED TAGGED O3 SIMULATION, which carries 13 tracers:

(1 ) O3      : Total O3
(2 ) O3Strt  : O3 from the Stratosphere      (tropopause - atm top   )
(3 ) O3Ut    : O3 produced in Upper Trop     (350 hPa    - tropopause) 
(4 ) O3Mt    : O3 produced in Middle Trop    (PBL top    - 350 hPa   )
(5 ) O3Row   : O3 produced in Rest of World  (surface    - PBL top   )
(6 ) O3PcBl  : O3 produced in Pacific BL     (surface    - PBL top   )
(7 ) O3NaBl  : O3 produced in N. American BL (surface    - PBL top   )
(8 ) O3AtBl  : O3 produced in Atlantic BL    (surface    - PBL top   )
(9 ) O3EuBl  : O3 produced in European BL    (surface    - PBL top   )
(10) O3AfBl  : O3 produced in N. African BL  (surface    - PBL top   )
(11) O3AsBl  : O3 produced in Asian          (surface    - PBL top   )
(12) O3Init  : O3 initial conditions         (all levels             )
(13) O3USA   : O3 produced over the USA      (all levels             )

Note that in both the simple and extended simulations, the stratospheric O3 tracer is listed as tracer #2. This is necessary for compatibility with the Linoz stratospheric ozone chemistry mechanism.

--Bob Y. 11:38, 11 April 2014 (EDT)

Modifications to the tagged O3 simulation

Write initial and final Ox mass to file

This update was included in v11-02b and approved on 16 Jun 2017.

The budget computation code needs the initial and final mass of O3 in [kg] in order to compute the difference in accumulated mass. The accumulated mass is then compared to the budget computed with all terms (emissions, deposition, chemical, etc) for validation. Although in theory we can get the initial and final masses from the restart file and convert units, it is more straightforward to print the values to a file.

--Lizzie Lundgren (talk) 19:36, 18 May 2017 (UTC)

PBLTOP and MTTOP bug fix

This update was validated in the 1-month benchmark simulation v11-01d (approved 12 Dec 2015).

PBLTOP and MTTOP are set to levels corresponding to 750 hPa and 350 hPa respectively in routine get_regional_pox in tagged_ox_mod.F. Values used for GEOS-4 are erroneously used for GEOS-5, GEOS-FP, MERRA, and MERRA-2, resulting in PBL and middle troposphere levels that are below what they should be for those MET fields. This bug fix updates the levels for GEOS-5, GEOS-FP, MERRA, and MERRA-2 from 5 to 16 for PBLTOP and from 10 to 27 for MTTOP.

--Lizzie Lundgren (talk) 14:58, 22 October 2015 (UTC)
--Bob Yantosca (talk) 17:49, 8 January 2016 (UTC)

European grid definition bug fix

This update was validated in the 1-month benchmark simulation v11-01d (approved 12 Dec 2015).

There is an error in the European grid definition within tagged_ox_mod.F in v10-01. Please see the tagged ozone simulation wiki page for details.The longitude definition of ITS_A_EUR is incorrectly defined as greater than -15.0 and greater than or equal to 55.0, and thus has no upper bound. It should instead have an upper bound of 55.0. If running the extended tagged ozone simulation, you should amend the definition to be as follows:

ITS_IN_EUR = ( Y >=  36.0 .and. ( X >  -15.0 .and. X <=   55.0 ) )

--Lizzie Lundgren (talk) 14:12, 28 September 2015 (UTC)
--Bob Yantosca (talk) 17:50, 8 January 2016 (UTC)

Make it easier to switch between the simple and extended simulations

This update was validated in the 1-month benchmark simulation v10-01c (approved on 29 May 2014).

If you wish to activate the EXTENDED TAGGED O3 SIMULATION, then look for this line of code at the top of module GeosCore/tagged_ox_mod.F:

!#define USE_ALL_TAGO3_TRACERS 1

and remove the ! comment character, so that the line now looks like this:

#define USE_ALL_TAGO3_TRACERS 1

Then make realclean and recompile GEOS-Chem, etc.

--Bob Y. 17:04, 30 May 2014 (EDT)
--Bob Yantosca (talk) 17:51, 8 January 2016 (UTC)



Previous | Next | Guide to GEOS-Chem simulations