From 5ea27d0045ca66e2b610058f5e0d290ff32ec6e8 Mon Sep 17 00:00:00 2001 From: "Eoin Whelan (Irish Met Service)" Date: Tue, 26 Nov 2024 07:13:12 +0000 Subject: [PATCH] Add MetCoOp L90 defintion and some updated doc (#54) * Add MetCoOp L90 defintion and some updated doc * doc markdown correct. Thanks Roel! --- doc/src/jbdiagnose.md | 15 ++++++++++++++- share/levdef/mcL90.def | 37 +++++++++++++++++++++++++++++++++++++ 2 files changed, 51 insertions(+), 1 deletion(-) create mode 100644 share/levdef/mcL90.def diff --git a/doc/src/jbdiagnose.md b/doc/src/jbdiagnose.md index b70d64c..5f11e92 100644 --- a/doc/src/jbdiagnose.md +++ b/doc/src/jbdiagnose.md @@ -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: diff --git a/share/levdef/mcL90.def b/share/levdef/mcL90.def new file mode 100644 index 0000000..e2175c0 --- /dev/null +++ b/share/levdef/mcL90.def @@ -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,