Skip to content

Commit

Permalink
Add MetCoOp L90 defintion and some updated doc (#54)
Browse files Browse the repository at this point in the history 
* Add MetCoOp L90 defintion and some updated doc

* doc markdown correct. Thanks Roel!
  • Loading branch information
@ewhelan
ewhelan authored Nov 26, 2024
1 parent cf81a71 commit 5ea27d0
Show file tree
Hide file tree
Showing 2 changed files with 51 additions and 1 deletion.
15 changes: 14 additions & 1 deletion doc/src/jbdiagnose.md
View file
Original file line number Diff line number Diff line change
Expand Up @@ -45,7 +45,20 @@ plotjbdiag -l 50 -t vercor -p QQ -r ./ -e IRELAND75L65
The development and implementation of the HARMONIE data assimilation are now in a quite advanced stage with pre-operational testing at most of the HIRLAMmember weather services. At the same time, the technical and scientific knowledge about various components of this data assimilation is not so widely spread and a joint effort is needed forthe validation of these locally installed system. One of the important components of the HARMONIE data assimilation is the background error statistics. A tool for diagnostics of the background error statistics directly from the background error statistics files has therefore been developed and applied to the recently developed statistics files for HARMONIE at mesoscale (AROME)resolution (2.5 km grid resolution) at SMHI, met.no, FMI, DMI, KNMI and AEMET. This note describes the diagnostic tool (jbdiagnose) and presents results from a comparison of the different statistics from the AROME implementations. Since this diagnostic software is quite fresh, there may certainly be coding errors that could affect the results. Any comments on the comparison and the graphs presented here are therefore most welcome.

## 2 The diagnostic tool
The first part of the diagnostic tool is a simple standalone fortran program (jbdiagnose.F90) that reads the two main background error statistics files (*.bal and *.cv), calculates various diagnostic quantities and writes these diagnostic quantitiesin simple ASCII files that can be used by, for example, GNUPLOT to produce various graphs. The input background error statistics files are opened with the local file names stabal96.bal and stabal96.cv and a namelist variable gsizein, describing the grid resolution, and namelist variables describing the vertical levels (AHALF, BHALF) are needed. Various information about the domain geometry is extracted from the input files. Oneexample of a script to run the diagnostic tool is provided in jbdiagnose.sh. The output ASCII files are stored in the
The first part of the diagnostic tool is a simple standalone fortran program (`jbdiagnose.F90`) that reads the two main background error statistics files (`jb_filename.bal` and `jb_filename.cv`), calculates various diagnostic quantities and writes these diagnostic quantitiesin simple ASCII files that are then used by `plotjbdiag` to produce various graphs. The input background error statistics files are opened with the local file names `stabal96.bal` and `stabal96.cv` and a namelist variable `gsizein`, describing the grid resolution, and namelist variables describing the vertical levels (`ahalf_nl` and `bhalf_nl`, ) are needed.

The following vertical level definitions are available for use:

| Name | Description |
| :----------- | :------------------------------- |
|`hirlL60.def` |HIRLAM 60 level definition |
|`harmL65.def` |Harmonie 65 level definition |
|`mcL90.def` |MetCoOp 90 level definition |
|`mfL70.def` |Meteo France 70 level definition |
|`mfL90.def` |Meteo France 90 level definition |


Various information about the domain geometry is extracted from the input files. One example of a script to run the diagnostic tool is provided in jbdiagnose.sh. The output ASCII files are stored in the
subdirectory ”diag”.

Output files of the following forms are produced:
Expand Down
37 changes: 37 additions & 0 deletions share/levdef/mcL90.def
View file
Original file line number Diff line number Diff line change
@@ -0,0 +1,37 @@
nlev_nl = 90,
ahalf_nl = 0., 2000.00000000000, 3642.00332501264, 5177.53961669522, 6662.09739706628,
8110.48222129512, 9530.30684407722, 10926.3069617276, 12301.7468057949, 13327.5368976706,
14073.1650967992, 14651.0137450350, 15101.0938691614, 15439.1727774029, 15676.5285223989,
15822.1774031477, 15883.8711062535, 15868.6014128583, 15782.8735081322, 15632.8584439763,
15424.4777645827, 15163.4481734636, 14855.3019956629, 14505.3928595053, 14118.8924958548,
13700.7824755393, 13255.8434207066, 12788.6433961634, 12303.5266304084, 11804.6033309086,
11295.7410845528, 10780.5581361136, 10262.4186924396, 9744.43029341695, 9229.44321271586,
8720.05179519538, 8218.59759866618, 7727.17418185409, 7247.63336503162, 6781.59278265137,
6330.44454660470, 5895.36484295635, 5477.32429293581, 5077.09891958129, 4695.28157388994,
4332.29368793664, 3988.39723662434, 3663.70680407303, 3358.20166479320, 3071.73780345379,
2804.05981005167, 2554.81259947939, 2323.55291578538, 2109.76059177861, 1912.84954403207,
1732.17849179690, 1567.06139587965, 1416.77762020195, 1280.58182460273, 1157.71360251619,
1047.40688151717, 948.899108425472, 861.440243751066, 781.582196394637, 706.316393310713,
635.553940344025, 569.203061913437, 507.169152912668, 449.354814063341, 395.659869445242,
345.981364548398, 300.213542705108, 258.247797130129, 219.972594963491, 185.273368594718,
154.032368019741, 126.128465853892, 101.436903586727, 79.8289632206295, 61.1715417388026,
45.3265954120208, 32.1504040294518, 21.4925763528995, 13.1946662044524, 7.08816797743150,
2.99144588773448, 0.704631299452974, 0.000000, 0.000000, 0.000000, 0.000000,
bhalf_nl = 0., 0.0000000000, 0.0000000000, 0.0000000000, 0.0000000000,
0.00000000000000, 0.00000000000000, 0.00000000000000, 0.00000000000000, 0.00331491174057500,
0.00926834903200809, 0.0169561232023512, 0.0261759225583379, 0.0367548025100129, 0.0485647110232709,
0.0615002136004181, 0.0754685218677021, 0.0903844465059903, 0.106167664155872, 0.122741192790625,
0.140030545460615, 0.157963283623013, 0.176468812519666, 0.195478324359616, 0.214924830326513,
0.234743243204367, 0.254870485263989, 0.275245604349600, 0.295809886668503, 0.316506958638316,
0.337282872882441, 0.358086175445440, 0.378867952751117, 0.399581857892898, 0.420184116626371,
0.440633513995234, 0.460891362913670, 0.480921456286709, 0.500690004403914, 0.520165559413050,
0.539318928687490, 0.558123078858878, 0.576553032207206, 0.594585756994382, 0.612200053202724,
0.629376435003770, 0.646097011140782, 0.662345364264847, 0.678106430123150, 0.693366377361291,
0.708112488571589, 0.722333043097408, 0.736017201990554, 0.749154895415224, 0.761736712697979,
0.773753795138595, 0.785197731621296, 0.796060456999156, 0.806334153166072, 0.816011152679994,
0.825083844757460, 0.833544583422562, 0.841385597562489, 0.848858814996254, 0.856220374783933,
0.863469599962813, 0.870605738764383, 0.877627952467973, 0.884535301531071, 0.891326729598322,
0.898001044880795, 0.904556898248369, 0.910992757174897, 0.917306874393843, 0.923497249724166,
0.929561582953762, 0.935497214826483, 0.941301051912107, 0.946969469178671, 0.952498180956417,
0.957882065797460, 0.963114921767607, 0.968189112378189, 0.973095031707568, 0.977820250871524,
0.982348053901130, 0.986654660599537, 0.990703096568891, 0.994431543129278, 0.997630260155212, 1.0000000000,

0 comments on commit 5ea27d0

Please sign in to comment.