Script non_annual

Script non_annual calls the following programs:

Program readgeia for GEIA NO$_{x}$ and SO$_{2}$ emissions

Input files are

• G_NOX_hi.1x1 GEIA NO$_{x}$ emissions above 100m
• G_NOX_lo.1x1 GEIA NO$_{x}$ emissions below 100m
• G_SO2_hi.1x1 GEIA SO$_{2}$ emissions above 100m
• G_SO2_lo.1x1 GEIA SO$_{2}$ emissions below 100m

Output files are GNOXs.d and GSO2s.d These are only used to describe the seasonal cycle and vertical distribution of the emissions, not the geographical distribution or total.

Program geia2tm3 for GEIA CO, isoprene and NMHC emissions

Script non_annual first preprocesses GEIA emissions compiled by Sander Houweling.

• CO biogenic emissions from vegetation (and ocean?), input file cobio90yr1.1a, output file cobio.d
• isoprene monthly emissions, input file isop90mn1.1a, output file fisop.d
• anthropogenic VOC yearly emissions translated to CBM-4 elements, input file evavocyr1.1a, output file feva.d
• biomass burning VOC yearly emissions translated to CBM-4 elements, input file bmbvocyr1.1a, output file fbmb.d
• anthropogenic VOC yearly other emissions translated to CBM-4 elements, input file othvocyr1.1a, output file foth.d

Fortran program geia2tm3 does the interpolation from source to target grid by calculating the fractional area of overlapping cells.
The emissions are valid for 1990.

Program geiach42tm3 for methane emissions

This program interpolates methane emissions compiled by Sander Houweling to the TM gird. Input files are:

• ch4rice_SH.d : Rice field emissions from Fung (Based on [Matthews et al., 1991])
• ch4stat_SH.d : Various stationary emissions (ruminants, termites, oceans, waste treatment, soil oxidation, a.o.)
• ch4wetl_SH.d : Wetland emissions (Based on [Walter, 1998])
• ch4bmbseas_SH.d : Seasonally dependent biomass burning
• ch4bmbnoseas_SH.d : Seasonally independent biomass burning

and the output files ch4rice.d, ch4stat.d, ch4wetl.d, ch4bmbs.d, and ch4bmbns.d resp. Software to compile these emissions needs to be included in this package. One may furthermore note that program geiach42tm3 is very similar to program geia2tm3.

The Fung data are available from
http://www.giss.nasa.gov/data/ch4fung/

program compvolc for volcanic emissions

Script non_annual uses compvolc to calculate volcanic emissions. Input files are heights_TM3.d which gives the heights of the TM levels and thus depends on the chosen number and location of TM levels. The input file cont_volc.1x1 has a number of descriptive lines followed by data lines in format (2(I3),1x,E10.5,1x): The unit in the input file is Mg/d. The input file loc_volc.1x1 is identical (except for a header line) to the file volcalt.txt that can be downloaded from the GEIA web site (checked 11 October 2003) at http://www.geiacenter.org/data/volcalt.txt. It contains the for 1509 volcanoes the name, latitude, longitude and altitude (m). The input file cont_volc.1x1 iis identical to the file COVavg1.1a that can be downloaded from the GEIA web site (checked 11 October 2003) at
http://www.geiacenter.org/data/COVavg1.1a. It contains the annual maximum !, non-TOMS measured, SO2 fluxes (Mg/d) from 43 continuously erupting volcanoes.

Program ancat for NO$_{x}$ emissions by aircraft

There are 5 ANCAT/EC2 files for 1991/1992:

• ancat_jul91 : civil aircraft emissions for July 1991
• ancat_oct91 : civil aircraft emissions for October 1991
• ancat_jan92 : civil aircraft emissions for January 1992
• ancat_apr92 : civil aircraft emissions for April 1992
• ancat_mil91-92 : miltary aircraft emissions for one year

Fuel and emissions have been disaggregated on a global grid of 1 x 1 degree x 1 km. Each file consists of a number of lines defining the fuel use and NO$_{2}$ emission in a cell. The data are organised in 5 columns: Latitude, Longitude, Level, Fuel, Emission. The data have written in fortran format: (3I5,2E13.5), see the accompanying readme file.

• Latitude: A negative sign indicates a Southern latitude. Range -90 - -1, +1 - +90.
• Longitude: A negative sign indicates a Western longitude. Range -180 - -1, +1 - +180.
• Level: Indicates vertical level. Range -1 - +16.
  -1 = taxiing; 0 = 0-1km altitude;
  1 = 1-2km altitude; 2 = 2-3km altitude;
  ...
  15 = 15-16km altitude; 16 = 16+km altitude.
  
• fuel used in units of kg/month/cell (except military which is in kg/year)
• NO$_{x}$ emitted in units of kg NO$_{2}$/month/cell (except military which is in kg NO$_{2}$/year)

There are also forecasts available for 2015 (not used).

Program zch4 for methane surface concentrations

input file ch43tm1.dat output file zch4.d

Program ch4relds for methane loss rate factors

Script non_annual calls program ch4relds to process the RIVM 2-D model methane loss rates due to reactions with OH, Cl ad O1D from input file ch4des_rivm2d.d into scaling factors for the TM methane loss rates due to the reaction with OH that are stored in the TM file ch4rds.d.

Program zsch4 for the methane boundary condition at 10hPa

Script non_annual calls program zsch4 to .. input file ch4_10hPa.data output TM file zsch4.d.

Program uars_ratio for the HNO3:O3 upper boundary condition

Ths program calculates the output TM file uars_ratio.d by interpolating and extrapolating! latitudinally by Hermite cubic interpolation with muherm. The original data in input file uars_hno3_o3_ratio_10hPa.d are given for 12 months at 39 mean latitudes at 4 degree intervals, i.e. at -76, -72 .. 76 degrees N.

Program interp_dms for the DMS ocean surface water concentrations

This program interp_dms interpolates the DMS monthly mean sea surface concentrations from the input file dms_array.dat to the the TM grid using Hermite interpolation (subroutine muherm) and writes the output to TM file DMSconc.d.

Program rspiro for the DMS and H2S land emissions

This program regrids the land emissions from [Spiro et al., 1992]. Input files are:

• inv.sol.dms : DMS emission by soils
• inv.sol.h2s : H2S emission by soils
• inv.veg.dms : DMS emission by vegetation
• inv.veg.h2s : H2S emission by vegetation

The output TM file DMSland.d contains the sum of DMS and H2S emissions.

Program fk for precipitation

input file precvtm1.dat output TM file fk.d. These are very old data ! We could use ECMWF data for this now. The program assumes the mean surface pressure is 984 hPa and uses the TM3 hybrid level structure to regrid.

Program pland for pland fraction

input file pland720_360.d output TM file pland.d. Should perhaps not be used anymore (if it is used at all) since land fraction is also avaiable from ECMWF (PvV).

Program biomseas_tm3 for the seasonal distribution of biomass burning

Default year is 1985. input file output TM file bmb_cycle.d.

Parameters from ECHAM for calculation of dry deposition

Program intersoil for soil parameters

input file soilpH_T42 output TM file soilph.d.

Program interLAI for leaf area indices

12 monthly input files LAI_T42_NN (with NN = 01 to 12) from ECHAM at T42 resolution are interpolated to output TM file lsmlai.d. Program is similar to intersoil.

Program interp_nh3 for

The natural emissions of ammonia from the input file nh3nat.1x1 are interpolated from the source grid to the TM target grid by calculating the fractional overlap of the grid cells. The output TM file is NATNH3.d. The program has been based upon ancat.