[ncl-talk] Help on reading IGRA 2 Sounding
Appo derbetini
appopson4 at gmail.com
Mon May 1 05:02:00 MDT 2017
Thank you very much Denis,
You did exactly what i fighting with.
Now, I looking how to calculate mean over selected days.
First, I am trying to fill unavailable level with missing value.
So that each day will have the same number of level.
In this case I can easily perform my analysis.
Best regards,
Appo
2017-04-29 21:29 GMT+01:00 Dennis Shea <shea at ucar.edu>:
> See attached
>
> For example:
>
> %> cd /Users/appo/Data/IGRA2
>
>
> drwxr-xr-x 4 2200 136 Apr 29 13:43 RAOB/
> -rw-r--r--@ 1 2200 3008 Apr 29 12:06 igra2-country-list.txt
> -rw-r--r--@ 1 2200 245729 Apr 29 12:07 igra2-station-list.txt
> -rw-r--r-- 1 2200 8489 Apr 29 14:21 rdIGRA2.ncl
>
> The RAOB directory could have many files but for test just one.
>
> -rw-rw-r--@ 1 2200 12956441 Apr 27 16:42 GBM00064500-data.txt
>
> %> ncl rdIGRA2.ncl
>
> [The actual print out looks nicer. ]
>
>
> 0)
> (0) ------------- nr=0 id=GBM00064500 lat=0.45 lon=9.4167 ------------
> (0) yyyy=1973 mm=1 dd=1 hh=12 rt=9999 nlev=6
> (0)
> (0) l1 l2 etim pres pf zgph zf temp tf rhum dpdp
> wspd wdir
> (0) 1 0 -9999 85000 -9999 -9999.0 -9999.0 -9999.0 4
> 60.0
> (1) 1 0 -9999 70000 -9999 -9999.0 -9999.0 -9999.0 14
> 60.0
> (2) 1 0 -9999 50000 -9999 -9999.0 -9999.0 -9999.0 16
> 80.0
> (3) 1 0 -9999 40000 -9999 -9999.0 -9999.0 -9999.0 14
> 80.0
> (4) 1 0 -9999 30000 -9999 -9999.0 -9999.0 -9999.0 8
> 70.0
> (5) 1 0 -9999 25000 -9999 -9999.0 -9999.0 -9999.0 5
> 60.0
> (0)
> (0) ------------- nr=1 id=GBM00064500 lat=0.45 lon=9.4167 ------------
> (0) yyyy=1973 mm=1 dd=2 hh=12 rt=9999 nlev=10
> (0)
> (0) l1 l2 etim pres pf zgph zf temp tf rhum dpdp wspd
> wdir
> (0) 1 0 -9999 85000 -9999 -9999.0 -9999.0 -9999.0 3
> 100.0
> (1) 1 0 -9999 70000 -9999 -9999.0 -9999.0 -9999.0 13
> 60.0
> (2) 1 0 -9999 50000 -9999 -9999.0 -9999.0 -9999.0 14
> 80.0
> (3) 1 0 -9999 40000 -9999 -9999.0 -9999.0 -9999.0 9
> 50.0
> (4) 1 0 -9999 30000 -9999 -9999.0 -9999.0 -9999.0 3
> 0.0
> (5) 1 0 -9999 25000 -9999 -9999.0 -9999.0 -9999.0 7
> 60.0
> (6) 1 0 -9999 20000 -9999 -9999.0 -9999.0 -9999.0 10
> 80.0
> (7) 1 0 -9999 15000 -9999 -9999.0 -9999.0 -9999.0 18
> 100.0
> (8) 1 0 -9999 10000 -9999 -9999.0 -9999.0 -9999.0 17
> 90.0
> (9) 1 0 -9999 7000 -9999 -9999.0 -9999.0 -9999.0 22
> 70.0
>
> [SNIP]
>
> 0)
> (0) ------------- nr=12141 id=GBM00064500 lat=0.45 lon=9.4167
> ------------
> (0) yyyy=2017 mm=4 dd=26 hh=12 rt=1046 nlev=74
> (0)
> (0) l1 l2 etim pres pf zgph zf temp tf rhum dpdp wspd
> wdir
> (0) 2 1 -9999 101100 B -9999 30.2 B -9999.0 10.0 5
> 180.0
> (1) 1 0 -9999 100000 110 B 27.8 B -9999.0 10.0 5
> 150.0
> (2) 1 0 -9999 92500 796 B 23.4 B -9999.0 10.0 9
> 125.0
> (3) 1 0 -9999 85000 1530 B 20.0 B -9999.0 13.0 7
> 150.0
> (4) 1 0 -9999 70000 3173 B 9.6 B -9999.0 6.0 15
> 80.0
> (5) 1 0 -9999 50000 5900 B -6.1 B -9999.0 3.9 5
> 65.0
> (6) 1 0 -9999 40000 7610 B -15.7 B -9999.0 6.0 7
> 75.0
> (7) 1 0 -9999 30000 9730 B -30.5 B -9999.0 18.0 8
> 40.0
> (8) 1 0 -9999 25000 10990 B -40.9 B -9999.0 13.0 7
> 25.0
> (9) 1 0 -9999 20000 12470 B -53.7 B -9999.0 31.0 9
> 70.0
> (10) 1 0 -9999 15000 14250 B -68.7 B -9999.0 27.0 16
> 75.0
> (11) 2 2 -9999 10100 -9999 -77.5 B -9999.0 24.0 7
> 30.0
> (12) 1 0 -9999 10000 16630 B -77.5 B -9999.0 24.0 6
> 35.0
> (13) 2 0 -9999 9320 -9999 -75.5 B -9999.0 25.0 -9999
> -9999.0
> (14) 2 0 -9999 9180 -9999 -9999.0 -9999.0 -9999.0 12
> 310.0
> (15) 2 0 -9999 8520 -9999 -9999.0 -9999.0 -9999.0 18
> 280.0
> (16) 2 2 -9999 8500 -9999 -78.3 A -9999.0 24.0 18
> 275.0
> (17) 2 0 -9999 7970 -9999 -9999.0 -9999.0 -9999.0 11
> 260.0
> (18) 2 0 -9999 7440 -9999 -76.1 A -9999.0 25.0 -9999
> -9999.0
> (19) 2 0 -9999 7230 -9999 -9999.0 -9999.0 -9999.0 10
> 320.0
> (20) 1 0 -9999 7000 18660 B -78.7 A -9999.0 24.0 8
> 325.0
> (21) 2 0 -9999 6740 -9999 -80.1 A -9999.0 24.0 -9999
> -9999.0
> (22) 2 0 -9999 6500 -9999 -9999.0 -9999.0 -9999.0 3
> 235.0
> (23) 2 0 -9999 6150 -9999 -9999.0 -9999.0 -9999.0 7
> 190.0
> (24) 2 0 -9999 5890 -9999 -70.7 A -9999.0 26.0 -9999
> -9999.0
> (25) 2 0 -9999 5820 -9999 -9999.0 -9999.0 -9999.0 12
> 240.0
> (26) 2 0 -9999 5060 -9999 -9999.0 -9999.0 -9999.0 28
> 240.0
> (27) 1 0 -9999 5000 20630 B -70.9 A -9999.0 26.0 28
> 240.0
> (28) 2 0 -9999 4870 -9999 -70.7 A -9999.0 26.0 -9999
> -9999.0
> (29) 2 0 -9999 4700 -9999 -9999.0 -9999.0 -9999.0 10
> 280.0
> (30) 2 0 -9999 4610 -9999 -62.3 A -9999.0 29.0 -9999
> -9999.0
> (31) 2 0 -9999 4410 -9999 -63.9 A -9999.0 28.0 -9999
> -9999.0
> (32) 2 0 -9999 4230 -9999 -9999.0 -9999.0 -9999.0 22
> 305.0
> (33) 2 0 -9999 3520 -9999 -59.9 A -9999.0 30.0 -9999
> -9999.0
> (34) 2 0 -9999 3160 -9999 -54.7 A -9999.0 31.0 -9999
> -9999.0
> (35) 1 0 -9999 3000 23810 B -55.1 A -9999.0 31.0 18
> 265.0
> (36) 2 0 -9999 2740 -9999 -56.7 A -9999.0 30.0 -9999
> -9999.0
> (37) 2 0 -9999 2580 -9999 -53.3 A -9999.0 32.0 -9999
> -9999.0
> (38) 2 0 -9999 2240 -9999 -9999.0 -9999.0 -9999.0 15
> 230.0
> (39) 2 0 -9999 2200 -9999 -55.7 A -9999.0 31.0 -9999
> -9999.0
> (40) 2 0 -9999 2040 -9999 -51.5 A -9999.0 32.0 -9999
> -9999.0
> (41) 1 0 -9999 2000 26400 B -51.3 A -9999.0 32.0 3
> 225.0
> (42) 2 0 -9999 1970 -9999 -9999.0 -9999.0 -9999.0 0
> 0.0
> (43) 2 0 -9999 1800 -9999 -50.7 A -9999.0 32.0 9
> 65.0
> (44) 3 0 -9999 -9999 600 -9999.0 -9999.0 -9999.0 5
> 95.0
> (45) 3 0 -9999 -9999 900 -9999.0 -9999.0 -9999.0 6
> 120.0
> (46) 3 0 -9999 -9999 1500 -9999.0 -9999.0 -9999.0 7
> 150.0
> (47) 3 0 -9999 -9999 2100 -9999.0 -9999.0 -9999.0 10
> 130.0
> (48) 3 0 -9999 -9999 3300 -9999.0 -9999.0 -9999.0 15
> 80.0
> (49) 3 0 -9999 -9999 3900 -9999.0 -9999.0 -9999.0 13
> 65.0
> (50) 3 0 -9999 -9999 4200 -9999.0 -9999.0 -9999.0 12
> 60.0
> (51) 3 0 -9999 -9999 4500 -9999.0 -9999.0 -9999.0 10
> 65.0
> (52) 3 0 -9999 -9999 5100 -9999.0 -9999.0 -9999.0 8
> 70.0
> (53) 3 0 -9999 -9999 6900 -9999.0 -9999.0 -9999.0 2
> 90.0
> (54) 3 0 -9999 -9999 8700 -9999.0 -9999.0 -9999.0 13
> 50.0
> (55) 3 0 -9999 -9999 11400 -9999.0 -9999.0 -9999.0 6
> 45.0
> (56) 3 0 -9999 -9999 14400 -9999.0 -9999.0 -9999.0 17
> 70.0
> (57) 3 0 -9999 -9999 15600 -9999.0 -9999.0 -9999.0 15
> 125.0
> (58) 3 0 -9999 -9999 16200 -9999.0 -9999.0 -9999.0 2
> 100.0
> (59) 3 0 -9999 -9999 16500 -9999.0 -9999.0 -9999.0 7
> 25.0
> (60) 3 0 -9999 -9999 17100 -9999.0 -9999.0 -9999.0 12
> 310.0
> (61) 3 0 -9999 -9999 17400 -9999.0 -9999.0 -9999.0 18
> 280.0
>
>
>
>
>
>
> On Fri, Apr 28, 2017 at 1:49 AM, Appo derbetini <appopson4 at gmail.com>
> wrote:
>
>> I am trying to read IGRA 2 Sounding Data Files.
>> For example i want to read this file ftp://ftp.ncdc.noaa.gov/p
>> ub/data/igra/data/data-por/GBM00064500-data.txt.zip
>> Description of the data format are available here
>> ftp://ftp.ncdc.noaa.gov/pub/data/igra/data/igra2-data-format.txt
>>
>> But I am have a lot of difficulties to do it.
>>
>> There is someone with a script to do it?
>>
>> Any help is welcome.
>>
>> Thank you very much in advance.
>>
>> Appo
>>
>>
>> _______________________________________________
>> ncl-talk mailing list
>> ncl-talk at ucar.edu
>> List instructions, subscriber options, unsubscribe:
>> http://mailman.ucar.edu/mailman/listinfo/ncl-talk
>>
>>
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20170501/934513bb/attachment.html
More information about the ncl-talk
mailing list