GAMAP bugs and fixes: Difference between revisions
(fix to combine ND49 timeseries) |
(update on using PS-TOP in gamap) |
||
| Line 1: | Line 1: | ||
== Using PS-TOP in GAMAP with GEOS-Chem v8+ outputs == | |||
GEOS-Chem does not write the category PS-TOP into the tracerinfo.dat anymore (see [[Outstanding_issues_yet_to_be_resolved#Pressure_diagnostic_:_PS-PTOP_vs._PEDGE|PS-PTOP vs. PEDGE]]). However some routines still try to get the pressure surface with that category name. You can simply replace PS-TOP by PEDGE-$ in your program for them to work. To be able to use older GEOS-Chem output and new ones, you may use a coding similar to the one I did for | |||
'''ctm_column_du.pro''', which you can find [http://www.as.harvard.edu/ctm/geos/wiki_docs/ctm_column_du.pro here]. | |||
Note that to read pressure surface from time series diagnostic (ND48-51), you need to keep PS-TOP in your IDL routines and replace PEDGE-$ by PS-TOP in the tracerinfo.dat files created by GEOS-Chem. This is a bug, and we expect to get rid of all references to PS-TOP in a future release of GAMAP. | |||
--[[User:Plesager|phs]] 16:25, 2 September 2008 (EDT) | |||
== Bug in GC_COMBINE_ND49 == | == Bug in GC_COMBINE_ND49 == | ||
a minor bug could give incorrect results when either /DMAX or /DAVG was used, and timeseries from more than one month were processed at once. The output OUTTIME vector (if requested) was not correct, and if a smaller TIME span than the original series was requested, you could end up with a wrong series. | a minor bug could give incorrect results when either /DMAX or /DAVG was used, and timeseries from more than one month were processed at once. The output OUTTIME vector (if requested) was not correct, and if a smaller TIME span than the original series was requested, you could end up with a wrong series. | ||
Revision as of 20:25, 2 September 2008
Using PS-TOP in GAMAP with GEOS-Chem v8+ outputs
GEOS-Chem does not write the category PS-TOP into the tracerinfo.dat anymore (see PS-PTOP vs. PEDGE). However some routines still try to get the pressure surface with that category name. You can simply replace PS-TOP by PEDGE-$ in your program for them to work. To be able to use older GEOS-Chem output and new ones, you may use a coding similar to the one I did for ctm_column_du.pro, which you can find here.
Note that to read pressure surface from time series diagnostic (ND48-51), you need to keep PS-TOP in your IDL routines and replace PEDGE-$ by PS-TOP in the tracerinfo.dat files created by GEOS-Chem. This is a bug, and we expect to get rid of all references to PS-TOP in a future release of GAMAP.
--phs 16:25, 2 September 2008 (EDT)
Bug in GC_COMBINE_ND49
a minor bug could give incorrect results when either /DMAX or /DAVG was used, and timeseries from more than one month were processed at once. The output OUTTIME vector (if requested) was not correct, and if a smaller TIME span than the original series was requested, you could end up with a wrong series. The fix will be available in the next GAMAP release (v2.13). Meanwhile you can get it here.
--phs 14:55, 15 August 2008 (EDT)