The program pos_chem
The program pos_chem modifies chemical fields (O3, HCL etc) in init files on lower and upper bound. The values are interpolated from climatology data sets.
The options for the execution of pos_chem can be set in the configuration file pos_chem.inp:
init_BNO_01100112.nc ! init file eqlat ! coordinate1 (lat, eqlat) theta ! coordinate2 (theta, zeta) /usr/nfs/sat_data/haloe ! directory for clim files 3 ! number of species ! species Clams, species Clim, lower bound, upper bound, clim filename HCl HCl 0. 500. haloe_twod_climat_acpt.nc O3 O3 0. 500. haloe_twod_climat_acpt.nc TRACER O3 0. 500. haloe_twod_climat_acpt.nc
In this example the species HCL, O3 and TRACER are replaced with values from climatology data sets below theta=0 and above theta=500 K. The interpolation is made in eqlat-theta-space (coordinate 1 and 2). The init file must contain the specified coordinates eqlat and theta (e.g. added with pos_dyn).
The used climatology data set is /usr/nfs/sat_data/haloe/haloe_twod_climat_acpt.nc.