External, initial and
boundary data
Rodica Claudia DUMITRACHE
Amalia IRIZA
Necessary Data to run the System
COSMO
topography
ICON/ GME
analysis/forecast
int2lm
ICON/ GME
topography
COSMO
analysis/lbc data
Namelist Input for the INT2LM
The execution of INT2LM can be controlled by 6 NAMELIST groups:
•CONTRL parameters for the model run
•GRID_IN specifying the domain and the size of the driving model grid
•LMGRID specifying the domain and the size of the LM grid
•DATA controlling the grib input and output
•PRICTR controlling grid point output
Obs:
- All NAMELIST groups have to appear in the input file INPUT in the order given above.
- Default values are set for all parameters, you only have to specify values, that have to be
different from the default.
External Data
For this Training we provide special data sets for Europe, Africa, Asia and South America. Only
for Europe we provide several resolutions, for Africa, Asia and South America we only provide
data sets with 7 km resolution.
The Model Domain
LMGRID - specify the domain of your choice, you have to choose the following parameters:
1. The lower left grid point of the domain in rotated coordinates.
2. The horizontal resolution in longitudinal and latitudinal direction in degrees.
3. The size of the domain in grid points.
4. The geographical coordinates of the rotated north pole, to specify the rotation
(or the angles of rotation for the domain).
necessary
necessary
5. The vertical resolution by specifying the vertical coordinate parameters.
optional
FE13 / A. Smalla
The Model Domain
The Horizontal Grid:
The COSMO-Model uses a rotated spherical coordinate system, where the poles are moved and
can be positioned such that the equator runs through the center of the model domain.
necessary
necessary
The Vertical Grid:
the number of the vertical levels must be given together
with their position in the atmosphere.
optional
FE13 / A. Smalla
The Model Domain
Example of 50 vertical grid for extra-tropical regions
necessary
necessary
Example of 50 vertical grid for tropical regions:
optional
FE13 / A. Smalla
Coarse Grid Model Data
Possible Driving Models
With the INT2LM it is possible to interpolate data from several driving models to the COSMO
grid.necessary
Up to now the following input models are supported in NWP mode:
• ICON: New non-hydrostatic global model from DWD with an icosahedral model grid
• GME: Old hydrostatic global model from DWD with an icosahedral model grid
• IFS: Global spectral model from ECMWF; data is delivered on a Gaussian grid
• COSMO: Data from a coarser COSMO-Model run can be usednecessary
to drive higher resolution runs
Specifying Regular Grids: GRID_IN
•startlat_in_tot latitude of lower left grid point of total input domain
•startlon_in_tot longitude of lower left grid point of total input domain
•dlon_in longitudinal resolution of the input grid
•dlat_in latitudinal resolution of the input grid optional
•ie_in_tot longitudinal grid point size
•je_in_tot latitudinal grid point size
•ke_in_tot vertical grid point size
•pollat_in geographical latitude of north pole (if applicable; or angle of rotation)
•pollon_in geographical longitude of north pole (if applicable; or angle of rotation)
FE13 / A. Smalla
Coarse Grid Model Data
Specifying
necessarythe ICON Grid:
The parameters that have to be specified for running with ICON data are:
•yicon_grid_cat directory of the file describing the horizontal ICON grid.
•yicon_grid_lfn name of the file describing the horizontal ICON grid.
•ke_in_tot ke (vertical dimension) for input grid.
necessary
•nlevskip number of levels that should be skipped for the input model.
At the moment DWD runs a resolution of about 13 km for the global ICON grid. If the global data
are used, the parameter yicon_grid_lfn has to be specified to
yicon_grid_lfn = icon_grid_0026_R03B07_G.nc,
where the "G" marks the global grid.
optional
For NWP partners, who run the COSMO-Model for daily operational productions, DWD provide
special grid files for the corresponding regional domain. The file names of these grid files contain
the name of the region and look like:
yicon_grid_lfn = icon_grid_<name-of-region>_R03B07.nc
FE13 / A. Smalla
Specifying Data Characteristics
the Namelist group DATA.
- some
additional specifications of the data have to be specified.
necessary
- every data set you can specify the format and the library to use for I/O by setting the
appropriate Namelist variables to:
’grb1’
’apix’
’api1’
’api2’
’ncdf’
to read or write GRIB1 data with DWD GRIB1 library.
to read GRIB1 or GRIB2 data with ECMWF grib_api.
necessary
to write GRIB1 data with ECMWF grib_api.
to write GRIB2 data with ECMWF grib_api.
to read or write NetCDF data with the NetCDF library.
External Parameters
The following variables have to be specified for the dataset
of the external parameters for the
optional
COSMO-Model domain:
ie_ext=965, je_ext=773 size of external parameters
ylmext_lfn=’cosmo_d5_07000_965x773.g1_2013111400’ name of external parameter file
ylmext_cat=’/tmp/data/ext1/’ directory of this file
ylmext_form_read=’grb1’ to specify the data format (or: ’ncdf’, ’apix’)
FE13 / A. Smalla
Specifying Data Characteristics
External
Parameters for the Coarse Grid Input Model:
necessary
The INT2LM also needs the external parameters for the coarse grid driving model. The
following variables have to be specified for this dataset:
•yinext_lfn=’icon_extpar_0026_R03B07_G_20141202.nc’, name
of external parameter file
necessary
for global ICON domain
•yinext_lfn=’icon_extpar_<name-of-region>_R03B07_20141202.nc’,
parameter file for a region
name
•yinext_cat=’/tmp/data/ext2/’, directory of this file
•yinext_form_read=’grb1’, to specify the data optional
format ( or ’apix’ or ’ncdf’)
FE13 / A. Smalla
of
external
Specifying Data Characteristics
Coarse Grid Input Data:
The necessary
following variables have to be specified for the input data sets of the driving model:
• yin_cat=’/tmp/data/input/’, directory of input data
• yin_hhl=’/tmp/data/input/’, name of the vertical grid HHL-file
(necessary only for COSMO and ICON GRIB2)
necessary
• yin_form_read=’grb1’, to specify the data format (or ’apix’ or ’ncdf’)
• ytunit_in=’f’, to specify the time unit and hence the file name of the input data
• yinput_type=’forecast’, to specify the time unit and hence the file name of the input data
COSMO Grid Output Data:
optional
The following variables have to be specified for the output data sets of the COSMO-Model:
• ylm_cat= ’/tmp/data/out/’, directory of the results
• ylm_form_write=’grb1’, data format for output data (or ’api1’, ’ncdf’)
• ytunit_out=’f’, to specify the time unit and hence the file name of the output data
FE13 / A. Smalla
Parameters to Control the INT2LM Run
the Namelist group CONTRL
necessary
Forecast Time, Domain Decomposition and Driving Model
To specify the date and the length of the forecast
• ydate_ini date when the forecast begins
• hstart start hour relative to the start of the forecast
necessary
• hstop stop hour relative to the start of the forecast
• hincbound increment for boundary update given in hours
To specify which data should be computed
• linitial compute initial data
• lboundaries compute boundary data
optional
To specify the driving model
• yinput_model the name of the driving model is specified as a character string. The
following names are possible: ’ICON’, ’GME’, ’IFS’, ’COSMO’, ’CM’.
To specify the number of processors
• nprocx number of processors in X-direction (longitudinal)
• nprocy number of processors in Y-direction (latitudinal)
• nprocio number of processor for asynchronous IO (usually 0)
FE13 / A. Smalla
Parameters to Control the INT2LM Run
•
lforest provide additional external parameters for the distribution of evergree and
deciduous forest
necessary
•
lsso provide additional external parameters to run the subgrid scale orography scheme
•
lprog_qi interpolate cloud ice from the driving model
•
lprog_qr_qs interpolate rain and snow from the driving model
•
lmulti_layer_in use multi-layer input soil levels
•
lmulti_layer_lm compute multi-layer output soil levels
•
llake provide additional external parameters for the FLake model
•
lprog_rho_snow interpolate prognostic snow density from the driving model
•
lfilter_oro to apply a filter for the orography
•
lfilter_pp to apply a filter for the pressure deviation (after vertical interpolation)
•
necessary
optional
lbalance_pp to compute a hydrostatic balanced pressure deviation (after vertical
interpolation)
FE13
/ A. Smalla
Parameters to Control the INT2LM Run
Example of Namelist Variables for /CONTRL/
necessary
necessary
optional
FE13 / A. Smalla
Interpolation using the ICON data
On Thursday, January 15, 2015 2:31 PM, Ulrich Schaettler (email information)
necessary
Changes to data sets distributed earlier:
•There are 3 constant files necessary now for running the interpolation ICON-2-COSMO
with INT2LM:
•icon_extpar_<your-domain-name>_R03B07_20141202.nc:
external
ICON
necessary
parameters based on the used ICON grid (from 02nd December 2014) and your domain
(NetCDF file). Note that these parameters have been updated and are really different to the
ones you got before!
•icon_grid_<your-domain-name>_R03B07.nc: external grid file to specify the ICON
grid points covering your domain (NetCDF file)
•icon_hhl_<your-domain-name>_R03B07.g2: external file to specify the heights of
the half levels (hhl) for the ICON grid points covering your domain (GRIB2 file)
optional
Note, that we are NOT distributing these 3 files daily, so you have to take them from the test
data set and store them for your operational runs (best in the CONST-directory for the input
model)
FE13 / A. Smalla
Interpolation using the ICON data
On Thursday, January 15, 2015 2:31 PM, Ulrich Schaettler (email information)
necessary
•
We are NOT distributing the full 90 ICON levels, but only the lower levels, which are
necessary to run the COSMO-Model, which is not going as high in the atmosphere, as
ICON. For that we introduced a Namelist variable in INT2LM: nlevskip. For your
applications you have to set nlevskip=20 (in group /GRID_IN/).
necessary
•
At the moment we cannot write proper UUIDs (unique universal identifiers) for the ICON
data, therefore the namelist variable lcheck_uuidOfHGrid has to be set to .FALSE.
(in group /CONTRL/). Later, these UUIDs will be used to check that the ICON forecast
data corresponds properly to the 3 ICON constant data sets.
optional
GME_2_ICON.pdf -> detailed information
FE13 / A. Smalla
Example of INPUT file
&CONTRL
ydate_ini='2016021500', ydate_bd='2016021500',
hstart=0.0, hstop=78.0, hincbound=3.0,
linitial=.TRUE., lboundaries =.TRUE.,
nprocx=4, nprocy=4, nprocio=0, lreorder=.FALSE.,
yinput_model='ICON',
ltime_mean=.TRUE.,
lfilter_oro=.TRUE., eps_filter=0.1,
ilow_pass_oro=1, ilow_pass_xso=0, rxso_mask=0.0,
lfilter_pp=.FALSE., lbalance_pp=.FALSE., norder_filter=5,
luvcor=.TRUE.,
lmulti_layer_in=.TRUE., lmulti_layer_lm=.TRUE.,
lprog_qi=.TRUE., lprog_qr_qs=.TRUE.,
lsso=.false., lforest=.TRUE., lbdclim=.FALSE., lt_cl_corr=.TRUE.,
llake=.FALSE., llake_coldstart=.FALSE.
lseaice=.TRUE., lprog_rho_snow=.TRUE.,
itype_albedo=1,
itype_ndvi=0, idbg_level=2,
leps_bc=.FALSE.,
lcheck_uuidOfHGrid=.false.,
lasync_io=.FALSE.,
lvertwind_bd=.TRUE.,
/
&GRID_IN
ni_gme = 384, i3e_gme = 91, kcontrol_fi =15, ke_soil_in=7,
ke_in_tot=90,
nrootdiv_icon=03, nbisect_icon=07,
yicon_grid_cat = '/export/home/cosmo/COSMO_OPERATIV/COSMO_7km/CONST',
yicon_grid_lfn = 'icon_grid_romania_R03B07.nc'
nlevskip=20,
FE13 // A. Smalla
Example of INPUT file
&LMGRID
startlat_tot = -16.5, startlon_tot = 4.0,
pollat=32.5,
pollon=-170.0,
dlon=0.0625,
dlat=0.0625,
ielm_tot=201,
jelm_tot=177,
kelm_tot=40,
ke_soil_lm=7,
lanalyt_calc_t0p0=.TRUE.,
necessary
ivctype=2, irefatm=2, delta_t=75.0, h_scal=10000.0,
vcflat=11430.0,
vcoord_d=22700.0, 20800.0, 19100.0, 17550.0, 16150.0, 14900.0, 13800.0, 12785.0,
11875.0, 11020.0, 10205.0, 9440.0, 8710.0, 8015.0, 7355.0, 6725.0, 6130.0,
5565.0, 5035.0, 4530.0, 4060.0, 3615.0, 3200.0, 2815.0, 2455.0, 2125.0,
1820.0, 1545.0, 1295.0, 1070.0, 870.0, 695.0, 542.0, 412.0, 303.0,
214.0, 143.0, 89.0, 49.0, 20.0, 0.0,
lnewVGrid=.true.,
/
&DATABASE
optional
/
FE13 / A. Smalla
Example of INPUT file
&DATA
ie_ext=961, je_ext=769,
ylmext_lfn='lm_d1_07000_961x769_new.g1',
ylmext_cat='/export/home/cosmo/COSMO_OPERATIV/COSMO_7km/CONST',
ylmext_form_read='apix',
ylm_hhl='COSMO_EU_HHL_002_041_001.g2'
yinext_lfn='icon_extpar_romania_R03B07_20141202.nc',
yin_hhl='icon_hhl_romania_R03B07.g2'
yinext_form_read='ncdf',
yinext_cat='/export/home/cosmo/COSMO_OPERATIV/COSMO_7km/CONST',
necessary
yin_form_read='apix',
yin_cat='/export/home/cosmo/COSMO_OPERATIV/COSMO_7km/DATA00',
ylm_form_write='grb1',
ylm_cat='/export/home/cosmo/COSMO_OPERATIV/COSMO_7km/DATA00'
nprocess_ini = 134, nprocess_bd = 135,
nl_soil_in=2, nl_soil_lm=2,
l_ke_in_gds=.TRUE.,
/
&PRICTR
igp_tot = 36, 40, 48, 44, 48, 85, 77
optional
jgp_tot = 30, 94, 38, 26, 26, 96, 12
lchkin=.TRUE., lchkout=.TRUE.,
lprps=.FALSE.,
/
&EPSCTL
iepsmem_bc=3,
iepstyp_bc=202,
iepstot_bc=8,
/
FE13 / A. Smalla
Thank you!
necessary
optional
FE13 / A. Smalla