• pos
• pos_dyn.f90

Program pos_dyn

This program completes a position file (lat, lon, theta/zeta) with some dynamical fields (PV, EQLAT, TEMP, PRESS etc.). These fields are interpolated from isentropic files and added to the original pos file.

The options for the execution of pos_dyn can be set in the configuration file pos_dyn.inp:

pos_XXX_97100412.nc      ! input pos-file
testdir                  ! directory for input/output files
isen_ECN_97100412.nc     ! isentropic file
isendir                  ! its directory
X                        ! eqlat=N/S/A/X for northern/southern/all/none hemisphere
4                        ! number of parameters to add
SH
O3
THETA
PV
2                        ! number of parameters to change
H2O  380 4000
O3  0 500

There are two parameterlists in the configuration file:

1. The parameters are added or replaced in the init file. If O3 or SH are added/replaced they will be named O3_ECMWF and SH_ECMWF.
2. The init file must already contain these parameters. For each parameter two levels must be specified. Below the first and above the second level the parameters are interpolated from isentropic file. The species SH and O3 in isentropic files are stored as mass mixing ratio, but H2O and O3 in init files are stored as volume mixing ratios. Therefor these species are converted before they are written to the init file:

   • SH__{init_file} = SH__{isen_file}

     H2O__{init_file} = (28.9644/18.015) * SH__{isen_file}

     O3__{init_file} = (28.9644/48.0) * O3__{isen_file}