External, initial and boundary data. Rodica Claudia DUMITRACHE Amalia IRIZA

Transcription

1 External, initial and boundary data Rodica Claudia DUMITRACHE Amalia IRIZA

2 Necessary Data to run the System COSMO topography ICON/ GME analysis/forecast int2lm COSMO analysis/lbc data ICON/ GME topography

3 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.

4 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). 5. The vertical resolution by specifying the vertical coordinate parameters.

5 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. The Vertical Grid: the number of the vertical levels must be given together with their position in the atmosphere.

6 The Model Domain Example of 50 vertical grid for extra-tropical regions Example of 50 vertical grid for tropical regions:

7 Coarse Grid Model Data Possible Driving Models With the INT2LM it is possible to interpolate data from several driving models to the COSMO grid. 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 used 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 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)

8 Coarse Grid Model Data Specifying the 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. 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. 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

9 Specifying Data Characteristics the Namelist group DATA. - some additional specifications of the data have to be specified. - every data set you can specify the format and the library to use for I/O by setting the appropriate Namelist variables to: grb1 to read or write GRIB1 data with DWD GRIB1 library. apix to read GRIB1 or GRIB2 data with ECMWF grib_api. api1 to write GRIB1 data with ECMWF grib_api. api2 to write GRIB2 data with ECMWF grib_api. ncdf 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 COSMO-Model domain: ie_ext=965, je_ext=773 size of external parameters ylmext_lfn= cosmo_d5_07000_965x773.g1_ 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 )

10 Specifying Data Characteristics External Parameters for the Coarse Grid Input Model: 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_ nc, name of external parameter file for global ICON domain yinext_lfn= icon_extpar_<name-of-region>_r03b07_ nc, name of external parameter file for a region yinext_cat= /tmp/data/ext2/, directory of this file yinext_form_read= grb1, to specify the data format ( or apix or ncdf )

11 Specifying Data Characteristics Coarse Grid Input Data: The 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 ( only for COSMO and ICON GRIB2) 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: 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

12 Parameters to Control the INT2LM Run the Namelist group CONTRL 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 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 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)

13 Parameters to Control the INT2LM Run lforest provide additional external parameters for the distribution of evergree and deciduous forest 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) lbalance_pp to compute a hydrostatic balanced pressure deviation (after vertical interpolation)

14 Parameters to Control the INT2LM Run Example of Namelist Variables for /CONTRL/

15 Interpolation using the ICON data On Thursday, January 15, :31 PM, Ulrich Schaettler ( information) Changes to data sets distributed earlier: There are 3 constant files now for running the interpolation ICON-2-COSMO with INT2LM: icon_extpar_<your-domain-name>_r03b07_ nc: external ICON 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) 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)

16 Interpolation using the ICON data On Thursday, January 15, :31 PM, Ulrich Schaettler ( information) We are NOT distributing the full 90 ICON levels, but only the lower levels, which are 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/). 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. GME_2_ICON.pdf -> detailed information

17 Example of INPUT file &CONTRL ydate_ini=' ', ydate_bd=' ', 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, /

18 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., ivctype=2, irefatm=2, delta_t=75.0, h_scal= , vcflat= , vcoord_d= , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , 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 /

19 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_ nc', yin_hhl='icon_hhl_romania_r03b07.g2' yinext_form_read='ncdf', yinext_cat='/export/home/cosmo/cosmo_operativ/cosmo_7km/const', 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 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, /

20 Thank you!