[ncl-talk] Need a copy of rdIGRA2.ncl

Dennis Shea shea at ucar.edu
Mon Feb 18 20:03:55 MST 2019


Hi Michael

I have placed some IGRA2 and (older)  IGRA scripts and data onto ftp.


Of course, you will have to change the following paths.

   dirIGRA2 = "/Users/shea/Data/IGRA2/" ;---Directory with country and
station files
   dirRAOB  = dirIGRA2+"RAOB/"             ;---Directory with RAOB files

Some files will need to be unzipped:

   %> gzip -d .....

ftp ftp.cgd.ucar.edu
anonymous
your_email
cd pub/shea/IGRA2
ls
mget *
cd RAOB
mget *

cd ../../IGRA
ls
mget *
quit

==========ftp/pub/shea/IGRA2/RAOB> cd ..
/ftp/pub/shea/*IGRA2*

-rw-r--r-- 1 shea cgdcas   3008 Feb 18 19:38 igra2-country-list.txt
-rw-r--r-- 1 shea cgdcas 245729 Feb 18 19:38 igra2-station-list.txt
drwxr-xr-x 2 shea cgdcas   3896 Feb 18 19:43 *RAOB/*
-rw-r--r-- 1 shea cgdcas   8489 Feb 18 19:38 *rdIGRA2.ncl*

/ftp/pub/shea/IGRA2/RAOB

-rw-r--r-- 1 shea cgdcas 1649980 Feb 18 19:41 GBM00064500-data.txt.gz
-rw-r--r-- 1 shea cgdcas     357 Feb 18 19:40 GBM-test.txt.gz

/ftp/pub/shea*/IGRA*

-rw-r--r-- 1 shea cgdcas 7702995 Feb 18 19:09 98646.dat
-rw-r--r-- 1 shea cgdcas 5090164 Feb 18 19:09 98646.nc
-rw-r--r-- 1 shea cgdcas    8139 Feb 18 19:09 derived2netCDF.ncl
-rw-r--r-- 1 shea cgdcas    3625 Feb 18 19:09 derived-countries.txt
-rw-r--r-- 1 shea cgdcas   90774 Feb 18 19:09 derived-stations.txt
-rw-r--r-- 1 shea cgdcas    3650 Feb 18 19:09 foo.txt
-rw-r--r-- 1 shea cgdcas  118156 Feb 18 19:09 IGRA.INDEX.nc
-rw-r--r-- 1 shea cgdcas   18709 Feb 18 19:09 rdIGRA_sonde.ncl
-rw-r--r-- 1 shea cgdcas    8638 Feb 18 19:09 readme.txt


======

 I think you can just run the sample to get a feel for it.

%> ncl *rdIGRA2.ncl*

HTH

Good Luck
D

On Mon, Feb 18, 2019 at 4:33 PM Michael Notaro <mnotaro at wisc.edu> wrote:

> Dennis, Appo, or other NCL users,
>
> Does anyone have a copy of this rdIGRA2.ncl script,
>
> alluded to below?
>
> Thanks, Michael
>
>
>
> [ncl-talk] Help on reading IGRA 2 Sounding *Dennis Shea* shea at ucar.edu
> <ncl-talk%40ucar.edu?Subject=Re:%20%5Bncl-talk%5D%20Help%20on%20reading%20IGRA%202%20Sounding&In-Reply-To=%3CCAOF1d_7A9kaSrwUFZ0qBLEnMyEnW8swM415QhX5jpYGuZUuHQw%40mail.gmail.com%3E>
> *Sat Apr 29 14:29:54 MDT 2017*
>
>
>    - Previous message: [ncl-talk] Help on reading IGRA 2 Sounding
>    <http://mailman.ucar.edu/pipermail/ncl-talk/2017-April/008790.html>
>    - Next message: [ncl-talk] Error with Pressure vs latitude plot
>    <http://mailman.ucar.edu/pipermail/ncl-talk/2017-April/008797.html>
>    - *Messages sorted by:* [ date ]
>    <http://mailman.ucar.edu/pipermail/ncl-talk/2017-April/date.html#8801> [
>    thread ]
>    <http://mailman.ucar.edu/pipermail/ncl-talk/2017-April/thread.html#8801>
>     [ subject ]
>    <http://mailman.ucar.edu/pipermail/ncl-talk/2017-April/subject.html#8801>
>     [ author ]
>    <http://mailman.ucar.edu/pipermail/ncl-talk/2017-April/author.html#8801>
>
> ------------------------------
>
> 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 <http://mailman.ucar.edu/mailman/listinfo/ncl-talk>> wrote:
>
> >* I am trying to read IGRA 2 Sounding Data Files.
> *>* For example i want to read this file ftp://ftp.ncdc.noaa.gov/ <ftp://ftp.ncdc.noaa.gov/>
> *>* pub/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 <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*
>
>
>
> Michael Notaro
> Associate Director
> Nelson Institute Center for Climatic Research
> University of Wisconsin-Madison
> Phone: (608) 261-1503
> Email: mnotaro at wisc.edu
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20190218/bc7b66df/attachment.html>


More information about the ncl-talk mailing list