Ilan Levy
2013-Mar-05 20:11 UTC
[R] Reading Wyoming radiosonde data files with RadioSonde package
Hi,
I need to do some analysis on historic daily radiosonde data I download
from the Wyoming Univ. web page (
http://weather.uwyo.edu/upperair/sounding.html).
I am trying to use the RadioSonde package (V 1.3), but the format of the
files from Wyoming don't match what RadioSonde is expecting.
Has anyone used the Radiosonde package on the Wyoming data?
Here is a sample of the Wyoming file format:
40179 Bet Dagan Observations at 00Z 22 Feb 2013
-----------------------------------------------------------------------------
PRES HGHT TEMP DWPT RELH MIXR DRCT SKNT THTA THTE
THTV
hPa m C C % g/kg deg knot K K
K
-----------------------------------------------------------------------------
70.0 18430 -64.5 265 91 446.1
446.1
59.6 19417 -60.9 265 83 475.1
475.1
50.0 20500 -62.5 265 75 495.8
495.8
46.7 20920 -64.3 501.2
501.2
38.3 22137 -63.5 532.5
532.5
33.3 23012 -56.3 573.2
573.2
30.0 23670 -59.9 580.8
580.8
28.1 24078 -60.5 590.0
590.0
20.5 26056 -57.1 656.0
656.0
Station information and sounding indices
Station number: 40179
Observation time: 130222/0000
Station latitude: 32.00
Station longitude: 34.81
Station elevation: 35.0
Mean mixed layer potential temperature: 0.00
Mean mixed layer mixing ratio: 0.00
Here is the code I tried:
filename <- 'D:\\Data\\sounding_test3.txt'
datakey <- "------"
varkey <- " PRES"
unitkey <- " hPa"
sample.sonde <- getsonde(filename, datakey, varkey, unitkey)
Error in getsonde(filename, datakey, varkey, unitkey) :
(getsonde): could not find a unique match for the data string
Thanks,
Ilik
Win7 OS 64-bit, R version 2.13.0.
[[alternative HTML version deleted]]
Clint Bowman
2013-Mar-05 21:45 UTC
[R] Reading Wyoming radiosonde data files with RadioSonde package
Sure looks as if that second "------------------..." line is causing
another attempt at parsing for varkey and unitkey which for 70.0 and 59.6
just won't make sense to getsonde.
It's a pain but I'd experiment by removing that second
"--------..." line.
Clint Bowman INTERNET: clint at ecy.wa.gov
Air Quality Modeler INTERNET: clint at math.utah.edu
Department of Ecology VOICE: (360) 407-6815
PO Box 47600 FAX: (360) 407-7534
Olympia, WA 98504-7600
USPS: PO Box 47600, Olympia, WA 98504-7600
Parcels: 300 Desmond Drive, Lacey, WA 98503-1274
On Tue, 5 Mar 2013, Ilan Levy wrote:
> Hi,
> I need to do some analysis on historic daily radiosonde data I download
> from the Wyoming Univ. web page (
> http://weather.uwyo.edu/upperair/sounding.html).
> I am trying to use the RadioSonde package (V 1.3), but the format of the
> files from Wyoming don't match what RadioSonde is expecting.
>
> Has anyone used the Radiosonde package on the Wyoming data?
> Here is a sample of the Wyoming file format:
>
>
> 40179 Bet Dagan Observations at 00Z 22 Feb 2013
>
>
-----------------------------------------------------------------------------
> PRES HGHT TEMP DWPT RELH MIXR DRCT SKNT THTA THTE
> THTV
> hPa m C C % g/kg deg knot K K
> K
>
-----------------------------------------------------------------------------
> 70.0 18430 -64.5 265 91 446.1
> 446.1
> 59.6 19417 -60.9 265 83 475.1
> 475.1
> 50.0 20500 -62.5 265 75 495.8
> 495.8
> 46.7 20920 -64.3 501.2
> 501.2
> 38.3 22137 -63.5 532.5
> 532.5
> 33.3 23012 -56.3 573.2
> 573.2
> 30.0 23670 -59.9 580.8
> 580.8
> 28.1 24078 -60.5 590.0
> 590.0
> 20.5 26056 -57.1 656.0
> 656.0
>
>
> Station information and sounding indices
>
> Station number: 40179
> Observation time: 130222/0000
> Station latitude: 32.00
> Station longitude: 34.81
> Station elevation: 35.0
> Mean mixed layer potential temperature: 0.00
> Mean mixed layer mixing ratio: 0.00
>
> Here is the code I tried:
>
> filename <- 'D:\\Data\\sounding_test3.txt'
> datakey <- "------"
> varkey <- " PRES"
> unitkey <- " hPa"
> sample.sonde <- getsonde(filename, datakey, varkey, unitkey)
>
> Error in getsonde(filename, datakey, varkey, unitkey) :
> (getsonde): could not find a unique match for the data string
>
> Thanks,
> Ilik
>
> Win7 OS 64-bit, R version 2.13.0.
>
> [[alternative HTML version deleted]]
>
> ______________________________________________
> R-help at r-project.org mailing list
> https://stat.ethz.ch/mailman/listinfo/r-help
> PLEASE do read the posting guide
http://www.R-project.org/posting-guide.html
> and provide commented, minimal, self-contained, reproducible code.
>