source: trunk/ToBeReviewed/LECTURE/GRIB/scan_grib_nbits.pro @ 67

Last change on this file since 67 was 67, checked in by pinsard, 18 years ago

miscellaneous modifications according to cerbere.lodyc.jussieu.fr: /usr/home/smasson/SAXO_RD/

  • Property svn:executable set to *
File size: 932 bytes
Line 
1FUNCTION scan_grib_nbits, num, recstart
2
3;
4  nrec =  n_elements(recstart)
5  nbits = bytarr(nrec)
6;
7  FOR i = 0L, nrec-1 DO BEGIN
8    offset = recstart[i]+8
9    a = assoc(num, bytarr(1, /nozero), offset-1)
10    sizepds = bit2int([binary(a[1]), binary(a[2]), binary(a[3])])
11    flag = binary(a[8])
12    gdsnotomitted = flag[0]
13    bmsnotomitted = flag[1]
14    ddd =  bit2int([binary(a[27]), binary(a[28])], /checkneg)
15;
16    offset = offset+sizepds
17;
18    IF gdsnotomitted THEN BEGIN
19      a = assoc(num, bytarr(1, /nozero), offset-1)
20      sizedds = bit2int([binary(a[1]), binary(a[2]), binary(a[3])])
21      offset = offset+sizedds
22    ENDIF
23;
24    IF bmsnotomitted THEN BEGIN
25      a = assoc(num, bytarr(1, /nozero), offset-1)
26      sizebms = bit2int([binary(a[1]), binary(a[2]), binary(a[3])])
27      offset = offset+sizebms
28    ENDIF
29
30    a = assoc(num, bytarr(1, /nozero), offset-1)
31    nbits[i] = a[11]
32
33  ENDFOR
34;
35  RETURN, nbits
36END
Note: See TracBrowser for help on using the repository browser.