Difference between revisions of "Sandbox"

From Geos-chem
Jump to: navigation, search
(GEOS-Chem Working Groups)
 
(169 intermediate revisions by 3 users not shown)
Line 1: Line 1:
Dear GEOS-Chem Users,
+
== Instructions for adding tracers to GEOS-Chem ==
  
Here is the June 2009 GEOS-Chem Newsletter.
+
'''''By [mailto:mpayer@seas.harvard.edu Melissa Payer]'''''
  
''Bob Yantosca, Philippe Le Sager, and Claire Carouge''<br>
+
These instructions describe how to add advected tracers to GEOS-Chem.  Advected tracers are carried in the <tt>STT</tt> array (which is declared in <tt>GeosCore/tracer_mod.F</tt>), and are transported by the winds.
''geos-chem-support@as.harvard.edu''
+
  
== GEOS-Chem Update ==
+
NOTE: Some advected tracers are also chemical species, which comprise GEOS-Chem's NOx-Ox-HC-aerosol chemistry mechanism.
  
=== GEOS-Chem v8-02-02 ===
+
#<p>Add tracers to input.geos under Tracer Menu</p>
 +
##<p>Increase “Number of tracers” accordingly</p>
 +
##<p>Add tracer number, name, molecular weight, and emitted species to end of tracer list</p>
 +
#<p>Add tracer(s) to restart file – see separate instructions</p>
 +
#<p>Modify the following source code routines:</p>
 +
##<p><tt>Headers/CMN_SIZE_mod.F</tt></p>
 +
###<p>Increase <tt>NNPAR</tt> by number of tracers being added</p>
 +
###<p>If the species have anthropogenic emissions, increase NEMPARA accordingly</p>
 +
###<p>If the species have biogenic emissions, increase NEMPARB accordingly</p>
 +
###<p>If the species have biomass burning emissions, increase NBIOMAX accordingly</p>
 +
##<p><tt>GeosCore/tracerid_mod.F</tt></p>
 +
###<p>Increase <tt>NNNTRID</tt> by the number of tracers being added NOTE: <tt>NNNTRID</tt> must be equal to <tt>NNPAR</tt> in <tt>CMN_SIZE_mod.F</tt></p>
 +
###<p>Define tracer IDs as needed:</p>
 +
####<p>IDxxxx – ID for SMVGEAR species</p>
 +
####<p>IDTxxxx – GEOS-Chem tracer ID</p>
 +
####<p>IDExxxx – GEOS-Chem emission ID</p>
 +
####<p>IDBFxxxx – GEOS-Chem biofuel ID</p>
 +
####<p>IDBxxxx – GEOS-Chem biomass ID</p>
 +
###<p>In subroutine <tt>TRACERID</tt>,</p>
 +
####<p>Under “Assign tracer, biomass, biofuel, and anthro emission ID’s,” add case for added tracers to initialize IDTxxxx, IDBFxxxx, IDBxxxx</p>
 +
####<p>Under “Initialize the IDExxx flags,” add IF statement for added tracer IDs and initialize IDExxxx</p>
 +
####<p>Under “Fill IDEMS with appropriate tracer ID #’s,” add IF statement for added tracers</p>
 +
####<p>Under “Print additional information,” add WRITE statements for added tracers</p>
 +
###<p>In subroutine SETTRACE  (NOTE: Only do the following if species is included in SMVGEAR chemical mechanism<p>
 +
####<p>Add IF statement, to initialize IDxxxx</p>
 +
####<p>In subroutine INIT_TRACERID,</p>
 +
####<p>Zero all IDs added (IDxxxx, IDTxxxx, IDExxxx, etc.)</p>
 +
##<p>Other modules that may need to be updated include:</p>
 +
###<p>wetscav_mod.F</p>
 +
###<p>drydep_mod.F</p>
 +
###<p>emissions_mod.F</p>
 +
###<p>biomass_mod.F</p>
 +
###<p>biofuel_mod.F</p> etc.
  
[[GEOS-Chem v8-02-02]] "beta" has now been released.  This version includes:
 
  
*new boundary layer mixing scheme (developed by Jintai Lin)
 
*updated Bond emission inventories for BC & OC (prepared by Eric Leibensperger). 
 
  
This version may be downloaded from:
+
Instructions for adding species to GEOS-Chem
  
  ftp ftp.as.harvard.edu
+
Modify run directory files
  cd pub/geos-chem/beta_releases
+
  get GEOS-Chem.v8-02-02-beta.stdrun.tar.gz
+
  
'''''NOTE: We currently working on enabling access to our CVS repository to outside users so that you will be able to obtain the GEOS-Chem source code without having to download a large TAR file.  Stay tuned...'''
+
1) Add tracers to input.geos under Tracer Menu
 +
c) Increase “Number of tracers” accordingly
 +
d) Add tracer number, name, molecular weight, and emitted species to end of tracer list
  
=== GEOS-Chem Steering Committee ===
+
2) Add tracer(s) to restart file – see separate instructions
  
The [http://www.as.harvard.edu/ctm/geos/geos_steering_cmte.html GEOS-Chem Steering Committee] held a telecon on June 11, 2009. Click [http://www.as.harvard.edu/ctm/geos/word_pdf_docs/steering_cmte/gc_sc_minutes_061109.pdf HERE] to view the meeting minutes.
+
3) Modify globchem.dat to include additional species, kinetic reactions (including emissions and drydep), or photolysis reactions
 +
NOTE: If you modify globchem.dat and plan to run GEOS-Chem using the KPP chemical solver, you will need to generate new gckpp*.F90 files. For more information see the following wiki pages:
 +
http://wiki.seas.harvard.edu/geos-chem/index.php/KPP_solvers_FAQ#What_are_the_cons_of_using_KPP.3F
 +
http://wiki.seas.harvard.edu/geos-chem/index.php/Interfacing_GEOS-Chem_with_KPP#Generating_KPP_input_files_from_GEOS-Chem_globchem.dat
  
Since last April's GEOS-Chem meeting, a number of [[GEOS-Chem model development priorities]] that were identified by the Steering Committee and Working Groups have now been implemented into GEOS-Chem:
+
4) If necessary, modify:
 
+
a) ratj.d – contains species names and branching ratios for FAST–JX photolysis species
{| border=1 cellspacing=0 cellpadding=5
+
For more information: http://acmg.seas.harvard.edu/geos/doc/man/chapter_5.html#5.4.1
|- bgcolor="#cccccc"
+
b) jv_spec.dat – contains cross-sections and quantum yields for FAST–JX photolysis species
!Item
+
For more information: http://acmg.seas.harvard.edu/geos/doc/man/chapter_5.html#5.4.3
!Author(s)
+
!Completion Date
+
|-
+
|Improved BL mixing scheme
+
|Jintai Lin
+
|Now included in [[GEOS-Chem v8-02-02]];<br>released 08 Jun 2009
+
|-
+
|BC/OC aerosol emissions from Bond inventory
+
|Eric Leibensperger
+
|Now included in [[GEOS-Chem v8-02-02]];<br>released 08 Jun 2009
+
|-
+
|Chemical updates
+
|Mat Evans<br>Fabien Paulot<br>Jingqiu Mao<br>Paul Palmer
+
|Updated chemical mechanism and photolysis files<br><tt>globchem.dat</tt>, <tt>jv_spec.dat</tt>, <tt>ratj.d</tt><br>(see [http://www.as.harvard.edu/ctm/geos/wiki_docs/chemistry/chemistry_updates_v4.pdf documentation])<br>Now included in [[GEOS-Chem v8-02-01]];<br>released 26 May 2009
+
|}
+
 
+
At the June 10, 2009 telecon, the GEOS-Chem Steering Committee identified the following remaining [[GEOS-Chem model development priorities|model development items]] as top priorities (here listed in order):
+
 
+
{| border=1 cellspacing=0 cellpadding=5
+
|- bgcolor="#cccccc"
+
!#
+
!Development priority
+
!Point person(s)
+
!Status
+
|-
+
|1.
+
|GEOS-5 ozone columns
+
|GC support team
+
|Validation is ongoing
+
|-
+
|2
+
|KPP implementation
+
|GC support team<br>Virginia Tech
+
|In beta testing
+
|-
+
|3.
+
|TOMAS implementation
+
|Win Trivitayanurak<br>GC Support Team
+
|Telecon to discuss implementation<br>scheduled for June 16th
+
|-
+
|4.
+
|Linoz
+
|Dylan Jones and<br>U. Toronto group
+
|In beta testing
+
|-
+
|5.
+
|Updated methane simulation
+
|Christopher Pickett-Heaps<br>Kevin Wecht
+
|Debugging and code cleanup is ongoing
+
|-
+
|6.
+
|Updated CO2 Simulation
+
|Ray Nassar
+
|&nbsp;
+
|-
+
|7.
+
|Updated MEGAN emission from terpenes
+
|Dylan Millet
+
|&nbsp;
+
|-
+
|8.
+
|COADS ship emissions
+
|Chulkyu Lee
+
|&nbsp;
+
|-
+
|9.
+
|Include updated aerosol optics 
+
|Colette Heald<br>Randall Martin<br>and aerosols Working Group
+
|&nbsp;
+
|-
+
|10.
+
|Include aromatic SOA
+
|Daven Henze
+
|&nbsp;
+
|-
+
|11.
+
|Archive OH fields from benchmark<br>for use in off-line simulations
+
|GC support team
+
|&nbsp}
+
|}
+
 
+
=== GEOS-Chem Working Groups ===
+
 
+
A couple of the [http://www.as.harvard.edu/ctm/geos/geos_steering_cmte.html GEOS-Chem Working Groups] have been renamed:
+
* "Carbon Gases" and "Organics" is now "Carbon Gases & Organics" (co-chairs Dylan Jones, Dylan Millet)
+
* "Regional Modeling" is now "Regional Air Quality" (co-chairs Rokjin Park, Yuxuan Wang)
+
 
+
We ask that each GEOS-Chem user join the [http://www.as.harvard.edu/ctm/geos/geos_working_groups.html Working Group] that is most relevant to his/her area of research.  The Working Groups are meant to foster communication and information sharing between GEOS-Chem users, as well as to identify [[GEOS-Chem_model_development_priorities|priorities for model development]] to the [http://www.as.harvard.edu/ctm/geos/geos_steering_cmte.html GEOS-Chem Steering Committee].
+
 
+
For more information about a particular working group, please contact the Working Group Chairperson(s).
+
 
+
=== Subscribing to the GEOS-Chem email lists ===
+
 
+
There are now several GEOS-Chem email lists:
+
 
+
{| border=1 cellspacing=0 cellpadding=5
+
|- bgcolor="#cccccc"
+
!Email list address
+
!Intended audience
+
|-
+
|<tt>geos-chem@seas.harvard.edu</tt>
+
|All GEOS-Chem users and developers<br>'''''NOTE: This is the email list to which information about <br>new model versions, bugs & fixes, etc. will be sent'''''
+
|-
+
|<tt>geos-chem-adjoint@seas.harvard.edu</tt>
+
|Members of the [[GEOS-Chem Adjoint|Adjoint Model and Data Assimilation Working Group]]
+
|-
+
|<tt>geos-chem-aerosols@seas.harvard.edu</tt>
+
|Members of the [[Aerosols_overview:_projects_and_developments|Aerosols Working Group]]
+
|-
+
||<tt>geos-chem-carbon@seas.harvard.edu</tt>
+
|Members of the [[Carbon Gases and Organics Working Group]]
+
|-
+
|<tt>geos-chem-emissions@seas.harvard.edu</tt>
+
|Members of the [[Emissions Working Group]]
+
|-
+
|<tt>geos-chem-oxidants@seas.harvard.edu</tt>
+
|Members of the [[Chemistry Issues|Oxidants and Chemistry Working Group]]
+
|-
+
|<tt>geos-chem-regional@seas.harvard.edu</tt>
+
|Members of the [[Regional Air Quality Working Group]]
+
|}
+
 
+
To subscribe to any of these lists, please visit the page: https://lists.seas.harvard.edu/mailman/listinfo and then click on the name of the list to which you wish to subscribe.  You will be taken to a page where you can specify your email address and a password and then click SUBMIT.  You will then receive a confirmation email from the system.  Follow the directions in the confirmation email and then you will be added to the list.
+
 
+
'''''PLEASE NOTE: If you wish to send an email to any of the above GEOS-Chem mailing lists, then please be sure to send your email from the same address with which you subscribed to the list.  Otherwise your email may bounce back to you. The GEOS-Chem mailing lists are set up to reject emails from non-list members.  This is done in order to prevent spam.'''''
+
 
+
=== Being a member of the GEOS-Chem community ===
+
 
+
We would like to encourage all new GEOS-Chem users to take a moment to read our [[GEOS-Chem welcome letter for new users]]. 
+
 
+
We expect that each GEOS-Chem user will possess at least some familiarity with the following:
+
 
+
* The Fortran programming language (F77 OK, F90 preferred)
+
* IDL (preferred), or other graphics software (e.g. Splus, Matlab)
+
* The Unix environment (i.e. commands, file paths, etc.)
+
* Unix scripts (not strictly required but is very helpful to know!)
+
 
+
The GEOS-Chem support team cannot provide you with general instruction in programming languages or answer general programming questions.  We would like to invite you to view to view the list of tutorials and other relevant resources that are listed on the GEOS-Chem web page at [http://www.as.harvard.edu/ctm/geos/geos_resources.html '''GEOS-Chem Resources Page'''.]
+
 
+
We invite you to read our document on 
+
[http://www.as.harvard.edu/ctm/geos/word_pdf_docs/resources_guidelines_tips.pdf '''User Resources, Code Development Guidelines, and Debugging Tips'''].
+
 
+
Each GEOS-Chem user is also expected to adhere to our list of [http://www.as.harvard.edu/ctm/geos/geos_best_practices.html '''best practices'''].
+
 
+
We also would like to invite all GEOS-Chem users outside Harvard who have run the 1-month benchmark simulation on their own systems to update the table of '''[[GEOS-Chem supported platforms and compilers]]''' with your timing information.
+
 
+
=== Adjoint Model Update ===
+
 
+
A successful meeting was held at Virginia Tech to develop a standard version of the GC adjoint merging the capabilities developed by individual groups: participants were Kumaresh Singh, Adrian Sandu, Monika Kopacz, Changshub Shim, and Daven Henze. The result was an integrated adjoint code based on v8-02-01 with GEOS-5 and including updates from the different groups, improved file formats and input files, etc.
+
 
+
== GEOS-5 met field reprocessing ==
+
 
+
=== Status ===
+
As described in the [[GEOS-Chem Newsletter April-May 2009|April-May 2009 newsletter]], we are reprocessing the entire archive of GEOS-5 met field data.  As of May 6, 2009 the following fields have been completed:
+
 
+
* 2003: December only
+
* 2004: entire year
+
* 2005: entire year
+
* 2006: entire year
+
* 2007: thru May 2007
+
* 2008: entire year
+
 
+
for each of the following geographic domains:
+
 
+
* China nested 0.5 x 0.667
+
* Global 2 x 2.5
+
* Global 4 x 5
+
 
+
MOTE:
+
 
+
The files have been TAR'd and Gzip'd for your convienience (as this dramatically reduces the size of the data). In each directory there are files <tt>01.tar.gz</tt>, <tt>02.tar.gz</tt> .. <tt>12.tar.gz</tt>, each of which contains a single month of met field data.  GEOS-Chem users outside Harvard may obtain these data files via anonymous FTP:
+
 
+
  ftp ftp.as.harvard.edu
+
  cd pub/geos-chem/data/GEOS_0.5x0.666_CH.d/GEOS_5/YYYY/MM/
+
  mget *
+
  cd pub/geos-chem/data/GEOS_2x2.5.d/GEOS_5/YYYY/MM/
+
  mget *
+
  cd pub/geos-chem/data/GEOS_4x5.d/GEOS_5/YYYY/MM/
+
  mget *
+
 
+
where <tt>YYYY</tt> and <tt>MM</tt> are the year and month that you wish to download.
+
 
+
+
However, to simplify the FTP process you might want to use a 3rd party FTP software program with a GUI interface, such as [http://fireftp.mozdev.org/ Mozilla FireFTP].
+
 
+
We are currently working on reprocessing GEOS-5 met data for 2006; however, this effort has been slowed down by some filesystem problems at Harvard.
+
 
+
== ESMF update ==
+
 
+
Development of the GEOS-Chem column chemistry module continues to proceed swiftly.  The SMVGEAR solver has been columnized and is now in testing to ensure that it gives similar results to the existing GEOS-Chem.  Also, the dust, sea salt, and carbon aerosol chemistry routines have been columnized and are working correctly. 
+
 
+
--[[User:Bmy|Bob Y.]] 11:33, 15 June 2009 (EDT)
+
 
+
 
+
 
+
== feel free to experiment here ==
+

Latest revision as of 21:25, 10 May 2023

Instructions for adding tracers to GEOS-Chem

By Melissa Payer

These instructions describe how to add advected tracers to GEOS-Chem. Advected tracers are carried in the STT array (which is declared in GeosCore/tracer_mod.F), and are transported by the winds.

NOTE: Some advected tracers are also chemical species, which comprise GEOS-Chem's NOx-Ox-HC-aerosol chemistry mechanism.

  1. Add tracers to input.geos under Tracer Menu

    1. Increase “Number of tracers” accordingly

    2. Add tracer number, name, molecular weight, and emitted species to end of tracer list

  2. Add tracer(s) to restart file – see separate instructions

  3. Modify the following source code routines:

    1. Headers/CMN_SIZE_mod.F

      1. Increase NNPAR by number of tracers being added

      2. If the species have anthropogenic emissions, increase NEMPARA accordingly

      3. If the species have biogenic emissions, increase NEMPARB accordingly

      4. If the species have biomass burning emissions, increase NBIOMAX accordingly

    2. GeosCore/tracerid_mod.F

      1. Increase NNNTRID by the number of tracers being added NOTE: NNNTRID must be equal to NNPAR in CMN_SIZE_mod.F

      2. Define tracer IDs as needed:

        1. IDxxxx – ID for SMVGEAR species

        2. IDTxxxx – GEOS-Chem tracer ID

        3. IDExxxx – GEOS-Chem emission ID

        4. IDBFxxxx – GEOS-Chem biofuel ID

        5. IDBxxxx – GEOS-Chem biomass ID

      3. In subroutine TRACERID,

        1. Under “Assign tracer, biomass, biofuel, and anthro emission ID’s,” add case for added tracers to initialize IDTxxxx, IDBFxxxx, IDBxxxx

        2. Under “Initialize the IDExxx flags,” add IF statement for added tracer IDs and initialize IDExxxx

        3. Under “Fill IDEMS with appropriate tracer ID #’s,” add IF statement for added tracers

        4. Under “Print additional information,” add WRITE statements for added tracers

      4. In subroutine SETTRACE (NOTE: Only do the following if species is included in SMVGEAR chemical mechanism<p>

        1. <p>Add IF statement, to initialize IDxxxx

        2. In subroutine INIT_TRACERID,

        3. Zero all IDs added (IDxxxx, IDTxxxx, IDExxxx, etc.)

    3. Other modules that may need to be updated include:

      1. wetscav_mod.F

      2. drydep_mod.F

      3. emissions_mod.F

      4. biomass_mod.F

      5. biofuel_mod.F

        etc.


Instructions for adding species to GEOS-Chem

Modify run directory files

1) Add tracers to input.geos under Tracer Menu c) Increase “Number of tracers” accordingly d) Add tracer number, name, molecular weight, and emitted species to end of tracer list

2) Add tracer(s) to restart file – see separate instructions

3) Modify globchem.dat to include additional species, kinetic reactions (including emissions and drydep), or photolysis reactions NOTE: If you modify globchem.dat and plan to run GEOS-Chem using the KPP chemical solver, you will need to generate new gckpp*.F90 files. For more information see the following wiki pages: http://wiki.seas.harvard.edu/geos-chem/index.php/KPP_solvers_FAQ#What_are_the_cons_of_using_KPP.3F http://wiki.seas.harvard.edu/geos-chem/index.php/Interfacing_GEOS-Chem_with_KPP#Generating_KPP_input_files_from_GEOS-Chem_globchem.dat

4) If necessary, modify: a) ratj.d – contains species names and branching ratios for FAST–JX photolysis species For more information: http://acmg.seas.harvard.edu/geos/doc/man/chapter_5.html#5.4.1 b) jv_spec.dat – contains cross-sections and quantum yields for FAST–JX photolysis species For more information: http://acmg.seas.harvard.edu/geos/doc/man/chapter_5.html#5.4.3