A step to the local use of orography in higher resolution : testing configuration e923 on DEC workstation

(more details neva.pristov++at++rzs-hm.si )

The complete version is on delage: ~mrpe691/report/c923/report_923WS.ps

There is an interest to make possible preparation of ALADIN climatic files at remotely located centres and in next step to include local data in higher resolution.

In this article first some basic information about configuration 923 is given, then all required input files for each step are listed. Modifications for running this configuration on workstation are described. Finally the results obtained on workstation are presented.

We can conclude that configuration e923 can run on other platforms. For this few changes in the code are maybe necessary and enough disk space is needed for all required input files.

1.   Basic information about configuration e923¹

Configuration e923 is a procedure, which interpolates climatic information from global fields into an arbitrarily chosen limited area. It creates the so-called clim files, containing the following 2D grid point fields:

• constants describing orography
• constants describing surface, soil, vegetation, which may be monthly dependent fields
• monthly climatological values for soil or surface variables

A clim file is created in 3 to 6 steps:

• step 1 - definition of orography
• step 2 - definition of surface, soil and vegetation characteristics without annual cycle
• step 3 - definition of monthly climatological values, modification of albedo and emissivity according to the climatology sea-ice limit
• step 4 - definition/modification of the vegetation and surface characteristics
• step 5 - modification of fields created by step 2 or 4 over land from high resolution datasets (for each month)
• step 6 - modification of climatological values

The first 4 steps are required for the ISBA surface scheme. For the non-ISBA case the first 3 steps are required and the input files for step 2 differ. Steps 5 and 6 are usually performed.

2.   How to run configuration e923

For each step a special namelist is needed, they are almost the same except for the value of N923 in NAMMCC, which is changing according to step, and in namelist for step 5 information about local database is added in NAMCLI. Each step different files with corresponding data are read. The output fields written in file Const.Clim can be in any projection or in lat./lon. grid. This is determined with the definition of geometry in namelists. The master program is run without options. Steps 4 to 6 have to be run for each month.

For each step different files with corresponding data are required. They are listed in table 1, here follow information from which base data are prepared and where the files on delage can be found.

Table 1: List of all input global fields with some additional information about size of file, number of bytes per data, number of points and resolution of global data. In one box are files for one step for ISBA case.

• Step 1

Datasets are prepared from GLOBE25 data in resolution 2'30''

Input data are on directory:

/home/m/mrpe/mrpe603/RELIEF_G/GLOB95

• Step 2

Datasets are prepared from different sources (Meteosat, NOAA-4, CLIMAP, ISLSCP, climatology from Webb, tables, classifications Wilson &Henderson-Sellers, Mahfouf)

Input data are on directory:

/home/m/mrpe/mrpe603/SURFACE_G/version1/i3e

• Step 3

Datasets are prepared from US Navy data, CLIMAP, climatology from NCAR, AMIP

Input data are on directory:

/home/m/mrpe/mrpe603/N108/i3e

• Step 4

Datasets are prepared from Henderson-Sellers classification + tables

Input data are on directory:

/home/m/mrpe/mrpe603/SURFACE_G/version1/i3e

• Step 5

Datasets are result from the combination of a AVHRR-derived vegetation mapping over Europe and the ESA forest mask. The initial data have a resolution of 2 km but have been interpolated onto a regular 0.1°x0.1° grid. The area is ELATSW=30°, ELONSW=-25°, ELATNE=72°, ELONNE=61°

Input data are on directory:

/home/m/mrpe/mrpe603/SURFACE_L/EUROPEb_v1/i3e

• Step 6

Datasets were created from one-year assimilation experiment performed with ARPEGE and truncation T79. This was necessary to obtain climatological values consistent with ISBA and the new vegetation.

Input data are on directory

/home/m/mrpe/mrpe603/CLIM_G/version1/i3e

3.   Code for DEC workstation (cycle AL09.07)

Changes in code

Besides other changes which are already done for configurations ee927, 001 and fullpos (see reports from Radnoti, Jerman or Spaniel), some additional changes related to reading of binary files are needed. The only problem appears in reading input files for step 1 in eincli1.F because on DEC record length in direct access files is given in units of 4 bytes.

eincli1.F : IARRAY32 is integer*4

by opening binary files RECL should be divided by 4

OPEN(31, FILE='Water_Percentage', ... , RECL=NDATX/4)

Local domain

The first step of configuration e923 can be run also on data set for local domain (LGLOBE=.FALSE.), not just for the data for whole globe (LGLOBE=.TRUE.). In order to be able to run it on local domain with the version Al09 routines gtoptx2.F, gtoptxy.F, gtopty2.F, locmaxi.F,eincli1.F, einter1.F prepared by E. Cordoneanu (December 1996) should be included (phased).

This is very useful for workstation version. Files for step one with data over whole globe are very huge. From them data for domain (a little bigger than for clima file) can be extracted. For detailed descriptions see the complete version.

Spectral fit of orography

In the export version of the code routines for minimizations in variational method are not included. They are needed if LKEYF=.TRUE. when variational method is used for spectral representation of the earth topography. For non commercial use following routines which are needed can be used: m1qn3.F m1qn3a.F ystbl.F mlis0.F dds.F dd.F.

Results

Orography fields obtained on Fujitsu or on workstation are the same in the case without spectral fit of orography. When this method is used there are some differences in values, because of different calculations. On Slovenian domain the highest difference is around 150 m (see figure 1). Other fields in clima files, which are connected to the orography, also show some differences.

Figure 1: Surface topography for the Slovenian domain;

a) from operational clim file created on Fujitsu with spectral fit LKEYF=.T.,

b) from clim file created on workstation without spectral fit (LKEYF=.F.),

c) from clim file created on workstation with spectral fit (LKEYF=.T.),

d) difference between fields c and a

Sources of modified routines for workstation version (gtoptx2.F, gtoptxy.F, gtopty2.F, locmaxi.F, eincli1.F, einter1.F, m1qn3.F m1qn3a.F ystbl.F mlis0.F dds.F dd.F), with example of Makefile and namelists for each step of configuration e923 are on delage: ~mrpe691/report/c923/e923workstation.tar.gz.

4.   How to run configuration e923 on workstation

Transfer all needed input files from delage. If you have problems with space extract subdomain from global data for files needed for step 1.

Prepare namelists (examples can be found on delage in the script ~mrpe691/report/c923/job923_n or inside ~mrpe691/report/c923/e923workstation.tar.gz )

Try to run the model, be careful: steps 4 to 6 have to be done for each month.

If there are problems some changes should be made. You can take subroutines from delage ~mrpe691/report/c923/e923workstation.tar.gz and compile ALADIN program again.

5.   Changes in cycle AL11

Here are listed some information about configuration e923 in cycle AL11.

• NAMCLI need not to be filled since cycle AL11
• Routines needed for minimization are included in export version.
• The changes needed for local domain computation have been introduced in new cycle.
• Some smaller input files are on dino also.

Example of the script with namelist for Fujitsu (from D.Giard) can be found on delage ~mrpe691/report/c923/job923_al11.

Home