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.