Re: reading ascii vertical sounding data (Upper air data from Univ. of Wyoming)

From: Mary Haley <haley_at_nyahnyahspammersnyahnyah>
Date: Fri Jan 04 2013 - 13:44:40 MST

Basit,

This is not a trivial file to read if there's no "consistency" with how the lines are written.

I wrote a "lazy" script that reads the file. I say "lazy" because it assumes that once it reaches a line with 11 values on it, then it stores the data values to a temporary variable.

It doesn't try to find any specific header or footer lines. If you need this data to be split up depending on which section it is written in,
then hopefully you can get enough from this script to modify it yourself. It also writes the data to a file.

--Mary

On Jan 3, 2013, at 11:10 PM, BasitAli Khan wrote:

> Hi Dave,
>
> Its only one ascii file with one years data of 365 days with two soundings
> each day. Every sounding has both head and footer just like the snippet of
> the data I copied in my first email. The number of rows (pressure levels)
> in each sounding also varies on different days. Since the header and
> footer repeats (month, day, time also varies in some part of the header)
> for every sounding, none of the example including asc4.txt applies to my
> case.
>
> I will greatly appreciate if you could give some ideas, clues, advise to
> handle this dataset.
>
> Best regards,
>
> ----
> Basit A. Khan, Ph.D.
> Postdoctoral Research Fellow
> Division of Physical Sciences & Engineering
> Office# 3204, Level 3, Building 1,
> King Abdullah University of Science & Technology
> 4700 King Abdullah Blvd, Box 2753, Thuwal 23955 6900,
> Kingdom of Saudi Arabia.
>
> Office: +966(0)2 808 0276, Mobile: +966(0)5 0860 3617
> E-mail: basitali.khan@kaust.edu.sa
> Skype name: basit.a.khan
>
>
>
>
>
>
>
>
>
> On 1/4/13 7:01 AM, "Dave Allured - NOAA Affiliate" <dave.allured@noaa.gov>
> wrote:
>
>> Basit,
>>
>> It is important to know what exactly is in a single file. Is there
>> just one header and one footer per file? It is not clear whether you
>> first posted a single ascii file, or two or three files pasted
>> together.
>>
>> If there is only one header and one footer, are the number of header
>> lines and footer lines always the same? If the answer is "yes" to all
>> of these, then I think you can use a variation of example asc4.txt on
>> the Read ASCII example page. Will this work?
>>
>> http://www.ncl.ucar.edu/Applications/read_ascii.shtml#asc4
>>
>> --Dave
>> (Please reply to list only!)
>>
>> On Wed, Jan 2, 2013 at 5:49 AM, BasitAli Khan
>> <BasitAli.Khan@kaust.edu.sa> wrote:
>>> Dear all,
>>>
>>> I have vertical sounding data of the following type (upper air - Wyoming
>>> Univ.) that i want to load and process in ncl. Its one year data and i
>>> have
>>> two soundings for each day. For every sounding there is header and
>>> footer on
>>> top and bottom of the numeric data. I am just wondering as how can i
>>> load
>>> ONLY numeric data fields into variables while ignoring all the header
>>> and
>>> footer info. I am using asciiread function to load this file as string
>>> datatype but don't know how to tell ncl to exclude any redundant info.
>>> Unfortunately, I could not find any relevant examples on
>>> http://www.ncl.ucar.edu/Applications/read_ascii.shtml and would greatly
>>> appreciate if anybody could advise a solution.
>>>
>>> Best regards,
>>> Basit
>>>
>>>
>>> 41024 OEJN Jeddah (King Abdul Aziz) Observations at 00Z 02 Jan 2012
>>>
>>>
>>> -------------------------------------------------------------------------
>>> ----
>>> PRES HGHT TEMP DWPT RELH MIXR DRCT SKNT THTA THTE
>>> THTV
>>> hPa m C C % g/kg deg knot K K
>>> K
>>>
>>> -------------------------------------------------------------------------
>>> ----
>>> 1013.0 17 22.8 16.8 69 12.02 20 10 294.9 329.5
>>> 297.0
>>> 1012.0 26 22.8 16.9 70 12.15 350 14 294.9 330.0
>>> 297.1
>>> 1000.0 136 22.4 18.8 80 13.85 345 16 295.6 335.5
>>> 298.0
>>> 984.0 275 21.3 18.4 83 13.73 345 20 295.8 335.5
>>> 298.3
>>> 930.0 763 17.6 17.0 97 13.31 350 14 296.8 335.4
>>> 299.2
>>>
>>> Station information and sounding indices
>>>
>>> Station identifier: OEJN
>>> Station number: 41024
>>> Observation time: 120102/0000
>>> Station latitude: 21.70
>>> Station longitude: 39.18
>>> Station elevation: 17.0
>>> 41024 OEJN Jeddah (King Abdul Aziz) Observations at 12Z 02 Jan 2012
>>>
>>>
>>> -------------------------------------------------------------------------
>>> ----
>>> PRES HGHT TEMP DWPT RELH MIXR DRCT SKNT THTA THTE
>>> THTV
>>> hPa m C C % g/kg deg knot K K
>>> K
>>>
>>> -------------------------------------------------------------------------
>>> ----
>>> 1012.0 17 27.0 14.0 45 10.02 320 14 299.1 328.7
>>> 300.9
>>> 1011.0 26 24.6 12.6 47 9.14 322 15 296.8 323.6
>>> 298.5
>>> 1004.0 90 24.2 13.5 51 9.77 340 19 297.0 325.6
>>> 298.8
>>> 1000.0 126 24.0 14.0 54 10.14 335 19 297.1 326.8
>>> 299.0
>>> 943.0 633 18.7 12.7 68 9.89 325 16 296.8 325.7
>>> 298.6
>>> 927.0 781 17.2 12.3 73 9.81 340 18 296.7 325.4
>>> 298.4
>>> 925.0 800 17.0 12.3 74 9.80 345 18 296.7 325.3
>>> 298.4
>>> 907.0 967 15.7 11.8 78 9.71 345 14 297.0 325.4
>>> 298.7
>>> 893.0 1100 14.6 11.5 82 9.63 20 10 297.2 325.4
>>> 298.9
>>> 887.0 1157 14.8 10.5 75 9.06 9 5 298.0 324.7
>>> 299.6
>>> Station information and sounding indices
>>>
>>> Station identifier: OEJN
>>> Station number: 41024
>>> Observation time: 120102/1200
>>> Station latitude: 21.70
>>> Station longitude: 39.18
>>> Station elevation: 17.0
>>> Showalter index: 10.08
>>> Lifted index: 5.13
>>> LIFT computed using virtual temperature: 4.91
>>>
>>> 41024 OEJN Jeddah (King Abdul Aziz) Observations at 00Z 03 Jan 2012
>>>
>>>
>>> -------------------------------------------------------------------------
>>> ----
>>> PRES HGHT TEMP DWPT RELH MIXR DRCT SKNT THTA THTE
>>> THTV
>>> hPa m C C % g/kg deg knot K K
>>> K
>>>
>>> -------------------------------------------------------------------------
>>> ----
>>> 1013.0 17 20.0 13.0 64 9.37 50 5 292.1 318.9
>>> 293.7
>>> 1011.0 35 21.6 15.6 69 11.14 50 5 293.8 325.9
>>> 295.8
>>> 1000.0 134 21.4 14.4 64 10.41 50 4 294.6 324.6
>>> 296.4
>>> 964.0 450 18.9 14.5 76 10.91 40 2 295.1 326.6
>>> 297.0
>>> 925.0 806 16.0 14.7 92 11.49 5 4 295.7 328.9
>>> 297.7
>>> 915.0 898 15.0 14.3 96 11.32 335 5 295.6 328.3
>>> 297.6
>>> 893.0 1105 14.6 13.1 91 10.72 268 8 297.2 328.5
>>> 299.1
>>> 876.0 1268 15.7 13.1 85 10.93 215 10 300.0 332.2
>>> 301.9
>>> 868.0 1346 16.2 13.1 82 11.03 197 15 301.3 334.0
>>> 303.3
>>>
>>>
>>> ________________________________
>
>
> ________________________________
>
> This message and its contents including attachments are intended solely for the original recipient. If you are not the intended recipient or have received this message in error, please notify me immediately and delete this message from your computer system. Any unauthorized use or distribution is prohibited. Please consider the environment before printing this email.
> _______________________________________________
> ncl-talk mailing list
> List instructions, subscriber options, unsubscribe:
> http://mailman.ucar.edu/mailman/listinfo/ncl-talk

_______________________________________________
ncl-talk mailing list
List instructions, subscriber options, unsubscribe:
http://mailman.ucar.edu/mailman/listinfo/ncl-talk

Received on Fri Jan 4 13:44:57 2013

This archive was generated by hypermail 2.1.8 : Fri Jan 04 2013 - 15:34:10 MST