|
|
(52 intermediate revisions by 3 users not shown) |
Line 1: |
Line 1: |
− | This page contains information about files produced by GEOS-Chem simulations, including diagnostic data and restart files used for initial conditions. For information the input files that ship with the GEOS-Chem run directories, please see our [[GEOS-Chem Input Files|''GEOS-Chem Input Files'' wiki page]].
| + | __FORCETOC__ |
| + | '''''[[GEOS-Chem log files|Previous]] | [[GEOS-Chem diagnostic output files|Next]] | [[Getting Started with GEOS-Chem]]''''' |
| + | #[[Minimum system requirements for GEOS-Chem|Minimum system requirements (and software installation)]] |
| + | #[[Configuring your computational environment]] |
| + | #[[Downloading GEOS-Chem source code|Downloading source code]] |
| + | #[[Downloading GEOS-Chem data directories|Downloading data directories]] |
| + | #[[Creating GEOS-Chem run directories|Creating run directories]] |
| + | #[[GEOS-Chem configuration files|Configuring runs]] |
| + | #[[Compiling GEOS-Chem|Compiling]] |
| + | #[[Running GEOS-Chem|Running]] |
| + | #[[GEOS-Chem output files|Output files]] |
| + | #*[[GEOS-Chem log files|Log files]] |
| + | #*<span style="color:blue">'''Restart files'''</span> |
| + | #*[[GEOS-Chem diagnostic output files|Diagnostic output files]] |
| + | #[[Python tools for use with GEOS-Chem]] |
| + | #[[GEOS-Chem_coding_and_debugging|Coding and debugging]] |
| + | #[[GEOS-Chem_overview#Further_reading|Further reading]] |
| | | |
− | == Restart files in GEOS-Chem 12 ==
| |
| | | |
− | <span style="color:green">'''''These updates were included in [[GEOS-Chem 12#12.1.0|GEOS-Chem 12.1.0]], which was released on 26 Nov 2018.'''''</span>
| + | This content has been migrated to the [https://geos-chem.readthedocs.io/en/latest/gcc-guide/04-data/restart-files.html '''Restart files''' chapter of <tt>geos-chem.readthedocs.io</tt>]. |
| | | |
− | In [[GEOS-Chem 12#12.1.0|GEOS-Chem 12.1.0]], several restart file updates were introduced. These include:
| |
| | | |
− | === New fields in GEOS-Chem restart file ===
| + | ----- |
− | | + | '''''[[GEOS-Chem log files|Previous]] | [[GEOS-Chem diagnostic output files|Next]] | [[Getting Started with GEOS-Chem]]''''' |
− | Several new fields have been added to the GEOS-Chem restart file in an attempt to remove differences between single and multi-segmented GEOS-Chem simulations.
| + | |
− | | + | |
− | #Module-level variables <tt>H2O2s</tt> and <tt>SO2s</tt> from <tt>wetscav_mod.F</tt> have been added to State_Chm (as <tt>State_Chm%H2O2AfterChem, StateChm%SO2AfterChem</tt>). These fields will be output to the GEOS-Chem restart file and then initialized to the values saved in the restart file at the start of the next simulation. Prior to this update, both H2O2s and SO2s were initialized to the H2O2 and SO2 tracer concentrations at the start of every simulation. '''''This change will impact multi-segmented runs only'''''.
| + | |
− | #Module-level variables <tt>DRY_TOTN</tt> and <tt>WET_TOTN</tt> from <tt>get_ndep_mod.F</tt> have been added to State_Chm (as <tt>State_Chm%DryDepNitrogen, StateChm%WetDepNitrogen</tt>). These fields will be output to the GEOS-Chem restart file and then initialized to the values saved in the restart file at the start of the next simulation. Prior to this update, both variables were initialized to zero at the start of every simulation. Storing them in the restart file may improve accuracy of soil NOx emissions over multi-segmented runs. '''''This change will impact multi-segmented runs only'''''.
| + | |
− | #Move <tt>State_PSC</tt> from the HEMCO restart file to the GEOS-Chem restart file.
| + | |
− | #Save out instantaneous met fields <tt>TMPU1, SPHU1, PS1DRY, PS1WET, DELPDRY</tt> to the GEOS-Chem restart file. These will be used to initialize the met fields at the start of the timestep, otherwise they will be set to the values of those fields at the end of the timestep.
| + | |
− | | + | |
− | While this update was added to GEOS-Chem Classic in [[GEOS-Chem_12#12.1.0|12.1.0]], it will be added to GCHP in 12.2.0.
| + | |
− | | + | |
− | === Restart collection in History.rc ===
| + | |
− | | + | |
− | GEOS-Chem restart files are now saved out via the history component. A new [[List_of_diagnostics_archived_to_netCDF_format#The_Restart_collection|Restart collection]] has been defined in <tt>HISTORY.rc</tt> and fields saved out to the restart file can be modified in that file. | + | |
− | | + | |
− | === Read restart file via HEMCO ===
| + | |
− | | + | |
− | GEOS-Chem restart files are now read in via HEMCO. The entries listed below have been added to <tt>HEMCO_Config.rc</tt> (and may vary slightly for different simulation types). These fields are obtained from HEMCO and copied to the appropriate <tt>State_Chm</tt> and <tt>State_Met</tt> fields in the new routine <tt>Get_GC_Restart</tt> (found in <tt>GeosCore/hcoi_gc_main_mod.F90</tt>).
| + | |
− | | + | |
− | #==============================================================================
| + | |
− | # --- GEOS-Chem restart file ---
| + | |
− | #
| + | |
− | # PSC state only needed for UCX
| + | |
− | #==============================================================================
| + | |
− | (((GC_RESTART
| + | |
− | * SPC_ ./GEOSChem.Restart.$YYYY$MM$DD_$HH$MNz.nc4 SpeciesRst_?ALL? $YYYY/$MM/$DD/$HH CS xyz 1 * - 1 1
| + | |
− | * TMPU1 ./GEOSChem.Restart.$YYYY$MM$DD_$HH$MNz.nc4 Met_TMPU1 $YYYY/$MM/$DD/$HH E xyz 1 * - 1 1
| + | |
− | * SPHU1 ./GEOSChem.Restart.$YYYY$MM$DD_$HH$MNz.nc4 Met_SPHU1 $YYYY/$MM/$DD/$HH E xyz 1 * - 1 1
| + | |
− | * PS1DRY ./GEOSChem.Restart.$YYYY$MM$DD_$HH$MNz.nc4 Met_PS1DRY $YYYY/$MM/$DD/$HH E xy 1 * - 1 1
| + | |
− | * PS1WET ./GEOSChem.Restart.$YYYY$MM$DD_$HH$MNz.nc4 Met_PS1WET $YYYY/$MM/$DD/$HH E xy 1 * - 1 1
| + | |
− | * DELPDRY ./GEOSChem.Restart.$YYYY$MM$DD_$HH$MNz.nc4 Met_DELPDRY $YYYY/$MM/$DD/$HH E xyz 1 * - 1 1
| + | |
− | * KPP_HVALUE ./GEOSChem.Restart.$YYYY$MM$DD_$HH$MNz.nc4 Chem_KPPHvalue $YYYY/$MM/$DD/$HH E xyz 1 * - 1 1
| + | |
− | * WETDEP_N ./GEOSChem.Restart.$YYYY$MM$DD_$HH$MNz.nc4 Chem_WetDepNitrogen $YYYY/$MM/$DD/$HH E xy 1 * - 1 1
| + | |
− | * DRYDEP_N ./GEOSChem.Restart.$YYYY$MM$DD_$HH$MNz.nc4 Chem_DryDepNitrogen $YYYY/$MM/$DD/$HH E xy 1 * - 1 1
| + | |
− | * SO2_AFTERCHEM ./GEOSChem.Restart.$YYYY$MM$DD_$HH$MNz.nc4 Chem_SO2AfterChem $YYYY/$MM/$DD/$HH E xyz 1 * - 1 1
| + | |
− | * H2O2_AFTERCHEM ./GEOSChem.Restart.$YYYY$MM$DD_$HH$MNz.nc4 Chem_H2O2AfterChem $YYYY/$MM/$DD/$HH E xyz 1 * - 1 1
| + | |
− | (((+STATE_PSC+
| + | |
− | * STATE_PSC ./GEOSChem.Restart.$YYYY$MM$DD_$HH$MNz.nc4 Chem_StatePSC $YYYY/$MM/$DD/$HH E xyz count * - 1 1
| + | |
− | )))+STATE_PSC+
| + | |
− | )))GC_RESTART
| + | |
− | | + | |
− | The <tt>CS</tt> cycle flag was added as an option to HEMCO in [[GEOS-Chem 12#12.1.0|GEOS-Chem 12.1.0]] to tell HEMCO to skip fields that aren't found in the provided file. This is useful when certain species aren't found in the restart file and, in that case, GEOS-Chem will initialize that species to the background concentrations specified in the species database.
| + | |
− | | + | |
− | --[[User:Melissa Payer|Melissa Sulprizio]] ([[User talk:Melissa Payer|talk]]) 16:30, 7 November 2018 (UTC)
| + | |
− | | + | |
− | == Viewing and manipulating restart files in netCDF format ==
| + | |
− | | + | |
− | There are many free and open-source software packages readily available for visualizing and manipulating netCDF files. These tools will reduce the need for the GEOS-Chem user community to rely on IDL (and GAMAP), which can be prohibitively expensive for some user groups. Some recommend tools are listed below.
| + | |
− | | + | |
− | :1. '''[http://www.unidata.ucar.edu/software/netcdf/old_docs/docs_4_1/netcdf/ncdump.html ncdump]''': This command-line tool generates a text representation of netCDF data and can be used to quickly view the variables contained in a netCDF file. For example:
| + | |
− | | + | |
− | # Display header information (dimensions, variables, attributes) only
| + | |
− | ncdump -h GEOSChem_restart.201308010000 | less
| + | |
− |
| + | |
− | # Display header information followed by data values for variable(s) specified
| + | |
− | ncdump -v SPC_NO GEOSChem_restart.201308010000 | less
| + | |
− | | + | |
− | :2. '''[http://meteora.ucsd.edu/~pierce/ncview_home_page.html ncview]''': Visualization package for netCDF files, recommended for a quick and easy look at netCDF files.
| + | |
− | | + | |
− | :3. '''[https://www.giss.nasa.gov/tools/panoply/ Panoply]:''' Data viewer for netCDF files. This package offers an alternative to ncview. From our experience, Panoply works nicely when installed on the desktop, but is slow to respond in the Linux environment.
| + | |
− | | + | |
− | :4. '''[http://nco.sourceforge.net/ NCO]''' and '''[https://code.zmaw.de/projects/cdo CDO]''': Command-line tools for manipulating and analyzing netCDF files.
| + | |
− |
| + | |
− | :5. Other programming languages may be used for viewing or manipulating netCDF files, including but not limited to:
| + | |
− | :*[https://www.ncl.ucar.edu/index.shtml NCL]
| + | |
− | :*[http://www.python.org/ Python]
| + | |
− | :*[https://www.mathworks.com/products/matlab.html Matlab]
| + | |
− | | + | |
− | Some of the tools listed above, such as ncdump and ncview, may come pre-installed on your system. Others like NCO, CDO, and Panoply may need to be installed or loaded (e.g. via the <tt>module load</tt> command). Check with your system administrator or IT staff to see what is available on your system.
| + | |
− | | + | |
− | --[[User:Melissa Payer|Melissa Sulprizio]] ([[User talk:Melissa Payer|talk]]) 15:58, 17 January 2017 (UTC)
| + | |
− | | + | |
− | === Regridding netCDF restart files ===
| + | |
− | | + | |
− | The same netCDF tools listed above can also be used to regrid netCDF restart files produced by GEOS-Chem.
| + | |
− | | + | |
− | :1. '''[https://code.zmaw.de/projects/cdo CDO]''': The Climate Data Operators include tools for regridding netCDF files. For example, the following command will apply distance-weighted regridding:
| + | |
− |
| + | |
− | cdo remapdis,gridfile infile.nc outfile.nc
| + | |
− | | + | |
− | :For <tt>gridfile</tt>, you can use the files in ftp://ftp.as.harvard.edu/gcgrid/data/ExtData/HEMCO/grids/. See the [https://www.rsmas.miami.edu/users/rajib/cdo.pdf#page=106&zoom=auto,70.8662,789.2717 Interpolation section in the CDO Guide] for more information.
| + | |
− | | + | |
− | <blockquote>'''''Bram Maasakkers wrote:'''''
| + | |
− | :I have noticed a problem regridding a 4x5 restart file to 2x2.5 using cdo 1.9.4. When I use:
| + | |
− | | + | |
− | ::<code>cdo remapdis,geos.2x25.grid GEOSChem_restart.200901010000.nc GEOSChem_restart.2009_2x25.nc</code>
| + | |
− | | + | |
− | :The last latitudinal band (-89.5) remains empty and gets filled with the standard missing value of cdo, which is really large. This leads to immediate problems in the methane simulation as enormous concentrations enter the domain from the South Pole. For now I’ve solved this problem by just using bicubic interpolation:
| + | |
− | | + | |
− | ::<code>cdo remapbic,geos.2x25.grid GEOSChem_restart.200901010000.nc GEOSChem_restart.2009_2x25_bicub.nc</code>
| + | |
− | </blockquote>
| + | |
− | | + | |
− | :2. '''[http://nco.sourceforge.net/ NCO]''': The netCDF Operators also include tools for regridding. See the [http://nco.sourceforge.net/nco.html#Regridding Regridding section of the NCO User Guide] for more information.
| + | |
− | | + | |
− | :3. Jenny Fisher provided an IDL script <tt>regridh_restart_nc.pro</tt> which is now included with [http://acmg.seas.harvard.edu/gamap/doc/ GAMAP]. Jenny wrote:
| + | |
− | | + | |
− | :<blockquote>The one caveat that I cannot figure out is something weird with the units. I have done this in such a way that it reads in whatever attributes (including units) are in the original file and writes the same units back to the regridded file. But when GEOS-Chem reads the new file it gives an “incompatible units” error — although printing them in the code shows something identical to what the case statement is searching for, as far as I can tell. If I use ncatted to overwrite the units (with identical name) they are fine. So it must be something to do with the way IDL encodes the units string, but I have no idea what. The relevant nco command is <tt>ncatted -O -a units,,m,c,"mol mol-1" INFILENAME.nc</tt></blockquote>
| + | |
− | | + | |
− | :4. Regridding routines written in [http://www.ncl.ucar.edu/ NCL] are also available in the [[NCL_tools_for_GEOS-Chem|NCL4GC package]] available on Bitbucket.
| + | |
− | | + | |
− | --[[User:Melissa Payer|Melissa Sulprizio]] ([[User talk:Melissa Payer|talk]]) 22:06, 19 January 2017 (UTC)
| + | |
− | | + | |
− | ==== Creating a netCDF restart file or adding new species to a netCDF restart file ====
| + | |
− | | + | |
− | You have a few options for adding new species to a netCDF restart file:
| + | |
− | | + | |
− | :1. In [[GEOS-Chem v11-01]], if the model cannot find a species in the netCDF restart file, it should set the initial concentration to the background concentration. The default background concentration is 1e-20 v/v. If you want to change that to a different value, you can add <tt>BackgroundVV = {VALUE}_fp</tt> to the call to <tt>Spc_Create</tt> for that species in <tt>Headers/species_database_mod.F90</tt>. To confirm your species is getting set to the background value you can turn on ND70 debug output, which will print the initial species concentrations (min/max) read from the restart file to the log file. At the end of your simulation, a new restart file will be saved out with the new species included. We recommend spinning up GEOS-Chem for an appropriate duration to achieve reasonable concentrations for the new species.
| + | |
− | | + | |
− | :2. You can use CDO and NCO to copy the restart field for one species to a new species. For example:
| + | |
− | | + | |
− | module load nco
| + | |
− | module load cdo
| + | |
− |
| + | |
− | # Extract field SPC_PMN from the original restart file
| + | |
− | cdo selvar,SPC_PMN initial_GEOSChem_rst.4x5_standard.nc GEOSChem_restart.PMN.nc
| + | |
− |
| + | |
− | # Rename selected field to SPC_NPMN
| + | |
− | ncrename -h -v SPC_PMN,SPC_NPMN GEOSChem_restart.PMN.nc GEOSChem_restart.NPMN.nc
| + | |
− |
| + | |
− | # Append new species to existing restart file
| + | |
− | ncks -h -A -M GEOSChem_restart.NPMN.nc initial_GEOSChem_rst.4x5_standard.nc
| + | |
− | | + | |
− | :3. Sal Farina wrote a simple Python script for adding a new species to the netCDF restart file:
| + | |
− | | + | |
− | #!/usr/bin/env python
| + | |
− | import netCDF4 as nc
| + | |
− | import sys
| + | |
− | import os
| + | |
− |
| + | |
− | for nam in sys.argv[1:]:
| + | |
− | f = nc.Dataset(nam,mode='a')
| + | |
− | try:
| + | |
− | o = f['SPC_OCPI']
| + | |
− | except:
| + | |
− | print "SPC_OCPI not defined"
| + | |
− |
| + | |
− | f.createVariable('SPC_SOAP',o.datatype,dimensions=o.dimensions,fill_value=o._FillValue)
| + | |
− | soap = f['SPC_SOAP']
| + | |
− | soap[:] = 0.0
| + | |
− | soap.long_name= 'SOAP tracer'
| + | |
− | soap.units = o.units
| + | |
− | soap.add_offset = 0.0
| + | |
− | soap.scale_factor = 1.0
| + | |
− | soap.missing_value = 1.0e30
| + | |
− |
| + | |
− | f.close()
| + | |
− | | + | |
− | --[[User:Melissa Payer|Melissa Sulprizio]] ([[User talk:Melissa Payer|talk]]) 23:33, 7 November 2017 (UTC)
| + | |
− | | + | |
− | ==== Cropping netCDF restart files to nested domains ====
| + | |
− | | + | |
− | If needed, regrid a coarse netCDF restart file to the nested grid resolution using GAMAP routine [http://acmg.seas.harvard.edu/gamap/doc/by_alphabet/gamap_r.html#REGRIDH_RESTART_NC regridh_restart_nc.pro].
| + | |
− | | + | |
− | Cropping netCDF files can be easily achieved with tools such as [https://code.zmaw.de/projects/cdo CDO] or [http://nco.sourceforge.net/ NCO]. For example, CDO has a SELBOX operator for selecting a box by specifying the lat/lon bounds:
| + | |
− | | + | |
− | cdo sellonlatbox,lon1,lon2,lat1,lat2 infile outfile
| + | |
− | | + | |
− | See page 44 of the [https://code.zmaw.de/projects/cdo/embedded/cdo.pdf CDO guide] for more information.
| + | |
− | | + | |
− | --[[User:Melissa Payer|Melissa Sulprizio]] ([[User talk:Melissa Payer|talk]]) 19:04, 3 May 2017 (UTC)
| + | |
− | | + | |
− | === Vertical coordinates in netCDF files produced by GEOS-Chem ===
| + | |
− | | + | |
− | All netCDF files produced by GEOS-Chem (i.e. diagnostic files and restart files) adhere to the [[Preparing_data_files_for_use_with_HEMCO#The_COARDS_netCDF_standard|the COARDS netCDF convention]] for the lon, lat, and time dimensions.
| + | |
− | | + | |
− | For the vertical dimension, we have chosen to use the following coordinate variables, emulating the file format of the NCAR Community Earth System Model (CESM):
| + | |
− | | + | |
− | variables:
| + | |
− | double lev(lev) ;
| + | |
− | lev:long_name = "hybrid level at midpoints (1000*(A+B))" ;
| + | |
− | lev:units = "level" ;
| + | |
− | lev:positive = "down" ;\
| + | |
− | lev:standard_name = "atmosphere_hybrid_sigma_pressure_coordinate" ;
| + | |
− | lev:formula_terms = "a: hyam b: hybm p0: P0 ps: PS" ;
| + | |
− | double hyam(lev) ;
| + | |
− | hyam:long_name = "hybrid A coefficient at layer midpoints" ;
| + | |
− | double hybm(lev) ;
| + | |
− | hybm:long_name = "hybrid B coefficient at layer midpoints" ;
| + | |
− | double ilev(ilev) ;
| + | |
− | ilev:long_name = "hybrid level at interfaces (1000*(A+B))" ;
| + | |
− | ilev:units = "level" ;
| + | |
− | ilev:positive = "down" ;
| + | |
− | ilev:standard_name = "atmosphere_hybrid_sigma_pressure_coordinate" ;
| + | |
− | ilev:formula_terms = "a: hyai b: hybi p0: P0 ps: PS" ;
| + | |
− | double hyai(ilev) ;
| + | |
− | hyai:long_name = "hybrid A coefficient at layer interfaces" ;
| + | |
− | double hybi(ilev) ;
| + | |
− | hybi:long_name = "hybrid B coefficient at layer interfaces" ;
| + | |
− | double P0 ;
| + | |
− | P0:long_name = "reference pressure" ;
| + | |
− | | + | |
− | The <tt>lev</tt> variable is used for data that is placed on the midpoints between vertical levels. This is an "approximate" eta coordinate, which is close to 1 at the surface and close to zero at the atmosphere top.
| + | |
− | | + | |
− | lev = 0.99250002413, 0.97749990013, 0.962499776, 0.947499955, 0.93250006,
| + | |
− | 0.91749991, 0.90249991, 0.88749996, 0.87249996, 0.85750006, 0.842500125,
| + | |
− | 0.82750016, 0.8100002, 0.78750002, 0.762499965, 0.737500105, 0.7125001,
| + | |
− | 0.6875001, 0.65625015, 0.6187502, 0.58125015, 0.5437501, 0.5062501,
| + | |
− | 0.4687501, 0.4312501, 0.3937501, 0.3562501, 0.31279158, 0.26647905,
| + | |
− | 0.2265135325, 0.192541016587707, 0.163661504087706, 0.139115, 0.11825,
| + | |
− | 0.10051436, 0.085439015, 0.07255786, 0.06149566, 0.05201591, 0.04390966,
| + | |
− | 0.03699271, 0.03108891, 0.02604911, 0.021761005, 0.01812435, 0.01505025,
| + | |
− | 0.01246015, 0.010284921, 0.008456392, 0.0069183215, 0.005631801,
| + | |
− | 0.004561686, 0.003676501, 0.002948321, 0.0023525905, 0.00186788,
| + | |
− | 0.00147565, 0.001159975, 0.00090728705, 0.0007059566, 0.0005462926,
| + | |
− | 0.0004204236, 0.0003217836, 0.00024493755, 0.000185422, 0.000139599,
| + | |
− | 0.00010452401, 7.7672515e-05, 5.679251e-05, 4.0142505e-05, 2.635e-05,
| + | |
− | 1.5e-05 ;
| + | |
− | | + | |
− | The <tt>lev</tt> variable may be used for quick plotting. To compute the actual pressure at the midpoint of the grid box (I,J,L), you will need to supply your own 2-D surface pressure field (e.g. saved from another diagnostic file):
| + | |
− | | + | |
− | Pmid = ( hyam(L) * PS(I,J) ) + hybm(L)
| + | |
− | | + | |
− | The <tt>ilev</tt> variable is used for data that is placed on vertical level edges or "interfaces" (hence the "i" in <tt>ilev</tt>). This is also an "approximate" eta coordinate.
| + | |
− | | + | |
− | ilev = 1, 0.98500004826, 0.969999752, 0.9549998, 0.94000011, 0.92500001,
| + | |
− | 0.90999981, 0.89500001, 0.87999991, 0.86500001, 0.85000011, 0.83500014,
| + | |
− | 0.82000018, 0.80000022, 0.77499982, 0.75000011, 0.7250001, 0.7000001,
| + | |
− | 0.6750001, 0.6375002, 0.6000002, 0.5625001, 0.5250001, 0.4875001,
| + | |
− | 0.4500001, 0.4125001, 0.3750001, 0.3375001, 0.28808306, 0.24487504,
| + | |
− | 0.208152025, 0.176930008175413, 0.150393, 0.127837, 0.108663, 0.09236572,
| + | |
− | 0.07851231, 0.06660341, 0.05638791, 0.04764391, 0.04017541, 0.03381001,
| + | |
− | 0.02836781, 0.02373041, 0.0197916, 0.0164571, 0.0136434, 0.0112769,
| + | |
− | 0.009292942, 0.007619842, 0.006216801, 0.005046801, 0.004076571,
| + | |
− | 0.003276431, 0.002620211, 0.00208497, 0.00165079, 0.00130051, 0.00101944,
| + | |
− | 0.0007951341, 0.0006167791, 0.0004758061, 0.0003650411, 0.0002785261,
| + | |
− | 0.000211349, 0.000159495, 0.000119703, 8.934502e-05, 6.600001e-05,
| + | |
− | 4.758501e-05, 3.27e-05, 2e-05, 1e-05 ;
| + | |
− | | + | |
− | To compute the actual pressure at the bottom and top edges of the grid box (I,J,L), you will need to supply your own 2-D surface pressure field (e.g. saved from another diagnostic file):
| + | |
− | | + | |
− | Pbot = ( hyai(L ) * PS(I,J) ) + hybi(L )
| + | |
− | Ptop = ( hyai(L+1) * PS(I,J) ) + hybi(L+1)
| + | |
− | | + | |
− | '''''[[User:Sebastian D. Eastham|Seb Eastham]] wrote:'''''
| + | |
− | | + | |
− | <blockquote>Daniel Rothenberg has been in touch regarding how vertical co-ordinates are set up in CESM. It's mostly the same (as COARDS), with the exception that EVERY file which includes a vertical co-ordinate also includes both <tt>lev</tt> and <tt>ilev</tt> dimensions. In this case, <tt>lev</tt> is the cell mid-points, and <tt>ilev</tt> is cell edges ("i"nterfaces). Both <tt>lev</tt> and ilev are set up just as we discussed before, although CESM doesn't include the <tt>formula_terms</tt> attribute, which I think we should use. In this case <tt>formula_terms</tt> would be <tt>ap: hyai b: hybi p0: P0 ps: ps</tt> for the <tt>ilev</tt> dimension and <tt>ap: hyam b: hybm p0: P0 ps: ps</tt> for the <tt>lev</tt> dimension.
| + | |
− | | + | |
− | #For <tt>lev</tt> and <tt>ilev</tt>, they should be set to <tt>A/1000 + B</tt>, where <tt>P0</tt> = 1000.
| + | |
− | #Both <tt>ilev</tt> and <tt>lev</tt> should be in all 3D files, even though only one of the two dimensions is used in a given file).
| + | |
− | #I think we can leave <tt>PS</tt> out of the files. It will make them non-compliant but I think it's one of those "accepted" non-compliances. If it turns out to be a problem we can revisit it later, but for plotting the expectation is that it's up to the user to figure out where they're getting <tt>PS</tt> from, rather than automatically archiving it in every file.
| + | |
− | </blockquote>
| + | |
− | | + | |
− | --[[User:Bmy|Bob Yantosca]] ([[User talk:Bmy|talk]]) 19:27, 5 September 2018 (UTC)
| + | |
− | | + | |
− | == GEOS-Chem v11-01 ==
| + | |
− | | + | |
− | === Impact of Output File Format Change ===
| + | |
− | | + | |
− | While migrating from binary punch to netCDF file formats is a structural change, several changes are also being introduced which will result in output differences. These changes include the following:
| + | |
− | | + | |
− | #Starting in v11-01f, species will always be written to the GEOS-Chem netCDF restart file and will always be initialized to restart file values if present in the file. In contrast, when using binary punch format restart files, species are initialized to background values stored in <tt>globchem.dat</tt> by default unless the option to read/write the species restart file is turned on in <tt>input.geos</tt>. ''This change will only impact multi-segmented full chemistry runs where users previously did not enable the species restart file''.
| + | |
− | #Starting in v11-01f, species units will be [mol/mol] in the netCDF restart file while they are in [molec/cm3/box] in the binary punch restart file. This change will facilitate the future replacement of tracers with species. However, due to the meteorology state-dependency of the unit conversion from [molec/cm3] to [mol/mol], there will be small differences in <tt>[molec/mc3]</tt> species values between using binary punch and Netcdf in multi-segmented runs. This is because the meteorology used at restart write is the timestep ''before'' the simulation end time while the meteorology used at restart read is the simulation start time which is the same as previous end time. ''This change will only impact multi-segmented full chemistry runs where users previously enabled the species restart file''.
| + | |
− | | + | |
− | === Output files in v11-01 ===
| + | |
− | | + | |
− | Output files in [[GEOS-Chem v11-01]] are listed below.
| + | |
− | | + | |
− | {| border=1 cellspacing=0 cellpadding=5
| + | |
− | |-bgcolor="#CCCCCC"
| + | |
− | !width="325px"|File
| + | |
− | !width="175px"|Data Format
| + | |
− | !width="450px"|Description
| + | |
− | | + | |
− | |-valign="top"
| + | |
− | |<tt>trac_avg.{met}_{grid}_{sim}.YYYYMMDDhh</tt>
| + | |
− | |BPCH
| + | |
− | |GEOS-Chem output diagnostic data
| + | |
− | | + | |
− | |-valign="top"
| + | |
− | |<tt>GEOSChem_restart.YYYYMMDDhhmm</tt>
| + | |
− | |NetCDF
| + | |
− | |GEOS-Chem output restart file containing advected tracer concentrations in units of [mol/mol]. If running a full chemistry simulation, then the restart file will also contain species concentrations, also in units of [mol/mol]. If running a mercury simulation, then the restart file will also contain ocean mercury data.
| + | |
− | | + | |
− | |-valign="top"
| + | |
− | |<tt>HEMCO_restart.YYYYMMDDhhmm</tt>
| + | |
− | |NetCDF
| + | |
− | |HEMCO output restart file
| + | |
− | | + | |
− | |-valign="top"
| + | |
− | |<tt>HEMCO_diagnostics.YYYYMMDDhhmm</tt>
| + | |
− | |NetCDF
| + | |
− | |HEMCO output diagnostics data
| + | |
− | | + | |
− | |-bgcolor="#cccccc"
| + | |
− | !colspan="21"|Optional output files:
| + | |
− | | + | |
− | |-valign="top" | + | |
− | |<tt>stations.YYYYMMDDhh</tt>
| + | |
− | |BPCH
| + | |
− | |Diagnostic output file that contains ND48 station timeseries output as specified in <tt>input.geos</tt>.
| + | |
− | | + | |
− | |-valign="top"--[[User:Lizzie Lundgren|Lizzie Lundgren]] ([[User talk:Lizzie Lundgren|talk]]) 16:04, 10 November 2016 (UTC)
| + | |
− | |<tt>tsYYYYMMDDhh.bpch</tt>
| + | |
− | |BPCH
| + | |
− | |Diagnostic output file that contains ND49 instantaneous timeseries output as specified in <tt>input.geos</tt>.
| + | |
− | | + | |
− | |-valign="top"
| + | |
− | |<tt>ts_24hr_avg.YYYYMMDDhh.bpch</tt>
| + | |
− | |BPCH
| + | |
− | |Diagnostic output file that contains ND50 24-hour average timeseries output as specified in <tt>input.geos</tt>.
| + | |
− | | + | |
− | |-valign="top"
| + | |
− | |<tt>ts_satellite.YYYYMMDDhh.bpch</tt>
| + | |
− | |BPCH
| + | |
− | |Diagnostic output file that contains ND51 local-time-average (satellite) timeseries output as specified in <tt>input.geos</tt>.
| + | |
− | | + | |
− | |-valign="top"
| + | |
− | |<tt>paranox_ts.YYYYMMDDhh.bpch</tt>
| + | |
− | |BPCH
| + | |
− | |Diagnostic output file that contains ND63 ship timeseries output as specified in <tt>input.geos</tt>.
| + | |
− | | + | |
− | |}
| + | |
− | | + | |
− | == Previous issues that are now resolved ==
| + | |
− | | + | |
− | === Enable compression in netCDF-4 output files ===
| + | |
− | | + | |
− | <span style="color:green">'''''This update was included in [[GEOS-Chem v11-02#v11-02a|v11-02a]] (approved 12 May 2017).'''''</span>
| + | |
− | | + | |
− | For more information about this issue, please see [[The_NcdfUtilities_package#Enable_compression_in_netCDF-4_output_files|this post on our ''The NcdfUtilities package'' wiki page]].
| + | |
− | | + | |
− | --[[User:Bmy|Bob Yantosca]] ([[User talk:Bmy|talk]]) 17:50, 24 May 2018 (UTC)
| + | |
− | | + | |
− | === Improve write speed of netCDF output files ===
| + | |
− | | + | |
− | <span style="color:green">'''''This update was included in [[GEOS-Chem v11-01#v11-01 public release|GEOS-Chem v11-01 public release]]'''''</span>
| + | |
− | | + | |
− | For more information about this issue, please see [[The_NcdfUtilities_package#Improve_write_speed_of_netCDF_output_files|this post on our ''The NcdfUtilities package'' wiki page]].
| + | |
− | | + | |
− | --[[User:Bmy|Bob Yantosca]] ([[User talk:Bmy|talk]]) 19:05, 8 March 2017 (UTC)
| + | |
− | | + | |
− | === GAMAP can now read GEOS-Chem restart files in netCDF format ===
| + | |
− | | + | |
− | <span style="color:green">'''''This update was included in [[GEOS-Chem v11-01#v11-01 public release|GEOS-Chem v11-01 public release]]'''''</span>
| + | |
− | | + | |
− | Starting in [[GEOS-Chem v11-01]], all restart files are now saved in [[Preparing_data_files_for_use_with_HEMCO#The_COARDS_netCDF_standard|COARDS-compliant netCDF format]]. We have had to make some minor modifications to both GEOS-Chem and GAMAP in order to allow GAMAP to read these files. The table below gives a summary of these modifications.
| + | |
− | | + | |
− | {| border=1 cellspacing=0 cellpadding=5
| + | |
− | |-valign="top" bgcolor="#CCCCCC"
| + | |
− | !width="100px"|GEOS-Chem or GAMAP?
| + | |
− | !width="250px"|File
| + | |
− | !width="650px"|Modification
| + | |
− | | + | |
− | |-valign="top"
| + | |
− | |GEOS-Chem
| + | |
− | |<tt>GeosCore/gamap_mod.F</tt>
| + | |
− | |In routine <tt>INIT_TRACERINFO</tt>, we now write metadata for all species (advected or not) to the tracerinfo.dat file under the ND45 tracer concentration diagnostic section. Because the netCDF restart file contains concentrations for both advected and non-advected species, we need to make sure that the <tt>tracerinfo.dat</tt> file created by GEOS-Chem contains metadata for all species.
| + | |
− | | + | |
− | |-valign="top"
| + | |
− | |GAMAP
| + | |
− | |<tt>internals/ctm_open_file.pro</tt>
| + | |
− | |The prior algorithm always assumed that a netCDF file would end in either <tt>.nc</tt> or <tt>.nc4</tt>. We now have removed this restriction. We now split the filename string on <tt>.</tt> and then examine the substrings for <tt>nc</tt> or <tt>nc4</tt> (case-insensitive).
| + | |
− | | + | |
− | |-valign="top"
| + | |
− | |GAMAP
| + | |
− | |<tt>internals/ctm_read_coards.pro</tt>
| + | |
− | |Added some minor modifications to read netCDF restart files:
| + | |
− |
| + | |
− | #Determine vertical grid from the number of layers (if the <tt>Model</tt> global attribute is not specified).
| + | |
− | #Assign category <tt>DXYP</tt> to variable <tt>AREA</tt>, which is included in the netCDF restart file.
| + | |
− | #Assign category <tt>IJ-AVG-$</tt> to variables beginning with either <tt>SPC_</tt> or <tt>TRC_</tt>. Also remove the <tt>SPC_</tt> and <tt>TRC_</tt> from the variable name internally so that the variable name will match the metadata in <tt>tracerinfo.dat</tt>.
| + | |
− | | + | |
− | |}
| + | |
− | | + | |
− | --[[User:Bmy|Bob Yantosca]] ([[User talk:Bmy|talk]]) 19:29, 23 January 2017 (UTC)
| + | |
− | | + | |
− | == Obsolete information ==
| + | |
− | | + | |
− | These sections pertain to data that is no longer relevant in the latest GEOS-Chem versions. We shall preserve this information here for reference.
| + | |
− | | + | |
− | === Legacy files made obsolete in v11-01 ===
| + | |
− | | + | |
− | [[Image:Obsolete.jpg]]
| + | |
− | | + | |
− | The following files were rendered obsolete starting in [[GEOS-Chem v11-01]].
| + | |
− | | + | |
− | {| border=1 cellspacing=0 cellpadding=5
| + | |
− | |-bgcolor="#CCCCCC"
| + | |
− | !width="350px"|File
| + | |
− | !width="325px"|Description
| + | |
− | !width="375px"|Reason for removal
| + | |
− | | + | |
− | |-valign="top"
| + | |
− | |<tt>smv2.log</tt>
| + | |
− | |SMVGEAR II log file containing information about reactions and species used by the ND65 prod-loss diagnostic.
| + | |
− | |This was removed when SMVGEAR was replaced by [[FlexChem]].
| + | |
− | | + | |
− | |-valign="top"
| + | |
− | |<tt>trac_rst.geosfp_4x5_fullchem.YYYYMMDDhhmm.mp</tt> or<br><tt>trac_rst.geosfp_4x5_fullchem.YYYYMMDDhhmm.sp</tt>
| + | |
− | |Tracer restart file that contains instantaneous tracer concentrations at simulation end time.
| + | |
− | |The new netCDF restart file format (prefix <tt>GEOSCHEM_restart_</tt>) contains concentrations for all species (whether advected or not).
| + | |
− | | + | |
− | |-valign="top"
| + | |
− | |<tt>spec_rst.geosfp_4x5_fullchem.YYYYMMDDhh</tt>
| + | |
− | |Species restart file that contains instantaneous species concentrations at simulation end time. This file saves instantaneous concentrations of all chemical species (listed in globchem.dat) on all levels and may be used to initialize another GEOS-Chem simulation.
| + | |
− | |The new netCDF restart file format (prefix <tt>GEOSCHEM_restart_</tt>) contains concentrations for all species (whether advected or not).
| + | |
− | | + | |
− | |-valign="top"
| + | |
− | |<tt>ocean_rst.YYYYMMDDhhmm.nc</tt>
| + | |
− | |Ocean mercury restart file containing several quantities pertaining to the ocean mercury module saved at the end of the GEOS-Chem simulation. Like the tracer restart file, this file may be used to initialize another GEOS-Chem [[Mercury|mercury simulation]].
| + | |
− | |The new netCDF restart file format (prefix <tt>GEOSCHEM_restart_</tt>) contains concentrations for all species (whether advected or not).
| + | |
− | | + | |
− | |}
| + | |
− | | + | |
− | | + | |
− | === GEOS-Chem v10-01 and earlier versions ===
| + | |
− | | + | |
− | [[Image:Obsolete.jpg]]
| + | |
− | | + | |
− | After a GEOS-Chem simulation has completed successfully, you will see several new files in the run directory. Some of these are generated for all simulations while others may be simulation-dependent.
| + | |
− | | + | |
− | For full-chemistry simulations, the following files are saved out:
| + | |
− | | + | |
− | {| border=1 cellspacing=0 cellpadding=5
| + | |
− | |-bgcolor="#CCCCCC"
| + | |
− | !width="350px"|File
| + | |
− | !width="750px"|Description
| + | |
− | | + | |
− | |-valign="top"
| + | |
− | |<tt>trac_avg.geosfp_4x5_fullchem.YYYYMMDDhh.mp</tt> or<br><tt>trac_avg.geosfp_4x5_fullchem.YYYYMMDDhh.sp</tt>
| + | |
− | |Diagnostic output file that contains time-averaged output for diagnostics enabled in <tt>input.geos</tt>.
| + | |
− | | + | |
− | <b>Data in this file are in [http://acmg.seas.harvard.edu/gamap/doc/Chapter_6.html#6.2 "binary punch" format] and may be viewed and manipulated using [http://acmg.seas.harvard.edu/gamap GAMAP].</b>
| + | |
− | | + | |
− | |-valign="top"
| + | |
− | |<tt>trac_rst.geosfp_4x5_fullchem.YYYYMMDDhhmm.mp</tt> or<br><tt>trac_rst.geosfp_4x5_fullchem.YYYYMMDDhhmm.sp</tt>
| + | |
− | |Tracer restart file that contains instantaneous tracer concentrations at simulation end time. This file saves instantaneous concentrations of all transported tracers (listed in the [http://acmg.seas.harvard.edu/geos/doc/man/chapter_5.html#TracerMenu Tracer Menu] of <tt>input.geos</tt>) on all levels and may be used to initialize another GEOS-Chem simulation.
| + | |
− | | + | |
− | <b>Data in this file are in [http://acmg.seas.harvard.edu/gamap/doc/Chapter_6.html#6.2 "binary punch" format] and may be viewed and manipulated using [http://acmg.seas.harvard.edu/gamap/ GAMAP].</b>
| + | |
− | | + | |
− | |-valign="top"
| + | |
− | |<tt>HEMCO_restart.YYYYMMDDhh.nc.mp</tt> or<br><tt>HEMCO_restart.YYYYMMDDhh.nc.sp</tt>
| + | |
− | |HEMCO restart file containing several quantities tracked by HEMCO for use in initializing another GEOS-Chem simulation. This file includes data used to calculate [[Hudman_et_al_2012_soil_NOx_emissions_algorithm|soil NOx emissions]], [[Ship_emissions#PARANOX_ship_plume_model|PARANOX ship plume chemistry]], and [[MEGAN_v2.1_plus_Guenther_2012_biogenic_emissions|MEGAN biogenic emissions]].
| + | |
− | | + | |
− | <b>Data in this file are in [http://www.unidata.ucar.edu/software/netcdf/docs/faq.html#whatisit netcdf format].</b>
| + | |
− | | + | |
− | |-valign="top"
| + | |
− | |<tt>v10-01.geosfp_4x5_fullchem.mp</tt> or<br><tt>v10-01.geosfp_4x5_fullchem.sp</tt>
| + | |
− | |GEOS-Chem simulation log file.
| + | |
− | | + | |
− | <b>Open this file in a text editor to diagnose compile and run errors.</b>
| + | |
− | | + | |
− | |-valign="top"
| + | |
− | |<tt>HEMCO.log.mp</tt> or<br><tt>HEMCO.log.sp</tt>
| + | |
− | |HEMCO log file containing detailed information about the emissions read from disk into HEMCO.
| + | |
− | | + | |
− | <b>Open this file in a text editor to diagnose HEMCO errors.</b>
| + | |
− | | + | |
− | |-valign="top"
| + | |
− | |<tt>smv2.log</tt>
| + | |
− | |SMVGEAR II log file containing information about reactions and species used by the ND65 prod-loss diagnostic. This file is only produced when SMVGEAR, and thus it will not be created for the [[GEOS-Chem_chemistry_mechanisms#Specialty_simulations|specialty simulations]].
| + | |
− | | + | |
− | <b>Open this file in a text editor to see if SMVGEAR read the [http://acmg.seas.harvard.edu/geos/doc/man/chapter_5.html#Globchem <tt>globchem.dat</tt>] file properly.</b>
| + | |
− | | + | |
− | |-valign="top"
| + | |
− | |<tt>diaginfo.dat</tt>
| + | |
− | |File generated by GEOS-Chem that contains diagnostic quantities for use with [http://acmg.seas.harvard.edu/gamap/ GAMAP].
| + | |
− | | + | |
− | |-valign="top"
| + | |
− | |<tt>tracerinfo.dat</tt>
| + | |
− | |File generated by GEOS-Chem that contains tracer name metadata for use with [http://acmg.seas.harvard.edu/gamap/ GAMAP].
| + | |
− | |}
| + | |
− | | + | |
− | | + | |
− | Additional files may be created by GEOS-Chem, depending on your simulation type and the options specified in the <tt>input.geos</tt> file. Other optional GEOS-Chem output files include:
| + | |
− | | + | |
− | {| border=1 cellspacing=0 cellpadding=5
| + | |
− | |-bgcolor="#CCCCCC"
| + | |
− | !width="350px"|File
| + | |
− | !width="750px"|Description
| + | |
− | | + | |
− | |-valign="top"
| + | |
− | |<tt>spec_rst.geosfp_4x5_fullchem.YYYYMMDDhh</tt>
| + | |
− | |Species restart file that contains instantaneous species concentrations at simulation end time. This file saves instantaneous concentrations of all chemical species (listed in [http://acmg.seas.harvard.edu/geos/doc/man/chapter_5.html#Globchem <tt>globchem.dat</tt>]) on all levels and may be used to initialize another GEOS-Chem simulation.
| + | |
− | | + | |
− | NOTE: If you are going to be running a very long GEOS–Chem simulation and must split your simulation into several stages (i.e. in order to stay within the computational time limits of your system), then you should set LSVCSPEC to T in the [http://acmg.seas.harvard.edu/geos/doc/man/chapter_5.html#ChemMenu Chemistry Menu] of <tt>input.geos</tt>. That will make sure that the chemical species concentrations are preserved when the next run stage starts. Otherwise, GEOS-Chem use the default species concentrations specified in globchem.dat at the beginning of each run.
| + | |
− | | + | |
− | <b>Data in this file are in [http://acmg.seas.harvard.edu/gamap/doc/Chapter_6.html#6.2 "binary punch" format] and may be viewed and manipulated using [http://acmg.seas.harvard.edu/gamap/ GAMAP].</b>
| + | |
− | | + | |
− | |-valign="top"
| + | |
− | |<tt>stations.YYYYMMDDhh</tt>
| + | |
− | |Diagnostic output file that contains ND48 station timeseries output as specified in <tt>input.geos</tt>.
| + | |
− | | + | |
− | <b>Data in this file are in [http://acmg.seas.harvard.edu/gamap/doc/Chapter_6.html#6.2 "binary punch" format] and may be viewed and manipulated using [http://acmg.seas.harvard.edu/gamap/ GAMAP]. For more information on working with ND48 timeseries output, see our [http://acmg.seas.harvard.edu/gamap/doc/timeseries_brief_tutorial.pdf timeseries tutorial].</b>
| + | |
− | | + | |
− | |-valign="top"
| + | |
− | |<tt>tsYYYYMMDDhh.bpch</tt>
| + | |
− | |Diagnostic output file that contains ND49 instantaneous timeseries output as specified in <tt>input.geos</tt>.
| + | |
− | | + | |
− | <b>Data in this file are in [http://acmg.seas.harvard.edu/gamap/doc/Chapter_6.html#6.2 "binary punch" format] and may be viewed and manipulated using [http://acmg.seas.harvard.edu/gamap/ GAMAP]. For more information on working with ND49 timeseries output, see our [http://acmg.seas.harvard.edu/gamap/doc/timeseries_brief_tutorial.pdf timeseries tutorial].</b>
| + | |
− | | + | |
− | |-valign="top"
| + | |
− | |<tt>ts_24hr_avg.YYYYMMDDhh.bpch</tt>
| + | |
− | |Diagnostic output file that contains ND50 24-hour average timeseries output as specified in <tt>input.geos</tt>.
| + | |
− | | + | |
− | <b>Data in this file are in [http://acmg.seas.harvard.edu/gamap/doc/Chapter_6.html#6.2 "binary punch" format] and may be viewed and manipulated using [http://acmg.seas.harvard.edu/gamap/ GAMAP].</b>
| + | |
− | | + | |
− | |-valign="top"
| + | |
− | |<tt>ts_satellite.YYYYMMDDhh.bpch</tt>
| + | |
− | |Diagnostic output file that contains ND51 local-time-average (satellite) timeseries output as specified in <tt>input.geos</tt>.
| + | |
− | | + | |
− | <b>Data in this file are in [http://acmg.seas.harvard.edu/gamap/doc/Chapter_6.html#6.2 "binary punch" format] and may be viewed and manipulated using [http://acmg.seas.harvard.edu/gamap/ GAMAP].</b>
| + | |
− | | + | |
− | |-valign="top"
| + | |
− | |<tt>paranox_ts.YYYYMMDDhh.bpch</tt>
| + | |
− | |Diagnostic output file that contains ND63 ship timeseries output as specified in <tt>input.geos</tt>.
| + | |
− | | + | |
− | <b>Data in this file are in [http://acmg.seas.harvard.edu/gamap/doc/Chapter_6.html#6.2 "binary punch" format] and may be viewed and manipulated using [http://acmg.seas.harvard.edu/gamap/ GAMAP].</b>
| + | |
− | | + | |
− | |-valign="top"
| + | |
− | |<tt>ocean_rst.YYYYMMDDhhmm.nc</tt>
| + | |
− | |Ocean mercury restart file containing several quantities pertaining to the ocean mercury module saved at the end of the GEOS-Chem simulation. Like the tracer restart file, this file may be used to initialize another GEOS-Chem [[Mercury|mercury simulation]].
| + | |
− | | + | |
− | <b>Data in this file are in the in [http://www.unidata.ucar.edu/software/netcdf/docs/faq.html#whatisit netcdf format].</b>
| + | |
− | | + | |
− | |-valign="top"
| + | |
− | |<tt>plane.log.YYYYMMDD</tt>
| + | |
− | |Diagnostic output file containing planeflight information scheduled via the [[Planeflight diagnostic|<tt>Planeflight.dat</tt>]] file. This diagnostic is turned on in the [[GEOS-Chem_Input_Files#Planeflight_Menu|Planeflight Menu]] of <tt>input.geos</tt>.
| + | |
− | | + | |
− | <b>The data in this text file can be read and plotted using GAMAP routines [http://acmg.seas.harvard.edu/gamap/doc/by_alphabet/gamap_c.html#CTM_READ_PLANEFLIGHT <tt>CTM_READ_PLANEFLIGHT</tt>] and [http://acmg.seas.harvard.edu/gamap/doc/by_alphabet/gamap_p.html#PLANE_PLOT <tt>PLANE_PLOT</tt>].
| + | |
− | | + | |
− | |}
| + | |
− | | + | |
− | --[[User:Melissa Payer|Melissa Sulprizio]] 15:42, 3 April 2015 (EDT)
| + | |