This directory contains output from STILT footprint computations using
GEOS input files, for sample receptors at the WLEF tower, during
January and July 2013.  The sample receptors are the 6 above ground
levels (AGL, in m: 11, 30, 76, 122, 244, 396) at the WLEF tower
(45.9451 N, 90.2732 W), at 08 UTC and 20 UTC on each day of the
month.

Footprint plots:

 WLEF-AvgJanFootprints.pdf: Maps of footprints, plotted separately by
                            level and time of day, averaged over all
                            receptors during the month of January
                            2013.
                            day1, day2: footprints summed over all
                             hourly periods during day1 and day 2.
                            10day: footprints summed over all
                             hourly periods of the 10-day back trajectory.

 WLEF-AvgJulFootprints.pdf: The same for July 2013.

Footprint tar files:

 FootPlots_plJan2013 - footprints for back trajectory computations for
                       WLEF January 2013, using STILT coupled to GEOS_572
 FootPlots_plJul2013 - footprints for back trajectory computations for
                       WLEF July 2013, using STILT coupled to GEOS_FP

Each tar file contains netcdf files, one for each receptor.
Sample files names are listed below for the first and last day of
January:

 foot2013x01x01x08x45.95Nx090.27Wx00011.nc
 foot2013x01x01x08x45.95Nx090.27Wx00030.nc
 foot2013x01x01x08x45.95Nx090.27Wx00076.nc
 foot2013x01x01x08x45.95Nx090.27Wx00122.nc
 foot2013x01x01x08x45.95Nx090.27Wx00244.nc
 foot2013x01x01x08x45.95Nx090.27Wx00396.nc
 foot2013x01x01x20x45.95Nx090.27Wx00011.nc
 foot2013x01x01x20x45.95Nx090.27Wx00030.nc
 foot2013x01x01x20x45.95Nx090.27Wx00076.nc
 foot2013x01x01x20x45.95Nx090.27Wx00122.nc
 foot2013x01x01x20x45.95Nx090.27Wx00244.nc
 foot2013x01x01x20x45.95Nx090.27Wx00396.nc

 (...)

 foot2013x01x31x08x45.95Nx090.27Wx00011.nc
 foot2013x01x31x08x45.95Nx090.27Wx00030.nc
 foot2013x01x31x08x45.95Nx090.27Wx00076.nc
 foot2013x01x31x08x45.95Nx090.27Wx00122.nc
 foot2013x01x31x08x45.95Nx090.27Wx00244.nc
 foot2013x01x31x08x45.95Nx090.27Wx00396.nc
 foot2013x01x31x20x45.95Nx090.27Wx00011.nc
 foot2013x01x31x20x45.95Nx090.27Wx00030.nc
 foot2013x01x31x20x45.95Nx090.27Wx00076.nc
 foot2013x01x31x20x45.95Nx090.27Wx00122.nc
 foot2013x01x31x20x45.95Nx090.27Wx00244.nc
 foot2013x01x31x20x45.95Nx090.27Wx00396.nc

File contents:

 Each file contains gridded footprints, computed in hourly intervals
 for two separate grids:

  foot1 - a global 1 deg x 1 deg grid (360x180)
  footnearfield1 - a limited-area 0.1 x 0.1 deg grid centered on the
                   receptor (310x210)

 Global footprints are computed for all 10 days of the back
 trajectories (240 hourly values), near field footprint only for the
 first 5 days of the back trajectories.

 An ncdump -h output for a sample file is given below.

funnel[78]% ncdump -h  /scratch/white/day14/trn/RData.plJan2013/foot2013x01x01x08x45.95Nx090.27Wx00011.nc
netcdf foot2013x01x01x08x45.95Nx090.27Wx00011 {
dimensions:
	nchar = 500 ;
	nident = 1 ;
	norig = 1 ;
	foot1lon = 360 ;
	foot1lat = 180 ;
	foot1date = 240 ;
	footnearfield1lon = 310 ;
	footnearfield1lat = 210 ;
	footnearfield1date = 120 ;
variables:
	int nchar(nchar) ;
		nchar:long_name = "nchar" ;
	char ident(nident, nchar) ;
		ident:units = "string" ;
		ident:long_name = "identifier string" ;
	float origagl(norig) ;
		origagl:units = "m" ;
		origagl:missing_value = -1.e+34f ;
		origagl:long_name = "original receptor height above ground" ;
		origagl:standard_name = "height" ;
		origagl:description = "before rounding for STILT" ;
	float origlat(norig) ;
		origlat:units = "degrees_north" ;
		origlat:missing_value = -1.e+34f ;
		origlat:long_name = "original receptor latitude" ;
		origlat:standard_name = "latitude" ;
		origlat:description = "before rounding for STILT" ;
	float origlon(norig) ;
		origlon:units = "degrees_east" ;
		origlon:missing_value = -1.e+34f ;
		origlon:long_name = "original receptor longitude" ;
		origlon:standard_name = "longitude" ;
		origlon:description = "before rounding for STILT" ;
	char origutctime(norig, nchar) ;
		origutctime:units = "UTC" ;
		origutctime:long_name = "original receptor time" ;
		origutctime:description = "in UTC, before rounding for STILT" ;
	char origutctimeformat(norig, nchar) ;
		origutctimeformat:units = "none" ;
		origutctimeformat:long_name = "original receptor time (origutctime) format" ;
	double foot1lon(foot1lon) ;
		foot1lon:units = "degrees_east" ;
		foot1lon:long_name = "foot1lon" ;
		foot1lon:standard_name = "longitude" ;
		foot1lon:description = "degrees longitude of center of grid boxes" ;
	double foot1lat(foot1lat) ;
		foot1lat:units = "degrees_north" ;
		foot1lat:long_name = "foot1lat" ;
		foot1lat:standard_name = "latitude" ;
		foot1lat:description = "degrees latitude of center of grid boxes" ;
	double foot1date(foot1date) ;
		foot1date:units = "days since 2000-01-01 00:00:00 UTC" ;
		foot1date:long_name = "foot1date" ;
	float foot1(foot1date, foot1lat, foot1lon) ;
		foot1:units = "ppm per (micromol m-2 s-1)" ;
		foot1:missing_value = -1.e+34f ;
		foot1:long_name = "gridded stilt footprint" ;
		foot1:description = "aggregates particle footprints on lon/lat/time grid starting at stilt start time (see foothr description), (i.e. foot[,,1] is a lat/lon grid of footprints aggregated in grid boxes and in time prior to the stilt start time up to and including footdate[1], which is foothr[1] hours prior to the stilt start time; foot[,,2] is a lat/lon grid of footprints aggregated in grid boxes and in time prior to footdate[1] up to and including footdate[2], which is foothr[2] hours prior to the stilt start time; etc)" ;
	float foot1hr(foot1date) ;
		foot1hr:units = "hours" ;
		foot1hr:missing_value = -1.e+34f ;
		foot1hr:long_name = "stilt footprint hours back from stilt start time" ;
		foot1hr:description = "stilt start time = hourly ceiling(receptor initialization time in ident) i.e. if ident time=17:01, stilt start time = 18:00" ;
	double footnearfield1lon(footnearfield1lon) ;
		footnearfield1lon:units = "degrees_east" ;
		footnearfield1lon:long_name = "footnearfield1lon" ;
		footnearfield1lon:standard_name = "longitude" ;
		footnearfield1lon:description = "degrees longitude of center of grid boxes" ;
	double footnearfield1lat(footnearfield1lat) ;
		footnearfield1lat:units = "degrees_north" ;
		footnearfield1lat:long_name = "footnearfield1lat" ;
		footnearfield1lat:standard_name = "latitude" ;
		footnearfield1lat:description = "degrees latitude of center of grid boxes" ;
	double footnearfield1date(footnearfield1date) ;
		footnearfield1date:units = "days since 2000-01-01 00:00:00 UTC" ;
		footnearfield1date:long_name = "footnearfield1date" ;
	float footnearfield1(footnearfield1date, footnearfield1lat, footnearfield1lon) ;
		footnearfield1:units = "ppm per (micromol m-2 s-1)" ;
		footnearfield1:missing_value = -1.e+34f ;
		footnearfield1:long_name = "gridded stilt footprint" ;
		footnearfield1:description = "aggregates particle footprints on lon/lat/time grid starting at stilt start time (see foothr description), (i.e. foot[,,1] is a lat/lon grid of footprints aggregated in grid boxes and in time prior to the stilt start time up to and including footdate[1], which is foothr[1] hours prior to the stilt start time; foot[,,2] is a lat/lon grid of footprints aggregated in grid boxes and in time prior to footdate[1] up to and including footdate[2], which is foothr[2] hours prior to the stilt start time; etc)" ;
	float footnearfield1hr(footnearfield1date) ;
		footnearfield1hr:units = "hours" ;
		footnearfield1hr:missing_value = -1.e+34f ;
		footnearfield1hr:long_name = "stilt footprint hours back from stilt start time" ;
		footnearfield1hr:description = "stilt start time = hourly ceiling(receptor initialization time in ident) i.e. if ident time=17:01, stilt start time = 18:00" ;

// global attributes:
		:Title = "Output from STILT" ;
		:Conventions = "CF-1.6" ;
}