## page was renamed from pos_chem.f90
=== 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''.