The package traj

The package traj allows to determine the trajectories of an ensemble of air parcels. Using assimilated winds and temperatures derived from meteorological data (currently UKMO or ECMWF), both forward and backward can be calculated. The numerical integration is based on the 4th order Runge-Kutta schema.

How to get the source code

The source code of the trajectory program is under control of the CVS (Concurrent Version System) utility. To get the source code the environment variable CVSROOT must be set to /usr/local/icg/icg1/archive. The command

cvs co -P traj

creates a directory traj in the current working directory and copies the current version of the trajectory package to this directory.

The directory traj includes the following subdirectories:

• doc: latex file for documentation

• mod: module files

• obj: object files

• source: source code

• templates: templates for configuration file

Compiling

The program can be compiled on different platforms with the makefile stored in the directory traj:

gmake traj

The program will be compiled for the used platform and an executable traj will be created.

Configuration file

The options for the execution of traj can be set in the configuration file traj.inp:

1)  example1/pos_EXP_97032512.nc  ! Name of initial positions file
2)  /dat_icg1/rad/example         ! Input directory
3)  I3U                           ! Prefix of wind files
4)  outdir                        ! Output directory
5)  TRU               ! Prefix for output file
6)  24                ! 1/2/3/4/6/8/12/24 -> Data every n hours
7)  y                 ! y/n - Backward/Forward trajectories
8)  y                 ! y/n - 3d(with theta_dot)/2d(without theta_dot) run
9)  THETA_DOT         ! short name of variable containing THETA-DOT information
10) 1                 ! eq. starts (a=0) diff. starts: a=1/2 same end/length
11) 12 23 03 1997     ! Endtime(hh dd mm yy) for a=0/1, lenght (hrs) for a=2
12) 1800              ! First timestep in seconds
13) 60                ! Second timestep in seconds
14) 0 1 0             ! Writing frequency (minutes hours days)
15) N.Thomas          ! user name
16) 4                 ! number of output parameters
17) TEMP              ! output parameters (one short name per line)
    PRESS             ! possible parameters are: TEMP,PRESS,U,V,PV,SZA
    PV
    SZA

Remarks:

1. Name of initial positions file (if it is not in current directory the path must be specified)
2. Directory with wind files

3. Prefix of wind files consisting of three characters:
   The filenames of isentropic files must correspond to the format isen_xxx_yymmddhh.nc. The prefix xxx must be specified in this line.
   The filenames of isobaric files must correspond to the format ecmwf_yymmddhh.nc. The prefix in this line is ignored.

4. Output directory

5. Prefix for output file consisting of three characters:
   The output filename corresponds to the format traj_xxx_3d_yymmddhh_yymmddhh.nc. The characters xxx must be specified in this line.

6. Interval between wind files in hours

7. Specification, if forward or backward trajectories should be calculated:
   y = backward trajectories
   n = forward trajectories

8. 2d or 3d run
   y = 3d run (with theta_dot)
   n = 2d run (without theta_dot)

9. Short name of variable containing THETA-DOT information

10. If the trajectories start at different times, in this line is specified if they should all end at the same time or if they should have the same length:
    1 = same end time for all trajectories
    2 = same length for all trajectories
    If all trajectories start at the same time, this specification is ignored.

11. If all trajectories start at the same time or if the trajectories start at different times but should end at the same time, the end time must be specified (hour, day, month, year). If the trajectories start at different times and should have the same length, the lenght of trajectories must be specified (in hours).
12. Timestep for calculation in seconds
13. This specification is only considered if the trajectories start at different times. The timestep is used in the period between the first starting and the last starting trajectory.

14. Writing frequency in the format mm hh dd, e. g. if 0 0 1 is specified, the data will be written out every 24 hours.

15. User name
16. Number of output parameters
17. List of output parameters:

    Possible output variables are TEMP, U, V, PV, SZA and for isentropic files also PRESS.

Initial positions of air parcels

This file includes the starting positions of the trajectories (longitudes, latitudes and levels) and the corresponding starting times. To create it the package pos can be used. It has to be a file in NetCDF data format with the following variables:

• TIME_INIT: time in Julian Seconds

• LAT: latitudes in deg N, valid range [-90,90]

• LON: longitudes in deg E, valid range [0,360]

• THETA|ZETA: potential temperature in K or
  PRESS: pressure in hPa

The filename must correspond to the following format:

• pos_xxx_yymmddhh.nc
  xxx : 1-3 characters
  yymmddhh : starting time
  e.g.: pos_ABC_99022812.nc

If the vertical coordinate ZETA is used, there must be a global attribute exp_VERTCOOR_name which is set to "zeta".

Example for initial positions file

Meteorological data

The data files with meteorological data used by the trajectory program must contain the coordinate variables:

• time: (one) time in julian seconds

• LAT: latitudes in deg N, valid range [-90,90]

• LON: longitudes in deg E, valid range [0,360]

• THETA|ZETA: potential temperature in K (in ascending order) or
  PRESS: pressure levels in hPa (in descending order)

LAT and LON must build a regular grid.

Besides the files must include at least the four-dimensional data sets for the variables:

• U: westerly wind

• V: southerly wind

• TEMP: temperature

• PV(isentropic files)|PV_ISOBAR(isobaric files): potential vorticity

• PRESS(isentropic files): pressure

• for 3d runs: THETA_DOT|ZETA_DOT: vertical velocity in K/day or
  OMEGA: vertical velocity in Pa/sec
  The name of this variable is specified in line 9 of traj.inp.

The filenames must correspond to the format:

• isen_xxx_yymmddhh.nc
  xxx : 3 characters
  yymmddhh : valid time
  e. g.: isen_ABC_99022812.nc

Example for isentropic file

Output file

The output file is created in NetCDF data format and is named:

• traj_xxx_3d_yymmddhh_yymmddhh.nc

The characters xxx are read from the configuration file (line 5) and the two times are the start and end time of trajectory calculation.

The dimensions in the trajectory output file are:

• time: number of written timesteps

• NPARTS: number of start positions/trajectories

The file contains the following variables:

• time(time): output times in julian seconds

• DATE_TIME(time,7): readable time (year, month, day, hour, minute, second, millisecond)

• LAT(time,NPARTS): latitudes in deg N

• LON(time,NPARTS): longitudes in deg E

• THETA|ZETA(time,NPARTS): potential temperature in K or
  PRESS: pressure in hPa

It also contains the variables specified in the configuration file (line 17). Possible output variables are TEMP, U, V, PV, SZA and for isentropic files also PRESS.

Further variables can be added to the trajectory output file using the program traj_add.

If there are different starting times for the trajectories the data set also contains variables with the prefix start_ (start_time, start_LAT, start_LON etc.). These are arrays with the dimension NPARTS containing the values at the starting times of the trajectories.

If the trajectories end at different times, the values at the end times will be written on variables with the prefix end_ (end_time, end_LAT etc.).

Example for trajectory output file

• History of Changes