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:       cl...@ecy.wa.gov
Air Quality Modeler             INTERNET:       cl...@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@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.


______________________________________________
R-help@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.

Reply via email to