search

ARPA-SIMC / libsim

Command-line utilities and fortran modules for meteorological data processing
GNU General Public License v2.0
7 stars 1 forks source link

v7d_transform daily average error with B14198 #65

Closed edigiacomo closed 5 years ago

edigiacomo commented 5 years ago

With libsim-6.3.0-1 and dballe-8-0.14.20181213gitd19d9303.fc28 (e.g. ventotto.metarpa), the daily average is not correct:

$ v7d_transform --variable-list=B14198 --comp-stat-proc=0:0 --comp-step='1 00' --comp-frac-valid=0.1 --comp-start="2018-12-08 00:00" --comp-stop="2018-12-10 00:00"  --input-format=BUFR --output-format=BUFR camse-rad1havg.bufr camse-rad24havg.bufr
$ dbamsg dump --interpreted camse-rad24havg.bufr 
#0[0] generic message, 44.60016,12.07738, dt: 2018-12-09T00:00:00, 2 contexts:
Level -,-,-,-, tr -,-,-,
001019 LONG STATION OR SITE NAME(CCITTIA5): Camse
001194 [SIM] Report mnemonic(CCITTIA5): locali
004001 YEAR(YEAR): 2018
004002 MONTH(MONTH): 12
004003 DAY(DAY): 9
004004 HOUR(HOUR): 0
004005 MINUTE(MINUTE): 0
004006 SECOND(SECOND): 0
005001 LATITUDE (HIGH ACCURACY)(DEGREE): 44.60016
006001 LONGITUDE (HIGH ACCURACY)(DEGREE): 12.07738
007030 HEIGHT OF STATION GROUND ABOVE MEAN SEA LEVEL (SEE NOTE 3)(M): 1.0
007031 HEIGHT OF BAROMETER ABOVE MEAN SEA LEVEL (SEE NOTE 4)(M): 0.0
Level 1,-,-,-, tr 0,0,86400
014198 [SIM] Visible radiation (downward)(W/M**2): 124.483

With libsim-6.2.8-1 and dballe-7.35-1:

$ v7d_transform --variable-list=B14198 --comp-stat-proc=0:0 --comp-step='1 00' --comp-frac-valid=0.1 --comp-start="2018-12-08 00:00" --comp-stop="2018-12-10 00:00"  --input-format=BUFR --output-format=BUFR camse-rad1havg.bufr camse-rad24havg.bufr
$ dbamsg dump --interpreted camse-rad24havg.bufr 
#0[0] generic message, rep_memo: , 44.60016,12.07738, ident: , dt: 2018-12-09T00:00:00, 2 contexts:
Level 1,-,-,-, tr 0,0,86400 1 vars:
014198 [SIM] Visible radiation (downward)(W/M**2): 56.201
Level -,-,-,-, tr -,-,-, 6 vars:
001019 LONG STATION OR SITE NAME(CCITTIA5): Camse
001194 [SIM] Report mnemonic(CCITTIA5): locali
005001 LATITUDE (HIGH ACCURACY)(DEGREE): 44.60016
006001 LONGITUDE (HIGH ACCURACY)(DEGREE): 12.07738
007030 HEIGHT OF STATION GROUND ABOVE MEAN SEA LEVEL (SEE NOTE 3)(M): -1.0
007031 HEIGHT OF BAROMETER ABOVE MEAN SEA LEVEL (SEE NOTE 4)(M): 0.0
edigiacomo commented 5 years ago

I forgot the attachment camse-rad1havg.bufr.zip

edigiacomo commented 5 years ago

Sorry, it seems that it's a dballe bug: with libsim-6.3.0-1 and dballe-7.36-1 the result is correct.