Documentation of the Changes in INT2LM
Version 1.20

31.08.2012

This is the next step towards the next big release INT2LM 2.0.

1. Changes to subroutines computing the Reference Atmospheres

(by Michael Baldauf)

The same changes that have been done in COSMO-Model 4.24 have now been implemented in the INT2LM.

They have been extended to also compute the temperature t0hl of the reference state on half levels.
The computation of the k-indices for special pressure levels has been removed from these subroutines and have been put to a new subroutine k_index_of_pressure_levels (also in meteo_utilities).
For irefatm=1 the pressure and temperature of the reference state on full levels can now be computed analytically, instead of just averaging the values from the half levels. For that purpose, a new Namelist switch lanalyt_calc_t0p0 has been introduced in the group /LMGRID/ (with Default=.FALSE.).
If the COSMO-Model should run with the new fast-waves solver, it is necessary to explicitly set this variable lanalyt_calc_t0p0=.TRUE. This corresponds to the setting of itype_fast_waves=2 in the COSMO-Model!
For irefatm=2 (subroutine reference_atmosphere_2) t0 is now calculated directly as a function of the height (instead of using p0), which is a bit more intuitive.

The interfaces in external_data.f90 and src_read_coarse_grid.f90 have been adapted accordingly.

Back to Contents

2. Unification of Utility Modules with latest version of COSMO-Model (4.25)

(by Ulrich Schätter)

io_utilities.f90
There have been changes in the subroutine make_fn to allow for a 10- or a 14-digit date variable. There have also been some corrections to update the string for date appropriately.
All interfaces of make_fn in src_gme_interpol.f90, src_read_coarse_grid.f90, and src_lm_output.f90 have been adapted accordingly.
meteo_utilities.f90
Changes in the subroutines reference_atmosphere_xx (see above).
utilities.f90
Adapted the length of character strings in subroutine get_utc_date and compute also the minutes and seconds. There have been necessary changes and adaptations of all date-variables in the modules:
- clm_utilities.f90
- int2lm_org.f90
- src_2d_fields.f90
- src_gme_interpol.f90
- src_read_coarse_grid.f90
- src_lm_output.f90

Back to Contents

3. Adaptations for sub-hourly analysis update with LETKF

(by Ulrich Schättler)

When using the LETKF method for data assimilation in the future, there will be the need to have an assimilation cycle of only few minutes, and not of hours. Therefore the COSMO-Model has been adapted to read and write files, where also minutes and seconds are included. This has now also been implemented in the INT2LM. Now laf-files or lbf-files (when using ytunit_in, ytunit_out) can have the form

lbfd20130504121500 for 04th of May 2013, 12 UTC + 15 Minutes and 00 seconds
instead of the hitherto existing lbfd2013050412

The main changes are:

Subroutine make_fn from io_utilities.f90
The file names, that are created by this subroutine, can have 4 additional digits.
This is controlled by the Namelist switch ydate_ini for the output data of INT2LM (i.e. for the lbf-files):
- If ydate_ini is specified with 10 digits (as it was up to now), also the file names are as they were before.
- If ydate_ini is specified with 14 digits (including minutes and seconds), the file names for the output files contain 4 additional digits (laf-files and lbf-files with ytunit_out='d').
For the input-files of INT2LM this is controlled by the Namelist switch ydate_bd:
- If ydate_bd is specified with 10 digits (as it was up to now), also the file names for the input are as they were before.
- If ydate_bd is specified with 14 digits (including minutes and seconds), the file names for the input files contain 4 additional digits.
Subroutine get_utc_date: to also compute the minutes and seconds in a precise manner, the computation of the return parameter acthour has been modified, and therefore the results are slightly different.

The date variables that are returned by the routine get_utc_date from module utilities have been extended:

Variable	Before	After
`yactdate1`	`2013050412`	`20130504121532`
`yactdate2`	`SAT 04.05.2013 12 UTC`	`SAT 04.05.2013 12:15:32 UTC`

This also changes some formats in the ASCII output files.

NOTE:
Most modules have been adapted to the new format of the date variables.
BUT all routines dealing with meta data of NetCDF I/O have not been adapted!!!

Back to Contents

4. Changes for NetCDF I/O

(by Burkhardt Rockel)

Global attributes for NetCDF I/O have been added according to the COSMO-Model, which can be set via namelist in the group /DATA/:

Group	Name	Meaning	Default
`/DATA/`	`yncglob_institution`	originating center name	`'-'`
	`yncglob_title`	title string for the output	`'-'`
	`yncglob_source`	program name and version	`'-'`
	`yncglob_project_id`	identification of the project of simulation	`'-'`
	`yncglob_experiment_id`	identification of the experiment of simulation	`'-'`
	`yncglob_contact`	contact e.g. email address	`'-'`
	`yncglob_references`	URL, report etc.	`'-'`
	`ncglob_realization`	number of the realization of the experiment	`'-999'`

Back to Contents

5. Technical Changes

data_int2lm_io.f90:
Increased nmaxlist to 60 to be able to allow all possible new interpolation features.
external_data.f90:
Adapted interfaces to subroutines reference_atmosphere_xx;
Allocation of zhsurfs_lm in any case, because it is used in subroutine calls.
io_utilities.f90:
Correction in case of east_add_in, west_add_in, south_add_in, north_add_in is not 0.
src_cleanup.f90:
Allocation of ztotal_times on all PEs (not only on PE 0), because it is used in all PEs.
src_coarse_interpol.f90:
Removed unused variables from USE-lists and added some more debug prints.
src_gribtabs.f90:
Corrected standard name for W;
Changes to read T_S for both climate and forecast versions in case of netCDF input.
src_lm_output.f90:
Introduction of additional global attributes in case of netCDF output, which can be set via namelist (analog the attributes in the namelist in COSMO);
Several corrections regarding netCDF output.
src_memory.f90:
Allocation of lmask_lm in any case, because it is used as an argument to subroutine calls.
src_namelists.f90:
Introduction of additional global attributes in case of netCDF output which can be set via namelist (analog the attributes in the namelist in COSMO);
Change print statements to account for new namelist parameter yinput_model
Changes to write T_S for both climate and forecast versions in case of netCDF output.
src_read_coarse_grid.f90:
Correction in case of east_add_in, west_add_in, south_add_in, north_add_in is not 0.
Added a grib_set 'missing_value' to allow for successive retrieve of a field with the desired value for missing data, necessary in case of frames.
src_read_ext.f90:
Correction in case of east_add_in, west_add_in, south_add_in, north_add_in is not 0.

Back to Contents

6. Changes of the Namelists

There were the following changes for the Namelist variables:

Group Name Meaning Default

/LMGRID/

lanalyt_calc_t0p0

If it is set to .TRUE., the values for the reference state of t0 and p0 are computed analytically. If it is set to .FALSE., they are only averaged between the half levels.
There is the following dependency between INT2LM and the COSMO-Model regarding this switch:

COSMO-Model	INT2LM
Running with Leapfrog	`lanalyt_calc_t0p0 = .FALSE.` is a necessary condition.
Running with old Runge-Kutta fast-waves (`itype_fast_waves=1`	`lanalyt_calc_t0p0 = .FALSE.` is a necessary condition.
Running with new Runge-Kutta fast-waves (`itype_fast_waves=2`	`lanalyt_calc_t0p0 = .TRUE.` is a necessary condition.

'.FALSE.'

/DATA/

yinext_form_read yin_form_read

The value to choose the grib_api for packing / unpacking GRIB data has been changed to be conform with other models at DWD. The old value was 'grax', now it is 'apix'.

'grb1': to choose the DWDLIB for Grib1
'apix': to choose the grib_api for Grib1 or Grib2
'ncdf': to choose NetCDF input

'grb1'

yncglob_institution originating center name '-'

yncglob_title title string for the output '-'

yncglob_source program name and version '-'

yncglob_project_id identification of the project of simulation '-'

yncglob_experiment_id identification of the experiment of simulation '-'

yncglob_contact contact e.g. email address '-'

yncglob_references URL, report etc. '-'

ncglob_realization number of the realization of the experiment '-999'

Back to Contents

7. Changes of the Results

There are no changes of the results in the NWP applications!

Back to Contents