Difference between revisions of "GEOS-Chem Adjoint v35"

From Geos-chem
Jump to: navigation, search
(Implemented new ship emissions scheme)
(Overview)
Line 67: Line 67:
 
|Forward Model Bug Fix
 
|Forward Model Bug Fix
 
|Christoph Keller
 
|Christoph Keller
 +
|Yanko Davila (CU Boulder)
 +
|v35a
 +
|Approved 30 Oct 2013
 +
|-bgcolor="#CCFFFF"
 +
 +
|-
 +
|[[#Bug_in_ship_CO_emissions| Bug in ship CO emissions ]]
 +
|Forward Model Bug Fix
 +
|Chris holmes
 +
|Yanko Davila (CU Boulder)
 +
|v35a
 +
|Approved 30 Oct 2013
 +
|-bgcolor="#CCFFFF"
 +
 +
|-
 +
|[[#Bugs_in_ship_emissions_and_non-local_PBL_mixing| Bugs in ship emissions and non-local PBL mixing ]]
 +
|Forward Model Bug Fix
 +
|Chris holmes
 
|Yanko Davila (CU Boulder)
 
|Yanko Davila (CU Boulder)
 
|v35a
 
|v35a
Line 80: Line 98:
 
|Approved 30 Oct 2013
 
|Approved 30 Oct 2013
 
|-bgcolor="#CCFFFF"
 
|-bgcolor="#CCFFFF"
 +
  
 
|-
 
|-
Line 117: Line 136:
  
 
--[[User:Yanko|Yanko]] 15:28, 24 July 2013 (EDT)
 
--[[User:Yanko|Yanko]] 15:28, 24 July 2013 (EDT)
 +
 +
===== Bug in ship CO emissions =====
 +
 +
EMEP ship emissions inventory was previously overwriting all ship CO emissions globally.  [[Ship_emissions#Bugs_in_ship_emissions_and_non-local_PBL_mixing | See forward wiki. ]]
 +
 +
--[[User:Yanko|Yanko]] 19:02, 30 October 2013 (EDT)
 +
 +
===== Bugs in ship emissions and non-local PBL mixing =====
 +
 +
Ship NOx emissions with the ICOADS-AMVER inventory were 3.3 times too low, due to a NO2 vs. N molecular weight error. Emissions should be 5.0Tg(N)/yr. As a result, O3 and HNO3 concentrations were underestimated, especially over the N. Atlantic and N. Pacific. [[Ship_emissions#Bugs_in_ship_emissions_and_non-local_PBL_mixing | See forward wiki. ]]
 +
 +
--[[User:Yanko|Yanko]] 19:02, 30 October 2013 (EDT)
  
 
===== Bug fix for 3-hourly GFED3 emissions =====
 
===== Bug fix for 3-hourly GFED3 emissions =====

Revision as of 23:02, 30 October 2013

Overview

BETA RELEASE

What's new in this version

GEOS-Chem Adjoint v35 contains the following major updates and improvements.

Feature Type Submitted by Committed by Version Status
Add support for NEI2008 Emissions Forward Model Update Daven Henze (CU Boulder) Yanko Davila (CU Boulder) v35 Approved 24 Jul 2013
Fixed a bug in Henry's constant Forward Model Bug Fix Fabien Paulot (Harvard) Yanko Davila (CU Boulder) v35 Approved 24 Jul 2013
Prevent negative emissions over Canada Forward Model Bug Fix Fabien Paulot (Harvard), Christoph Keller Yanko Davila (CU Boulder) v35 Approved 24 Jul 2013
Add support for new deposition based cost functions Adjoint Model Update Fabien Paulot (Harvard), Daven Henze (CU Boulder) Yanko Davila (CU Boulder) v35 Approved 24 Jul 2013
Bug fix for 3-hourly GFED3 emissions Forward Model Bug Fix Christoph Keller Yanko Davila (CU Boulder) v35a Approved 30 Oct 2013
Bug fix in grid_mod.F90 Forward Model Bug Fix Christoph Keller Yanko Davila (CU Boulder) v35a Approved 30 Oct 2013
Bug in ship CO emissions Forward Model Bug Fix Chris holmes Yanko Davila (CU Boulder) v35a Approved 30 Oct 2013
Bugs in ship emissions and non-local PBL mixing Forward Model Bug Fix Chris holmes Yanko Davila (CU Boulder) v35a Approved 30 Oct 2013
Eliminate several more array temporaries Forward Model Bug Fix GEOS-Chem Support Team Yanko Davila (CU Boulder) v35a Approved 30 Oct 2013
Implemented new ship emissions scheme Adjoint Model Update Daven Henze (CU Boulder), Yanko Davila (CU Boulder) Yanko Davila (CU Boulder) v35a Approved 30 Oct 2013


- Updates to forward model

Add NEI2008 Emissions

Implemented NEI2008 Emissions Inventory

See benchmark results gcadj_std_v35 tests 003 and 023

--Yanko 15:28, 24 July 2013 (EDT)


- Bug fixes in forward model

Fixed a bug in Henry's constant

Small error in the definition of the Henry's constant, which is used in the wet deposition routine. See forward wiki.

--Yanko 15:28, 24 July 2013 (EDT)

Prevent negative emissions over Canada

Found a nasty bug in CAC emissions. We should set the emissions to 0 when 2002 exceeds 2005 by more than a factor of three. See forward wiki.

--Yanko 15:28, 24 July 2013 (EDT)

Bug in ship CO emissions

EMEP ship emissions inventory was previously overwriting all ship CO emissions globally. See forward wiki.

--Yanko 19:02, 30 October 2013 (EDT)

Bugs in ship emissions and non-local PBL mixing

Ship NOx emissions with the ICOADS-AMVER inventory were 3.3 times too low, due to a NO2 vs. N molecular weight error. Emissions should be 5.0Tg(N)/yr. As a result, O3 and HNO3 concentrations were underestimated, especially over the N. Atlantic and N. Pacific. See forward wiki.

--Yanko 19:02, 30 October 2013 (EDT)

Bug fix for 3-hourly GFED3 emissions

Running with 3-hourly emissions from the first day of a month seems to work fine, but trying to start on any other day appears to cause a problem. See forward wiki.

--Yanko 14:38, 30 October 2013 (EDT)

Bug fix in grid_mod.F90

The sine of the last latitude edge (i.e. the north pole) is not calculated, which can cause problems for the regridding. See forward wiki.

--Yanko 14:38, 30 October 2013 (EDT)

Eliminate several more array temporaries

In many areas of GEOS-Chem, we pass arrays as arguments from one routine to another. But if this is not done properly, it can cause GEOS-Chem to use an excessive amount of memory and take longer to run, especially at very fine resolutions. See forward wiki.

--Yanko 14:38, 30 October 2013 (EDT)


- Updates to adjoint model

Add support for new deposition based cost functions

Added support for new deposition based cost functions.

See benchmark results gcadj_std_v35 tests 012 and 023

--Yanko 15:28, 24 July 2013 (EDT)

Implemented new ship emissions scheme

Implemented new ship emissions scheme matching forward model paranox module. Changes to the adjoint include, but are not limited to writing the adjoint of paranox and updating setemis_adj module.

Ship_emissions#PARANOX_ship_plume_model

Ship_emissions#PARANOX_ship_plume_model

See benchmark results gcadj_std_v35a tests

--Yanko 14:38, 30 October 2013 (EDT)

- Bug fixes in adjoint model

Small bug fix in wetscav_adj_mod.f

Fixed a small bug related to the ALPHA2 variable on WETDEP_ADJ subroutine in wetscav_adj_mod.f

--Yanko 15:28, 24 July 2013 (EDT)