source: trunk/SRC/ToBeReviewed/LECTURE/GRIB/read_grib_pds.pro @ 214

Last change on this file since 214 was 157, checked in by navarro, 18 years ago

header improvements + xxx doc

  • Property svn:keywords set to Id
File size: 1.1 KB
Line 
1;+
2; @file_comments
3;
4;
5; @categories
6;
7;
8; @param NUM
9;
10;
11; @param RECSTART
12;
13;
14; @returns
15;
16;
17; @restrictions
18;
19;
20; @examples
21;
22;
23; @history
24;
25;
26; @version
27; $Id$
28;-
29FUNCTION read_grib_pds, num, recstart
30;
31;
32  compile_opt idl2, strictarrsubs
33;
34  offset = recstart+8
35  a = assoc(num, bytarr(1, /nozero), offset-1)
36;
37  sizepds = bit2int([binary(a[1]), binary(a[2]), binary(a[3])])
38 ;
39  paramtableversion = (a[4])[0]
40 ;
41  centerid = (a[5])[0]
42 ;
43  procid = (a[6])[0]
44 ;
45  gridid = (a[7])[0]
46;
47  flag = binary(a[8])
48  gdsnotomitted = flag[0]
49  bmsnotomitted = flag[1]
50;
51  paramunitid = (a[9])[0]
52;
53  levtype = (a[10])[0] 
54  levalue1 = (a[11])[0]
55  levalue2 = (a[12])[0]
56;
57  year = (a[13])[0]
58  month = (a[14])[0]
59  day = (a[15])[0]
60  hour = (a[16])[0]
61  minute = (a[17])[0]
62  timeunit = (a[18])[0]
63  p1 = (a[19])[0]
64  p2 = (a[20])[0]
65  timerange = (a[21])[0]
66  n1 = (a[22])[0]
67  n2 = (a[23])[0]
68  nbmiss = (a[24])[0]
69  century = (a[25])[0]
70;
71  subcenterid = (a[26])[0]
72;
73  d =  bit2int([binary(a[27]), binary(a[28])], /checkneg)
74;.........
75
76
77  RETURN, {size:sizepds, gdsnotomitted:gdsnotomitted, bmsnotomitted:bmsnotomitted, d:d}
78END
Note: See TracBrowser for help on using the repository browser.