Initial condition inputs¶
The INI_CONCS.0.nc file contains the initial concentration fields in molecules/cm3 on the model’s grid at the beginning of the first hour. Among the species in INI_CONCS.0.nc, only those which are listed in the ACTIVE_SPECIES file of the chemical scheme plus the prescribed species are used. The concentration fields include the “normal” concentrations and the “conco” concentrations (for use in twostep) in case of a chained simulation i.e when the initial conditions are actually the end.nc file of a previous CHIMERE run. If the simulation is not to be chained, the “conco” are ignored. If no concentrations are provided for a given species, its initial field is set at 1d-17 (see
check that no species is set at null initial concentrations due to forgetting it in INI_CONCS.0.nc since no error is then raised.
The initial concentrations provided in dry air mole fractions can be converted into total air mole fractions for use by CHIMERE, following the specifications of the user in the ACTIVE_SPECIES file of the chemical scheme.
See also the fields plugins for CHIMERE for documentation on this file and the information to provide to the CIF to generate it.
The format of INI_CONCS.0.nc is as follows, with NVERTI, NZONAL and NMERID given the relevant values, NS the number of species in the file (the “normal” and “conco” concentrations do not use the same name), UNITS is either ppb, in which case the values are converted into molec/cm3, or anything else, in which case it is assumed that they already are molec/cm3. XX PUT A CHECK OF THIS IN CHIMEREXX
dimensions: Time = 1 ; DateStrLen = 19 ; SpStrLen = 23 ; west_east = NZONAL ; south_north = NMERID ; bottom_top = NVERTI ; Species = NS ; variables: char Times(Time, DateStrLen) ; char species(Species, SpStrLen) ; float lon(south_north, west_east) ; float lat(south_north, west_east) ; double S1(bottom_top, south_north, west_east) ; S1:units = "UNITS" ; . . . double SNS(bottom_top, south_north, west_east) ; SNS:units = "UNITS" ; double S1_o(bottom_top, south_north, west_east) ; S1_o:units = "UNITS" ; . . . double SNS_o(bottom_top, south_north, west_east) ; SNS_o:units = "UNITS" ;
- Important points:
there is only one record so that Time is not a record dimension.
the species names can be as long as SpStrLen = 23 characters including the _o
the concentrations must be written as
floatto ensure enough precision for the adjoint.