observingclouds / eurec4a_snd Goto Github PK

I would like to reuse this code to convert radiosondes from the same different platforms (Vaisala, MeteoModem) for a field experiment called SOFOG3D (SOuth westFOGs 3D experiment for processes study) conducted by Meteo France/CNRM on fog properties in France during last winter.

Maybe a GNU/GPL license is a good one: this allows the use of the code (all-purpose), its modification and the re-publication of the modified code. If the code is modified, the source must be cited and the modifications must be listed.

For know I planned to cite the software in the NetCDF file in the attribute created_with, for example created_with = "L1_mwx.py 10.5281/zenodo.4391098" with adding the Licence.

Thanks a lot for all the work already done!

The time unit of "flight_time" is currently written to be "seconds since 1970-01-01" although it is actually "seconds since launch_time".

This should be fixed to be consistent with the other L2 datasets.

One simulated sounding made it into the sounding dataset of the Ron Brown, namely '20200117_2202'.

This needs to be deleted.

According to the Operational Guideline quicklooks should have a filename as follows:

File nomenclature
The file name convention is like this PlatformName_Instrument_Variable_YYYYMMDD_HHmmSS_SpaceCov_SpaceRes_Lev el_ForcastNB.ext

• PlatformName (mandatory) : Name of the platform
• Instrument and Variable (at least of this information)
• YYYYMMDD_HHmmSS (mandatory) : Datetime of the quicklook formatted like
  YYYYMMDD_HHmmSS, seconds are optional. ex: 20190702_2345
• SpaceCov (optional): Space resolution ex. : 180W-180E-40S-40N (this is used to plot
  the quicklook on the interactive map)
• SpaceRes (optional): Space resolution, ex : 025 for 0,25° ; 001 for 0,01° ; 1km for 1 km
• Level (optional): pressure level or altitude, ex : 50hPa, 200m
• ForcastNB (optional) : Used for forecast quicklooks. It indicates the forecast delay
  after the forecasts production. Ex: 06 (forecast 6hours after the run), 12 (forecast
  12hours after the run)

The soundings "20200114_0409" and "20200114_0242" are split, although they are essentially one sounding.

Check whether the used formulas are correct/precise enough.

Since the sonde measures relative humidity directly, it should be checked, which formula is used in order to derive the dewpoint to be consistent.

Currently this package includes the software to convert BUFR, MWX, COR files to netCDF files, but also specific scripts and configurations that were specific to the EUREC4A campaign.

It would be good to separate these and have a software that can be used for future campaigns without having too much 'clutter'.

Currently only files with one continues profile are supported but the radiosonde sounding files should include also the descending branch. This should automatically be detected and written in a separate output file.

Time reference should be the 2020-01-01 to omit any issues arising due to the disregard of leap-seconds.

Describe the bug
I get two unhelpful warnings when running 'sounding_converter', even though it does convert my BUFR file. One warning is about not being a git repository, the second is that 'No version of the converter could be found.'

To Reproduce
Steps to reproduce the behavior:
Run

test_eurec4a_sounding_converter sdeszoek$ sounding_converter -i ./BUFR/bufr309052_all_20190328_1813_0.bfr -c ./meta_information.ini`

Result:
These two error/warning messages:

fatal: not a git repository (or any of the parent directories): .git
WARNING - root - main - No version of the converter could be found! Please consider the installation via conda or if this is not working clone the git repository

Expected behavior
I think the INFO lines, though verbose, are the expected correct behavior.

INFO - root - main - Version of script: 0.0.43 (conda package), -- (git version)
INFO - root - main - Files to process ['bufr309052_all_20190328_1813_0.bfr']
INFO - root - main - Reading file number 0
INFO - root - main - DONE: BUFR/bufr309052_all_20190328_1813_0.bfr converted to NC/RHB_SoundingDescentProfile_15N54W_20190328_1813.nc

Error message
Paste error message if available (make sure to remove personal details)

Environment (please complete the following information):

• OS: MacOS 'Darwin ... 18.7.0 Darwin Kernel Version 18.7.0: Sat Oct 12 00:02:19 PDT 2019; root:xnu-4903.278.12~1/RELEASE_X86_64 x86_64'
• python-version 3.7.3
• eurec4a-version 0.0.43 (conda list | grep eurec4a)

Create reader for meteomodem raw data to directly convert the data to netCDF without the need to use the BUFR files and add uncertainty to the humidity measurements ( BUFR outputs Td and not the measured RH).

The pressure unit in L1 should be also Pa and not hPa

The calculation of ascent rate is not yet implemented