Difference between revisions of "Advection scheme TPCORE"

From Geos-chem
Jump to: navigation, search
(Overview)
(Overview)
Line 2: Line 2:
  
 
== Overview ==
 
== Overview ==
Several versions of the TPCORE algorithm are used in GEOS-Chem. There are two versions of the algorithm, the version originally implemented in GEOS-Chem and a version based on GMI model. Then, there are modified versions for the nested grid simulations. The modifications for the nested grids do not change the algorithm itself but only the definition of the region where transport occurs. See the table to know which version of TPCORE is used in different simulations.
 
 
 
{| border=1 cellspacing=0 cellpadding=5
 
|- bgcolor="#CCCCCC"
 
! Simulation type
 
! TPCORE version
 
|-align="center"
 
| GEOS-3
 
| original TPCORE with M. Prather's pressure-fixer
 
|-align="center"
 
| GEOS-4, GEOS-5 and GCAP<br>before v8-01-03
 
| modified original TPCORE with P. Cameron-Smith's pressure-fixer<br>(shouldn't be used, see below)
 
|-align="center"
 
| GEOS-4, GEOS-5 and GCAP<br>after v8-01-03
 
| based on GMI model with P. Cameron-Smith's pressure-fixer
 
|-align="center"
 
| nested 1x1
 
| modified original TPCORE with M. Prather's pressure-fixer
 
|-align="center"
 
| nested 0.5x0.666
 
| modified original TPCORE with P. Cameron-Smith's pressure-fixer<br>(do not use for polar regions)
 
|}
 
 
 
 
There are a couple of different versions of TPCORE currently implemented into GEOS-Chem.
 
There are a couple of different versions of TPCORE currently implemented into GEOS-Chem.
  
 
;TPCORE v7.1: This version is used with the GEOS-3 met fields.  The source code is contained in F90 module <tt>tpcore_mod.f</tt>.
 
;TPCORE v7.1: This version is used with the GEOS-3 met fields.  The source code is contained in F90 module <tt>tpcore_mod.f</tt>.
;GMI TPCORE: This version is used with the GEOS-4 and GEOS-5 met fields.  The source code is contained in the F90 module <tt>tpcore_fvdas_mod.f90</tt>.  This version of TPCORE was taken from the GMI source code.
+
;GMI TPCORE: This version is used with the GEOS-4 and GEOS-5 met fields and for GCAP simulations.  The source code is contained in the F90 module <tt>tpcore_fvdas_mod.f90</tt>.  This version of TPCORE was taken from the GMI source code.
;TPCORE for GEOS-3 nested grids: This version contains the various modifications for the GEOS-3 nested grids (i.e. how to deal with boundary conditions, etc.).  The source code is in F90 module <tt>tpcore_window_mod.f</tt>.  This version was based on TPCORE v7.1 and was modified by Yuxuan Wang..  
+
;TPCORE for GEOS-3 nested grids (1x1): This version contains the various modifications for the GEOS-3 nested grids (i.e. how to deal with boundary conditions, etc.).  The source code is in F90 module <tt>tpcore_window_mod.f</tt>.  This version was based on TPCORE v7.1 and was modified by Yuxuan Wang..  
;TPCORE for GEOS-5 nested grids: This version contains the various modifications for the GEOS-5 nested grids (i.e. how to deal with boundary conditions, etc.).  The source code is in F90 module <tt>tpcore_geos5_window_mod.f90</tt>.  This version was based on a previous implementation of TPCORE by S-J Lin and Kevin Yeh (but was NOT the GMI implementation), and was modified by Yuxuan Wang and Dan Chen.
+
;TPCORE for GEOS-5 nested grids (0.5x0.666): This version contains the various modifications for the GEOS-5 nested grids (i.e. how to deal with boundary conditions, etc.).  The source code is in F90 module <tt>tpcore_geos5_window_mod.f90</tt>.  This version was based on a previous implementation of TPCORE by S-J Lin and Kevin Yeh (but was NOT the GMI implementation), and was modified by Yuxuan Wang and Dan Chen. It shouldn't be used over polar regions.
 +
 
  
 
=== TPCORE replaced with version based on GMI model ===
 
=== TPCORE replaced with version based on GMI model ===

Revision as of 17:37, 17 February 2010

NOTE: The page is in construction.

Overview

There are a couple of different versions of TPCORE currently implemented into GEOS-Chem.

TPCORE v7.1
This version is used with the GEOS-3 met fields. The source code is contained in F90 module tpcore_mod.f.
GMI TPCORE
This version is used with the GEOS-4 and GEOS-5 met fields and for GCAP simulations. The source code is contained in the F90 module tpcore_fvdas_mod.f90. This version of TPCORE was taken from the GMI source code.
TPCORE for GEOS-3 nested grids (1x1)
This version contains the various modifications for the GEOS-3 nested grids (i.e. how to deal with boundary conditions, etc.). The source code is in F90 module tpcore_window_mod.f. This version was based on TPCORE v7.1 and was modified by Yuxuan Wang..
TPCORE for GEOS-5 nested grids (0.5x0.666)
This version contains the various modifications for the GEOS-5 nested grids (i.e. how to deal with boundary conditions, etc.). The source code is in F90 module tpcore_geos5_window_mod.f90. This version was based on a previous implementation of TPCORE by S-J Lin and Kevin Yeh (but was NOT the GMI implementation), and was modified by Yuxuan Wang and Dan Chen. It shouldn't be used over polar regions.


TPCORE replaced with version based on GMI model

NOTE: The GMI version of TPCORE was implemented into GEOS-Chem v8-01-03.

Dylan Jones and Hongyu Liu both independently found that the existing TPCORE transport code used to perform the advection for GEOS-4 and GEOS-5 meteorology (tpcore_fvdas_mod.f90, by S-J Lin and Kevin Yeh) causes an overshooting in the polar stratopsheric regions. Claire Carouge has supplied a fix for this issue.

Hongyu Liu wrote:

[Here are] the radon plots for all GEOS series as well as the info on GEOS-Chem version and met field version. This is much less of a problem (or not a problem at all) in GEOS-Chem/GEOS-STRAT and GEOS-Chem/GEOS-3. Actually GEOS-Chem/GEOS-STRAT is very close to GMI/GEOS-STRAT (not shown). Compare plots below. All these simulations use same options for tpcore (IORD=3, JORD=3, KORD=7).
It seems that GMI does not have this problem. You might want to compare the plots on Pages 7 & 10 of this PDF file.

Dylan Jones wrote:

Attached are the PowerPoint slides. We added slide 9 which compares the mixed layer depth for GEOS-Chem (with GEOS-5) and GMI (with GEOS-4); we see the same problem in GEOS-Chem with GEOS-5, but we do not see it in GMI.
All the GEOS-4 GEOS-Chem runs were with version v7-02-04. The GEOS-5 GEOS-Chem runs were with v8-01-01.

Claire Carouge replied:

I think I finally got tpcore to work. You have some plots attached for Radon. So far, I ran 5 months of simulation and each month looks much better: no polar spike and much much lower stratospheric tracer. The concentration of Radon is never null in the stratosphere but I don't think it is in GMI neither, I think the range on the plot was cut off. So I did the plot with the same cut off. Let me know if you see something strange I haven't seen.
I don't know for sure what was the problem but I think it was a problem with the definition of the polar cap. In GEOS-Chem, we only averaged values for cells at the poles on 1 band of latitude. In GMI, they use what they call an enlarged polar cap and thus they average values at the poles on 2 bands of latitude.
The pressure fixer was taken from GMI and thus was written for an enlarged polar cap. The parts of the code that were explicitly labelled for the enlarged polar cap had been removed but not other parts that were not explicitly labelled. So we ended up with a hybrid pressure fixer used with a tpcore with a not enlarged olar cap. I'm guessing this was creating the polar spikes.
But, the tpcore in GEOS-Chem and GMI are not easy to compare. I'm not sure there was not a difference in the transport itself and I can't be totally sure the pressure fixer was well introduced in tpcore. I have some doubts about this.
So as GMI was an example of a working (and clean!) algorithm with the enlarged polar cap, I decided to introduce the enlarged polar cap in GEOS-Chem. For this:
  1. I slightly modified the pressure-fixer to return to the exact version from GMI
  2. I changed tpcore to the one from GMI as I was sure the pressure fixer was well introduced in (and the algorithm is cleaner).

Dylan Jones wrote:

Claire,
Attached are some comparisons of the CO/O3 correlations in the UTLS with the new tpcore.
The files called *.gmi_geos.png are with the old version of tpcore and those called *.gmi_geos4_tpcore are with the corrected tpcore. As you can see, there is much less scatter in the correlations, which suggest less anomalous mixing with the new code. The scatter is more similar to what we see in GMI. There is still an offset in the stratosphere (for high ozone in the plots), but I think this may be due to linoz rather than to the transport scheme. This is a huge improvement for GEOS-Chem. Thanks very much for fixing it.

Hongyu Liu replied:

Thanks Claire! Well done. It's clear that the polar overshooting problem is now fixed.

Also see this PDF document for more information about the new TPCORE version as installed in GEOS-Chem v8-01-03.

References

  1. Lin, S.-J., and R. B. Rood, 1996: Multidimensional flux form semi-Lagrangian transport schemes, Mon. Wea. Rev., 124, 2046-2070.
  2. Lin, S.-J., W. C. Chao, Y. C. Sud, and G. K. Walker, 1994: A class of the van Leer-type transport schemes and its applications to the moisture transport in a General Circulation Model, Mon. Wea. Rev., 122, 1575-1593.

--Bob Y. 09:49, 17 February 2010 (EST)