Output files

Most of the output files are netcdf. Some choices are available to manage them.

dumpnctype

The precision of the float in the output files. Called NetCDF_outtype in CHIMERE.

dumpncoutput

Choice to dump or not concentrations in out.nc. Called NetCDF_output in CHIMERE.

dumpncpar

Choice to dump or not various variables in par.nc. Called NetCDF_parout in CHIMERE.

nparammax:

Maximum number of miscellaneous variables to be written in par.nc. Experienced users can increase it to add more parameters (see iniparams, outprint).

Deposition output

If required, the dry and wet depositions are written, cumulated every nsavedepos hours, in dep.nc. This file is in the following the format with NH, NZONAL, NMERID, NVERTI and NS given the relevant values, NDEPOS the total number of species to be deposited (counts two for a species that is both dry and wet deposited), NDEPDRY/WET the number of dry/wet deposited species:

dimensions:
	Time = UNLIMITED ; // (NH currently)
	DateStrLen = 19 ;
	SpStrLen = 20 ;
	Species = NDEPOS ;
	west_east = NZONAL ;
	south_north = NMERID ;
variables:
	float lon(south_north, west_east) ;
		lon:units = "degrees_east" ;
		lon:long_name = "Longitude" ;
	float lat(south_north, west_east) ;
		lat:units = "degrees_north" ;
		lat:long_name = "Latitude" ;
	char species(Species, SpStrLen) ;
	char Times(Time, DateStrLen) ;
	float S1_DRY(Time, south_north, west_east) ;
		S1_DRY:units = "molecules/cm2" ;
        float S2_DRY(Time, south_north, west_east) ;
		S2_DRY:units = "molecules/cm2" ;
        .
        .
        .
	double SNDEPDRY(Time, south_north, west_east) ;
		SNDEPDRY:units = "molecule/cm2" ;
	float S1_WET(Time, south_north, west_east) ;
		S1_WET:units = "molecules/cm2" ;
        float S2_WET(Time, south_north, west_east) ;
		S2_WET:units = "molecules/cm2" ;
        .
        .
        .
	double SNDEPWET(Time, south_north, west_east) ;
		SNDEPWET:units = "molecule/cm2" ;

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" ;

Ouput of miscellaneous variables

If required, a list of variables other than the concentrations are written in par.nc. This file is in the following the format with PP, NZONAL and NMERID given the relevant values, the list of parameters can be modified by experienced users in iniparams and outprint, increasing nparammax if necessary:

dimensions:
	Time = UNLIMITED ; // (PP+1 currently)
	DateStrLen = 19 ;
	west_east = NZONAL ;
	south_north = NMERID ;
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" ;
	char Times(Time, DateStrLen) ;
	float winhs(Time, south_north, west_east) ;
		winhs:units = "m/s" ;
		winhs:long_name = "Lower layer horizontal wind" ;
	float winh2(Time, south_north, west_east) ;
		winh2:units = "m/s" ;
		winh2:long_name = "Second layer horizontal wind" ;
	float winvs(Time, south_north, west_east) ;
		winvs:units = "m/s" ;
		winvs:long_name = "Lower layer vertical wind" ;
	float winvt(Time, south_north, west_east) ;
		winvt:units = "m/s" ;
		winvt:long_name = "Top layer vertical wind" ;
	float sphu(Time, south_north, west_east) ;
		sphu:units = "Kg/Kg" ;
		sphu:long_name = "Specific humidity in lower layer" ;
	float airm(Time, south_north, west_east) ;
		airm:units = "molecule/cm**3" ;
		airm:long_name = "Air density in lower layer" ;
	float temp(Time, south_north, west_east) ;
		temp:units = "degC" ;
		temp:long_name = "Air temperature in lower layer" ;
	float winxs(Time, south_north, west_east) ;
		winxs:units = "m/s" ;
		winxs:long_name = "Velocity-equivalent of mixing ex" ;
	float winxt(Time, south_north, west_east) ;
		winxt:units = "m/s" ;
		winxt:long_name = "Velocity-equivalent of mixing ex" ;
	float thlay(Time, south_north, west_east) ;
		thlay:units = "m" ;
		thlay:long_name = "Lower layer thickness" ;
	float hght(Time, south_north, west_east) ;
		hght:units = "m" ;
		hght:long_name = "Boundary layer height" ;
	float atte(Time, south_north, west_east) ;
		atte:units = "0-1" ;
		atte:long_name = "Radiation attenuation factor" ;
	float zeni(Time, south_north, west_east) ;
		zeni:units = "0-1" ;
		zeni:long_name = "Cosine of zenith angle" ;
	float tem2(Time, south_north, west_east) ;
		tem2:units = "degC" ;
		tem2:long_name = "2m air temperature" ;
	float usta(Time, south_north, west_east) ;
		usta:units = "m/s" ;
		usta:long_name = "Friction velocity" ;
	float aerr(Time, south_north, west_east) ;
		aerr:units = "s/m" ;
		aerr:long_name = "Aerodynamic resistance" ;
	float wsta(Time, south_north, west_east) ;
		wsta:units = "m/s" ;
		wsta:long_name = "Convective velocity" ;
	float depo(Time, south_north, west_east) ;
		depo:units = "cm/s" ;
		depo:long_name = "Deposition rate" ;
	float topc(Time, south_north, west_east) ;
		topc:units = "mm/hour" ;
		topc:long_name = "Total precipitation" ;
	float clwc(Time, south_north, west_east) ;
		clwc:units = "kg/Kg" ;
		clwc:long_name = "Cloud liquid water + ice content" ;
	float humi(Time, south_north, west_east) ;
		humi:units = "%" ;
		humi:long_name = "Relative humidity" ;