Concentration field output

If required, the 4D concentration fields are written at an hourly resolution. The values written are the concentrations obtained at the end of the specified hour i.e. they are not the average over this hour. This file is in the following the format with PP, NZONAL, NMERID, NIVOUT and NVERTI given the relevant values, NS the total number of species to write with their UNITS (default is ppb vol) and reference (WD = dry or wet) as indicated in OUTPUT_SPECIES.

dimensions:
	Time = UNLIMITED ; // (PP+1 currently)
	DateStrLen = 19 ;
	west_east = NZONAL ;
	south_north = NMERID ;
	bottom_top = NIVOUT ;
	vcoord_dim = NVERTI ;
variables:
	double lon(south_north, west_east) ;
		lon:units = "degrees_east" ;
		lon:long_name = "Longitude" ;
	double lat(south_north, west_east) ;
		lat:units = "degrees_north" ;
		lat:long_name = "Latitude" ;
	float a_vcoord(vcoord_dim) ;
		a_vcoord:units = "no_units" ;
		a_vcoord:long_name = "A_sigma_coefficient" ;
	float b_vcoord(vcoord_dim) ;
		b_vcoord:units = "no_units" ;
		b_vcoord:long_name = "B_sigma_coefficient" ;
	char Times(Time, DateStrLen) ;
	float S1(Time, bottom_top, south_north, west_east) ;
		S1:units = "UNITS WD" ;
		S1:long_name = "S1 Concentration" ;
        .
        .
        .
	float SNS(Time, bottom_top, south_north, west_east) ;
		SNS:units = "UNITS WD" ;
		SNS:long_name = "SNS Concentration" ;