[R] Reading Wyoming radiosonde data files with RadioSonde package

Clint Bowman clint at ecy.wa.gov
Tue Mar 5 22:45:09 CET 2013


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.
>



More information about the R-help mailing list