source: trunk/ToBeReviewed/PLOTS/LABEL/label_gmt.pro @ 39

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

upgrade of PLOTS/LABEL according to cerbere.lodyc.jussieu.fr: /usr/home/smasson/SAXO_RD/ : files

  • Property svn:eol-style set to native
  • Property svn:executable set to *
  • Property svn:keywords set to Author Date Id Revision
File size: 1.5 KB
Line 
1;
2; Apply GMT palette into IDL color/intervals system
3;
4PRO label_gmt, min, max, intervalle, ncoul, ncontour, level_z2d, coul
5@common
6@com_eg
7
8   IF pal_type NE '2dom' THEN BEGIN
9
10;  color defined in lec_pal_gmt.pro
11
12      ncontour = ncont_gmt
13      level_z2d = levels_gmt
14      coul = coul_gmt
15      max = max_gmt
16      intervalle =  -1
17
18   ENDIF ELSE BEGIN
19
20      ; grey_shade palette (case 1 in label)
21      IF finite(min) EQ 0 THEN read, ' Grey shade needs a min,max :', min, max
22      ncontour  = fix((max-min)/intervalle)
23      level_z2d = min + intervalle*findgen(Ncontour)
24      max=level_z2d(Ncontour-1)+intervalle
25
26      print, '     Number of contour intervals, plotting min & max ', ncontour, min, max
27      print, ' '
28
29      ; color index
30
31      IF idx_pal EQ 0 THEN BEGIN
32
33         ; build palette
34         red = lonarr(99)
35         red[*] = 255
36         red[50:98] = long((100.-float(grey_shade))/100.*255.)
37
38         IF field.origin EQ 'diff' THEN BEGIN
39     
40           ; difference plot : lighter below first negative interval
41            red[51:98] = long((100.-float(grey_shade_2))/100.*255.)
42            red(1:48) = long((100.-float(grey_shade))/100.*255.)
43            red[50] = 255
44
45         ENDIF
46
47         ; first color=black
48         ; last = white
49         red = [0, red]
50         red[99] = 255
51         ; gray scale
52         green = red
53         blue = red
54
55         tvlct, red, green, blue
56
57      ENDIF 
58 
59      mid_index = max(where(level_z2d LE fldatt.mid))
60      coul = findgen(Ncontour)+49-mid_index +2
61
62   ENDELSE
63
64
65END
Note: See TracBrowser for help on using the repository browser.