Documentation of the Changes in INT2LM
Version 2.08

30.10.2020

1. Renamed S_ORO to S_ORO_MAX

(by Ulrich Schättler)

The new hydrology scheme with the modified ground water runoff now has been implemented to the official COSMO code. In this connection an official GRIB code entry has been defined for the slope of the orography and an official shortname has been given, which now is S_ORO_MAX. Therefore the name S_ORO, used up to now in INT2LM, has been changed to S_ORO_MAX.

Back to Contents

2. Corrections for computing difference between initial date and boundary start date

(by Denis Blinov)

There were problems when running INT2LM with driving data that have an older reference date than the actual start date, if the difference between these dates is less than an hour (e.g. 15 minutes).

In src_read_coarse_grid.f90, the difference in minutes is calculated by the routine diff_minutes. Then the difference is computed in units of 15 minutes (shortest time interval that could be used in practice in GRIB1) with the command:

ndiff_ini_bd = - NINT( (imindif * 60 + (ibdsec - inisec)) / dt ) but a difference in seconds will be "deleted" by the integer arithmetic used. So it is sufficient to have ndiff_ini_bd = - NINT( (imindif * 60) / dt ).

Then the dubious correction was added "to account for seconds"
ndiff_ini_bdref = NINT((inimin * 60) / dt) ndiff_ini_bd = ndiff_ini_bd - ndiff_ini_bdref

Reason for that once was, to be able to have an offset between initial date and reference date of the boundary fields of less than 15 minutes. But the way it is done here is not correct:

For a reference date for the boundary fields, inimin is always 0, so ndiff_ini_bdref is 0, and ndiff_ini_bd is not changed. But for a shorter offset (e.g. only 15 minute), wrong values are computed here. These two lines just have to be removed.

Note, that this kind of counting in steps of 15 minutes does not allow to have an offset between initial date and reference date of the boundary fields of less than 15 minutes.

The same corrections are done in files src_icon_interpol.f90 and src_gme_interpol.f90.

Back to Contents

3. Corrections for using HadGEM data and new vertical interpolation

(by Burkhardt Rockel)

INT2LM crashed when using HadGEM global data as input in case of lcm2lm=.TRUE. and itype_profiles_vert_interp=2. The reason was that wrong values of two parameters for setting height levels have been used for vertical gradient computations in the call to vert_interp in src_vert_inter_lm.f90. This has been corrected now.

The resetting of itype_profiles_vert_interp to 1 for lcm2lm in src_namelists.f90 has been removed.

Back to Contents

4. Technical Changes and Bug Fixes

external_data.f90: (by Denis Blinov)
Print outs for Kinne variables only for itype_aerosol=3
src_2d_fields.f90: (by Denis Blinov)
Splitted itype_aerosol=3 from a block with itype_aerosol=2/5
src_cleanup.f90: (by Denis Blinov)
Deallocate Kinne variables only for itype_aerosol=3
src_gribtabs.f90: (by Burkhardt Rockel)
In case lcm2lm=.TRUE. the input of T_SNOW is optional, however, the input of W_SNOW is not. This stops the execution of INT2LM. Since W_SNOW is not urgently needed in climate mode, as it is build by COSMO implicitly, W_SNOW can get the "O" for optional reading.
src_lm_fields.f90: (by Lucio Torrisi)
Calculation of full pressure also if lgrh=.False.
src_lm_output.f90: (by Denis Blinov)
Corrected the actual forecast date ydate in SR write_nc_gdefs (lmmss_ini and .NOT. lmmss_ini had been swapped erroneously)
src_lm_output.f90 and data_int2lm_io.f90: (by Guy de Morsier)
To prevent differences between HSURF and HHL of + and -0.125 m use the packrate variable nrbit_phhl=24 to also pack hsurf.
src_read_hhl.f90: (by Lucio Torrisi)
Less strict test for lat/lon when reading HHL file: The allowed difference in lat or lon has been increased from 1e-5 to 1e-4.

Back to Contents

5. Changes of the Namelists

NONE

Back to Contents

6. Changes of the Results

NONE. If not using affected namelist settings described above.

Back to Contents