The package i3d

The program i3d.f90 can be used for interpolation of the CLaMS output (usually stored on irregular grid in init*.nc) on an arbitrary different irregular grid (e.g. experimental data). In particular, it can be applied for interpolation on a spatially highly resolved regular grid that, on the one side, does not destroy the fine structures in the spatial distributions calculated with CLaMS and, on the other side, makes easier the handling of the CLaMS output (e.g. plotting).

How to get the source code

Latitude range of the (external) grid where interpolations have to be determined. The source code of this package 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 {{{Start triangle for the walking triangle algorithm: y/n -- north/south hemisphere (does not work. Now, the n/2 (n-number of all triangles) is chosen as a start triangle). cvs co -P i3d }}} creates a directory i3d in the current working directory and copies the current version of the i3d package to this directory.

The directory i3d includes the following subdirectories:

Compiling

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

gmake i3d

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

Configuration file

The options for the execution of i3d are set in the configuration file i3d.inp: 

1)  y            ! north_hemisphere (y/n)
2)  30., 90.     ! lat_min, lat_max
3)  3            ! 0=nearest, 1=weighted, 2=both, 3=weighted_new
4)  2        ! number of tags
    CH4 
    O3
5)  n                                          ! use one init-file
6)  /private/CNx_ECN_3d_24_1.5                 ! input directory (init-files)
    CNx                                        ! prefix of init-files (3 character)
7)  data/mipas_isen_noonpos_030103_030131.nc   ! input file (aircraft/satellite data)
8)  outdir                                     ! output directory
9)  testout                                    ! output prefix
10) y                                          ! write warnings  (y/n) 

1)  y            ! north_hemisphere (y/n)
2)  30., 90.     ! lat_min, lat_max
3)  3            ! 0=nearest, 1=weighted, 2=both, 3=weighted_new
4)  2            ! number of tags
    CH4 
    O3
5)  y                            ! use one init-file
6)  test1/init_BN0_02100112.nc   ! init-file
7)  test1/pos_A_0_01100112.nc    ! input file (aircraft/satellite data)
8)  test1                        ! output directory
9)  testout                      ! output prefix
10) n                            ! write warnings  (y/n)

Remarks:

  1. Start triangle for the walking triangle algorithm:

    y/n -- north/south hemisphere (does not work. Now, the n/2 (n-number of all triangles) is chosen as a start triangle).

  2. Latitude range of the (external) grid where interpolations have to be determined.
  3. Interpolation type
  4. Number of species and their names
  5. ...
  6. Name of init file | Directory with CLaMS output and prefix of these files
  7. Input file with the external grid
  8. Output directory
  9. Name of the NetCDF output file

The walking triangle algorithm